Dissertation
Slamet Widodo
Freiberg, 2013
ANALYSIS OF DYNAMIC LOADING BEHAVIOUR FOR PAVEMENT
ON SOFT SOIL
To the Faculty of
Geowissenschaften, Geotechnik und Bergbau
of the Technische Universität Bergakademie Freiberg
approved
THESIS
DOKTOR-INGENIEUR
DR.-ING.
Submitted
i
ACKNOWLEDGMENTS
The author wishes to express his sincere gratitude to Prof. Dr.-Ing. Herbert Klapperich, the
chairperson of his dissertation committee, for his very kind encouragement, guidance,
friendship and inspiration given throughout the graduate work. The author is also grateful to
the other dissertation committee members, Prof. Dr.rer.nat. Rafig Azzam and Prof. Dr.
Masyhur Irsyam for their very precious time for reading this dissertation.
The author is very grateful to Dr.-Ing. Habil. Nandor Tamaskovics for his kind
encouragement, competent guidance, care and always inspires him. Thanks to all of the
author’s colleagues and staff at Soil Mechanics laboratory, Institut für Geotechnik for their
help, co-operation and valuable discussions and especially to Dr.-Ing. Ernst-Dieter Hornig
and Frau Helga Vanselow and also Hong Shen, Abdelazim Ibrahim, Jinyang Fu, Phan Luu
Minh Phuang, Irina, Anastasia Alferova, Ahmed Mohamed, and Sebastian Althoff.
Special thanks is addressed to some generous people who permit him using their data in part
of this thesis namely Suzanne van Eekelen, Prof. Marcio S.S. Almeida, Abdullah C. Hassandi
and Cao Wei-ping.
The financial support of the Overseas Scholarship Program of Master/PhD for Academic
Staff of Tanjungpura University with contract No. 1725.35/D4.4/2009 during the author’s
studies at the Technische Universität Bergakademie Freiberg is gratefully acknowledged. The
support by the GraFA (Graduierten-und Forschungsakademie) is also gratefully
acknowledged. Likewise, join working with DFG-Research Project is really helpful. The
author honored to thank to Prof. Dr. Chairil Effendy, the former Rector of Tanjungpura
University and Prof. Dr. Thamrin Usman as the Rector of Tanjungpura University.
Thanks warmly to the author’s Indonesian friends in Freiberg who made his years enjoyable:
families of Arief Wijaya, Herry Permana, Toni and Eki, Nazaruddin, Heru, Annisa, Anesia,
Miranti, Annas, Anto, and Linda. Also, special thanks goes to Adang Budiman that improves
his English and Prof. Zhuang for correcting the final manuscript.
The author is very grateful to his parents, special thanks to the author’s parents, Suhadi and
Sutrinah and also mother-in-law, Nurjanah who always pray for the success, happiness and
moral support, and encouragement.
Finally, the author is overwhelmed in thanking his wife, Juliana and his children, Basith,
Retnoayu, and Satrio, for their understanding and patience throughout the period of his study.
Also, special thanks goes to Fajar Indaharti, from her he learn deep patience and sincerity.
ii
ABSTRACT
The increasing need for regional development has led engineers to find safe ways to construct
the infrastructure of transportation on soft soils. Soft soil is not able to sustain external loads
without having large deformations. The geotechnical properties of soft soil which is known
for its low bearing capacity, high water content, high compressibility and long term
settlement as well.
In pavement engineering, either highway or runway as an infrastructure, a pavement
encompasses three important parts namely traffic load, pavement and subgrade. Traffic load
generated from tire pressure of vehicle and/or airplane wheels are usually around 550 kPa
even more on the surface of the pavement. Pavement generally comprises granular materials
with unbounded or bounded materials located between traffic load and subgrade, distributing
the load to surface of subgrade.
One of the promising soil improvement techniques is a piled embankment. When
geosynthetics layer is unrolled over piles, it is known as geosynthetics supported piled
embankment. Particularly in deep soft soil, when piles do not reach a hard stratum due to
large thickness of the soft soil, the construction is an embankment on floating piles.
Furthermore, because of different stiffness between piles and subsoil, soil arching effect
would be developed there.
By using Finite Element analysis, some findings resulted from experimental works and
several field tests around the world as field case studies are verified. Some important findings
are as follows: the stress concentration ratio is not a single value, but it would be changed
depending on the height of embankment, consolidation process of subsoil, surcharge of traffic
load, and tensile modulus of geosynthetics as well. Ratio height of embankment to clear piles
spacing (h/s) around 1.4 can be used as a critical value to distinguish between low
embankment and high embankment. When geosynthetics is applied to reinforce a
pavement/embankment, the vertical distance of geosynthetics layers and number of
geosynthetics layers depend on the quality of pavement material. The lower layer of
geosynthetics withstands a tensile stress higher than upper layer. Primary reinforcements for
geosynthetics in piled embankments are located at span between piles with maximum strains
at zones of adjacent piles. Traffic load that passes through on the surface of the pavement can
reduce the soil arching, but it can be restored during the off peak hours. Settlements of
embankments on floating piles can accurately be modelled using the consolidation
calculation type, whereas the end-bearing piles may be used the plastic calculation type.
Longer piles can be effectively applied to reduce a creep. By applying length of floating piles
more than 20% of soft soil depth, it would have a significant impact to reduce a creep on a
deep soft soil.
iii
TABLE OF CONTENTS
ACKNOWLEDGMENTS ii
ABSTRACT iii
TABLE OF CONTENTS iv
LIST OF FIGURES ix
LIST OF TABLES xiii
LIST OF NOTATIONS xv
CHAPTER 1. INTRODUCTION
1.1. Reinforced Base Course 2
1.2. Piles Supported Embankment 2
1.3. Motivation 3
1.4. Geological Indonesia 5
1.5. Floating Foundation 7
1.5.1. Wooden Mattress Foundation 7
1.5.2. Wooden Piled Foundation 8
1.5.3. Cementitious Piled Foundation 8
1.6. Research Aims 9
1.7. Layout of Thesis 9
iv
2.3.2. Stress Distribution on Unbounded Material 19
2.3.3. Stress Distribution on Bounded Material 20
2.3.4. Magnitude of Loading 20
2.4. Geosynthetics Reinforcement 21
2.4.1. Type of Geosynthetics 21
2.4.2. Geosynthetics Reinforcement Mechanism 21
2.4.2.1. Restraint and Confinement 22
2.4.2.2. Membrane Mechanism 22
2.4.2.3. Local Reinforcement 22
2.4.3. Geosynthetics Reinforcement Method 22
2.4.3.1. Giroud-Noiray Method (1981) 23
2.4.3.2. Giroud-Han Method (2004) 23
2.4.3.3. US Army Corps of Engineer Method (2003) 24
2.4.3.4. DuPont Method (2008) 24
2.5. Piled Embankments 26
2.5.1. Conventional Piled Embankment 26
2.5.2. Geosynthetics-Reinforced Piled Embankment 27
2.5.3. Soil Arching Concept 28
2.5.3.1. Rectangular Prism: Terzaghi (1943) 28
2.5.3.2. Rectangular Pyramid: Guido (1987) 29
2.5.3.3. Semicircular Arch: Hewlett & Randolph (1988) 30
2.5.3.4. Positive Projecting Subsurface Conduits: BS8006-1 (2010) 31
2.5.3.5. Multi Vaulted-Dome: German Standard (EBGEO 2010) 32
2.5.3.6. Arching Evolution 34
2.5.4. Definition of Terms 35
2.5.4.1. Stress Concentration Ratio 35
2.5.4.2. Column Stress Ratio 36
2.5.4.3. Efficacy 36
2.5.4.4. Stress Reduction Ratio 36
2.5.5. Load Transfer Mechanism 37
2.5.5.1. BS 8006 Method 38
2.5.5.2. Adapted Guido Method 38
2.5.5.3. Adapted Nordic Method 39
2.5.5.4. Adapted Terzaghi Method 39
2.5.5.5. Hewlett and Randolph Method 40
2.5.5.6. German Method (EBGEO 2010) 41
2.5.6. Column Design 41
2.5.7. Tension in Geosynthetics Reinforcement due to Vertical Stress 41
2.5.7.1. BS 8006 Method 43
2.5.7.2. EBGEO (2010) Method 43
2.5.8. Soil Resistance 44
2.5.9. Tension in Geosynthetics Reinforcement due to Lateral Sliding 44
2.5.10. Settlement Analysis 45
2.5.10.1. Different Settlement on the Surface of Embankment 45
2.5.10.2. Settlement on the Bottom of Embankment 46
2.5.10.2.1. Interaction Factor Method 46
2.5.10.2.2. Equivalent Raft Method 48
2.5.10.2.3. Equivalent Pier Method 48
2.5.10.2.4. Piled Raft Method 49
2.5.10.2.5. Japanese Method 49
v
2.5.10.2.6. Public Works Center Research Method 50
2.5.10.2.7. Scandinavian Method 51
2.5.10.3. Relative Settlement Reduction (RSR) 54
REFERENCES
viii
LIST OF FIGURES
xii
LIST OF TABLES
xiv
LIST OF NOTATIONS
xvi
Fs,G+Q,k vertical load at static and variable-dynamic loading
Ft tensile strength for parallel fiber
Fv shear strength
g plastic potential function
G shear modulus of material
G specific gravity
G specific gravity of wooden material
G0 initial shear modulus
G0ref elastic small strain shear modulus at reference pressure pref
Gmax maximum shear modulus of material
Go initial shear modulus of material
h height of embankment
H thickness of soil above the point
h required base course thickness
H thickness of soil above the point
h thickness of the first layer
H thickness of fill granular material
h thickness of first layer
H thickness of soil layer
h/s ratio of embankment height to clear spacing
Hc critical height of embankments without ground treatment
hg height of arching
hg = h for h < s / 2
hg = s / 2 for h > s / 2
h'o thickness of unreinforced aggregate base
Hr thickness of reinforced granular material
Hu thickness of unreinforced granular material
I moment inertia of beam
Io Bessel function of the first kind and order of zero
Ip placticity index
J geogrids aperture stability modulus
Jk tensile stiffness of geosynthetics
K earth pressure coefficient
k modulus of subgrade reaction
k foundation stiffness
k reaction modulus of subgrade
k stiffness of geosynthetics
k1, k2,k3 regression coefficients obtained from regression analysis
Ka active earth pressure coefficient (Ka = tan2 (45+/2)
Kagh active lateral earth pressure coefficient according to DIN 4085
Kb empirical bearing capacity factor
ki permeability
Kkrit critical ratio of major stress = tan2 [45o + 'k / 2]
KoNC Ko-value for normal consolidation
Kp coefficient of passive earth pressure
kv permeability
Kw coefficient instead of K
L length of span
xvii
L length of the column
L pile length
L column length
ℒ fourth order constitutive tensor
l radius of relative stiffness
L/d to optimum column aspect ratio
LL liquid limit
Lw clear spacing
m power for stress-level dependency of stiffness
m dimensionless parameter
M slope of critical state line
m mass per unit length of the beam
ME compressibility modulus
MV rigidity modulus
Mr resilient modulus
Msoil oedometer compression modulus of surrounding soil
mv coefficient of compressibility of untreated soft clay
M fourth-order tangent stiffness tensor of the material
n number of piles
N second order constitutive tensor
N number of axle passes
N number of blows from Standard Penetration Test (SPT)
N number of cycles
N number of passes of axle load
n porosity
n stress concentration ratio
Nc factor of bearing capacity
OC organic content
OCR state of preconsolidation stress
P constant concentrated vertical with moving load
P wheel load
p total applied pressure of the embankment
p mean principle stress
p support pressure from roof above underground structure
p tire inflation pressure
p* normalized loading (p/po)
PA axle load
Pa reference pressure (atmospheric pressure)
Pav average load on a pile within group
Paxle axle load
pc vertical stress on the column
peq constant on the ellipse in the p-q plane
pG,k static loading
pG+Q,k static and variable loading
Pi actual axle load
PI plasticity index
Po equivalent standard axle load of 80 kN
xviii
po nominal overburden total stress at the elevation of underground
po pressure on the circular footing
pz vertical pressure on the circular footing
q distributed surcharge load
q surcharge acting at the surface of the soil
q ultimate vertical stress in a stone column
q vertical stress at the subgrade surface
Q design dynamic load
qa quantity
qc cone tip resistance
qf ultimate deviatoric stress
Qr required design vertical load
qu unconfined compressive strength
qu unconfined compressive strength
r column radius
r horizontal distance from centerline
R radius of equivalent tire contact area
R aspect ratio of the group
R radius of circular area
R,r,,mR,mT five addtional parameters for iIntergranular strain concept
Rf failure ratio
RSR relative settlement reduction
s diagonal pile spacing
s allowable rutting depth
s column spacing
s pile spacing
s rut depth
S settlement of embankment supported by piles
S stiffness coefficient
s/d ratio of pile spacing to pile cap diameter
S0 embankment settlement constructed on soft soil without the support of piles
S1 compression of the stabilized volume for case 1
S1 settlement of a single pile under unit load
S2 compression of stabilized volume for case 2
s-a width of strip
Sc settlement of the columns
So settlement of untreated soil subjected to the actual load of embankment
Sp primary settlement
Sr degree of saturation
SRR stress reduction ratio
Ss settlement of untreated soil subjected to reduced pressure
Ss secondary settlement
su undrained shear strength
t time
T corrected thickness of base course
t1 time in first step
t2 time in next step
Teff effective aggregate thickness
xix
To initial aggregate thickness
TRP tension in geosynthetics
Jaumann stress rate
V constant velocity
v Poisson's ratio
vcr critical velocity
vur un-/reloading Poisson’s ratio
Vs shear wave
w lateral deflection of plate strip
w water content
w deformation at interface between two layers
w(r) vertical displacement of the surface
wi settlement of a pile i within a group of n piles
Wn water content
WT equivalent vertical load
z thickness of the soil overlying the element
ϵacc v volumetric strain
dstab consolidation settlement of stabilized ground
d consolidation settlement of unstabilized ground
P different stress
Sr differential settlement between the columns and the untreated soil
state ’compaction’ = n /n0
ij interaction factor for pile i due to any other pile j within the group
influencing factor due to the presence of geosynthetics layer
settlement reduction ratio
maximum sag (vertical deflection) of geosynthetics
settlement of roof
intergranular strain tensor
* normalized displacement ( / B)
acc(N) average accumulation curve
ampl strain amplitude
elastic strain rate
accumulation of deviatoric strain
accumulation of volumetric strain
e elastic strain
p plastic strain
plastic strain rate
strain
v vertical strain
k strain in geosynthetics
’ effective internal friction angle
c, hs, n, ed0, ec0, ei0, , eight material parameters for hypoplastic model
max peak secant angle of shearing resistance
internal friction angle of the soil
xx
unit weight of the soil
partial safety coefficient for the consideration of possibility
k unit weight of embankment
unit weight of the soil
ampl shear strain amplitude
H nominal vertical stress at the base of embankment
p plastic shear strain
unit weight
av average stress ratio
creep factor
installation damage factor
biological and chemical degradation factor
* modified swelling index
plastic multiplied factor
* modified compression index
Poisson's ratio
* modified creep index
1, 2 Poisson’s ratio of first and second layer
density of soil
average stress of applied embankment plus surcharge at level of top piles
a allowable strength
3 confining stress
c stress on the columns
c vertical stress on the pile cap
d deviatoric stress
creep strength of column
h' effective horizontal pressure on the columns
h horizontal stress
i radial (vertical) stress acting immediately beneath the crown of the arch
pile vertical stress on piles
rz vertical stress at interface between two layers with radius = r
s reduced pressure on the untreated soil due to embankment
s stress on the subsoil
soil vertical stress on soil
u ultimate strength
v vertical stress
z vertical stress at interface between two layers with radius = 0
flexural flexural strength of wooden material
'ro free-field lateral effective stress
compressive compressive strength of wooden material
comp or tensile tensile strength of wooden material
θ bulk stress (σ1 + σ2 + σ3 )
shear strength
oct octahedral shear stress
u,col undrained shear strength of the columns
xxi
shear strength of wooden material
(x0 , z) fixed coordinates along the length of strip and in the vertical direction
xxii
CHAPTER 1
Introduction
The increasing need for infrastructure development has led engineers to find an alternative on
soft soils. Engineers are continually faced with maintaining and developing pavement
infrastructure with limited financial resources. Traditional pavement design and construction
practices require high-quality material for fulfilment of construction standards. In many areas
of the world, quality materials are unavailable or in short supply. Due to these restrictions,
engineers seek alternative designs or innovative design practices.
Weak foundation soils have always been a challenge to geotechnical engineers when
designing infrastructure. Low bearing capacity, slope stability, lateral pressure, movements
and differential settlements are some major concerns. However, a variety of techniques are
available to address these issues, including preloading, deep mixing columns, stone columns,
use of lightweight fill and soil replacement. Flexible pavement layered system is shown in
Fig.1-1. Meanwhile, vertical load distribution induced by a moving wheel load is illustrated
in Fig. 1-2.
(a) (b)
Fig. 1-2 Stress distribution with depth in a flexible pavement, (a) High stress area directly
under wheel load (b) Reduced load at subgrade level
1
1.1. Reinforced Base Course
Fig. 1-3 shows a typical section of reinforced base course, which consists of an aggregate
base layer, a subgrade and a reinforcement layer usually placed between the base course and
subgrade. The base course and geogrid transmit the traffic load to the top of the subgrade.
Under repeated load, the behaviour of the base-geogrid-subgrade system is complicated. The
overall behaviour depends on the properties of geosynthetics, soil characteristics, and the
interaction between the soil and the reinforcement.
(a) (b)
Fig. 1-3 Relative load magnitude at subgrade layer level (a) Unreinforced flexible pavement
(b) Geosynthetics-reinforced flexible pavement
2
Settlements and differential settlements are also significantly reduced once the technique is
used successfully.
Geosynthetics-reinforced piles supported embankment is a promising technique when it is
constructed over weak soil as depicted in Fig. 1-4. It is also valuable in two ways; if there is
adjacent structure that is sensitive to differential settlement and if there is a need for
accelerated construction. This complex foundation system consists of strong piles in soft soil
and compressible soil. A bridging layer is placed above the columns consisting of coarse
grained soil and one or more layers of geosynthetics is situated on the bridging layer.
(a) (b)
Fig. 1-4 Geosynthetics-reinforced piles supported embankment (after Satibi, 2009)
1.3. Motivation
Stiff soils are preferable as a foundation for buildings, roadways and other infrastructures.
However, these types of foundation soil are not always available on site. In tropic regions like
Indonesia, soft soil deposits (peat, organic and inorganic soft soil) occur in lowland and
highland areas and generally termed as basin and valley of peat. In some places, the depth of
soft soil can reache 30 meters or even more.
Soft soil exhibits very low bearing capacity and it is not suitable for constructing
embankment, highway, runway, railway, building or any other load bearing engineering
structure. Organic soft soil in its natural state consists of water and decomposing plant
fragments with virtually no measurable bearing strength. Peat and organic soil are considered
as soft soil because they have high settlement value and even under moderately applied load.
The most extensive areas of peat soil are located in the northern hemisphere. It is estimated
about 1 billion acres of the world are peat or equal to 4.5% of total land area. In United States
peat is found in 42 states, with a total area of 30 million hectares. Canada has 170 million
hectares and the former USSR has 150 million hectares. In Japan, peat is widely distributed
throughout Hokaido with an area approximately 200 thousand hectares (Alwi, 2007).
The world tropical peat land is located in Asia, South Africa and Latin America covering
around 30 million hectares. Around two thirds (66%) are in Southeast Asia (Global
Environtment Centre). While in Indonesia peat soil covers about 20 million hectares of the
country land. Most extensive areas are located across eastern Sumatra, western and southern
Kalimantan and southern Papua. The soft areas are mostly dominated by peat overlaying soft
clay soils, whereas there are mainly soft clay soils spreading out in other places as depicted in
Fig. 1-5.
3
Kalimantan
Sumatera
Maluku
Sulawesi Papua
Java
The total area of soft soils (peat and soft clay soil) in Indonesia is about 60 million hectares,
which contribute around 30 percent of the total land area (Satibi, 2009). Distribution of peat
soils in some regions or provinces in Indonesia is reported by Radjaguguk (1991) resulted
from a study conducted by Euroconsult in 1983 and the final report by BB Litbang (2008).
This is shown in Table 1-1. Meanwhile, Table 1-2 illustrates the depth distribution of peat
soils.
Table 1-1 Distribution of peat soils in Indonesia
Area of peat soils in Indonesia (hectares)
Islands / Provinces
After Radjaguguk, 1991 After BB Litbang, 2008
Sumatera Island
Nanngroe Aceh 270 -
North Sumatera 335,000 -
West Sumatera 31,000 -
Riau 1,704,000 4,043,600
Jambi 900,000 716,839
South Sumatera 990,000 1,483,662
Bengkulu 22,000 -
Lampung 24,000 -
Kalimantan Island
West Kalimantan 4,610,000 1,729,980
Central Kalimantan 2,162,000 3,010,640
South Kalimantan 1,484,000 331,629
East Kalimantan 1,053,000 -
Sulawesi Island
Central Sulawesi 15,000 -
South Sulawesi 1,000 -
Southeast Sulawesi 18,000 -
Maluku Archipelago 20,000 -
Papua Island 4,600,000 7,001,239
Java Island 25,000 -
17,994,270 18,317,589
Source: Radjaguguk, 1999; BB Litbang, 2008
4
Table 1-2 Distribution of thickness for peat soils in Indonesia
Distribution based on thickness of peat (%)
Total area
Province Shallow Intermediate Deep
(hectares)
0-150 cm 100-200 cm >200 cm
Riau 8.6 10.7 80.7 486,339
Jambi 33.4 9.3 57.3 168,163
South Sumatera 63.0 11.5 25.5 317,784
West Kalimantan 39.5 34.6 25.9 100,754
South & Central Kalimantan 62.6 19.6 17.8 190,145
Total 1,263,185
Source: Radjaguguk, 1991
Fig. 1-6 Geological time scale showing the appearance of life forms and the occurrence of
major geological events (after MacKinnon et al., 1996)
5
As recently as of 25 million years ago, very recently on the geological time scale as in Fig. 1-
6, the Indonesian archipelago as we know it today simply did not exist, but the story began
much earlier than that.
The continental land masses are by no means permanent and the earth is in a dynamic state.
The outer solid part of the earth, the crust, is quite thin, like the rind of an orange. The are
two kinds of crust: oceanic and continental. Oceanic crust is usually young (0-200 millions
years), thin (5-15 km) and composed mostly of dense volcanic rocks. Continental crust often
has a core of older rocks (200 to 3,500 million years), is thicker (20-50 km) and is less dense
than oceanic crust, composed of rocks such as sandstones and granites. Western Indonesia,
comprised of Kalimantan, Sumatra, and West and Central Java, is composed predominantly
of continental crust, as is much of the shallow sea floor between these islands. Below the
crust of the earth is a zone where the rock is hotter and more plastic. Continental and oceanic
plates float on a fluid, underlying material (Simandjuntak, 1993).
Plate movements are very slow, only a few centimeters per year, but over a 60-million year
period (the time since the dinosaurs became extinct) a plate drifting only 1 cm per year would
have moved 600 km. The Indonesian region is dominated by three major plates, namely
Southeast Asian plate, the Indo-Australian plate and the Pacific plate, as well as several
smaller platelets that have sheared off Irian Jaya. The process is still continuing today
(Gorsel, 2012).
Soil conditions are very important in affecting the distribution of vegetation. There are five
factors in soil formation: lithology, climate, topography, biological organisms and time. In
Kalimantan, the majority of its soils have developed on the rolling plains and dissected hills
6
on sedimentary and old igneous rock. These soils range from strongly weathered and acid
ultisols to young inceptisols. In the south, extensive alluvial plains and peat soils extend into
the Java Sea. Weathering is strong in the humid tropics, favoured by both warmth and
moisture. Because of the high rainfall, soils are constantly wet and soluble constituents are
removed. This process is called leaching. High levels of weathering, leaching and biological
activity (degradation of organic matter) are characteristics of many Bornean soils (van
Gorsel, 2012).
9
Chapter 4 presents some experimental works in laboratories and field case studies for
various situations over the world in handling the soft soil. This chapter describes some
problems and simultaneously solves on the pertaining realm in soft soil such as excessive
settlement, horizontal thrust, vertical distance of geosynthetics layer, various kinds of piles,
static and cyclic loading.
Chapter 5 provides some analytical approaches of rutting on the surface of the pavement as
material response under static and moving load as well settlement on foundation. Moreover,
some constitutive models are introduced that can be applied in finite element analysis.
Moreover, geometrical idealization is needed to transform the real geometric shapes to finite
element analysis.
Chapter 6 discusses some interesting findings in Chapter 4 and then verified with numerical
analysis using Plaxis software package. Some important findings is presented such as optimal
vertical distance of geosynthetics layer in base course, influence of very soft soil against soil
arching and/or stress concentration ratio, determining critical height of embankment, various
types of load transfer platform, tensile strain of geosynthetics in piled embankment and creep
phenomenon in soft soil.
Chapter 7 elaborates excessive settlements of pavement construction supported by floating
piles over deep soft soil of two case studies in Supadio Airport Pontianak regarding with
runway reconstruction work and apron widening project.
Chapter 8 presents the main findings and concludes this study. Later on, recommendation
for further works is highlighted.
10
CHAPTER 2
Literature Review
Fig. 2-1 Atterberg limit for organic and inorganic soils (after Puslitbang, 2001)
There is also fine grained soil, in which Indonesian Geotechnical Guide-1, 2001 divides into
three groups based on their organic content, as shown in Table 2-1. Organic soil (O) is the
soil which has the organic content 25 to 75%. Organic soil is categorised into OL and OH
appropriate level of plasticity. Peat is a soil that has organic content more than 75%. Based on
its fiber content, peat soil is grouped into two: amorphous with fiber content less than 20%
and fibrous with the fiber content exceeds of 20%.
11
Table 2-1 Soil classification based on organic content
Organic content Soil group
>75% Peat
25 - 75% Organic soil
< 25% Low organic soil
Source: Indonesian Geotechnical Guide-1, 2001)
Peat as used today includes a vast range of peat, peaty organic soils, organic soils and soils
with organic content (Landva et al., 1983). According to the ASTM standard D 2487-00,
organic clay/silt with sufficient organic content will influence the soil properties. By
classification, an organic soil is a soil that would be classified as a clay/silt except that its
liquid limit value after oven drying is less than 75% of its liquid value before.
In geotechnical engineering, soft terms specify to clay with shear strength ranging from 12.5
to 25 kPa, whereas the very soft clay is below 12.5 kPa. Table 2-3 gives some clues of an
indication of the shear strength when it was identified in the field (Indonesian Geotechnical
Guide-4, 2001).
Table 2-3 Field indicator of undrained shear strength for soft clay
Consistency Field indication
Indonesian Public Works (1999) defines the soft soil as the soil that can be penetrated with
the thumb a minimum 25 mm or has the undrained shear strength less than 40 kPa on the
basis of field vane shear test. Soft soil may compose of inorganic and organic soil. Inorganic
soft soil generally consists of clay or silt and it has organic content between 0 and 25% or ash
content ranging from 75 to 100%. Whereas organic soft soil consists of clay or silt and
organic content between 25 and 75%, or ash content ranging from 25 to 75%.
12
Peat is a soil-forming consisting of the primary of the remaining plants and the organic
content exceeding 75%. Basically, all types of soil are Resen old in geological terms that are
less than 10000 years old. The geological period is also commonly known as the Holocene.
Undrained shear strength is an important parameter. This parameter for soft clay around the
surface in Indonesia ranges from 10 to 20 kPa. Undrained shear strength of 10 kPa is only
able to support an embankment about 2 m high.
2.1.7.3. Compressibility
Farrell et al. (1994) showed that the compression index Ireland peat associated with liquid
limit according to equation:
cc = k (WL – 10) (2-1)
where: k = 0.007 to 0.009
For fibrous peat the equation above cannot be applied. Consolidation tests on fibrous peat in
Barengbengkel show the value of cc up to 20. The vertical compression index is almost twice
of horizontal compression index values. Keep in mind that the value of high compression
index (cc) cannot be applied to the conventional calculation in small strain. Peat compression
index will decrease with increasing stress.
2.1.7.4. Permeability
Barry et al. (1992) performed the pumping test to support permeability in the forests of Riau
and revealed that the permeability between 10-2 and 10-4 m/s. They also compared it with
other research, as indicated in Table 2-5.
Table 2-5 Permeability values for peat soils
Discription of peat Permeability (m/s) Sources
-1
At the surface > 10 Hobbs (1986)
On the bottom 3 x 10-5 -idem-
Fen Acrotelm in Russia
near surface 3 x 10-5 -idem-
near bottom 6 x 10-7 -idem-
Peat soils in Irlandia 3 x 10-8 to 10-7 -idem-
Sphagnum peats
H8 to H10 6 x 10-8 -idem-
H3 10-5 -idem-
Sedge peat, H3 to H5 10-5 -idem-
Brushwood, H3 to H6 10-5 -idem-
Fibrous acidic Malaysia peat 2 x 10 to 6 x 10-8
-5
Toh et al. (1990)
After Barry et al., 1992 and Hobbs, 1986
15
development in this region. Some properties of peat soil both islands of Sumatera and
Kalimantan as reported Soepandji (1996, 1998) are depicted in Table 2-6.
Table 2-6 Properties of peat soils in Sumatera and Kalimantan islands
Kalimantan Sumatera
Properties
Pontianak Banjarmasin Duri Desa Tampan Musi
Ash Content (%) 1.2 3.29 21.96 25.2 50.7
Water Content (%) 632 198 235 338 235.4
Specific Gravity 1.42 1.47 1.6 1.55 1.82
Liquid Limit (%) 260 182 440 236 274
Plastic Limit (%) 196 148 377 309 194
Shrinkage Limit (%) - 28 - 59 -
pH value 4.8 6.5 3.9 3.6 3.3
Bulk Density - - 1.084 0.95 1.12
Compression Index, Cc - - 2.5-3.2 2.11 1.57
Recompression Index, Cr - - 0.07-0.13 0.107 0.05
Classification ,ASTM low ash, low ash, organic organic soil organic
D4427-92 (1997) acidic acidic soil soil
After Soepandji et al., 1996, 1998 cited Eka Priadi, 2008
16
Paxle
p = Contact Pressure
H Pavement or Embankment
Subgrade
Fig. 2-2 Contact area and tire pressure (after TENAX, 2001)
For single or dual wheel, contact pressure (p) is equal to their tire inflation pressure. The
contact pressure is assumed to be a circular area and its radius is calculated by:
(2-6)
(2-7)
(2-9)
It is similarly that the required thickness for a section reinforced using geosynthetics is:
(2-10)
Magnitudes of bearing capacity of the subgrade soft soil both unreinforced and reinforced
are different between methods each other as shown in Table 2-7.
17
Table 2-7 Comparison of bearing capacity in different design methods
Unreinforced Reinforced
Method Material
section section
Giroud-Noiray (1981) 3.14 cu 5.14 cu Geotextile
Giroud-Han (2004) 3.14 cu 5.14 cu Geotextile
5.71 cu Geogrid
USACE (2003) 2.8 cu 3.6 cu Geotextile
5.8 cu Geogrid
Barenberg (1975) 3.0 cu 6.0 cu Geotextile
DuPont Typar (2010) 3.14 cu 5.14 cu Geotextile
Philips (1987) 2.8 cu 5.0 cu Geotextile
Rodin (1965) 3.1 cu 6.2 cu Geogrid
Roadex III (2008) 4.0 cu - -
(2-11)
(2-12)
(2-13)
Secondary settlement:
(2-14)
18
Sources of dynamic actions are described in Fig. 2-3.
Dynamic Actions
The actions in terms of DIN 1054 can be differentiated into dynamic, cyclic, and shock-like
actions. Dynamic actions refer to high frequency. Inertia forces are not negligible and it can
critically influence system behavior. Cyclic actions refer to low frequency actions where the
inertia forces can generally be ignored (frequencies ≤1 to 2 Hz). The shock-like actions refer
to actions acting over a short period only. The time may be in the range of milliseconds up to
several second. Their upper bound is not fixed and inertia forces may also act.
Additional distinguishing criteria include load-time history characteristic, effective spatial
direction, source and frequency of occurrence. Load time history is shown in Fig. 2-4.
19
Fig. 2-5 Load distribution on granular material (after Giroud et al., 1981)
(2-15)
Standard axle load is different for each country. In Greece, Beskou et al. (2011) reports that
for a locomotive (or engine) is around 210 kN and 150 kN for carriage. In Indonesia, railway
infrastructure is subjected to a maximum axle load of 180 kN (Kepmenhub, 2000).
20
2.4. Geosynthetics Reinforcement
2.4.1. Type of Geosynthetics
Geosynthetics for earth reinforcement purpose is divided into two groups: permeable matter
(geotextile, geogrid, geocomposite and geolinier) and impermeable matter (geomembrane
and intermembrane). For ground improvement or earth reinforcement, it must have some
properties such as shear strength, tensile strength, and punch strength (Suryolelono, 1997).
Types of geosynthetics used in earth reinforcement are shown in Table 2-8.
Table 2-8 Type of geosynthetics used in reinforced earth (after Suryolelono, 1997)
Material Type Form
Geotextile Woven and Unwoven
Geogrid Uniaxial and Biaxial
Geosynthetics
Geolinier Bar
Geostrip Strip/tape
Geocell Assembled cell
1. Restraint + Confinement
2. Membrane mechanism
3. Local reinforcement
Fig. 2-7 Three stabilization mechanisms (after Perkins and Ismeik, 1997)
More attention has been given to the important practical application of geogrid reinforcement
incorporated at the base of a layer granular fill placed on a soft clay subgrade. This kind of
construction is commonly used for low-cost unpaved roads such as temporary site access
roads, low embankments, car parks and the working platform. The purpose of the fill is to
provide a suitable operating surface on which concentrated loads may be carried without the
subgrade failing or deforming excessively. It is now common practice to use a layer of
polymer geosynthetics at the base of the fill layer in order to separate the fill from the soft
soil beneath and to improve its load-carrying capacity.
The behaviour of such a system is complex and a number of procedures for the design have
been proposed, notably by Barenberg et al. (1975), Giroud and Noiray (1981), Sellmeijer et
al. (1982), Giroud et al. (1984). These procedures which are based on simplified deformation
22
mechanisms and some empiricism have provided the basis for satisfactory design. Some
methods work better than others for certain site conditions, material and traffic volume.
However, none is principles of plasticity and fully addresses the significance of shear forces
acting between the base and the subgrade. In 2004 Giroud and J. Han introduced a method
and improve the previous method. Also US Army Corps of Engineer in 2003 launched a
method particularly for low volume.
Where: h'o is thickness of unreinforced aggregate base (m); N is number of passes of axle
load; PA is axle load (Newton); s is rut depth (m) and cu is undrained cohesion of subgrade
soil (Pascal).
(2-18)
23
Nc = 5.71 and J=0.32 m-N/degree for Tensar BX 1100 reinforced base course
Nc = 5.71 and J=0.65 m-N/degree for Tensar BX 1200 reinforced base course
Fig. 2-8 Correlation between Cone index, CBR, and Shear strength (after Archer, 2008)
Next step is to determine the subgrade bearing capacity using equation below:
Subgrade bearing capacity = C. Nc (psi) (2-19)
Once the subgrade bearing capacity has been determined, the engineer or designer can refer
to one of the three relevant design charts (single wheel, dual wheel and tandem wheel).
Soil
CBR
2
cu [kN/m ]
cu [psi]
2
ME [MN/m ]
2
MV [MN/m ]
Fig. 2-9 Correlation chart for estimating the subgrade value (after Barenberg, 1975)
An unpaved road normally consists of unbound aggregate base. Inclusion geotextile between
subgrade and aggregate base allows for better aggregate compaction, subgrade consolidation,
reinforcement of the structure and to increase the ultimate bearing capacity of subgrade
around (2+) cu.
The first procedure is to determine initial aggregate thickness (TO) according to load and
subgrade conditions and then consider service life and aggregate efficiency. Fig. 2-10 below
describes this stage easily. The left side of the chart is the subgrade CBR and axle load Pi to
determine TO or alternatively using right one.
(a) (b)
Fig. 2-10 Design chart DuPont method (a) Compacted crushed stone thickness for 1000 axle
loads, (b) Factor to determine curve Pi
The second step is to make an adjustment of TO for service life or corrected thickness of base
course (T) using equation (2-20).
T = C . TO = [ 0.27 log ( Ni . ESAL) + 0.19 ] . TO (2-20)
ESAL = ( Pi / Po )3.95 (2-21)
The service life is expressed as the total number of 80 kN axle load application. The actual
axle load (Pi) is first converted to an equivalent standard axle load (Po=80 kN).
The last step is to adjust the use of different kind of material. This difference is accounted for
by using the aggregate efficiency . Angular crushed aggregate is the best because it
25
interlocks well and provides a high bearing capacity. Depending on availability, other
materials or blends can be used and Table 2-8 indicates typical thickness efficiency factors of
various surfacing and base materials.
Table 2-8 Aggregate efficiency
Material Efficiency
Paving Stone 2
Hot Mix / Dense-Macadam 2
Dense Surface Course 2
Soil-cement > 5 MPa compression 1.5
Soil bitumen 1.5
Hard crushed stone aggregate-‘standard’ 1.0
Medium crushed stone aggregate (CBR > 80%) 0.8
Hard round stone aggregate (CBR > 80%) 0.8
Medium round stone aggregate 0.5
Sandy gravel (CBR = 20-30%) 0.5
Crushed limestone 0.5
Loose gravel, compactable sand 0.4
Source: DuPont method, 2008
(a) (b)
Fig. 2-11 Typical piled embankments (a) Conventional piled embankment (b) Piled
embankment with basally reinforcement
26
The conventional pile-supported system or piled embankments without reinforcement
requires large pile caps and very close piles. This relies on soil arching in the embankment
fill and can lead to punching failure. It is essential thing to transfer the large embankment
loads to the piles and to avoid surface deformation due to large differential settlement
between caps. Consideration required to ensure the edges piles can take the lateral loads
imposed by lateral thrust, whereas piled embankment with reinforcement the reinforcement
take places the lateral embankment thrust loads so raking piles not required anymore.
A survey of various projects (Han, 1999) found that in conventional piled embankments the
coverage ratio (the ratio area of pile caps to total foundation area) is 60-70%, whereas in
geosynthetics reinforced piled supported embankments is reduced to about 10-20%.
27
2.5.3. Soil Arching Concept
Arching is defined by Mc Nulty (1965) as ’’the ability of a material to transfer loads from
one location to another in response to a relative displacement between the locations. A
system of shear stress is the mechanism by which the loads are transferred ’’. Consider soil
on a rigid base, there is no tendency for differential movement and hence no soil arching as in
Figure 2-12 (a) the stress acting at point of a is the overburden stress (H), where is the unit
weight of the soil and H is the height of the soil prism. When one of the local supports at the
point a is removed, the point a is in tension and a roof tension arch is formed. The true arch
collapses as the soil is not in equilibrium as depicted in Figure 2-12(b) meanwhile in the next
stage as in Figure 2-12(c) the soil settles in an inverted arch, the adjacent soil develops the
required shear strength and the soil reaches equilibrium state. The transfer of pressure from
the yielding portion to the stationary portion is so-called ’arching’.
The corresponding normal stress (horizontal stress) on the vertical surface of sliding ( h) is
given by:
h = K. v (2-24)
where: h = the horizontal stress (kN/m2)
K = the earth pressure coefficient (dimensionless parameter)
The shear strength of the soil (assuming the soil to be cohesionless) is determined by:
= h . tan (2-25)
where: = the shear strength (kN/m2)
= the internal friction angle of the soil (degree)
When the element is in equilibrium, the summation of the vertical forces must be equal to
zero. Therefore, the vertical equilibrium can be expressed as:
dv/dz = – K. v tan / B (2-26)
where: 2B = the width of strip (m)
z = the thickness of the soil overlying the element (m)
Using the boundary condition that v = 0 for z = 0, the partial differential equation can be
solved as follows (Terzaghi, 1943 and later McKelvey, 1994).
(2-27)
(2-28)
Handy (1985) proposed that the shape of the arched soil is catenary and suggested the use of
the coefficient Kw instead of K, by considering an arch of minor principal stress.
Kw = 1.06 (cos2 + Ka sin2 ) (2-29)
where: = 45 + /2
Russell et al. (2003) proposed that K could be conservatively taken as 0.5 and Potts &
Zdravkovic (2008) proposed that K = 1.0 gave good correspondence with the results of plane
strain finite element.
0.5(s-a) a
45o
The stress on the subsoil (s) resulted from self weight of the unsupported soil mass is equal
to the volume of triangle/pyramid multiplied by the soil unit weight (). Then, it is divided by
the area over which the soil prism acts. For the two dimensional situation, the stress acting on
the subsoil is:
s = (s-a) / 4 (2-30)
Meanwhile, for the three dimensional situation, the equation is modified to:
s = (s-a) / 3 √2 (2-31)
where: s-a = the width of strip (m)
From the equations above that the height of the embankment has no effect on the pressure
acting on the subsoil. Moreover, the friction angle of the fill material is not considered in this
case. Love&Milligan (2003) suggest that the Guido method may experience difficulties when
dealing with situations where support of the existing subsoil is very low. The Guido method
concentrates more on reinforcement rather than arching process.
a s
30
The 'arches of sand' transmit the majority of embankment load onto the pile caps, with the
subsoil carrying load predominantly from the 'infill' material below the arches. The arches are
assumed to be semi circular (in 2D) uniformly without overlapping.
The analysis considers equilibrium of an element at the 'crown' as in Fig. 2-16. Here the
tangential (horizontal) direction is the direction of major principal stress and the radial
(vertical) direction is the minor principal stress. Yielding is in the 'passive' condition since the
horizontal stress is the major principal stress.
a 0.5s
By using the boundary condition, considering vertical equilibrium, that the stress at the top of
arching layer is equal to the weight of material above acting on the outer radius of arch will
give a solution for the radial (vertical) stress acting immediately beneath the crown of the
arch ( i ). The vertical stress acting on the subsoil is then obtained by adding the stress due to
the infilling material beneath the arch, based on the maximum height of infill (s-a)/2.
s = i + (s-a) / 2 (2-32)
When considering the three dimensional solution the equation above is modified to:
s = i + (s-a) / √2 (2-33)
The vertical stress (s) is considered uniform here. Though, Low et al. (1994) introduced this
parameter to allow a possible non uniform vertical stress on the soft ground.
In case at the pile cap, the tangential (vertical) stress is the major principal stress and the
radial (horizontal) stress is the minor principal stress. In conjunction with overall vertical
equilibrium of the embankment, a value of s is obtained in the limit when the ratio of the
major and minor principal stress is Kp. In fact, yielding occurs in an ’active’ condition, since
the vertical stress is the major principal stress.
31
The equation proposed by Marston was derived from field tests at the Engineering
Experiment Station at Iowa State College in 1913. The equation is normally used to calculate
the reduced loads on buried pipes (see Fig. 2-17).
Fig. 2-17 Stress distribution of sub-surface conduit (after Marston et al., 1913)
(2-35)
32
Fig. 2-18 Theoretical arching model (after Zaeske et al., 2001)
EBGEO (2010) recommends the use of geosynthetics reinforcement but the arching effect
and the membrane tension are dissociated. After arching on embankments, it will develop a
redistribution of vertical stress at the surface between piles (zo,k) and vertical stress at
surface of piles ( zs,k) as depicted in Fig. 2-19.
Soft soil
Vertical stress at the surface between piles after arching on embankments is divided into two
groups, namely the vertical stress for static loading (zo,G,k) and vertical stress for static and
variable-dynamic loading ( zo,G+Q,k ).
1+ 2 . 2− (2-36)
1+ 2 . 2− (2-37)
Furthermore, the vertical stress at the surface of piles after arching on embankments resulted
from surcharge divided into two groups, namely the vertical stress for static loading (zs,G,k)
and vertical stress for static and variable-dynamic loading ( zs,G+Q,k ).
zs,G,k = [ ( k . h + pG,k ) – zo,G,k ] AE/As + zo,G,k (2-38)
zs,G+Q,k = [ ( k . h + pG+Q,k ) – zo,G+Q,k ] AE/As + zo,G+Q,k (2-39)
Where AE is the area of one cell pile embankment and As is a support surface for point and/or
linear bearing elements.
Vertical loads or forces on bearing element as:
Fs,G,k = zs,G,k . As (2-40)
Fs,G+Q,k = zs,G+Q,k . As (2-41)
Generally, the resultant force on bearing element conservatively is calculated using:
Fs,G,k = ( k . h + pG,k ). AE (2-42)
Fs,G+Q,k = ( k . h + pG+Q,k ). AE (2-43)
(2-45)
where : p = the support pressure from roof above underground structure (kN/m2)
po = the nominal overburden total stress at the elevation of underground (kN/m2)
= the settlement of roof (kN/m2)
B = the width of underground structure (m)
As shown in Fig. 2-21 that the GRC is divided into four parts namely the initial arching
phase, maximum arching, loading of the recovery stage, and the ultimate state.
Fig. 2-21 Generalized ground reaction curve (after Iglesia et al., 1999)
2.5.4.3. Efficacy
Hewlett &Randolph (1988); Low et al. (1994) defined efficacy, E, as the proportion of the
embankment weight carried by the piles rather than the subsoil. Efficacy increase (tending
towards 1.0) as the effect of arching increases.
(2-48)
(2-49)
(2-52)
(2-53)
(2-54)
36
foundation soil between the columns and hence the values of n, SRR, and CSR are equal to
one. If complete soil arching hypothetically develops, the entire applied embankment load
would be carried by the columns and no load would be carried by the soil between columns.
Relationship of the parameters is shown in Fig 2-22.
The design of the load transfer platform (Collin, 2004) based on the use of multiple layers is a
refinement of a method referred to Guido method. The primary assumptions for the beam
theory are:
1. a minimum of three layers of reinforced is applied to create the platform
2. spacing between layers is 200-450 mm
3. platform thickness is greater than or equal to one half the clear spacing piles
4. soil arching is fully developed within the depth of the platform
When one or more layers of geosynthetics reinforcement are placed in the fill above the
columns, the stress that would otherwise be applied to the foundation soil between columns is
37
assumed to be carried by geosynthetics. Thus, the SRR value can be used to represent the
portion of embankment load carried by geosynthetics reinforcement.
BS8006 defines a critical embankment height equal to 1.4(s-a). If the embankment height is
below the critical height, arching is not fully developed. The embankment and surcharge load
is converted to an equivalent vertical load, WT, between the piles which triangle load with WT
as maximum value at the centre of span, as determined by using equation:
(2-55)
If the embankment heights is greater than the critical height assumed that all load above the
critical height are transferred directly to columns as a result of arching in the embankment fill
(Kempton et al.1998).
(2-56)
(2-57)
(2-58)
HCritic
Fig. 2-25 Triangular soil arching of Nordic method (after Nordic guideline,2003)
The stress reduction ratio for the two dimensional approach, SRR2D, is given as:
(2-59)
Rogberck et al. (1998) provide a correction factor used to compute the geosynthetics tension
for an embankment supported by square piles caps. The net effect of this correction factor for
the three-dimensional stress reduction is the same as the two dimensional value and when a
surcharge load considered, the equation above is:
(2-60)
Equation 2-59 is applicable for embankment height greater than or equal to the critical height
and for lower embankments. The upper part of the triangle must be truncated to calculate the
stress reduction ratio.
(2-61)
where: K = coefficient of lateral earth pressure, and = internal friction angle of embankment
fill. The K value is assumed to be equal to one and the equation above is referred as Adapted
Terzaghi's method 1.
39
Several years later Russell et al. (2003) presented a modified version of the previous method.
They assume that the portion of the embankment fill that settles as a cruciform has a height of
n.H and the embankment fill above the settling cruciform is treated as surcharge. The stress
reduction ratio from latter method (Adapted Terzaghi Method 2) is expressed as:
(2-62)
(2-63)
Efficacy at the top of pile is:
(2-64)
(2-65)
The stress reduction ratio may be determined efficacy and its value at the arch crown is:
(2-66)
40
The stress reduction ratio at the top of the pile is:
(2-67)
1+ 2. 2 (2-68)
Kempfert et al. (2004) recommend using the equation to determine the load on the top of
geosynthetics reinforcement. Soil support from the subgrade beneath geosynthetics layer also
can be taken into account. The stress reduction ratio may be determined as follows:
(2-69)
The required design vertical load, Qr, in the column is determined according to the following
equation:
(2-70)
41
2.5.7. Tension in Geosynthetics Reinforcement due to Vertical Stress
Geosynthetics reinforcement is commonly used in soil by placing at the base of the
embankment. The tension will provide support between the pile caps (see Fig.2-28).
Geosynthetics reinforcement helps to transfer the weight of the embankment directly on to
the columns or piles (Lawson, 1992). At the edges of embankment it also prevents lateral
spreading (Hewlett&Rundolph, 1988). However, these two functions are normally considered
independently.
The effect of additional capacity to carry a vertical load could be added based on purely
tensile response. They proposed that assuming the geosynthetics was subjected to a uniform
vertical load and deforms as a parabola.
(2-71)
(2-72)
As can be seen in equation (2-72) that the strain, , increases as the square of . The tension
in the geosynthetics is assumed as linear response against strain as below:
(2-73)
where: k = the stiffness of geosynthetics and in another expression:
(2-74)
This can be re-arranged to express how the load which can be carried theoretically increase
with the sag:
(2-75)
Giroud, J.P. (1995) presented a correlation between deflection and strain in geosynthetics, it
is similar to Eq. 2-71, and re-writing it as follows:
42
(2-76)
(2-77)
The tension in the reinforcement is calculated taking into consideration the maximum
allowable strain in the reinforcement. Six percent of strain is considered the upper limit for
transferring the load to the piles. The upper limit should be reduced for shallow embankments
to prevent differential movement on the surface of the embankment. To avoid long term
localized deformations at the surface of the embankment, the long term strain should be kept
to a minimum and a maximum creep strain of 2% is permitted for permanent construction.
Fig. 2-29 Plan view of foundation and geosynthetics (after Kempfert et al., 2004)
43
The maximum strain in the geosynthetic reinforcement is dependent upon the tensile stiffness
of geosynthetic, Jk, the modulus of subgrade reaction, k, the vertical load, Fk, clear spacing,
Lw, and width of column, bErs . The value of strain in geosynthetics, k, can be determined
using the design chart in Fig. 2-30.
Fig. 2-30 Chart to determine strain in geosynthetics (after Kempfert et al., 2004)
The tensile force in the reinforcement can be calculated using the following equation:
TRP = k . Jk (2-78)
When applied two layers of geosynthetics, the calculated tensile force is divided with respect
to the ratio of their tensile moduli (Kempfert et al. 2004).
44
= the unit weight of embankment fill (kN/m3)
the height of embankment fill (m)
In EBGEO 2010, it is similar to the equation above, but any little bite different for coefficient
of active earth pressure (see Fig.2-31).
Fig. 2-32 Piled embankments showing arching soil, notations for geometry and settlement
(after Zhuang, 2009)
45
When a height of embankment is lower than the critical height, the differential settlement
potentially develops there and vice versa. Table 2-11 below shows the critical height of some
design methods.
Table 2-11 Critical height for different design methods
Method Critical Height, Hc
Terzaghi (1943) 2.5 (s-a)
Carlsson (1987) 1.87 (s-a)
Hewlett&Randolph (1988) 1.4 (s-a)
BS8006 (1995) 1.4 (s-a)
Horgan & Sarby (2000) 1.54 to 1.92 (s-a)
Kempfert et al. (2004) sg / 2
46
wi = Pav . S1 . ij (2-83)
where: Pav = average load on a pile within group; S1 = settlement of a single pile under unit
load; ij = interaction factor for pile i due to any other pile j within the group, corresponding
to the spacing sij between pile i and j.
Fig. 2-33 Superposition via the interaction factor method (after Poulos and Davis, 1980)
Simplified or closed-form expressions for the interaction factors have been developed, thus
enabling a simpler computer analysis. Mandolini and Viggiani (1997) have developed the
following expressions for the interaction factor, in one of the following forms:
= A ( s / d )B (2-84)
or = C + D ln ( s / d ) (2-85)
where: A,B,C,D = fitting parameters
s/d = ratio of pile spacing to pile cap diameter
The value of A ranged between 0.57 and 0.98, while the range of B was -0.60 to -1.20. The
value of C is equal to 1.0 and D = -0.26.
The original interaction factors published by Poulos (1968) were based on the assumption
that the soil was a homogenous elastic medium, having a constant modulus with depth. This
was clearly a great simplification of reality, and in subsequent years, some significant
improvements and extensions have been made including non-uniform soil modulus, influence
of bearing stratum and interaction between two dissimilar piles.
Influence of non-homogeneity (see Fig. 2-34) compares relationship between interaction
factor and s/d for three cases, namely a homogeneous soil layer with a constant modulus Es
with depth, a soil where the surface modulus (Es0) is 3 times that at the base (EsL), and a non-
homogeneous soil layer whose modulus varies linearly with depth from zero at the surface (a
Gibson soil) but which has the same average modulus as the uniform layer. Influence of
stiffness of the bearing stratum against interaction factor which is effect of the stiffness of
bearing stratum Es2 as a multiple of the overlying soil Es1.
47
(b)
(a)
(c)
Fig. 2-34 Interaction factor, (a) Simplified case (b) Influence factor for non-homogeneity for
soil layer (c) Influence factor for bearing stratum (after Poulos,1968)
48
The lower value is more relevant to predominantly bearing piles, while the larger one is more
applicable to predominantly friction or floating piles. Randolph (1994) has related the
accuracy of equivalent pier method to the aspect ratio R of the group, where:
R = ( n s / L )0.5 (2-87)
where: n = number of piles, s = pile spacing, L = pile length. The equivalent pier method
tends to overestimate stiffness for values of R less than about 3, for values of R of 1 or more
provide that the pile spacing is not greater than 5 diameters.
For shallow soft soils we have a possibility to build the piled embankment using end-bearing
piles which bottom of columns is located on the hard stratum, but for deep soft soil it is not so
effective and economic if we construct this kind of structure. Therefore, the floating piles are
the best choice to overcome this problem. Though, settlement on soft soil is too high because
of high compressibility and creep is also the main problem.
49
2.5.10.2.6. Public Works Research Center Method
Public Works Research Center method (2000) as cited Han (2003) has come up with a design
method for reinforced embankments on deep mixed column.
Settlement of the columns is given as:
Sc = c . L / Ec (2-91)
where: Sc = settlement of the columns
c = stress on the columns
L = length of the columns
Ec = modulus deformation of the columns
Modulus deformation has a correlation with unconfined compression strength (qu) of the
column depending on a kind of column.
Ec = 100 qu (2-92)
The settlement of the untreated soil is given by:
Ss = So. s / p (2-93)
where: Ss = settlement of untreated soil subjected to reduced pressure, s
So = settlement of untreated soil subjected to the actual load of embankment, p
s = reduced pressure on the untreated soil due to embankment
p = total applied pressure of the embankment
The differential settlement illustrated in Fig. 2-35 between the soil and the columns in the
absence of geosynthetics reinforcement is given by:
S = Ss - Sc (2-94)
Fig. 2-35 Settlement on the end-bearing piles (after PWRC method, 2000)
When there is an inclusion of geosynthetics layer, the differential settlement can be given into
account an influence factor due to the inclusion of the reinforcement.
Sr = Ss / [1 + 2 (Ss / p)] (2-95)
where: Sr = differential settlement between the columns and the untreated soil
= influence factor due to the presence of geosynthetics layer
This influence factor is related to the tensile stiffness of geosynthetics reinforcement. The
relation between the two factors can be seen in Fig.2-36.
50
b = pile spacing
Fig. 2-37 Stress distributions for columns of deep mixed foundation (after Smith,2005)
The stress concentration in the stiffer column must greater than stress in the surrounding soil,
and then re-writing Eq. 2-46 for the stress concentration ratio, n. and coverage ratio, as, from
Eq. 2-49.
(2-96)
The load applied on embankment is carried by both the column and the soft soil. The average
stress applied by embankment, , may be expressed as:
= col as + soil (1-as) (2-98)
51
The vertical stress is carried by column and the surrounding soil as illustrated in Fig. 2-37
above as follows:
(2-99)
(2-100)
Fig. 2-38 Stress-strain relationship in dry mixed column (after Broms, 1999)
(2-103)
52
Where: S1 = the compression of the stabilized volume for case 1, di = stratum (or sub-layer)
thickness within the reinforced depth, Msoil = oedometer compression modulus of
surrounding soil. From both equations above that settlement of the stabilized volume
decreases with increasing area replacement ratio and with increasing the column stiffness.
In the second case when the creep limit of the columns is reached. The columns cannot take
load anymore. Therefore, subsequent loads carried by the unstabilized soil between columns
govern the settlement. Compression of the stabilized volume is calculated using the following
equation:
(2-104)
where: S2 = the compression of stabilized volume for case 2, = the creep strength
Broms (2003) suggests that the observed settlement can often be larger than the calculated
settlement. The differences between estimated and observed settlements generally increase
with increasing column lengths.
For partly penetrating columns or floating piles, the compression of a stratum of thickness d*
below the reinforced depth can be estimated for both cases. For case 1, the applied load is
assumed to be transferred directly down through the reinforced depth, di, and then it is
distributed through the underlying layer, d*, with angle of 1H:2V as shown in Fig. 2-39. The
compression of the underlying stratum, d*, and the reinforced depth, di, are the total
settlement.
For case 2 when the column creep strength is reached, the load be carried by columns is
transferred directly down through the reinforced depth (di) and then it is distributed through
the underlying layer (d*) at angle of 1H:2V. Meanwhile, the applied load that exceeds the
creep strength of the columns is applied to the ground surface and distributed through the
underlying soil at angle of 1H: 2V.
di
For all methods mentioned above, settlement is a procress of water dissipation from soil body
in which it consists of lot of fine grained soil. Some methods can be applied to accelerate this
process such as vertical drains, pre-loading, vacuum technique, and electro-osmosis. Zhuang
53
et al. (2006) conducted electro-osmotic consolidation test for small-scale model. They
concluded that soil properties like c, , max exhibit increasing after treatment of electro-
kinetic geosynthetics (EKG).
Where S0 is embankment settlement constructed on soft soil without the support of piles and
S is the settlement of an embankment supported by piles.
54
CHAPTER 3
There are four kinds of material to be characterized regarding the topic research: subsoil,
embankment, geosynthetics, and piles. It is important to know well some properties of them.
Because stress and strain relationship in material is induced by loading either static or
dynamic loading, and the characteristic of loading, including the magnitude and movement of
loading on transportation infrastructure is necessary to be understood.
1. Unit weigth (gr/cm3) 1.63-1.76 1.68-1.75 1.67-1.74 1.75 – 1.86 1.68 – 1.73 1.63-1.89 1.61-1.79 1.73 1.59 – 1.71
2. Clay content (%) 22-43 25-40 24-34 30 - 52 30 -52 30-61 23-45 44 25 - 58
3. Liquid limit (%) 80-110 72-108 83-94 53 -107 53 - 107 72-130 62-90 81 82 -104
4. Plastic Index (%) 52-79 40-74 44-58 24 - 57 24 - 97 39-79 28-45 48 46 - 62
5. Linear shrinkage (%) 18-22 12-26 - - 18 15-27 - - -
6. Water content (%) 32-48 37-53 29-49 24 - 40 27 -32 40-55 38-53 34 34 - 52
7. Passing # 200 (%) 83-98 82-98 95-98 73 - 96 76-96 92-98 89-94 97 92 - 94
8. Classification of soil CH, clay CH,clay CH CH,Silty clay CH, clay CH, silty clay CH,clay CH,clay MH,Silt
9. Mineral of clay Montmori Montmori Montmori Montmori Montmori Montmori Montmori Montmori Montmori
lonite lonite lonite lonite lonite lonite, 60 % lonite 45 % lonite lonite, 10 %
10.Colour Browny grey Browny grey Blackish Greyey black Greyey black Blackish grey Blackish Blackish grey Blackish grey
grey grey
Location : - Depth (m) 2.0 -.3.0 1.0 -3.0 2.0 -5.0 1.0 - 4.0 2.0 – 5.0 1.0 -.4.0 1.0 -3.0 1.0 -2.0 1.0 - 3.0
- Site (KM) 36-43 13.8 38.8 37.5 – 58.25 23.0 -27.0 5..0 – 19.0 12.8 – 14.45 16.50 25.50 – 69.60
Source: Final Report of Guidance for Roadway Construction over Expansive Soil, 2003
Soft soils are widespread in a lot of locations in Indonesia. In Java island, the soft soils
mostly consist of clay and/or silt, whereas in Sumatra and Kalimantan they are not only soft
clay but some regions covered by peat soils. Java Island is most dense in population and the
55
island has many infrastructures particularly roadways, railways and runways constructed over
the soft soil.
Chen (1975) uses a single index based on plasticity index to identify expansive soil. For
Plasticity Index (PI) values ranging from 20 to 55 are high for swelling level and very high
for PI values more than 55. Whilst Seed et al. (1962) use equation Ac = PI / (CF-10) to
identify activity level of soil. For value of Ac that more than 1.25 the soil has high level of
activity.
1. Soil classification - Low organic Low organic Low organic Low organic
2. Water content, Wn % 58.22-169.98 35 – 91 81.6-109.4 81.6-91
3. Organic content, OC % 10 - 11.88 10-11.9
4. Specific gravity , Gs - 2.2-2.6 2.36 – 2.70 2.611 2.4-2.6
5. Unit weight, (kg/m3) kN/m3 12.3-16.49 13.9 – 18.5 13.94-14.31 13.9-14.3
6. Liquid Limit, LL % 20-70 17.10 - 62.46 84.52-99.25 62.5
7. Plastic Limit, PL % 17-35 14.62 – 38.44 39.6-50.22 35-39
8. Plasticity Index, PI % 5-35 2.48 – 29.67 34.3-59.6 29.7-35
9. Void ratio, e % 1.02-3.30 0.88 – 2.65 - 1-2.6
10. Cohesion (UD), c kN/m2 6-16.5 7.2-19.8 0.38-8.08 7.2-8
11. Friction angle (UD), o 1.0-20.3 3.03-13.66 13.61-27.26 13.6-13.7
12. Compression index, Cc - 0.16-0.34 0.278-1.663 - 0.28-0.34
13. Coeff. of compressibility - - 1.4E-3 to 7E-3 - 1.4E-3 to 7E-3
14. Permeability, kv m/day 4.4E-6 to 7.8E-4 - - 4.4E-6 to 7.8E-4
15. Unconf. comp. strength kN/m2 - 8.4-65.5 - 8.4-65.5
16. Oedometer modulus kN/m2 550 - - 550
17. Young’s modulus kN/m2 650-1166 - - 650-1166
Based on Chen (1975) approach, for soil in this region with PI value average 19.5 can be
classified as high for swelling level.
56
3.1.2. Mechanical Properties
Mechanical properties are a necessary thing when investigating strength and deformation of
soil during loading. Shear strength of cohesive soil can be determined by using shear tests
and/or triaxial tests in which parameters cohesion, c, and internal friction angle, , can be
measured. Whilst in predicting the rate of settlement for the soil can be carried out by using a
consolidation test or oedometer testing in the other to obtain consolidation index, cc.
3.1.2.2. Compression
When soil undergoes a loading, because of their relatively low permeability, their
compression is controlled by the rate at which water is squeezed out of the pores. The slope e
against log ’ plotted in normally consolidated soil is referred to as the compression index,
cc. The load increment ratio was uniform where the loading was from 25 to 800 kN/m2. The
compression index, cc, varies widely with the increasing depth, however, the depth does not
influence of cc.
Priadi (2008) characterized the Pontianak soft organic soil compressibility behavior that the
top layer (around 10 m deep) is highly compressible ranging from 0.5 to 1.38 with an average
value of about 0.8, whereas at below this layer ranging from 0.2 to 0.5 with an average value
is about 0.3. Meanwhile, the recompression index, cs, ranges widely from 0.03 to 0.25. Some
of the 1-D Oedometer test results are shown in Fig. 3-1.
Fig. 3-1 Oedometer test of Pontianak soft organic soil (after Priadi, 2008)
57
The over consolidated ratio, OCR, is defined as the ratio between the pre-consolidation stress
and the effective in-situ stress. OCR is a state parameter that indicates the amount of over-
consolidation of the soil (Brinkgreve, 2001). This value notably reduces with a depth.
Pontianak soft organic soils are heavily over consolidated from the ground surface to about 5
m depth due to the wetting and drying cycles during deposition. The over consolidation ratio
ranges from 2 to 11 at this layer, whereas OCR values range from 1.3 to 2 are found at 5 to 20
m depth.
58
Table 3-4 Various materials for embankments
Unit weight
No. Materials
(kN/m3)
1 Sand 18-22
2 Cohesive soil 16-19
3 Corduroy 7.0
4 Rubber slag 4.0-6.0
5 Pumice 10.9
6 Dreg saws < 10.0
7 Peat bales < 10.0
8 Expanded polystyrene (EPS) 0.2-0.4
Source: Indonesian Geotechnical Guidance-4 (2001)
The Indonesian Geotechnical Guidance-4 (2001) gives the design parameters when using
material as embankment fill as depicted in Table 3-5.
Table 3-5 Design parameters for embankment material
Geographical Zone
Parameters Unit
A B
Unit weight,
3
kN/m 18 20
Undrained shear strength, Su kN/m2 100 100
Cohesion, c ' kN/m2 10 5
Internal friction angle, ' [O] 35 30
Source: Indonesian Geotechnical Guidance-4(2001)
A Java island (vulcanic rocks)
B Sumatra, Kalimantan, Sulawesi. Papua island (sedimentary and metamorphic
rocks)
59
Webb et al. (1986) have reported a number test of cohesionless soils in repeated load triaxial
test following the AASHTO procedure as equations follows:
(3-3a)
(3-3b)
Sasongko (1996) also reported a correlation between CBR and Mr using repeated load triaxial
test in laboratory as shown equations follows:
(3-4a)
(3-4b)
(3-4c)
NCHRP (2004) provided a correlation between CBR and Mr for a number of cohesionless
soils in repeated load triaxial test following the AASHTO procedure. A typical equation for
medium clay sand is shown equation (3-5).
(3-5a)
(3-5b)
Several constitutive models have been proposed by many researchers for modeling resilient
moduli of soils and aggregates. Dunlap (1963) suggested the following relationship for
presenting resilient modulus:
(3-6)
(3-7)
(3-8)
(3-9)
Uzan (1992) introduced the octahedral shear stress in place of deviator stress in equation (3-
9), which provided a better explanation for the stress state of the material, in which the
60
normal and shear stress change during loading. The proposed model is known as the k1- k3
model. The universality of this model stems from its ability to conceptually represent all
types of soils from pure cohesive soils to non-cohesive soils.
(3-10)
where
The coefficients k1, k2, and k3 are constants, depending on the state and quality of unbound
granular materials. Since coefficient k1 is proportional to Young’s modulus, it should always
be positive as Mr can never be negative. The coefficient k2 should be positive, because
increasing the volumetric stress produces stiffening or hardening of the material, yielding
higher modulus. The coefficient k3 should be negative because an increase in the shear stress
softens the material, thereby yielding lower modulus. If nonlinear property coefficients k2 and
k3 are set to zero, then the model can be simplified as linear elastic. If k3 is zero, the behavior
could be non-linear hardening and if k2 is zero, the behavior is non-linear softening.
where is unit weight of embankment fill and cu is the undrained shear strength of subsoil
beneath embankments.
The behaviour of soil under cyclic loading is non-linear and dependent on some factors
including soil type, confining pressure, number of loading cycles and amplitude of loading.
Non linear hysteretic soil behaviour is commonly characterized by a viscous damping and
equivalent shear modulus (Seed and Idriss, 1970; Hardin and Drenevich, 1972). Definition of
damping is a measure of energy dissipation. It increases with increasing magnitude of cyclic
shear strain, whereas shear modulus decrease with increasing magnitude of cyclic shear
strain. It is also known that dynamic properties of soil are influenced by the plasticity index,
void ratio, relative density and number of cycles (Cabalar and Cevik, 2008).
61
Table 3-6 Correlation between Vs and N
Shear wave velocity, Vs
Author(s) Soil type
(meter/second)
Kanai, et al. (1966) All Vs = 19 N 0.6
Shibata (1970) Sand Vs = 32 N 0.5
Ohba & Toriuma (1970) Alluvial Vs = 85 N 0.31
Ohta, et al. (1972) Sand Vs = 87 N 0.36
Ohsaki & Iwasaki (1973) All Vs = 82 N 0.39
Cohesionless Vs = 59 N 0.47
Imai & Yoshimura (1975) All Vs = 92 N 0.329
Imai, et al. (1975) All Vs = 90 N 0.341
Imai (1977) All Vs = 91 N 0.337
Imai & Tonouchi (1982) All Vs = 97 N 0.314
Imai & Yoshimura (1990) All Vs = 76 N 0.33
Ohta & Goto (1978) All Vs = 85 N 0.348
Sands Vs = 88 N 0.34
Gravels Vs = 94 N 0.34
JRA (1980) Clays Vs = 100 N 0.333
Sand Vs = 80 N 0.333
Seed & Idriss (1981) All Vs = 61 N 0.5
Seed, et al. (1983) Sands Vs = 56 N 0.5
Sykora & Stokoe (1983) Granular Vs = 100 N 0.29
Okamota, et al. (1989) Dilluvial sands Vs = 125 N 0.3
Lee (1990) Sands Vs = 57 N 0.49
Clay Vs = 114 N 0.31
Silts Vs = 106 N 0.32
Yokota, et al. (1991) All Vs = 121 N 0.27
Jafary, et al. (1997) All Vs = 22 N 0.85
Clayey soils Vs = 27 N 0.73
Silty soils Vs = 22 N 0.77
Clayey & Silty Vs = 19 N 0.85
Source: Jafary et al., (2002)
Once shear wave velocity is determined, and then shear modulus of material, G, is obtained
using equation G= .Vs 2, which is the density of soil.
Shear modulus of air dry clean sands at a certain level of shear strain, , can be represented
approximately by the following empirical equation irrespectively of kinds of sands.
( = 10-6) (3-12)
( = 10-5) (3-13)
( = 10-4) (3-14)
62
Where G is shear modulus in kg/cm2, p is mean principle stress in kg/cm2 and e is the void
ratio. Eq. 3-5 is identical to the empirical equation for round Ottawa-sand proposed by
Hardin, et al. (1972).
All tests demonstrated the well known dependence of Gmax on effective confining pressure 'o
and density expressed in terms of the void ratio e. In other to an analytical expression for the
Gmax = Gmax ('o) relationship, the general equation suggested by Hardin (1972) is adopted.
(3-15)
Where: S is a stiffness coefficient. The experimental results are closely approximated by
setting S= 420 and n= 0.6. It should be noticed that the value of the exponent n is higher than
widely used for cohesive and cohesiveless soils which ranges between 0.4 and 0.5. The
variation of shear modulus with shear strain amplitude as expressed in the following
equation.
(3-16)
Damping ratio D was found to be essentially independent on confining pressure and density.
An average value Dmin=2% was determined from all tests. The increase of damping with
shear strain amplitude may be approximately expressed by:
(3-17)
(3-18)
3.3. Geosynthetics
3.3.1. Material Properties
The geosynthetics terminology may be based on the subdivision by PrEN ISO 10318.
According to this standard “Geosynthetics” is a generic term describing a product at least one
of whose components is made from a synthetic or natural polymer, in the form of a sheet, a
strip or a three dimensional structure, used in contact with soil and/or other materials in
geotechnical and civil engineering applications. As depicted in Fig. 3-2 that geosynthetics can
be differentiated into permeable and impermeable products.
63
Geosynthetics
Geocomposites
The use of a geosynthetics in pavement system reinforcement is to aid in support of the traffic
load. Traffic loads may be vehicular loads experienced over the life of the pavement. Base
(or subbase) reinforcement is a treatment using of a geosynthetics as a tensile element at the
bottom of base (or subbase) or within a base course and is designed to address the pavement
distress mode of pavement surface deformation or rutting and asphalt fatigue cracking.Whilst
subgrade restraint is the use of geosynthetics at the subgrade/subbase or subgrade/base
interface to increase the support of construction over a weak or low strength subgrade
(Barenberg, 1980; Steward et al., 1977; Giroud and Noiray, 1982; Holz et al., 1987).
The following benefits of using geosynthetics in roadways are identified (TenCate Mirafi,
2010):
1. Reducing the intensity of stress on subgrade (function: separation).
2. Preventing subgrade fines from pumping into the base (function: filtration).
3. Preventing contamination of the base materials allowing more open graded, free-
draining aggregates to be considered in the design (function : filtration).
4. Reducing the depth of excavation required for removal of unsuitable subgrade
materials (function: separation and reinforcement).
5. Reducing the thickness of aggregate required to stabilize the subgrade (function:
separation and reinforcement).
6. Minimizing disturbance of the subgrade during construction (function separation
and reinforcement).
7. Assisting the increase in subgrade strength over time (function: filtration).
8. Minimizing the differential settlement of roadway, which helps maintain
pavement integrity and uniformity (function: reinforcement).
9. Minimizing maintenance and extending the life of the pavement (function: all).
Others important findings from laboratory and/or field studies include the following:
64
1. An optimum benefits when the geosynthetics was placed at the bottom of a 200-
300 mm thick base layer.
2. For thicker base sections, the most beneficial reinforcement location appeared to
be in the middle of the base, where geogrids were found to perform best.
3. For thin bases (less than 200 mm), lack of separation was noted as a potential
problem for geogrids. Geogrid-geotextile composites tend to perform better for
thin bases, especially where subgrade strengths were below a CBR of 3%.
4. Reinforcement benefits were observed with subgrade strengths up to a CBR of
8%.
Benefits using geosynthetics as reinforcement can be defined by TBR (traffic benefit ratio)
and BCR (Base course reduction) as shown in Table 3-7. The TBR is defined as the ratio of
number of cycles necessary to reach the same rut depth for a test section containing
reinforcement to unreinforced section with the same section thickness and subgrade
properties. Furthermore, BCR is expressed as a percentage savings of the unreinforced base
course thickness.
Table 3-7 TBR and BCR resulted from laboratory and field test
Materials TBR BCR
Geotextiles:
Range 1 – 220 22 – 33 %
Typical value 1.5 – 10
Geogrids:
Range 0.8 – 670 30 – 50 %
Typical value 1.5 – 70
After TenCate Mirafi, 2010
Besides the ratio coefficients (TBR, BCR), the following properties are considered to
influence performance : tensile strength at 1%, 2% and 5% strain, coefficients of pullout and
direct shear, aperture size (grids) and percent open area (geotextiles) and stiffness properties
including the flexural rigidity and aperture stability. For subgrade restraint applications, the
properties of tensile strength at 2% and 5% strain are primarily related to geosynthetics
performance.
1. Weak Subgrade (CBR <3), For Thin (≤250 mm) Base Sections and Thick (>250 mm)
Base Sections.
When a weak subgrade exists, woven geotextile should be placed at the surface interface. In
addition, if the required base course is greater than 250 mm (10 in), a second layer of
reinforcement, biaxial geogrid, should be placed in the middle of the base course section.
2. Firm Subgrade (CBR >3), For Thin (≤250 mm) Base Sections and Thick (>250 mm)
Base Sections.
65
For a firm subgrade and relatively thin base course section is designed biaxial geogrid at the
subgrade interface. While geogrid can be placed in the middle of the base course when firm
subgrade and relatively thick base course section is designed.
(a) Weak Subgrade, Base course ≤250 mm (b) Weak Subgrade, Base course >250 mm
(c) Firm Subgrade, Base course ≤250 mm (d) Firm Subgrade, Base course >250 mm
Fig. 3-3 Position of geosynthetics in pavement design practice (after Mirafi, 2010)
(a) (b)
Fig. 3-4 Typical strain vs. Force behaviour of reinforcement (a) Exxon, 1989
(b) Carlson, 1987
66
Typical short term strength of geosynthetics is described in Table 3-8 for various raw
materials (EBGEO, 2010; Althoff, 2011).
Table 3-8 Tensile strengths of geosynthetics (after EBGEO, 2010)
Typical short term strengths Typical elongation
Raw material Product types [kN/m] at failure [%]
from to max. from to
AR Woven geogrids 40 1200 2200 2 4
(aramides) Woven geotextile 100 1400 300 2 4
PE Woven geogrids 20 150 300 15 20
(polyethylene) Extruded geogrids 40 150 200 10 15
Woven geotextile 30 200 400 15 20
PET Woven geogrids 20 800 1200 8 15
(polyester) Bonded geogrids 20 400 500 6 10
Woven geotextile 100 1000 1600 8 15
PP Woven geogrids 20 200 500 8 15
(polypropylene) Bonded geogrids 20 200 400 8 15
Extruded geogrids 20 50 - 8 20
Woven geotextile 20 200 600 8 20
PVA Woven geogrids 30 1000 1600 4 5
(polyvinyl alcohol) Woven geotextile 30 900 1800 4 5
When choosing geosynthetics for reinforcement, there are two factors that must be
considered well namely internal and external factor. Internal factors such as tensile strength,
creep properties, whereas external factors such as kind of embankment fill, endurance against
environment (ultra violet, acidic or alkaline matter, micro-organism. In other to cover all
conditions, strength of geosynthetics has to be adjusted using some partial factors. According
to Nordic guideline (2003) some conversion factors for design purpose are listed in Table 3-9
through Table 3-11.
Table 3-9 Conversion factors of geosynyhetic reinforcements
(after Nordic guideline, 2003)
Conversion parameters Conversion factor
Creep factor
Installation damage
Biological and chemical degradation
67
Table 3-11 Conversion factors for damage during installation
(after Nordic guideline, 2003)
Conversion parameters Conversion factor, Material factor, fd
Clay/silt 0.91 1.1
Sand 0.83 1.2
Gravel (Natural) 0.77 1.3
Gravel (Broken) 0.72 1.4
Chrused Rockfill 0.67 1.5
According to Swedish Road Administration publication 1992:10, the material factors for
biological and chemical degradation, fenv, may be assumed 1.1 as long as the pH value ranges
between 4 and 9, which gives a conversion factor of = 0.91.
Allowable tensile strength of geosynthetics for reinforcement design is defined as ultimate
tensile strength divided by the reduction factor (or partial factor).
(3-21)
There are three classifications of columnar foundation include (Han and Wayne, 2000)
namely:
flexible column (such as stone columns and lime columns)
semi-rigid columns (such as lime-cement and soil-cement columns)
rigid piles (such as concrete pile, timber piles, and vibro-concrete piles)
68
According to the new code SNI 2002 (Indonesian National Standardization, 2002), quality of
wooden material or quality code use mixed Letter and Number to declare Elastic Modulus as
in Table 3-13.
Table 3-13 Quality code for wooden material (after SNI, 2002)
Code of Young’s Flexural Tensile strength Tenslie strength Shear Compr. strength
quality modulus strength parallel fiber perpendic. fiber strength perpendic.
Fb Ft Fc Fv Fc
[MPa] [MPa] [MPa] [MPa] [MPa] [MPa]
E26 25000 66 60 46 6.6 24
E25 24000 62 58 45 6.5 23
E24 23000 59 56 45 6.4 22
E23 22000 56 53 43 6.2 21
E22 21000 54 50 41 6.1 20
E21 20000 50 47 40 5.9 19
E20 19000 47 44 39 5.8 18
E19 18000 44 42 37 5.6 17
E18 17000 42 39 35 5.4 16
E17 16000 38 36 34 5.4 15
E16 15000 35 33 33 5.2 14
E15 14000 32 31 31 5.1 13
E14 13000 30 28 30 4.9 12
E13 12000 27 25 28 4.8 11
E12 11000 23 22 27 4.6 11
E11 10000 20 19 25 4.5 10
E10 9000 18 17 24 4.3 9
Elastic modulus or Young’s modulus of wooden material can be estimated using equation
below (SNI, 2002).
(3-22)
Where: G is specific gravity of wooden material at water content 15%.
69
3.4.3. Stone Column
Hughes and Withers (1974) performed pioneering laboratory studies of sand columns within
a cylindrical chamber containing clay and used radiography to track the deformations
occurring within and outside the columns. They found that CCET (cylindrical cavity
expansion theory) represented the measured column behaviour very well and proposed that
the ultimate vertical stress (q) in a stone column could be predicted by:
(3-25)
where: ' is the friction angle of stone infill, 'ro is the free-field lateral effective stress and c
is the undrained shear strength.
The equation above is widely used in practice today. There are alternative approaches for
estimating the bearing capacity of a single column and column group, such as that recently
published by Etezad et al. (2006). The authors report an analytical treatment of bearing
capacity failure mechanisms. Failure mechanisms adopted are based upon the output from a
combination of finite element analysis and field trials.
Absolute and differential settlement restrictions usually govern the length and spacing of
columns, and the preferred method of estimating post-treatment in European practice was
developed by Priebe (1995). Although, this method is strictly applicable to infinite array of
columns and has some empiricism in its development.
Priebe's settlement improvement factor, n, defined as:
(3-26)
It is a function of the friction angle of stone ' , the soil's Poisson's ratio and an area
replacement ratio dictated by the column spacing. The area replacement ratio is defined as
Ac/A , where Ac = cross-sectional area of one column and A = total cross-sectional area of the
'unit cell' attributed to each column. Ac/A is related geometrically to the column radius, r, and
column spacing, s, according to:
(3-27)
Where: k is and 2/√3 for square and triangular column grids respectively.
70
Priebe's 'basic improvement factor' may be derived from the chart as shown in Fig. 3-6. Need
to be noted that the reciprocal area replacement ratio A/Ac is used on the chart.
A lower limit to the undrained strength of cu = 15 kPa is suggested for treatment with stone
column, although there have been situations where softer soils have been successfully
improved (Raju et al., 2004). In other hands, UK National House Building Council (NHBC,
1988) suggests that stone columns should not be used when Ip >40%. Wood et al., (2002)
conducted what is considered to be the most comprehensive laboratory model investigations
of large groups of columns. The results suggest that significant improvement in the bearing
capacity requires an area replacement ratio of 25 % or greater.
McKelvey et al. (2004) used a transparent medium with ’clay-like’ properties to allow visual
monitoring of the columns throughout the foundation loading. The main findings of this
research relate to optimum column aspect ratio L/d (L=column length, d=column diameter)
that in the case of ’short column’ (i.e. L/d=6), bulging took place over the entire length of
column. The ’long column’ (L/d=10) deformed significantly in the upper region whereas the
bottom portion remained undeformed. McKelvey et al. (2004) postulated a ’critical column
length’ of L/d=6, which is in keeping with earlier work (Hughes and Withers, 1974; Muir
Wood et al. 2004).
71
Table 3-14. Typical Japanese mixing installation parameters (after Kaiqiu, 2000)
Description Single drive shaft Double drive shaft Multibarrel drive shaft
Depth of stabilization 49 ft (15 m) > 49 ft (15 m) 98-131 ft (30 - 40 m)
Penetration velocity 2 - 3.3 ft/m 0.7 - 3.3 ft/m 3.3 - 6.6 ft/m
0.6 - 1.0 m/min 0.2 - 1.0 m/min 1.0 - 2.0 m/min
Withdrawal velocity 2 - 3.3 ft/m 0.7 - 3.3 ft/m 3.3 - 6.0 ft/m
0.6 - 1.0 m/min 0.2 - 1.0 m/min 1.0 - 1.5 m/min
Rotating speed 50 rpm 46 rpm 20 - 30 rpm (penetration)
40 - 60 rpm (withdrawal)
Deep mixing methods in the U.S. have been used on several projects either dry or wet
method. In general, dry mixed stabilization is appropriate for sites with relatively deep
deposits of very soft soil, and sufficient groundwater to hydrate both the lime and cement
(Esrig and Mac Kenna, 1999). Cohesive soils with moisture contents between 60% and 200%
are best suited for dry mixing.
While several different types of laboratory tests are used to evaluate the shear strength and
stiffness of deep mixed columns, the most frequently used is the unconfined compression test,
mainly because of the simplicity of the test. Many factors affect the unconfined compressive
strength because of a wide variety of soil types and binder mixes. The 28-day unconfined-
compressive strengths for soil treated by the wet method may range from 140 to 27000 kPa
(Haley&Aldrich, 2000; Kaiqiu, 2000; Tatsuoka&Kobayashi, 1983) whereas using the dry
method range from 14 to 2700 kPa (Hebib&Farrell, 2002; Jacobson et al, 2002; Kaiqiu,
2000). Unconfined compressive strengths, qu, for three projects in the U.S. are presented in
Table 3-15.
Table 3-15 Specified values of qu on deep mixing projects in the U.S.
Soil types / binder
Projects Specified qu Reference(s)
amount
Oakland Airport Wet method; Loose At 28 days, Average qu > Yang et al,2001
Roadway, California sandy fill and soft 1035 kPa, Minimum qu >
soil; 160-240 kg/m3 690 kPa
cement
Central Artery Project, Wet method; Fill At 56 days, Maximum qu Lambrechts et
Boston and organic soft > 26900 kPa, Minimum al,1998;
clay; 220-300 kg/m3 qu > 2100 kPa Maswoswe,2001
cement
I-95 Route 1, Wet method; Soft At 28 days, Average qu > Shiells et al,2003;
Alexandria organic clay; 300 1100 kPa, Minimum qu > Lambrechts et
kg/m3 cement 690 kPa al,1998
Source, M. Smith, 2008
Stabilization of soft organic soils with cement columns using the mix-in-place technique
(MIP) for a railway embankment at section of Büchen-Hamburg was upgrade in 2003 by the
German Railway company (Deutsche Bahn) to allow a train speed of 230 km/h
(Schwarz&Raithel, 2005). The cement columns (diameter 0.63 m and 5-8 m length) were
installed in a square 1.5×1.5 m grid, containing 2.5 to 3% cement, which can be characterized
72
as a wet deep mixing technique, the composition of binder (water, cement and bentonite) and
the water binder ratio (approx. 1.0). Each of 500 m3-treated soil, 6 unconfined compression
tests, was carried out after 28 days. According to the test, unconfined compressive strength
after 28 days of all samples exceeds the design criteria of qu ≥2.2 Mpa.
For cemented columns installed by the wet method, Takenaka (1995) reported that undrained
shear strength is equal to one-half of the unconfined compressive strength for those values
below several hundred kPa and become less than one-half when they are greater than several
hundred kPa. As a rule of thumb, Takenaka (1995) recommended that undrained shear
strength be taken as one-third of the unconfined compressive strength. Kivelo (1997) found
that the undrained shear strength can be less than one-half the unconfined compressive
strength at low confining pressures. However, when the total confining pressure exceeds 150-
250 kPa, the undrained shear strength becomes almost constant at a value equal to one half of
the unconfined compressive strength.
Fig. 3-7 Undrained shear strength of lime/cement column (after Kivelo, 1997)
The peak strength is typically reached at strains 0f 1% to 2% and decrease in strength once
the peak strength is exceeded (Kivelo, 1998). The residual strength of soil-cement is 65% to
90% of the unconfined compressive strength (Tatsuoka&Kobayashi, 1983).
The undrained secant modulus of elasticity, E50, which evaluated at 50% of the peak strength,
is a measure of soil-cement compressibility. Some researchers correlate E50 to unconfined
compressive strength for columns installed the dray method (Braker, 2000; Broms, 2003;
Jacobson et al., 2003; Navin&Filz, 2005). Whilst for cement treated soils using the wet
method also have been performed and presents relatively higher values of secant modulus of
elasticity than those using the dry method (Kawasaki et al., 1981, Navin&Filz, 2005, Fang et
al., 2001). The relationship between E50 and qu is provided in Table 3-16.
Table 3-16 Relationship between E50 and qu
Binder types E50 Reference(s)
50 - 180 qu Baker, 2000; Broms, 2003
Dry lime/cement
75 qu Jacobson et al., 2003
65 - 250 qu Baker, 2000; Broms, 2003
Dry cement
300 qu Navin and Filz, 2005
350 - 1000 qu Kawasaki et al., 1981
30 - 300 qu Fang et al., 2001
Wet cement
150 qu Mc Ginn and O'Rouke, 2003
300 qu Navin and Filz, 2005
Source:M. Smith, 2008
73
When the modulus of elasticity is used in design analysis, the secant modulus E50 is typically
used as the design value of column Ecol. The modulus of elasticity on samples prepared in the
laboratory is typically higher than modulus determined from coring test obtained in situ
actual columns (Broms, 2003).
The oedometer compression modulus Eoed is related to modulus of elasticity Ecol and
Poisson's ratio as follows:
(3-28)
Generally the Poisson's ratio of deep mixed treated soil is around 0.25 to 0.45 (Terashi,
2003). Therefore, Eoed is equal to 1.2 to 3 Ecol.
The total unit weight of treated soil using the dry method increases from 3% to 15% above
the untreated soil. Whilst tensile strength of soil improved by the wet method, it is 10% to
20% of unconfined compressive strength. Moreover, permeability of treated soil ranges from
10-7 to 10-8 m/s is routinely achievable.
Transport infrastructure such as roadway, railway and runway is mainly subjected by moving
load. Although at a certain situation they are static loading such as car parking at the parking
lot, airplane parking on a parking stand at the apron. When a vehicle passes through a
roadway, time loading at certain point on the surface of the roadway section depends on
velocity of the vehicle. Load repetition induced by wheels of the vehicles on the surface is
able to result in a rut depth during a service period of roadway.
3.5.2. Airplane
Compared to a wheel load generated by traffic on the roadway, the wheel load of airplane has
quite higher in magnitude than that of on the roadway. It depends on weight, tire pressure,
wheel configuration of airplane. Tire pressures of the airplane vary from 0.5 MPa to more
than 1.5 MPa.
Because weight of the airplane is quite heavy that will be transferred to the surface of the
pavement, wheel configuration plays an important role. For light aircraft, it uses a single
wheel configuration. Dual wheel configuration can be seen for a moderately weight of the
airplane. Furthermore, another configuration for heavy weight of the airplane is tandem and
also dual tandem.
In other to accommodate the various airplanes for operational movement, it needs to provide
aerodrome areal as shown in Table 3-17. The aerodrome reference code uses number and
letter codes to express class of airport (ICAO, 1999). Maximum allowable tire pressure
category consists of High (no pressure limit), Medium (pressure limited to 1.50 MPa), Low
(pressure limited to 1.00 Mpa) and Very low (pressure limited to 0.50 MPa. The higher tire
pressure indicates the heavier weight of the airplane.
74
Table 3-17 Aerodrome reference code (after ICAO, 1999)
Code element 1 Code element 2
Code Aeroplane reference field Code Wingspan Outer main gear wheel
number length letter span *
[1] [2] [3] [4] [5]
1 Less than 800 m A Up to but not including Up to but not including
15 m 4.5 m
2 800 m up to but not B 15 m up to but not 15 m up to but not
including 1200 m including 24 m including 24 m
3 1200 m up to but not C 24 m up to but not 6 m up to but not
including 1800 m including 36 m including 9 m
4 1800 m and over D 36 m up to but not 9 m up to but not
including 52 m including 14 m
E 52 m up to but not 9 m up to but not
including 65 m including 14 m
F 65 m up to but not 14 m up to but not
including 80 m including 16 m
* Distance between the outside edges of the main gear wheels
3.5.3. Trains
In Germany, according to Ril 836 for rail infrastructure, the subgrade or improved subgrade
has to be able to support a load above this surface layer with bearing capacity at least around
52 kPa. For high speed trains from 100 to 300 km/h, it needs the additional layer around 1 to
2 times of the superstructure thickness (Muncke et al., 1999; Kempfert et al., 1999).
It is different for each country with respect to a standard axle load. For instance in Greece,
Beskou et al., (2011) reported that locomotive (or engine) is around 210 kN and 150 kN for
carriage. In Indonesia, railway infrastructure is subjected to axle load maximum 180 kN
(Indonesian Railway Code, 2003).
75
CHAPTER 4
Experimental Works
and Field Case Studies
4.1. General
To understand well the pavement behaviour and/or embankment overlying soft soil, some
researchers have performed some experimental works either in laboratories using small scale
approach or monitoring full-scale directly in the field (Zaeske, 2001); Heitz, 2006; Wei-ping
et al., 2007; Hassandi et al., 2005; Schwarz et al., 2005; Djarwadi, 2006; Almeida et al.,
2008; Eekelen, 2009). Treatments for ground improvement in this study only focus on
horizontally inclusion of geosynthetics, vertically inclusion of piles and the combined
technique both of them. Meanwhile, some tests have been done to know an influence of
geosynthetics into granular soil (Bussert, 2006; Ruiken&Ziegler, 2008; Ruiken et al., 2010).
Furthermore, some intensive tests regarding with geosynthetics-soil interaction behaviour
have been carried out using interaction testing device (ITD) at Geotechnical Institute of
Technical University Bergakademie Freiberg (Aydogmus&Klapperich, 2004; Aydogmus,
2006; Tamaskovics&Klapperich, 2010; Althoff, 2010, 2011).
Fig. 4-2 Radial strains in laboratory test (after Ruiken and Ziegler, 2008)
Some studies to observe the horizontal deformation subjected to a vertical load has been
performed by researchers (Clayton et al., 1993; Soong & Koerner, 1997; Yang et al, 2009).
Ruiken and Ziegler (2010) also performed the deformation by inclusion geosynthetic as
shown in Figure 4-3. A sand sample in box sized (H×B×L) 1m×1m×0.45m compacted up to
100% standard proctor was used to do the experiment as shown in Fig. 4-3. A couple side
was made of steel and another couple side was glass with 106 mm thick. To reduce friction
between sand and glass, a thin latex membrane inserted in between and a uniform load up to
50 kN/m2 was applied at the surface of the sample box.
Fig. 4-3 Installation for measuring horizontal pressure (after Ruiken et al., 2010)
Inclusion of geogrid from one up to 5 layers has been applied and then the experiment result
was a curve correlating between horizontal deformation of vertical movable plane and
horizontal outward pressure on the vertical plane as depicted in Fig. 4-4. It was interesting
77
that for 2 and 3 inclusion of geogrids they almost coincided with each other. Moreover, using
of 5 layers did not give lower horizontal pressure than that of using 4 layers. Overall, more
and more numerous layers suggest that horizontal outward pressure was reduced.
1. Side frame
2. Base plate
3. Load plate with reinforcement
4. Threaded rods
5. Plug gauge with fine thread
6. Force measurement
7. Movable front plate
8. HDPE coating with PE membrance
9. Displacement transducer
10. Sand/gravel
11. Geosynthetic layer
Fig. 4-5 Test setup with movable front wall (after Bussert, 2006)
The test results showed that interaction of soil and geogrid mainly depends on the
geosynthetic layer spacing, soil grain size, geosynthetic aperture size as well as strength of
shape and extensional stiffness of geogrids. Load plate was subjected to a load (in kPa) and
strain of geosynthetics was measured based movement of movable wall. Magnitude of load
78
in kPa was recorded as the stress reduction at the front of the moving wall in horizontal
direction caused by different geogrid products as reinforcement and also without
reinforcement is shown in Fig. 4-6.
Fig. 4-6 Earth stress reduction at movable front wall with different geosynthetics
and without reinforcement of soil body (after Bussert, 2006)
Tests of friction and shear behaviour are costly and time-consuming. Pull-out test in
particular is one of the most expensive performance tests (Koerner, 2005). Thus, in previous
research a more economical technique for using the ITD (for IRV, ISV, and IPV) was
investigated in the institute (Tamaskovics and Klapperich, 2010; Althoff, 2010; 2011;
Widodo, 2013). Multi-Stage Large Shear-Frame test has some advantages. During
installation of sample, soil in the shear box is highly compacted and then the tests are carried
out at low velocity. These test series contain multi-stage tests (ISV, IRV, and IPV) using a
cohesive soil (Canitz-silt) in combination with twelve different geogrids from six different
producers which have some properties such as short time tensile strength (25-180 kN/m),
elongation, aperture size, junction strength, surface and thickness. According to the grain size
distribution, Canitz-silt is strongly sandy and slightly clayey silt.
My research activity under DFG-research project FY 2012/2013 in this laboratory was still
ongoing. Multi-stage tests for ISV, IRV, and IPV have been undertaken not only for Canitz-
silt soil, but also Hohen Bockär glassy sand and Kaolin soil. In order to reach maximum
density for three kinds of soil, Proctor test has been done in which the measured water
contents are 11.3%, 2.5%, and 28.07% for Canitz-silt soil, Hohen Bockär glassy sand and
Kaolin soil, respectively. All types of geosynthetics were tested in machine direction, with
and without cross-element in order to understand behaviour of geosynthetics-soil-interaction.
Fig. 4-8 Layout of test setup (a) Top view and (b) Side view (after Wei-ping et al., 2007)
The sand used in the model tests was come from Qiantang River Beach, China. Specific
gravity G = 2.64, the coefficient of uniformity Cu = 2.5 and the coefficient of curvature Cc =
0.96. Sand grains were subangular and predominantly quartz with D10 = 0.1 mm, D60 = 0.25
mm, Dmax = 2 mm, emax = 0.89 and emin = 0,54. The unit weight of the fill in the tank was
15.35-15.83 kN/m3 and corresponded to a relative density of 557%. Peak secant angle of
shearing resistance max was 44o. Three types of geosynthetics were used in the model tests
with biaxial tensile strength 0.35, 1.40 and 22.5 kN/m at 8% axial strain.
Based on 15 model tests having a ratio of embankment height to clear spacing, h/s, from 0.7
to 2.0 they obtained some interesting results. In Test 1 through Test 4, when h/s < 1.4, the
surface of the embankment was non-uniform. It implies that the differential settlement
occurred on the top of an embankment. In Test 5 through Test 7, when h/s > 1.6, the
settlements at the base of the embankment were non-uniform, but the embankment surface
remained almost horizontal. Deformations in embankments of model Tests 1 through Test 7
also suggest that the height of equal settlement plane is about 1.4 - 1.6 times the cap beam
clear spacing, i.e., he = (1.4-1.6)s. Model tests from 1 to 7 were dedicated to unreinforcement
embankment whereas model tests from 8 to 13 were intended to reinforcement embankment
using geosynthetic. By using geosynthetics either low embankment h/s = 0.7 or hight
embankment h/s = 1.8 produces the higher stress concentration ratio, because the
geosynthetics transfers the vertical load over geosinthetics to beam cap. In addition, from the
experimental test that the higher embankment indicated the higher stress concentration ratio
when using the same of tensile strength of geosinthetics.
In this experiment, four cross-sectional 16×16 cm concrete piles were installed to support an
embankment. The unit weight of piles = 24 kN/m3, modulus of elasticity E = 23,982
MN/m2, strain emax = 0.9%, Poisson's ratio = 0.2, Some properties of soft soil (peat) used in
the experiment consist of unit weight of piles = 8 kN/m3, modulus of elasticity E = 0.85
MN/m2 at ' = 100 kPa, water content w = 300-350%, cohesion c = 8.5 kPa, internal friction
angle ' = 24o, compression index cc = 2.48, organic content OC = 80.2%, permeability ki =
4.1E-7 m/s.
Geogrid used in the research were Polyester GW 60 PET and Polyvinyl GW 180 PVA.
Tensile strength for longitudinal and transversal direction are 60 kN/m and 180 kN/m
respectively. Elastic modulus of geogrids are 850 kN/m at a strain of 2% for GW 60 PET
and 3,800 kN/m at strain of 1.25%.
82
Some interesting findings were obtained from the study. Settlements at surface of
embankment induced by cyclic loading with the same magnitude of the load was deeper than
static loading. Numerous number of geosynthetics layers were able to reduce settlement at
the surface of the embankment.
Highest tensile force occurred at an interface between the pile and geosynthetics, and then the
lower tensile force was located on between piles in the orthogonal direction, the lowest
tensile force was at the middle of geosynthetics between piles in a diagonal direction.
Moreover, maximum strain max for low embankment (h = 0.35 m) is higher than high
embankment (h = 0.70 m). Meanwhile, settlement at the surface of the embankment for low
embankment was deeper than the high embankment over the increasing external load.
The influence of cyclic and dynamic loading could decrease the arching effect. Moreover,
vertical stress in the embankment fill will increase over the large number of repetitions.
83
The applied product Secugrid 120/40 R6 has a Federal Railway Agency's certification for
this application. A 0/45 mm crushed mix acted as covering material. The geogrid was
installed in 10 layers with a distance of 0.5 m. Within just a few weeks, whole structure had
been completed, so that the flow of regular traffic was quickly reinstated.
The Secugrid solution with geogrid width of 4.75 m allowed faster and more cost-efficient
installation. The deformations were remedied as further layers of Secugrid were installed to
reinforce the embankment. Measurements were taken from the upper edge of the
embankment to determine the degree of deformation. Any deformations were hardly noted,
which confirmed that Secugrid had allowed an existing road to be successfully and safely be
widened on an extremely soft subgrade at favourable costs.
According to Djarwadi (2006) for internal stability using the limit equilibrium method, the
construction needs a tensile strength of geotextile around 87 kN/m, whereas foundation
stability (external stability) needs a tensile strength of geotextile about 272,8 kN/m.
Compaction has to be done in each embankment thickness of 0.3 m using vibratory smooth
drum compactor 15 ton weight until compaction degree of proctor standard ≥ 97%, in which
are 8 passing with a constant speed 10 km/h. The compaction tests using the sand cone test
are for each 2500 m³ volume of embankment. The observation for settlements of 400 m long
road embankments is depicted in Fig. 4-14. We can see that settlements at the final elevation
can be more than 1.0 m during less than 3 weeks of execution.
Fig. 4-16 and 4-17 show a typical GRS RWs having a FHR facing constructed in center
Tokyo, in Yamanote Line near Shinjuku station. It had been built during 1995-2000. This
new type GRS RWs has been constructed in more than 600 sites in Japan and the total wall
length is recently more than 100 km as of March 2008 (after Palmeira, 2008).
Fig. 4-16 GRS-RWs having a FHR facing in Yamanote Line, a) Typical cross-section
b) Wall under construction c) Completed wall, (after Palmeira et al., 2008)
86
Fig. 4-17 Staged construction of GRS RWs with a FHR facing
(after Palmeira et al., 2008)
87
b)
a)
c)
Fig. 4-18 Embankment at Barra da Tijuca, a) General scheme of piled embankment, b) Pile
caps above the initial fill, c) Pile caps inside the initial fill (after Almeida et al., 2008)
Some findings are gained from the field test. The strains of geogrids at the face of the pile cap
are higher than those at half distance between caps. Meanwhile, strains of geosynthetics in
between pile caps at an orthogonal direction were higher than those at a diagonal direction.
Based on the damaged geogrids, a circular cap could be more effective than the square caps
in thoses cases. In addition, the result suggested that subgrade reaction should not be
considered for designing a geogrids-piled embankment.
Fig. 4-19 Measured strains in geogrid in points at: a) face of the pile cap, b) half distance
between caps (after Almeida et al., 2008)
88
4.3.2.2. Gebeng Bypass Highway field test
Hassandi and Edil (2005) reported a full scale test embankment to evaluate the performance
of different types of loading platform (LTP) supported on geopiers. Three types of LTP were
constructed, they include reinforced concrete LTP, geosynthetic-reinforced LTP with two
layers of geogrid (catenary beam LTP), and geosynthetic-reinforced LTP with three or more
layers of geogrids (beam LTP).
To evaluate the performance of the catenary and the beam LTP. A well instrumented full-
scale test embankment was constructed in Gebeng, Pahang State, Malaysia as shown in
Figure 4-20. The length of the embankment was approximately 90 m, 13.5 wide, and 3.5 m
high. The test embankment was divided into 4 major sections (section 1 through 4) and 2
controlled sections (C1 and C2) at the two ends of embankments.
Fig. 4-20 Layout of the major and the control sections of test embankment
(after Hassandi et al., 2005)
Site investigation carried out indicated a silty clay layer as deep as 15 m at some locations.
However, this soft layer generally ends at approximately 5-6 m below the original ground
level. The soils are composed of highly plastic clay and silt with natural water content w
between 35% and 61%. The field vane tests that the shear strength Su lies between 14 to 60
kPa with most of the values less than 25 kPa and sensitivity S varies from 3 to 11.
Geopiers were 75 cm in diameter with 5.5 m deep installed at all sections. Section 1 has the
beam LTP with four layers of geogrids and 3.25 m geopier spacing designed by Collin
Method. Section 2 has the beam LTP with three layers of geogrids with 2.5 m geopier
spacing designed by the Collin Method. Section 3 has the catenary LTP with two layers of
high strength with 2.5 m geopier spacing design by BS 8006. Section 4 has a continuous
reinforced concrete slab as LTP. Settlement plates were set directly above the geopier
elements and at the center between groups of geopier elements. Earth pressure cells were
positioned on the geopiers and in between them. Piezometers and extenometers were placed
at different depths. Vertical inclinometers were placed at the toe of the embankment in
89
between the geopiers, and horizontal inclinometers immediately under and across the
embankment in between the geopiers.
Geogrid biaxial Tensar SS 20 was utilized for reinforcement in Section 1 and 2. The distance
between each geogrid layer in this section was 300 mm. In Section 1, aggregate blanket was
laid between layers of geogrid to a thickness of 1.5 m. Meanwhile in Section 2 was 1.0 m.
Geogrid uniaxial Miragrid 24 XT was used for reinforcement in Section 3. The distance
between the geogrid layers was 75 mm and the selected aggregate blanket thickness was 0.3
m.
According to Nordic Hanbook (2003) it recommends that the differential settlement between
the subsoil and the columns shall not be more than 0.1 m and 0.2 m. Total and differential
settlements of embankment test from different section was indicated in Table 4-3. We can
look that the use of piled embankment would reduce total settlement. Section 3 as catenary
LTP suggests the highest differential settlement and the lowest total settlement compared to
other sections.
Table 4-3 Total and differential settlements at the base of the embankments
Another finding resulted from the vertical inclinometer readings indicates that the use of LTP
supported on columns could reduce the lateral displacement of the subsoil at the edge of the
embankment.
Fig. 4-21 Scheme of the full-scale test Kyoto road (after Eekelen et al., 2009)
Comparison of two and a half years of measurements on a full-scale test compared with the
calculation based on German EBGEO and the British BS 8006 showed that EBGEO gives
better predictions of the load distributed in the piled embankment than BS 8006. The
dynamic load influences the arching and tensile stress in the geosynthetics. However, the
arching can be restored during the rest period.
Peat soil has a water content of 80 to 330% and an organic content between 25 and 80%.
Underneath these soil layers, slightly silty sand layers with a thickness up to 8 m are present,
which are medium dense packed. At the base of the sand layers, boulder clay is present,
which has a soft to stiff consistency and a water content of 10 to 20%.
The cemented columns (diameter 0.63 m), totally 3,260 MIP columns of a length between 5
and 8 m, were installed in square 1.5 × 1.5 m grids using the MIP-technique. Using a single
auger, a cement slurry is injected continuously into the soil during penetration as well as
during retrieval of the auger. On the top of the MIP-columns two layers of Fortrac R PVA
91
geogrid type M 400/30-30 were unrolled. To obtained an uniform bearing platform for the
ballast bed, 2.5 to 3% cement was added to the filling material. A gap graded gravel-sand
mixture (soil group SI according to the German Standard DIN 18196) with a coefficient of
uniformity ≥6 was used. The filling material was placed in layers of maximum 30 cm
thickness in accordance to the Ril 836. Each layer was compacted to a degree of compaction
at least 98%.
Fig. 4-22 Foundation system at section of Büchen-Hamburg (after Schwarz et al., 2005)
For laboratory testing purpose, wet grab samples were extracted from 4.5% of column every
500 m3 of treated soil, 6 unconfined compressive tests were carried out after 28 days, to
determine the unconfined compressive strength qu. According to the tests, unconfined
compressive strength after 28 days of all samples exceeded the design criteria of qu ≥ 2.2
MN/m2.
Fig. 4-23 Installation of MIP-columns and placing of geogrids (after Schwarz et al., 2005)
92
The settlement behaviour of tracks was monitored by means of geodetic measurements of the
outer rail of both tracks (see Fig. 4-24). The measurements were conducted in 3 measurement
sections each 20 m in length, consisting of 5 measurement points with spacing 5 m. The
effectiveness of the executed improvement measures was proved by means of settlement
measurements. The measurements show that the track Hamburg-Berlin has settled up 7 mm
in a period of 6 months after reopening the track. This settlement can be considered as small
since usually a settlement of 10 mm to 15 mm will occur, due to compaction of ballast bed,
the protective layer and embankment, even if the soil condition is favourable.
Final erection and driving of piles have been executed in mid November 2009. Rod and plat
settlement gauges above some of the pile caps with the aim of monitoring the settlement of
the pile group were installed.
4.4. Summary
Some important findings regarding with experimental works in laboratories and case studies
in the fields are as follows:
Inclusion of geosynthetics into granular material can reduce horizontal outward
pressure. The more amount of geosynthetic layers increases, lateral pressure will be
significantly decreased. It depends on properties of geosynthetics such as tensile
modulus, tensile strength, and size of aperture. Inclusion of geosynthetics in cohesive
soil exhibits a reduced cohesion of compound material.
In the field, in case of high embankment, geosynthetics can be applied using multi-
layer ranging from 30 cm to 50 cm between layers. Moreover, reinforcement in soft
soil using geosynthetics cannot overcome excessive settlement.
Properties of soft soil can be easily recognized such as high water content, low
undrained shear strength, low friction angle, low elasticity modulus and low
compressibility.
Piled embankment is a promising method to solve problem when constructing
infrastructure over soft soil. In laboratory, it is mostly modelled as end-bearing piles
that refers to shallow soft soil. Meanwhile, there is no a model has been derived for
floating piles.
Critical height of embankment is an important parameter in which differential
settlement on surface of embankment is equal to zero or very few. This parameter can
94
be applied to determine minimum thickness of embakment in order to avoid excessive
roughness of surface pavement.
Differential settlement at surface of pavement can be reduced by utilizing a high
tensile strength of geosynthetics. Particularly, in case of low embakment in which
effective high of embakment is lower than critical height.
There are various types of piles made of such as concrete pile, wooden pile, stone
column or pier aggregate, stabilized column which can be used in a construction with
or without pile cap. Because of these various types of piles, it will provide a little bit
different final result for a construction.
In many cases of geosynthetic-reinforced piled embankment, layer of geosynthetics
functionate as a load transfer platform (LTP) which is directly laid down at base of
embankment to transfer a load to piles. Multi-layer of geosynthetics can also be
applied.
Strains of geosynthetics at the face pile are higher than those at half distance between
piles. Meanwhile strains at half distance between piles are higher than those at center
of four piles.
In cases of floating piles over soft soil, total settlements at surface of embankment are
high enough and they are still occurred because of creep. Otherwise, differential
settlement at base of embankment is small enough. Furthermore, influence of piles to
reduce total settlement is really obvious.
95
CHAPTER 5
For a surface load of mIo/mr, the vertical displacement of the surface is given as follow
(Milovic, 1992):
(5-1)
(5-2)
(5-3)
(5-4)
If the elastic properties (E and ) are equal in the two layers, the coefficient of K and L are
equal to zero and the above equations reduce to Boussinesq’s equations. The main
assumption in layered elastic theory that the two-layer system is linear elastic and there is no
relative displacement at the interface between two layers (perfectly rough interface).
Based on elastic analysis, Fox (1948) provided a solution to the vertical stress z on the top
of the second layer for a perfectly rough interface and perfectly smooth interface. Fig.5-2
provides the vertical stress on the axis for the case with a/h=1. Here a= radius of the circular
footing, h= thickness of the first layer, d= depth, pz= the vertical pressure on the circular
footing, po= the pressure on the circular footing. The vertical ratio of rough interface is
0.644, 0.292, 0.081 and vertical ratio of smooth interface is 0.722, 0.305, 0.082 for E1/E2= 1,
10, 100 respectively.
Fig. 5-2 Vertical stress distribution at the surface of second layer for two-layered system
(after Fox, 1948)
5.1.2. Analytical Method for Dynamic Response of Beam and Plate on Winkler Type
Elastic Foundation under Moving Loads
First of all, consider an infinitely extended elastic beam-like plate strip (modeling the
pavement) on an elastic foundation (Winkler springs and daspot) under a constant
97
concentrated vertical load P moving with a constant velocity V. The equation of the free
motion of this system is written as (Thompson, 1963):
(5-5)
where w = w(x0,t) is the lateral deflection of plate strip, k and c is the foundation stiffness and
damping, respectively. D=E·h3/12(1-v2) is the flexural rigidity of the plate strip with E and v
being the modulus of elasticity and Poisson's ratio, respectively. Mass density and the
thickness of plate strip h, while x0 and z are the fixed coordinates along the length of the strip
and in the vertical direction, respectively, whilst t denotes time.
Now, it is introduced a new coordinate system x, z ,which moves with the load P. Thus, one
has
(5-6)
with the new coordinate system, the previous equation becomes
(5-7)
indicating that w in the moving coordinate system is independent of time after all the
transient vibrations have disappeared and the motion of the plate is said to be steady-state.
By introducing some symbols
into Eq. 5-7 above, we may rewrite the equation into the form
(5-8)
Therefore, the characteristic equation of the equation is
(5-9)
with a discriminant , which may be positive, negative or zero. For instance, at case > 0,
after application of boundary conditions (w=dw/dx=0) at x = ± ∞ , compatible of deflection,
slope and bending moment of the plate (pavement) under loading and formulation of the
discontinuity in the shear load. Finally, we can obtain:
Now, consider an elastic beam of infinite extent (modeling the pavement) on elastic
foundation (Winkler springs and dashpot) under a vertical distributed (over a finite length)
and time-dependent load p (x0,t) moving with a constant velocity V. The equation of motion
of this system in fixed Cartesian coordinates (x0,z) is (Kim and Roesset, 2003):
(5-11)
98
where v = v (x0,t) is the lateral beam deflection, I the cross-sectional moment of inertia, m the
mass per unit length of the beam, k the foundation stiffness per unit length and c is the
viscous damping constant. By using coordinate system (x,y), which moving load during a
time x = x0 – V·t, we may rewrite the Eq. 5-11 in another form:
(5-12)
When a distributed load p moves harmonically with time, i.e., p(x,t) = p0(x)eit, where i= √-1
and the operational load frequency, one has that v(x,t) = v0(x) eit , the Eq. 5-12 becomes:
(5-13)
For the general case of time variation of the load p, one can apply onto the double Fourier
transform with respect to x and t and then we obtain:
(5-14)
where:
Finally, consider the three-dimensional extension for an elastic plate of infinite extent on an
elastic foundation (Winkler springs and dashpots) under a vertical distributed (over a finite
surface) and time-dependent load p(x0,y0,t) moving with a constant velocity V along the x-
direction. The equation of motion of this system in fixed Cartesian coordinates (x0,y0,z) is
(Kim and Roesset, 2003):
(5-15)
where D is the flexural rigidity of plate, w is deflection along the z-direction and x0,y is the
middle plane of the plate.
5.1.3. Analytical Method for Dynamic Response of Layered Half-space under Moving Loads
Firstly, consider an elastic layered three-dimensional half-space (whose top layer can model
the pavement) under a concentrated vertical time-dependent load moving on its surface with
constant velocity V (Grundmann et al., 1999). For a homogeneous elastic layer, the equations
of motion may be written:
(5-16)
where ui (i=1,2,3) is the displacement vector, and are the Lame elastic constants, is the
mass density, commas and overdots denote space and time differentiation, respectively. Axes
x and y corresponding to i=1 and i=2 denote the two horizontal directions, while axis z for i=3
denotes the vertical direction. The Eq. 5-16 above can be expressed using the aid of
Helmholtz’s decomposition as:
(5-17)
where the scalar and vector i functions satisfy the wave equations
(5-18)
with ijk being the alternating tensor and cp = √(l+2)/ and cs=√/ the dilatational and
shear wave velocities.
99
5.1.4. Critical Velocity
The velocity of high-speed train can approach or exceed the characteristic wave velocity of
the dynamic system comprising of the underlying soft ground, trackbed/embankment, and the
moving load. As the train’s velocity approaches some ’critical velocity’ large deformation
can occur. These motions could be dangerous for train and the integrity of the structure. The
velocity of Rayleigh surface in soft sandy soils may be vary low (90-130 m/s) even such
wave velocity in some soft soils in Netherlands can be as low as 29-47 m/s, which is
considerably lower than the railway's proposed design speed, the matter clearly required
investigation.
Kenney (1954) gave some insight into parameters that seem to be of importance to the critical
velocity, vcr, as the equation below:
(5-19)
100
(a) (b)
Fig. 5-3 Cyclic loading (a) Procedure of an explicit calculation of accumulation
(b) Evolution of total strain in a cyclic triaxial test (after Wichtmann, 2005)
The accumulation starting from the second cycle is calculated directly by means of an
equation of the shape:
(5-20)
where: : Jaumann stress rate, D: strain rate, Dacc: given accumulation rate, E: elastic
stiffness, without following the strain path during the particular cycles. The equation leads to
an accumulation of stress (e.g. = -E: Dacc at D=0) and/or strain (D= Dacc at =0). When =0
the strain follows the average accumulation curve acc(N).
The strain amplitude ampl is assumed constant for the explicit calculation. The explicit
calculation may be interrupted after definite numbers of cycles and ampl can be updated in
the an implicit so-called control cycle.
(5-23)
and the material constant C1 and C2. The tensor ampl contains the amplitudes of the particular
strain component. The latter transformation in Eq. (5-23) is valid for the case of cyclic simple
shear tests with a constant shear strain amplitude ampl. The compaction rate is
obtained from Eq. 5-24:
(5-24)
101
The model of Sawicki & Swidzinski was performed on tests with a relatively low number of
cycles (N < 103). A major drawback of the model is that only the volumetric accumulation,
whereas the deviatoric one not be considered.
The model of Bouckovalas et al. (1984) predicts the cyclic flow rule correctly: vacc vanishes
on the critical state line, while qacc becomes zero at av = 0. A power law is used for the
dependence of the accumulation rates on the number of cycles. In the model, a constant a = 3
and c = -1.5 are chosen and factor f remains vague in the mathematical definition. The model
uses a state variable for the historiotropy which considers also the amplitude of the cycles.
The model gives a too strong amplitude-dependence with approximate value acc ~ (ampl)3.
102
theory is so-called as associated plasticity. But, for Mohr-Coulomb yield function, this theory
will give overestimate dilatancy prediction. Therefore, besides yield function f, it is used a
plastic potential function g. At case g ≠ f, it is so-called as non-associated plasticity.
Generally, strain rate can be written as:
(5-30)
where is a plastic multiplied factor. For pure elastic, is equal to zero and for plastic
behaviour, is a positive value.
(5-31)
(5-32)
Fig. 5-4 Basic principle of elastic perfectly plasticity (after Hill, 1958)
The equations can be used to obtain a relationship between effective stress rate and strain rate
for elastoplastic model (Smith & Griffith, 1982; Vermeer, 1982):
(5-33)
with:
Parameter can be used as a switch. If material behaviour is elastic, value is equal to zero,
whilst value is one when material behaviour is plastic. Plastic theory above is limited only
for a continuous and smooth yield surface, and not included into a multi-surface of yield likes
MC model. For multi-yield surface, Koiter (1960) took into consider to flow vertices
involving two or more plastic potential functions:
(5-34)
103
(5-35d)
(5-35e)
(5-35f)
two parameters of plastic model emerge into yield function are the friction angle and
cohesion c well known in geotechnical engineering. These yield functions together will form
a hexagonal cone in main stress space as shown in Fig. 5-5.
Fig. 5-5 Mohr-Coulomb yield surface in main stress space (c=0) (after Smith et al., 1982)
Besides yield functions, there are six plastic potential functions for MC model:
(5-36a)
(5-36b)
(5-36c)
(5-36d)
(5-36e)
(5-36f)
Plastic potential function has a third plastic parameter, namely dilation angle . This
parameter is needed to model an increment of plastic volumetric strain positively as it
actually occurs on stiff soil.
For c > 0, Mohr-Coulomb criteria allow to tensile stress, but in fact, soil only reminds this
stress very small, even nothing. Plaxis software can model by using constraint of tensile
stress. This constraint incorporates with three additional yield function:
(5-37a)
(5-37b)
(5-37c)
The allowable tensile stress t can be setup (default) as zero. For three yield functions, these
can be used with an associated flow rule. The Mohr-Coulomb model needs five parameters.
Generally, these parameters can be obtained from laboratory tests as follows:
E : Young‘s modulus [kN/m2]
104
v : Poisson's ratio [-]
: friction angle [o]
c : cohesion [kN/m2]
: dilation angle [o]
(5-38)
The ultimate deviatoric stress qf and quantity qa in Eq.5-38 are defined as:
(5-39)
The above relationship for q is derived from the Mohr-Coulomb failure criteria, which
involves the strength parameters c and p. As soon as q = qf , the failure criterion is satisfied
and perfectly plastic yielding occurs. The ratio between qf and qa is given by the failure ratio
Rf , which should obviously be smaller than 1. Rf = 0,9 often is a suitable default setting.
(5-40)
E50ref is a reference stiffness modulus corresponding to the reference stress pref. The actual
stiffness depends on the minor principal stress 3', which is the effective confining pressure
in a triaxial test. The amount of stress dependency is given by the power m. Von Soos (2001)
showed that m values range between 0.4 and 1.0. In order to simulate a logarithmic stress
dependency, as observed for soft clay, the power should be taken equal to 1.0 and for m value
sands is 0.5. As a secant modulus E50ref is determined from a triaxial stress-strain-curve for a
mobilization of 50% of the maximum shear strength qf.
(5-41)
where Eurref is the reference Young's modulus for unloading and reloading, corresponding to
the reference pressure ref. The un-/reloading path is modeled as purely (non-linear) elastic.
The elastic components of strain e are calculated according to a Hookean type of elastic
relation.
(5-42)
For drained triaxial test stress paths with 2 = 3 = constant, the elastic Young's modulus Eur
remains constant and the elastic strain are given by equations:
(5-43)
(5-44)
(5-45)
Fig. 5-7 Successive yield loci for various values of the hardening
parameter p and failure surface (after Schanz et al., 1999)
In contrast to elastic perfectly Mohr-Coulomb (MC) model, in the Hard Soil (HS) model,
plastic strains may occur before the limit MC-failure stress reached. The HS model
incorporates two other yield surfaces, which are not fixed in principal stress space, but they
may expand and soil hardening is simulated due to plastic straining. As illustrated in Fig. 5-
8, distinction is made between two types of hardening, which are shear hardening and
compression hardening. For the shear hardening law, a yield function fs is introduced, which
is a function of the triaxial loading stiffness E50 and for the compression hardening a yield
function fc is formulated.
Fig. 5-8 Yield surface of HS model, a) Successive yield loci for shear hardening and
compression hardeing in p-q space, b) Total yield contour in principal stress space
(after Schanz et al., 1999)
107
where cv is the critical state friction angle, being a material constant independent of density
(Schanz & Vermeer, 1996), and m is the mobilized friction angle:
(5-49)
The above equations correspond to the well-known stress-dilatancy theory (Rowe 1962,
Rowe, 1971). The essential property of the stress-dilatancy theory is that the material
contracts for small stress ratio m < cv, whilst dilatancy occurs for high stress ratio m < cv.
At failure, when the mobilized friction angle is equal to the failure angle, p, it is found the
previous equation, that:
(5-50)
Hence, the critical state angle can be computed from failure angle p and p . The definition
of flow rule is equivalent to the definition of the plastic potential functions g12 and g13
according to:
(5-51a)
(5-51b)
By using the Koiter-rule (Koiter 1960) for yielding depending on two yield surface (Multi-
surface plasticity) one finds:
(5-52)
In fact, maximum soil stiffness is observed at very low strain levels, e.g. Strains smaller than
10-5 (Atkinson and Sällfors, 1991). Special devices is needed to identify stiffness at very
small strain. Biarez and Hicher (1994) gave a simple correlation for quarts sand as follows:
(5-53)
where e is the void ratio of the soil and p is the mean stress. Alphan (1970) also estimated
preliminary estimation of the E0. Hardin and Drnevich (1972) formulated the decay of
stiffness when strains increase.
(5-54)
where E is the actual secant modulus at the corresponding shear strain , E0 is the initial
stiffness of soil and 0.7 is the shear strain at 70 percent from the initial stiffness E0 .
The initial shear modulus G0 is determined from the relationship between E0 and Poisson's
ratio v as in the following equation
(5-55)
(5-56)
The stiffness degradation due to plastic straining is modelled by involving material
hardening. Therefore, before reaching plastic material behaviour, the formulation of small
strain stiffness is cut off at the unloading-reloading stiffness Eur. There are two additional
input parameters are required for the HS-Small model: the elastic small-strain shear modulus
G0ref at reference pressure pref and the curve-decay value 0.7 in primary loading.
(5-57)
where c is the strain up to the end of consolidation, t is the time measured from the
beginning of loading, tc is the time to end of primary consolidation and CB is a material
constant (see Fig. 5-10). Rewriting above the equation as:
(5-58)
(5-59)
where: C = CB (1 + eo)
Butterfield (1979) alsodescribed a slightly different possibility for secondary compression
as:
(5-60)
Originally, Hencky (1928) used the subscript '0' denoting initial values whilst the superscript
'H' for denotating logaritmic strain. For small strains it is possible to show that:
(5-62)
110
1D Creep model
The SSC model described here is on the basis of the work done by Vermeer and Neher
(1999) on elastic visco-plastic creep model. Vermeer and Neher (1999) adopted Bjerrum's
idea to find an analytical expression for quantity of c. For the one dimensional creep, two
strain components need to be modelled as shown in Fig. 5-11.
Fig. 5-11 Standard Oedometer Test, a) Stepwise loading in e-log vs. ’ plot
b) Void ratio vs. Time (After Vermeer and Neher, 1999)
First of them is the more or less elastic deformation, as directly observed in unloading and
reloading condition. The other component of strain is irreversible and time dependent.
Volumetric strain implies a change of void ratio and it is convenient to formulate the
deformation in terms of void ratio e, and the change of void can be expressed by a equation
below:
(5-63)
where the superscripts e and c refer to the elastic and creep component respectively. The
elastic change of void ratio is formulated as follows:
(5-64)
where Cs is the swelling index, which can be the unloading-reloading index Cur. Whilst the
creep deformation is represented using power law
(5-65)
where is a particular reference time, which can mostly be taken equal to one day. Cis the
secondary compression index that is also referred to as the creep index and Cc is the
compression index obtained from an oedometer test. From the equation above that the creep
rate depends on the OCR value. Some typical soil data give Cs= Cc/10 and C = Cc/30. This
will give of about 27.
The preconsolidation stress p increases during creep according to the differential equation
(5-66)
(5-67)
111
where p0 is the initial preconsolidation stress for ec = e0c. Subscript ‘0’ denotes the initial
value. Finally, the creep rate formulation can be obtained as the equation below:
(5-68)
The effective stress ’ may be either larger or smaller than p0. For a special case of a
constant effective stress, the differential creep formulation can be integrated analytically to
obtain the logarithmic creep law
(5-69)
where t = 0 for e = e0, and then
(5-70)
(5-71)
with
Fig. 5-12 shows that the stress measure peq is constant on the ellipse in the p-q plane. In fact,
the ellipses are from the Modified Camclay Model as introduced by Roscoe and Burland
(1968).
Fig. 5-12 Diagram of peq-ellipse in a p-q plane (after Vermeer and Neher, 1999)
Soil parameter M represents the slope of the so-called ’critical state line’ and it is defined as:
(5-72)
where cv is the constant volume friction angle, also referred to as critical-state friction angle.
The preconsolidation pressure changes during creep according to the law
(5-73)
where * and * are a modified compression index and a modified swelling index
respectively. In case of small strain, it gives
112
(5-74)
and
(5-75)
(5-76)
(5-77)
(5-78)
Parameter used in Soft Soil Creep model has a relationship with Camclay parameters model.
The modified compression index * is / (1+e) , the modified swelling index * is / (1+e) ,
and the modified creep index * is C / 2.3(1+e). The ratio of the unloading/primary loading
stress, * / *, cannot be smaller than 1 and should normally be between 2 and 10. Users
should be very wary of values outside this range' for most practical cases the value falls
within the range of 3 to 7. Secondly, there is the creep ratio, (*-*)/*, to consider. This
ratio can have a wide range of values, normally between 5 and 25, where high values
represent stiff soil with little creep and small values represent soft soils with a considerable
amount of creep. For most practical cases the ratio falls within the range of 10 to 20 and if the
creep ratio is over 25 one could reconsider the use of the creep model.
Table 5-1 Material parameters for the Soft Soil Creep model
(after Vermeer and Neher, 1999)
Symbols Name of parameters
* a modified compression index
* = * / 2-10 a modified swelling index
* = (*-*) / 10-20 a modified creep index
vur = 0.15 un-/reloading Poisson's ratio
c' effective cohesion
' effective friction angle
K0NC height of normal consolidation surface
OCR the state of preconsolidation stress
Fig. 5-13 Influence of n (a) and hs (b) on Oedometric curves (after Herle and Gudehus, 1999)
(5-81)
where mean stress ps1 and ps2 may be calculated from axial stresses using the Jaky formula K0
= 1 – sin c, and ep1 and ep2 are the void ratios corresponding to the stress ps1 and ps2. Tangent
compression indices corresponding to the limit values of interval ps1 and ps2 (Cc1 and Cc2) can
be approximated by secant moduli between loading steps preceding and following steps ps1
and ps2. The parameter hs can be expressed as:
(5-82)
where Cc is a secant compression index calculated from limiting values of the calibration
interval ps1 and ps2, in which ps and ep are averages of the limit values of p and e in this range.
Further model parameters are the reference void ratio ed0, ec0 and ei0, corresponding to the
densest, critical state and loosest particle packing at the zero mean stress. Bauer (1996)
114
formulated the reference void ratio ed, ec and ei corresponding to the non-zero stress
depending on the mean stress, as the Eq. 5-83 and then shown in Fig. 5-14.
(5-83)
Fig. 5-14 Dependency of the reference void ratio of ed0, ec0 and ei0 on the mean stress
(after Herle and Gudehus, 1999)
The initial void ratio emax can be considered equal to the critical state void ratio at zero
presure ec0. Void ratio ed0 and ei0 can approximately be obtained from empirical relations.
The physical meaning of ed0 is the void ratio at maximum density, void ratio ei0 represents
intercept of the isotropic normal compression line with p = 0 axis. The ratio ei0 /ec0 ≈ 1.2 was
derived from considering skeleton consisting of ideal spherical particles (Herle and Gudehus,
1999). Whereas the minimum void ratio ed0 should be obtained by densification of a granular
material by means of cyclic shearing with small amplitude under constant pressure. If there
are no data, it can be estimated using an empirical relation ed0 /ec0 ≈ 0.4.
The last two parameters and should be calibrated by means of single-element simulations
of the drained triaxial tests. Parameter controls independently the shear stiffness and
controls the peak friction angle.
Fig. 5-15 Geometrical idealizations of one cell piled embankment (a) Axisymmetric
idealization (b) Idealization for plane strain after Bergado (c) Idealization plane strain using
equivalent stiffness (d) 3D geometry (after Satibi, 2009)
For Axisymmetric, the three dimensional one cell of the piled embankment is transformed in
to a circular cell using the area of the pile and the soil the same. Fig. 5-15(a) shows the
transformation of the squared cell to a circular cell which one radian of the circular piled
embankment cell is used.
Plane strain after Bergado: The three dimensional grid of piles can be transformed into
continuous walls with an equivalent thickness teq in plane strain model as indicated in Fig. 5-
15(b). By keeping an improved area ratio (Ac/AE) is constant, the thickness of the continuous
116
wall is calculated based on the consideration of this ratio. Due to symmetrical geometry, only
a half of the plane strain geometry is used for the FE calculations.
(5-85)
Plane strain with equivalent pile stiffness: An alternative method to transform the three
dimensional grid of piles into a continuous wall in plane strain condition is by assuming
equivalent wall stiffness. The equivalent stiffness of the wall Eeq is taken as the proportional
average of the pile and soil stiffness. Hence,
(5-86)
where Ec, Es and Aw are pile stiffness, soft soil stiffness and wall area as illustrated in Fig.5-
15(c). It is worthy to note that when using this approach, the improved area ratio becomes
larger.
3D geometry: Three dimensional dimension of the actual case can be best analyzed using 3D
geometry. Here, half of the piled embankment cell is considered the 3D FE analysis as
described in Fig. 5-15(d).
According to Irsyam et al. (2008) he modeled bamboo piles in the Bamboo piles-Mattress
system as soil reinforcement for embankment on soft clay, and results of FE-analysis was in a
good agreement with field measurement. To model pile he applies ‘node-to-node anchor’ in
which it was commonly used and available in Plaxis. By modelling pile using elasto-plastic
‘node-to-node anchor’, soil in between wooden piles is not confined and then actual vertical
pile capacity can be obtained. The stiffness of elasto-plastic ‘node-to-node anchor’ is taken
from the estimated displacement of wooden piles in the vertical direction.
5.4. Summary
118
CHAPTER 6
6.1. Introduction
This chapter discusses some findings resulted from some experimental works and several
field case studies. Most of the results are directly obtained from observations through devices
installed both in the laboratories and in the field. By using FE-calculations, the results are
then verified with the main aim to do comparisons in order to be able to properly understand
the behavior of material.
D = 50 cm
Fig. 6-1 Schematic work procedure for FE analysis
119
Strength of granular material for the base course or subbase course can be expressed by CBR
(California Bearing Ratio) or resilient modulus (Mr). Table 6-1 provides a correlation
between CBR and Mr using Eq. 3-2 through Eq. 3-5, as explained in Chapter 3.
Table 6-1 Resilient moduli of base course and subbase course
Resilient moduli (Mr)
Source(s) Equations
CBR 40 CBR100 CBR150
Heukelom&Klomp Mr (psi) = 1500 CBR 60 ksi 150 ksi 225 ksi
(1962) Mr (MPa) = 10.34 CBR 413 MPa 1034 MPa 1551MPa
Webb&Campbell Mr (psi) = 3116 CBR 0.4779707 18 ksi 28 ksi 34 ksi
(1986) Mr (MPa) = 21.485 CBR 0.4779707 125 MPa 194 MPa 235 MPa
Sasongko (1996) Mr (MPa) = 10 CBR 400 MPa 1000 MPa 1500 MPa
Mr (MPa) = 4 CBR 160 MPa 400 MPa 600 MPa
NCHRP (2004) Mr (psi) = 2555 CBR 0.64 27 ksi 48 ksi 63 ksi
Mr (MPa) = 17.6 CBR 0.64 186 MPa 335 MPa 435 MPa
According to Sasongko (1996), resilient modulus correlates to 4 times CBR value when
CBR>28%. This is quite close with NCHRP’s equation. Whilst Heukelom’s equation
extensively developed for fine grained soil with maximum CBR of 10% provides too high
estimation, whereas Webb’s equation gives too small estimation. Based on data in Table 6.1,
the appropriate values of subbase course with CBR of 40% correlates to the resilient modulus
between 160 MPa and 190 MPa. Meanwhile, base course with CBR of 100% correlates to the
resilient modulus between 335 MPa and 400 MPa and base course very dense with CBR of
150% correlates to the modulus resilient between 435 MPa and 600 MPa.
In this FE analysis, three resilient moduli are used for representing resilient moduli of
pavement materials, namely 200 MPa, 400 MPa and 600 MPa respectively. Mohr-Coulomb
(MC) and the Hardening Soil (HS) model are used for analyzing this work. Besides elastic
modulus represented by resilient modulus, other parameters are also introduced for FA-
analysis, namely friction angle () and dilation angle (ψ). Dilation angle is an important
parameter because pavement material is a compacted soil. Material properties are shown in
Table 6-2.
Table 6-2 Material properties used in finite element analysis
Element and model
Properties Unit Granular material Geogrid
MC HS Elastic
Unit weight kN/m3 18 18 -
Friction angle [o] 38, 45 38, 45 -
Dilation angle [o] 8, 15 8, 15 -
Cohesion, c kPa 0.5 0.5 -
Young’s modulus, E MPa 200, 600 - -
Secant modulus MPa - 50, 100 , 150 -
Oedometer comp. modulusoed MPa - 50, 100 , 150 -
Un-/reloading modulusur MPa - 200, 400, 600 -
Poisson’s ratio, vur /v [-] -/0.3 0.2/- -
Tensile strength kN/m - - 30, 300
Drainage Drained Drained Non-porous
Calculation type Stage construction, plastic calculation
Pictures in the Table 6-3 through Table 6-6 show the patterns of failure which indicate zones
having the maximum horizontal displacements (red colour) when being subjected to a load
120
using MC-model for the elastic modulus of 200 MPa and 600 MPa. In Table 6-3 and Table 6-
4, some findings are resulted from FE-analysis using MC model, which the elastic modulus
of granular soil of 200 MPa and 600 MPa with tensile strength of geosynthetics of 30 kN/m.
It is clearly shown that granular material with high friction angle of 45o and dilation angle of
15o is able to support higher stress (around 37 kPa) than lower friction angle and dilation
angle, 38o and 8o respectively, which can support the vertical stress around 34 kPa. Moreover,
by using material with a higher friction angle and dilation angle can reduce horizontal
displacement smaller than lower one. Another important finding is that the number of
inclusion of geosynthetics in granular material contributes to reduce horizontal displacement.
Table 6-3 Horizontal displacements for MC-model, E=200 MPa and EA= 30 kN/m
Elastic modulus MC-model for E = 200 MPa , EA= 30 kN/m
Friction angle and
dilatant angle = 45 , = 15o
o
= 38o, = 8o
Number of
geosynthtetics layer
0 1 2 3 0 1 2 3
Patterns of horizontal
displacements
Table 6-4 Horizontal displacements for MC-model, E=600 MPa and EA= 30 kN/m
Elastic modulus MC-model for E = 600 MPa , EA= 30 kN/m
Friction angle and
dilatant angle = 45 , = 15o
o
= 38o, = 8o
Number of
geosynthtetics layer
0 1 2 3 0 1 2 3
Patterns of horizontal
displacements
121
Similar trend in previous findings, by inclusion of higher tensile strength of 300 kN/m in
granular material with elastic modulus of 200 kPa, it is able to support the vertical stress
enormously, particularly when granular material has the higher friction angle and dilation
angle, as indicated in Table 6-5. However, inclusion tensile strength of 300 kN/m in granular
material with an elastic modulus of 600 MPa only contributes slightly vertical stress, as
shown in Table 6-6.
Table 6-5 Horizontal displacements for MC-model, E=200 MPa and EA= 300 kN/m
Elastic modulus MC-model for E = 200 MPa , EA= 300 kN/m
Friction angle and
dilation angle = 45 , = 15o
o
= 38o, = 8o
Number of
geosynthtetics layer
0 1 2 3 0 1 2 3
Patterns of horizontal
displacements
Table 6-6 Horizontal displacements for MC-model, E=600 MPa and EA= 300 kN/m
Elastic modulus MC-model for E = 600 MPa , EA= 300 kN/m
Friction angle and
dilation angle = 45 , = 15o
o
= 38o, = 8o
Number of
geosynthtetics layer
0 1 2 3 0 1 2 3
Patterns of horizontal
displacements
We may resume some findings as presented in Table 6-3 through Table 6-6 that smaller
horizontal displacement can be obtained when we apply higher elastic modulus, higher
friction angle and dilation angle, higher tensile strength and number of geosynthetics.
122
Table 6-7 through Table 6-12 show some important results using HS-model for granular
material with an elastic modulus of 200, 400 and 600 MPa respectively. We may notice
results as shown in Table 6-7 through Table 6-9 that the higher values of elastic moduli,
friction angle and dilation angle provide higher vertical stress. When we apply the elastic
modulus of 200 MPa with tensile strength of geosynthetics around 30 kN/m, at least 3 layers
of geosynthetics is needed for granular materials having friction angle of 45o and dilation
angle of 15o whereas using lower friction angle and dilation angle (38o and 8o respectively),
sample soil fails even applied for 4 layers of geosynthetics.
Table 6-7 Horizontal displacements for HS-model, E=200 MPa and EA= 30 kN/m
Elastic modulus HS-model for E = 200 MPa , EA= 30 kN/m
Friction angle and
dilation angle = 45 , = 15o
o
= 38o, = 8o
Number of
geosynthtetics layer
0 1 2 3 0 1 2 3
Patterns of horizontal
displacements
Remarks Soil body Soil body More More Soil body Soil body Soil body Soil body
collapses collapses load step load step collapses collapses collapses collapses
Horizontal displacement 0.720 0.753 0.845 0.858 0.521 0.544 0.589 0.634
(mm)
Vertical displacement 1.228 1.303 1.375 1.422 1.111 1.121 1.268 1.250
(mm)
Force (kN) 0.633 0.634 0.637 0.64 0.634 0.634 0.633 0.636
Stress (kPa) 40.302 39.578 39.325 39.491 35.846 34.051 35.491 36.315
Table 6-8 Horizontal displacements for HS-model, E=400 MPa and EA= 30 kN/m
Elastic modulus HS-model for E = 400 MPa , EA= 30 kN/m
Friction angle and
dilation angle = 45 , = 15o
o
= 38o, = 8o
Number of
geosynthtetics layer
0 1 2 3 0 1 2 3
Patterns of horizontal
displacements
Remarks More More More More Soil body More More More
load step load step load step load step collapses load step load step load step
Horizontal displacement 0.339 0.323 0.361 0.348 0.357 0.330 0.346 0.324
(mm)
Vertical displacement 0.693 0.698 0.682 0.699 0.709 0.643 0.685 0.637
(mm)
Force (kN) 0.726 0.735 0.743 0.748 0.062 0.650 0.649 0.659
Stress (kPa) 39.614 39.925 40.427 42.273 30.267 35.026 36.017 35.798
123
At least a layer of geosynthetics has to be used when the elastic modulus of sample soil is 400
MPa with a friction angle of 38o and a dilation angle of 8o as shown in Table 6-8. Another
important finding is that use of the higher friction angle and dilation angle can slightly reduce
the horizontal displacement.
Table 6-9 Horizontal displacements for HS-model, E=600 MPa and EA= 30 kN/m
Elastic modulus HS-model for E = 600 MPa , EA= 30 kN/m
Friction angle and
dilation angle = 45 , = 15o
o
= 38o, = 8o
Number of
geosynthtetics layer
0 1 2 3 0 1 2 3
Patterns of horizontal
displacements
Table 6-10 through Table 6-12 are the results using HS-model for granular materials with
elastic modulus of 200, 400 and 600 MPa respectively and inclusion of geosynthetics with
tensile strength of 300 kN/m. Similar result as previous findings that higher values of elastic
modulus, friction angle, dilation angle and number of layers are able to reduce the horizontal
displacement and simultaneously increase the vertical stress.
Table 6-10 Horizontal displacements for HS-model, E=200 MPa and EA= 300 kN/m
Elastic modulus HS-model for E = 200 MPa , EA= 300 kN/m
Friction angle and
dilation angle = 45 , = 15o
o
= 38o, = 8o
Number of
geosynthtetics layer
0 1 2 3 0 1 2 3
Patterns of horizontal
displacements
Remarks Soil body Soil body More More Soil body Soil body More More
collapses collapses load step load step collapses collapses load step load step
Horizontal displacement 0.780 0.931 1.09 0.995 0.521 0.689 0.744 0.707
(mm)
Vertical displacement 1.549 1.756 1.882 1.914 1.111 1.267 1.459 1.451
(mm)
Force (kN) 0.630 0.654 0.722 0.817 0.634 0.635 0.678 0.725
Stress (kPa) 39.048 46.547 41.453 43.401 35.846 36.987 35.500 35.745
124
Table 6-11 Horizontal displacements for HS-model, E=400 MPa and EA= 300 kN/m
Elastic modulus HS-model for E = 400 MPa , EA= 300 kN/m
Friction angle and
dilation angle = 45 , = 15o
o
= 38o, = 8o
Number of
geosynthtetics layer
0 1 2 3 0 1 2 3
Patterns of horizontal
displacements
Remarks More More More More Soil body More More More
load step load step load step load step collapses load step load step load step
Horizontal displacement 0.337 0.364 0.355 0.340 0.357 0.326 0.308 0.278
(mm)
Vertical displacement 0.690 0.726 0.687 0.729 0.709 0.641 0.641 0.639
(mm)
Force (kN) 0.726 0.776 0.826 0.868 0.622 0.677 0.694 0.718
Stress (kPa) 39.551 44.254 47.231 50.274 30.267 37.434 33.885 35.517
Table 6-12 Horizontal displacements for HS-model, E=600 MPa and EA= 300 kN/m
Elastic modulus HS-model for E = 600 MPa , EA= 300 kN/m
Friction angle and
dilation angle = 45 , = 15o
o
= 38o, = 8o
Number of
geosynthtetics layer
0 1 2 3 0 1 2 4
Patterns of horizontal
displacements
125
1.2
0.8
Horizontal displacement (mm)
0.6
EA 30 kN/m Phi 45-Psi 15- E 600
0.4
EA 300 kN/m Phi 38-Psi 8- E 400
Fig. 6-2 suggests that at least two geosynthetics layer (vertical distance of 36.7 cm) of 300
kN/m tensile strengths on the elastic modulus of 200 MPa pavement material are needed, if
not soil sample will be collapsed. Whereas use of low tensile strength of 30 kN/m does
contribute stability of a soil sample. Meanwhile, at least a layer of geosynthetics (vertical
distance of geosynthetics around 55 cm) is needed when applying the elastic modulus of 400
MPa with friction angle of 38o and dilation angle of 8o. Another finding using material with
the elastic modulus of 600 MPa without the geosynthetics layer (soil sample height of 110
cm), soil sample does not fail. It is worthly to mention here that the best result is the use of
300 kN/m-tensile stiffness with the soil stiffness, phi, psi are 400 MPa, 38 degree, 8 degree
respectively. The reason is that the increased number of geosynthetics layers applied causes
horizontal displacements is smaller.
No doubt, the base course or subbase course material must be compacted on a construction.
Hence, this material has a dilation angle and also the material should have a friction angle at
least 30o. Fig. 6-3 presents some findings using MC-model. Generally, the higher values of
friction angle, dilation angle, tensile strength of geosynthetics and number of geosynthetics
layers will provide the higher vertical stress as shown in Fig. 6-3. In x-axis, the numbers for
126
the first, second and third are the friction angle, dilation angle, and number of geosynthetics
layers respectively.
60
(a) 30
20 EA=30 kN/m
EA=300 kN/m
10
0
38-8-0 38-8-1 38-8-2 38-8-3 45-15-0 45-15-1 45-15-2 45-15-3
Friction angle, Dilation angle, Number of layers
41
Max. stress achieved (kPa)
40
39
38 MC-model
37 E = 600 MPa
36
35
(b) 34
33 EA=30 kN/m
32
EA=300 kN/m
31
30
38-8-0 38-8-1 38-8-2 38-8-3 45-15-0 45-15-1 45-15-2 45-15-3
Friction angle, Dilation angle, Number of layers
Fig. 6-3 Maximum stress using MC-model at different elastic modulus of base course
(a) E=200 MPa (b) E= 600 MPa
Fig. 6-4(a) through Fig. 6-4(c) reveal some findings using HS-model as given in Table 6-7
through 6-12. When using MC-model, the model did not show stability of sample during
subjected by a load. Vertical stresses are influenced by the elastic modulus of granular
material, friction angle, dilation angle, tensile strength of geosynthetics and number of layer
inserted in granular soil as shown in Fig. 6-4. In x-axis, number of first, second and third is
the friction angle, dilation angle and number of geosynthetics layers respectively.
50
Max. stress achieved (kPa)
45
40
35
30
25
(a) 20
EA=30 kN/m
EA=300 kN/m
15 HS-model
10 E = 200 MPa Soil body collapses
5
0
38-8-0 38-8-1 38-8-2 38-8-3 45-15-0 45-15-1 45-15-2 45-15-3
Friction angle, Dilation angle, Number of layers
127
45
60
Max. stress achieved (kPa)
50
40
30
(c) EA=30 kN/m
20 HS-model
EA=300 kN/m E = 600 MPa
10
0
38-8-0 38-8-1 38-8-2 38-8-3 45-15-0 45-15-1 45-15-2 45-15-3
Friction angle, Dilation angle, Number of layers
Fig. 6-4 Maximum stress using HS-model at different elastic modulus of base course
Circles signs designated in Fig. 6-4 are the stresses achieved when the soil sample is
subjected to a vertical load. It can be seen that the higher friction angles can bear the higher
stresses as well as higher tensile strength of geosynthetics.Though it is not so obvious for the
elastic modulus of 400 MPa when geosynthetics with tensile strength of 300 kN/m is applied.
No doubt, the higher elastic modulus suggests higher stability when being subjected to a load.
Fig. 6-4 shows that the elastic modulus of 200 MPa needs at least 2 layers of geosynthetics
(36.7 cm-vertical distance of layers which the soil sample is 110 cm high) to withstand the
loading. Whilst, materials having the elastic moduli of 400 MPa with a friction angle of 38o
needs at least one layer of geosynthetics (55 cm vertical distance). Granular soils having the
elastic modulus of 600 MPa can withstand a load with vertical stress around of 35 kPa
without using a layer of geosynthetics.
35
ST2
ST3
30 ST4
ST2a
Vertical pressure (kPa)
25 ST3a
ST4a
20 ST1
MC top, phi 34 deg.
MC bot. phi 34 deg.
15 MC top, phi 44 deg.
MC bot. phi 44 deg.
10
0
0 200 400 600 800 1000 1200
Embankment height (mm)
(a) Vertical stress resulted from Test 7 and model (b) Vertical stress using FEM
Fig. 6-5 Influence of embankment height due to vertical stress
during embankment filling of Test 7
During the consolidation process, in this model, imitated by discharging water gradually from
the water bags, the surface of subsoil will go downward, and it means that the reduced
vertical support occurs over the surface of subsoil. In FE-analysis the process applies
prescribed displacement. It turns out that the stresses at the top pile increase until certain
maximum value and then gradually going down and more stable circumstance afterwards. It
means that the stress concentration ratio, n, would be changing during the consolidation
process. The stress concentration ratio for higher value of h/s will result in the stress
concentration ratio to be higher than those of smaller h/s (see Fig. 6-6). There is strongly
dependent on pile-subsoil relative displacement Sc, which the soil arching is developed at
maximum value when Sc=8-13 mm, whilst beyond this range the soil arching may be less.
129
80 12
ST2 h/s=0.7
60 h/2=1.2
ST4
50 8 h/s=1.4
ST2a
40 ST3a 6 h/s=1.6
30 ST4a
h/s=1.8
20 4 h/s=2
ST1
10 MC h/s=0.7
MC top 2
0 MC h/s=1.4
MC bottom
0 20 0 MC h/s=2.0
0 20
Pile-subsoil relative displacement (mm)
Pile-subsoil relative displacement (mm)
(a) (b)
Fig. 6-6 Variations of vertical stresses and stress concentration ratios, (a) Vertical stress
during water discharge Test 7, (b) Influence h/s on stress concentration ratios
130
A comparison for settlement at the surface of the embankment between measurement and
finite element analysis for Test 1 (C) and Test 7 (h/s=2) is shown in Table 6-14.
Table 6-14 Settlements at surface of embankment using FE-analysis
Settlements at surface of embankment (mm)
Finite element method
Measurement
Test number = 34 o = 44 o
Above Between Above Between Above Between
of pile piles of pile piles of pile piles
Test 1
11 37 13.73 35 13.5 34.5
(h/s=0.7)
Test 7
40 39 29.63 29.63 29.75 29.77
(h/s=2.0)
Experimental work and FE-analysis suggest that the height of the equal settlement plane, hc,
is about 1.4-1.6 times of the cap beam clear spacing or hc=(1.4-1.6)s. Furthermore, to ensure
that no differential settlement occurs at the surface of embankment, the embankment height
of 1.6s is necessary. Fig. 6-8 presents some deformation of embankment including mesh
deformation, the total vertical displacement, total vertical and horizontal strain for h/s=1.4
using FE analysis when subsoil is moved downward 35 mm.
(a) Deformation (b) Vertical strain (c) Horizontal strain (d) Vertical displacement
Fig. 6-8 Shapes of deformation for h/s=1.4
6.3.3. Comparison of Stress Concentration Ratio between Experimental Work and Several
Analytical Methods
It is useful to understand the parameter of stress concentration ratio on arching phenomenon
between experimental results and some analytical methods as well as the finite element
method. Several analytical methods have already been existed such as Method of Low et al.
(1994), Method of Terzaghi (1943), British Standard BS 8006 (1995).
There are seven model tests with different ratios of embankment height to clear spacing (h/s).
Test 1 has the smallest ratio of 0.7, whereas Test 7 has the highest ratio of 2. Meanwhile,
some values in between are 0.9, 1.2, 1.4, 1.6, and 1.8 respectively. Results of Test 1 through
Test 7 are then compared with analytical methods and also the finite element method.
Parameters of material fill are assumed for c=0 kPa, =44o and =15.5 kN/m3. Resume of
equations for 2D-analytical methods are shown in Table 6-15.
131
Table 6-15 Equations in analytical method of two-dimension
Method Stress on pile, p Stress on subsoil, s Remarks
Low et
al.(1994) =a/(s+a)
=0.8
a= width of pile cap
s = piles clear spacing
h=high of embankment
Terzaghi
(1943) K= 0.7
Valid for h/s < 2
BS8006
(1995) Cc=1.95 h/a- 0.18 for
end-bearing piles
Cc=1.55 h/a- 0.07 for
friction piles
hc=1.4s
Then, stress concentration ratio can be easily obtained which n=p /s. Herein, the
comparison is only intended for embankment without reinforcement namely from Test 1 with
h/s=0.7 to Test 7 with h/s=2.0 as depicted in Fig. 6-9.
4 5
Stress concentration ratio
Stress concentration ratio
3.5
4
3
2.5 3
2
Test 1, h/s=0.7 2 Test 2, h/s=0.9
1.5 Low et al. Low et al.
1 Terzaghi Terzaghi
1 BS 8006 (2D)
0.5 BS 8006 (2D)
MC model MC model
0 0
0 10 20 30 0 10 20 30
Pile-subsoil relative displacement (mm) Pile-subsoil relative displacement (mm)
(a) (b)
7 10
Stress concentration ratio
Stress concentration ratio
6
8
5
4 6
Test 3, h/s=1.2 Test 4, h/s=1.4
3 Low et al. 4 Low et al.
2 Terzaghi Terzaghi
BS 8006 (2D) 2 BS 8006 (2D)
1 MC model MC model
0 0
0 10 20 30 0 10 20 30
Pile-subsoil relative displacement (mm) Pile-subsoil relative displacement (mm)
(c) (d)
132
9 10
Stress concentration ratio
12
Stress concentration ratio
10
8
Test 7, h/s=2.0
(g)
6 Low et al.
4 Terzaghi
BS 8006 (2D)
2 MC model
0
0 10 20 30
Pile-subsoil relative displacement (mm)
Fig. 6-9 Comparison between test results, Analytical methods and Finite element analysis for
stress concentration ratio: (a) Test 1, (b) Test 2, (c) Test 3, (d) Test 4,
(e) Test 5, (f) Test 6, (g) Test 7
The figures show that Terzaghi method is always over-prediction in the stress concentration
ratio whereas the BS8006 method suggests strongly under-estimate. Meanwhile, Low et al.
method gives slightly larger result for low embankment h/s≤1.4 but this method good agrees
with high embankment h/s>1.4. Iglesias et al. (1999) have introduced a terminology namely
Ground Reaction Curve by means of normalising subsoil movement downwards to pile clear
spacing. This curve describes a curve of vertical stress over subsoil. When stress on top of the
piles is at a maximum value, the vertical stress over subsoil is at a minimum value.
Meanwhile, the curves above show that the critical values are between 8-13 mm. By using
pile clear spacing of 60 cm, it means that the critical values are from 1.3% to 2.2%.
133
Two control sections were also provided at the edge of the field test. These sections are
embankment without reinforcement. Layout and dimension of different LTP are shown in
Fig. 6-10.
The geosynthetic-reinforced LTPs consist of an aggregate layer with geogrid reinforcement
within the aggregate. Section 1 has a 1.5 m thick beam LTP with four layers of biaxial
extruded polypropylene geogrid (Tensar SS20) reinforcement spaced at 0.3 m apart vertically
within the LTP, and supported on geopiers spaced in a 3.25 m center-to-center square pattern.
Section 2 has a 1.0 m thick beam LTP spaced at 0.3 m apart vertically, and supported on
geopiers in spaced a 2.5 m center-to-center square patterns. Section 3 has a catenary LTP
with two layers of uniaxial woven polyester high-strength geogrid coated with polyvinyl
chloride (Miragrid 24XT) spaced at 75 mm apart vertically placed in a different direction
(longitudinal and transverse) of embankment. This section was also supported on geopiers
spaced at 2.5 m. Last section (section 4) was supported on 0.3 m thick continuous steel-
reinforced concrete LTP over geopiers spaced of 2.5 m.
The selected aggregate used in the LTPs consisted of well-graded crushed granitic rock with
fine material that was less than 3%. This aggregate is normally used as a sub-base layer for
road pavements. Then, over this aggregate blanket, the embankment was constructed using
gravelly sandy clay.
The rammed aggregate pier, called geopier, is a relative the new intermediate-depth columnar
foundation introduced in the construction industry (Fox and Cowell, 1998). Typically, drilled
holes were extended between 2 and 8 m below subsoil surface. In this field test, the initial
drilled diameter of the geopiers was 0.75 m and the initial depth of the drilled hole for the
geopier was 5.5 m.
Site investigation provides a soft silty clay/clayey silt layer as deep as 15 m at some location.
This layer is composed of highly plastic clay with natural content between 35% and 61%.
Field vane tests indicated that the shear strength ranges from 14 to 60 kPa, with most of the
values less than 25 kPa. Moreover, its value of sensitivity varies from 3 to 11.
135
The total settlement resulted from field test and calculation using finite element method is
compared. Treatment between compaction and no compaction using a load of 550 kPa
generated from wheels roller during the compaction process is presented. In addition, Mohr
Coulomb model for whole soil material (embankment and subsoil) and Mohr Coulomb model
for embankment combined with Soft Soil Creep (SSC) model for soft soil are also presented.
Figure 6-11(a) through Figure 6-11(d) show a comparison between field measurement and
finite element calculation for several types of LTP. Fig. 6-11(a) and Fig. 6-11(b) show that
the total settlement of geogrid from field measurement at the bottom layer would be
deformed deeper than at the upper layer of geogrid. Mohr Coulomb model using plastic
calculation with compaction load of 550 kPa can be used to predict the total settlement only
for the final step of compaction at the end of project execution. Although, it cannot follow the
creep phenomenon on soft soil after completion of project execution. Meanwhile, the SSC-
model for the soft soil is better and can follow the creep phenomenon when using the
consolidation calculation type and no compaction load.
4 200
Embankment height (m)
3.5 Compaction 180
160 Model MC Plastic fully compaction
3 load at layer 1 (bottom)
Embankment height (m)
Settlement (mm)
2.5 120 at layer 1 (bottom)
Model MC-SSC Plastic fully
2 100 compaction load at layer 1 (bottom)
0 0 Geogrid Layer 3
Geogrid Layer 4
(a) Section 1
4
140
Embankment height (m)
3.5 Compaction
120
Model MC Plastic fully compaction
Embankment height (m)
0.5 20
Geogrid Layer 2
0 0
Geogrid Layer 3
01-Apr
01-May
15-Jul
04-Feb
10-Feb
10-Feb
15-Feb
20-Feb
01-Mar
01-Jan
15 Agu
15-Jun
27-Jan
30-Jan
(b) Section 2
136
4 350 Embankment height (m)
Compaction
3.5 300 Model MC Plastic fully compaction
load at layer 1 (bottom), Sect. 3
3
250
Model MC-SSC Plastic no compaction
Settlement (mm)
2.5 load at layer 1 (bottom), Sect. 3
200
Embankment height (m)
15-Jul
10-Feb
13-Feb
16-Feb
20-Feb
01-May
05-Mar
25-Mar
7 Agu
15-Jan
15 Agu
15-Jun
01-Jan
19-Jan
4
Embankment height (m)
Compaction 700
3.5
Model MC Plastic fully compaction
600 load
3
Embenkment height (m)
300
1.5 Model MC Plastic adjusted
compaction load
1 200
Model MC-SSC Plastic no compaction
load
0.5 100
Control Section 1 (C1)
0 0
15-Apr
01-May
15-Jul
05-Mar
25-Mar
10-Feb
13-Feb
16-Feb
20-Feb
7 Agu
15 Agu
15-Jun
01-Jan
15-Jan
19-Jan
Similarly, settlement for a catenary LTP over the floating piles will be better predicted with
consolidation calculation type without compaction process regarding with creep as shown in
Fig. 6-11(c). Although, plastic calculation type can be used without considering creep with
compaction at the final of embankment hight. Meanwhile, Fig. 6-11(d) suggests that
settlement without piles will be deeper than those of using piles.
137
6.5. Low Embankment on the End-bearing Piles
In case, the use of end-bearing pile and load transfer platform of catenary (1 or 2 layers of
geosynthetics) Almeida et al. (2008) performed a field test of a low embankment (h/s <1.4) at
Barra da Tijuca district in Brazil. Two situations had already been done by means of
excavated and non-excavated zones. Purpose of the excavated zone is to know the influence
soil support on pile embankment. Lay out for the field test is described in Fig. 6-12.
Fig. 6-12 Layout of field test at Barra da Tijuca district in Brazil (after Almeida et al., 2008)
Properties of materials are shown in Table 6-18 and FE-calculation uses 2-D axisymmetric
and plane strain model. Pavement materials are approached with MC and HS-model and
subsoil with SSC-model. While, linear elastic model is applied for pile and geosynthetics.
138
Table 6-18 Material properties used in finite element analysis
at Barra da Tijuca district in Brazil
Element and model
Concrete Working Fill Subsoil Biaxial
Properties pile platform embankment Geogrid
Linear Linear
MC MC HS SSC
elastic elastic
Unit weight,
3
kN/m 24 17.5 18 18 12.5 -
Friction angle, [o] - 33.8 37 17 5 -
Dilation angle, [o] - 0 0 0 -
Cohesion, c kPa - 11 11 11 11 -
Young’s modulus, E MPa 24000 300 400 - 0.75 -
Secant modulus, E50 MPa - - - 130 - -
Oedometer comp. modulus, Eoed MPa - - - 150 - -
Un-/reloading modulus, Eur MPa - - - 400 - -
Modified compression index, [-] - - - - 0.226 -
Modified swelling index, [-] - - - - 0.012 -
Modified creep index, [-] - - - - 0.005 -
Poisson’s ratio, vur /v [-] -
-/0.15 -/0.3 0.2/- 0.15/- -
/0.3
Tensile strength, EA kN/m - - - - - 1400
Settlements resulted from field measurement and FE-calculation can be looked into in Table
6-19 both 2D and 3D layout for various positions. Meanwhile, Fig. 6-13 shows the settlement
using FE-calculation with the model using both the MC model and HS model for pavement
material and SSC model for the soft soil at the centre of 4 piles.
Table 6-19 Comparison of settlements between at base of embankment
from field measurement and FE-calculation
Unfortunately, FE-analysis finds the difficulties when modeling the excavated sector which
suggests huge deformation. Whereas, on non-excavated sector the settlements could be
approached using the plastic calculation type with involving loading during the compaction
process or consolidation calculation type . Fig. 6-13 shows settlements at the centre of 4 piles
between field measurement and FA-analysis. It can be seen that several weeks after the end
of construction, the settlement in the field was developed highly huge because of traffic and
weak anchorage.
139
140 60
Embankment height (m), Non-
120 50 excavated sector
Embankment height (cm)
Settlements (cm)
40 Excavated sector
80
Model MC-SSC Plastic calculation,
30
Non-excavated sector
60
20 Model HS-SSC Plastic calculation,
40 End of construction
Non-excavated sector
Strains of 3D layout resulted from field measurement particularly on the excavated zone and
compared with axisymmetric configuration are shown in Table 6-20.
Table 6-20 Comparison of tensile strains for geosynthetics between field measurement
and FE-analysis
Strain (%)
Deformation h/b Field FE-analysis
gauge Layout Position measure MC-SSC HS-SSC
ment
DG 01 3D Face of pile caps 1.5 2.05 1.14 1.52
DG 02 3D 1.5 1.73 1.14 1.52
DG 03 3D 1.5 1.51 1.14 1.52
DG 05 3D Half distance between pile 1.5 0.51 0.27 0.36
DG 09 3D caps, parallel to pile face 1.5 0.32 0.27 0.36
DG 06 3D 90o to DG 05 1.5 1.50 0.85 1.14
DG 10 3D 90o to DG 09 1.5 1.36 0.85 1.14
DG 07 3D Centre of four pile caps, 1.5 1.14 0.63 0.83
DG 08 3D align to pile array 1.5 0.97 0.63 0.83
It is clear that the maximum tensile strength of geosynthetics would be occurred at zone of
adjacent piles, and then the lower tensile strength would be found at half distance between
piles with parallel to pile face. This zone is so-called ’primary reinforcement’, whereas other
as ’secondary reinforcement’.
The Hegemann sludge mixture is a mixture of dredged material and additives containing
mainly clay and cement with the following properties: average unit weight =18.6 kN/m3, a
friction angle of =33.8o and a cohesion of 11.5 kPa. The Kyoto road was built over a 9 m
deep of soft soil with a reaction modulus k= 477 kN/m3. Table 6-21 shows soil properties
used in Kyoto road embankment.
Table 6-22 shows Young’s moduli of subsoil used for analysis. These values are found from
soil investigation on site and then tested in the laboratory.
141
Fig. 6-15 Isochrones of geogrids (after Eekelen, 2009)
By using Fig. 6-15, tensile stiffness of material used on the test of road embankment can be
estimated for long-term as indicated in Table 6-23.
Table 6-23 Tensile stiffnesses of geogrids (after Eekelen, 2009)
Direction on Time under Ultimate tensile Tensile stiffness J (kN/m2)
road load strength, UTS J=(% of UTS/strain) x UTS
construction (kN/m) Herein values at 2 % strain
Longitudinal 1 day 350 (25.0/2)x350=4375
Perpendicular 1 day 400 (25.0/2)x400=5000
Longitudinal 1 year 350 (22.8/2)x350=3990
Perpendicular 1 year 400 (22.8/2)x400=4560
Longitudinal 10 years 350 (22.1/2)x350=3868
Perpendicular 10 years 400 (22.1/2)x400=4420
142
Devices TPC t1, t2 and t3 were installed on top of reinforcement. These TPCs measured the
pressure imposed directly on the piles. The total pressure cells above reinforcement layer
(TPCt1, TPCt2, TPCt3) measure load A, whereas total pressure cell below reinforcement
layer (TPCb1) for measuring load (A+B). The vertical load B, which is the load carried by the
geosynthetics reinforcement (GR). By means of tensile forces in the geosynthetics, this load
is transferred to the piles. The curve presenting load B was determined by subtracting the
average measured pressure of TPC t1, TPC t2 and TPC t3 from TPC b1.
The transferred vertical load on each pile was calculated using the equation of
*H+p=18.6*1.15+p= 21.39+p kPa, where p is a surcharge. As soon as arching occurred, the
load was transferred laterally to the piles. Therefore, the pressure measured at TPCs would be
more than 21.39+p. The vertical distance between TPCs and the horizontal line at 21.39 kPa
was an indication of arching.
There are various kinds of material used in the model such as embankment fill, pavement,
subsoil, geosynthetics, pile cap, and pile. Therefore, the behaviour of the material also varies.
HS model is used to model granular and embankment fill. Linear elastic is applied for
modeling concrete pile, wooden pile, and hotmix asphalt. Elastoplastic is used to model
geogrid. Futhermore, SSC model is applied for subsoil. Properties of materials in
axisymmetric configuration of FE-analysis are shown in Table 6-24.
Table 6-24 Material properties used in finite element analysis of Kyoto road
Element and model
Concrete Wooden Hotmix Granular Fill Subsoil Biaxial
Properties Pile caps piles asphalt mixture embkmnt upper lower geogrid
Linear Linear Linear Elasto-
HS HS SSC SSC
elastic elastic elastic plastic
kN/m3 24 9.75 23 20 18.6 10.3 13.4 -
[o] - - - 44 33.8 20.3 17.43 -
[o] - - - 8 0 0 0 -
c kPa - - - 1 11.5 2.05 4.53 -
E50 MPa - - - 160 6.215 1.07 2.36 -
Eoed MPa - - - 130 4.975 - - -
E/Eur MPa 24000/- 12500/- 4000/- -/400 -/14.92 - - -
[-] - - - - - 0.176 0.138 -
[-] - - - - - 0.032 0.010 -
[-] - - - - - 0.011 0.006 -
vur /v [-] -/0.2 -/0.2 -/0.2 -/0.2 -/0.2 0.15/- -
EA kN/m - - - - - - 4100
kv m/s - - - 0.1 1.728E-4 6.86E-4 4.1E-4 -
Konc [-] - - - - - 0.653 0.700 -
Drainage Non Non Non Non
Drained Drained Undrained
Porous Porous Porous Porous
Distribution of load A, B, and C along 2 ½ years for both field measurement and finite
element analysis is shown in Fig. 6-17. Herein, there are three important parts of field
measurement. The first part is the fluctuation of vertical stress on the top pile over geogrid
along 2 ½ years. Secondly, it is the distribution of vertical stress on the top pile below
geogrid. This stress is additional stress due to load below the soil arching which geogrid bring
it to top pile. Last part is a fluctuation of vertical stress on the surface of subsoil below
geogrid. Each part is also compared with results using the analytical method and finite
element method.
143
At center of piles
(a)
At edge of piles
(b)
At center of span
At adge of span
At upper subsoil
(c)
At lower subsoil
Fig.6-17 Load Distribution for A, B and C of vertical stresses observed throughout 2 ½ years,
Prediction using analytical method and FE-calculation (a) Vertical stress (load A) directly on
pile above geosynthetics, (b) Vertical stress (load B) under geosynthetics, (c) Vertical stress
(load C) on subsoil
144
Fig. 6-17 shows that it took several months to develop the arching fully. Perhaps, this
phenomenon was caused by cementation and settlement process of the fill material. From
July 2006 onward, the arching measurements were relatively constant. The fluctuations were
mostly due to variation in the weather, moisture content and the alternating periods of traffic.
The prediction of the load acting directly on the piles using BS8006 is better than that of
EBGEO, which is much higher than the measured values. However, this value is not so
important because for the design purpose for piles the total load is conservatively assumed as
the applied load on piles.
When taking into account the support of the subsoil, the EBGEO gives a better approach for
predicting the vertical load on geosynthetics compared to BS8006. It is an utmost important
thing that the imposed load on geosynthetics directly determines the tensile forces on the
geosynthetics reinforcement.
By using HS model for fill material and SSC model for soft soil, also using plastic calculation
type with 550 kPa as traffic load to predict distribution of vertical stresses for load A, B, and
C, it shows that the vertical stresses resulted from field measurement is in range of FE-
analysis. Furthermore, contour images of stresses from axisymmetrical approach of FE-
analysis with1651 elements are depicted in Fig.6-18. As shown in the figure, critical stress on
pile is at the edge of the piles which punch take place.
145
passages of the day. However, during periods without traffic or rest period so the
geosynthetics reinforcement has an opportunity to restore again.
(a)
(b)
Fig. 6-19 Daily arching cycle on the Kyoto road, (a) Load A, (b) Load B
(after Eekelen et al, 2008)
The following Fig. 6-20 shows the increasing of load B after a passage of a rather heavy
truck, but the arching generally recovers during the rest of the day or weekend, although
other passage is still occurred.
Fig. 6-20 Passage of a truck of 397 kN with 2x2-axles (after Eekelen et al, 2008)
It is clear that higher loading implies the alteration of soil arching. When load at the surface
of the embankment is increased, geosynthetics will bear the additional load and then to be
transferred to pile.
6.7. Summary
Finite element analysis has huge ability and some advantages to solve a lot of problem in
geotechnical engineering. In Chapter 6, some interesting topics regarding with experimental
146
works and some field case studies in Chapter 4 were readily solved using 2D Plaxis software.
Some important findings can be explained as follows:
In case of high embakment, because of certain circumstance, utilizing multi-layer of
geosynthetics provides some benefits. Effectiveness of vertical distance of
geosynthetics depends on properties of geosynthetics mainly tensile strength and
tensile modulus. That is also influenced by elasticity modulus of granular soil, friction
angle and dilation angle. By inserting geosynthetics into granular soil, compound
material is able to bear a higher load and simultaneously reduces the horizontal
outward thrust. For a compacted soil, HS-model provides a good agreement with soil
behaviour. It can reveal the influence of number of layer and pattern of horizontal
displacement. By using this model, it is clearly described that maximum lateral
outward thrust will be occurred at the lower part of embankment for low elasticity
modulus of granular soil, and then at the middle and the upper part for medium and
very high elasticity modulus, respectively.
Vertical stress at top pile and on surface of subsoil will be increased in line with high
of embankment.
When subsoil is moving downward because of consolidation process or low bearing
capacity, vertical stress at top pile will be increased until reaching a maximum value
and then goes slow down and stable afterwards. Otherwise, on surface of subsoil it
will be decreased before getting a stable stage.
Stress concentration ratio (SCR) is influenced by a height of embankment. The more
higher of embankment, it will be increased. Maximum stress concentration ratio or
this critical value is reached when ratio relative displacement between pile and
subsoil to clear spacing of piles ranges from 1.3 to 2.2%.
Maximum SCR using Low et al.’s method suggests more close to the measured
findings in experimental works compared to others (Terzaghi and BS8006). Terzaghi
method provides upper limit, and otherwise BS8006 as lower limit.
Differential settlement at surface of embankment can be minimized or equal to be
nearly zero, when ratio height of embankment to clear spacing of piles is at least 1.4.
Shape of soil arching for a low relative density (55±7%) is isosceles triangle.
Load transfer platform (LTP) can be provided using catenary and/or beam system.
Placement of inclusion these materials must be inside of soil arching.
Using soft soil creep (SSC) model with consolidation calculation type agrees with
behaviour of soft soil compared to other model (MC) and other type of calculation
(plastic).
Influence of loading on piled embankment over undrained soft soil during a
compaction stage is slightly small.
In case of embankment over end-bearing piles, plastic calculation type can be applied
and behaviour of embankment material is better modelled using HS-model.
Meanwhile, for soft soil it can be applied the SSC-model.
Maximum strain of geosynthetics occurs at the face of the piles, and then the lower
one is at the half distance between piles. Meanwhile the lowest is at the centre of four
piles.
147
Finite element provides good results. For instance, in case study of Kyoto road, by
using the HS-model to model the embankment and the SSC-model to model the
subsoil which agree with the field results. Here, axisymmetric configuration with a
load of 550 kPa is applied and plastic calculation type is used for FE-analysis. In
comparison with another analytical method (BS8006), the EBGEO method suggests
better result which is close to the measured findings in the field.
When surcharge load induced by vehicular traffic passes through on surface of
pavement, shape of soil arching will be changed. Not only an increasing vertical stress
or increasing load at the top piles will be increased, but it also happens on surface of
subsoil.
For an embankment with stabilized material, for instance using cement, soil arching
effect fully starts after a period of time.
148
CHAPTER 7
This chapter presents analyses of two sites encountered in the field measurement on Supadio
airport projects, namely runway reconstruction work and Apron widening project. My works
was included in these sites which my role was as a pavement engineer during design stage
and intensively as supervision team when executing the project. The first one is a small part
of Runway Overlay Project in 2009 to level up PCN 50 of Supadio airport runway in order to
be able to serve heavier airplane. Latter is a project that was completed in 2006 to expand the
apron to anticipate high demand of air traffic in the future. Furthermore, finite element
method is applied and then compared with some results of field measurements. Main topic of
this chapter is to discuss settlement of pavement construction over soft soil.
7.1. Location
For purpose of case studies, Supadio airport is chosen regarding with settlement of pavement
construction over soft soil. There are two sites for these case studies. The first site is located
in the middle of the runway and the second one is situated on the new apron. Main problem
in the first site is the settlement of partial weak segment of runway, particularly on wheel
track of the airplane, because of the heavier kind of airplane. Whilst, discussion at the second
site is different settlement in adjacent segment between new construction and existing
construction as well as creep impact for pavement construction over deep soft soil. Aerial
view of location for case studies is presented in Fig. 7-1.
149
In the airport, there is a runway with 2500 m long and 30 m wide which direction of the
runway is 33-15. Whilst, the widening concrete apron with 80×100 m is located in northern
part of existing asphaltic apron. Everyday the airport serves air traffic started from 7 a.m to
10 p.m which the heaviest airplane is Boeing 737-400 with a maximum wheel load around 30
tonnes.
Asphalt Concrete 10 cm
Sand 10 cm
Fig. 7-2 Cross-section of Supadio airport runway (reconstruction pavement in the middle and
existing pavement on left-right side)
150
Wooden piles were driven downward as soon as possible after excavation was completed.
Base course as an important part of this work has to be stabilized. Stabilized base course
using cement content of 6% and then each more or less 30 cm deep of the base course layer
was compacted in this site to get maximum density until the whole depth of base course.
Above the stabilized layer, crushed stone with 15 cm thick was poured and then compacted.
This layer serves to prevent the poor impact of water contacting directly with the bottom
asphaltic surface on the top surface of stabilized base course due to an oxidation process
during curring time. Top layer of pavement is 10 cm thick asphaltic surface course with
Marshall stability 1000 kg. Fig. 7-3 shows the execution steps of runway reconstruction
work.
The aim of reconstruction work of Supadio airport runway in Indonesia is to provide a stable
surface for the operation of the airplane B737-400 which this kind of airplane is heavier than
B737-300 operated before. After about 4 months airplane B737-400 was operated, a
deflection appears around 2 cm that occured in this weak spot of runway. To overcome this
problem, wooden piles 12 cm in diameter and 4 m long using a square pattern of 40 cm piles
spacing was applied to support the weak section and also stabilize base course using 6%
cement.
Everyday the airport serves more or less 30 departures, and if during ramping of airplane to
apron the wheels of airplane touch an observed point around 0.5 second, hence during
airplanes passing through 4 months is equivalent to an hour of parking load and similarly for
36 months around a half day. It is relatively difficult to ensure a settlement caused by the
moving load using finite element method. Though, the settlement coming from an airplane is
very significant. No load means that there is no loading at the surface, whereas fully load is
similar to an airplane that parks at the parking stand all time. Hence, the moving load may
take a value in between.
There are various kinds of material used in the model such as embankment fill, pavement,
subsoil, and pile. Therefore, the behaviour of the material also varies. MC model is used to
model sand, cobble stone, and gravel. Linear elastic behaviour is applied to model wooden
pile, stabilized base course, and hotmix asphalt. Furthermore, SSC model is applied for
subsoil. Properties of material for FE-analysis is shown in Table 7-1. In the FE-analysis the
plane strain after Bergado (1994) was applied to transform from 3D to 2D.
151
Table 7-1 Material properties used in FE-analysis on Supadio airport runway reconstruction
work
Properties Element and model
Wooden Asphaltic Gravel Stabilized Subsoil Coble Sand
piles concrete base course stone
Symbol Unit Linear Linear MC Linear SSC MC MC
elastic elastic elastic
Unit weight, kN/m3 9.5 22.5 20 20 15.5 20 18
Friction angle,
[o] - - 35 44 13.7 33 30
Dilation angle,
[o] - - 0 - 0 0 0
Cohesion, c kPa - - 1 210 8 1 1
Young’s
MPa 12000 4000 400 660 - 300 120
modulus, E
Modified
compression [-] - - - - 0.1008 - -
index,
Modified
swelling index, [-] - - - - 0.0202 - -
Modified creep
[-] - - - - 0.00350 - -
index,
Poisson’s ratio,
[-] -/0.33 -/0.2 -/0.33 -/0.2 0.15/- -/0.33 -/0.33
vur /v
Drainage Non- Non-porous Drained Non-porous Undrained Drained Drained
porous
Calculation type Stage construction, Plastic calculation
Characteristic of soft soil is low bearing capacity and high compressibility when being
subjected by a load. Table 7-2 presents FE-analysis for settlement of the pavement that is
reinforced by wooden piles.
Table 7-2 Settlements and types of loading for pavement reinforced by wooden piles
6.60 mm 6.16 mm
17.70 mm 6.23 mm
152
Table 7-3 presents FE-analysis for settlement of pavement without wooden piles. Settlements
of pavement without wooden piles are deeper than those of with wooden piles. It means that
piles contribute to reduce the settlement of the pavement.
Table 7-3 Settlements and types of loading for pavement unreinforced by wooden piles
Settlement and types of loading
Description
Parking load Moving load
14.87 mm 13.93 mm
33.07 mm 14.06 mm
During the settlement process, soil body beneath a load will shove the outer of adjacent soil
to heave up. It can be seen some pictures in Table 7-2 and Table 7-3. Furthermore, Fig.7-4
describes clearly this phenomenon by looking into incremental displacement.
153
Fig. 7-5 provides settlements on the weak spot of runway both reinforced by wooden piles
and without piles. These results are compared to each other, as well as a comparison with the
settlement of existing pavement construction.
35
Settlements (mm)
20 20
18
16
15
Settlements (mm)
14
12
10 10
8
6
5
4
2
0 0
0 10 20 30 40 0 10 20 30 40
Months Months
It is clear that by using piles to support the pavement construction contributes smaller
settlements than those of no piles. Fig. 7-5 (b) indicates that a new type of construction for
the reconstruction work is not better than the type of the existing pavement construction when
the piles do not be applied to support it. For settlements induced by the moving airplanes, the
FE calculation suggests that after the reconstruction work is completed. The settlement is
more or less 6 mm and this value is stable until 36 months.
154
2
Concrete slab 30 cm, fc 400 kg/cm
2
Soil cement 10 cm, qu 50 kg/cm
120 cm Base course 25 cm, CBR 80%
Subgrade
In airport pavement engineering, there are two important classification values, the values for
airplane and pavement. Value for pavement is so-called PCN (Pavement Classification
Number). Meanwhile value for airplane is so-called ACN (Aircraft Classification Number).
The PCN value must be higher than the ACN value. It means that the pavement is able to
serve a certain airplane. Characteristic of airplane weight and Aircraft Characteristic Number
(ACN) for airplane of Boeing 737-400 is briefly listed in Table 7-4.
Table 7-4 Weight and ACN for airplane of Boeing 737-400 (after ICAO, 1999)
ACN for Rigid Pavement ACN for Flexible Pavement
Subgrade (MN/m3) Subgrade (CBR)
Aircraft All-up Load on Tire High Med Low Ultra High Med Low Ultra
Model mass/ mass one main pressure low low
empty gear leg (MPa) 150 80 40 20 15 10 6 3
(kg) (%)
From data above we know that the maximum weight of one main gear leg is 68,266×46.91%
(32,000 kg or 30 tonnes). For taking an assumption that the shape of the tire imprint is a
circle, the radius of the circle is around 28.32 cm.
Finite element analysis is used to estimate settlements on apron pavement. The consolidation
calculation type is chosen in estimating settlement for soft soil. The SSC model is applied to
model the soft soil. For soft soil which creep phenomenon is really obvious, the model is able
to imitate the creep behavior when being subjected by a load. The MC model is applied for
granular soil. Meanwhile, elastic linear model is applied for wooden piles and concrete slab.
Material properties of some elements for FE-analysis are shown in Table 7-5.
155
Table 7-5 Material properties used in finite element analysis of Supadio airport apron
widening project
Element and model
Wooden Portland Soil Base Sub- Sand Subsoil
Properties pile cement cement course base
concrete course
Symbol Unit Linear Linear MC MC MC MC SSC
elastic elatic
Unit weight, kN/m3 9.5 24 19 20 19.5 18 15.5
Friction angle, [o] - - 44 40 38 35 13.7
Dilation angle, [o] - - - 8 5 3 0
Cohesion, c kPa - - 210 1 1 1 8
Young’s modulus,
MPa 12000 30000 200 500 300 15 -
E
Modified
compression [-] - - - - - - 0.1008
index,
Modified swelling
[-] - - - - - - 0.0202
index,
Modified creep
[-] - - - - - - 0.0035
index,
Poisson’s ratio,
[-] -/0.33 -/0.2 -/0.2 -/0.33 -/0.33 -/0.33 0.15/-
vur /v
Permeability, ki m/day - - - 1 1 1 0.0001
OCR [-] - - - - - - 1.5
Drainage Non- Non- Non- Drained Drained Drained Un-
[-] porous porous porous drained
Calculation type Stage construction, Consolidation calculation
156
Settlements of apron pavement at the 5th year using 2D Plane strain of FE analysis is
described in the Fig.7-8. We can see that the estimated total settlement is around 25.0 cm and
there is no differential settlement (too small) at the base of the embankment between top piles
and surface of subsoil.
Comparison of settlements between the settlement at the first parking stand resulted from
field measurement and settlements from FE-analysis with different time and type of loading
can be seen in Fig. 7-9.
60
50
40
Settlements (cm)
30
No load of airplane
Load on parking stand (12 hours)
20
Load on parking stand (24 hours)
No load (Field measurement)
10
Actual load (Field measurement)
0
0 5 10 15 20 25 30 35
Years
157
Fig. 7-9 presents total settlements at surface of parking stand 1 during the 5 years. The
measured values from field measurement are located between settlements with no load of
airplane and those of 50% occupancy level of parking stand. Of course, actual settlements
from field measurement are higher than those of when no airplane that parks at the surface of
pavement. It is over-estimated when predicting the total settlement with airplanes fully park
on the apron which actually airplanes are coming, parking and then leaving alternately.
Fig. 7-10(a) through 7-10(d) show settlements at different length of piles either no load or
maximum load of airplane that parks on the surface of concrete pavement. From these
figures, we may notice settlements both on top piles and middle between piles at the base of
the pavement for two cases, namely no airplane parks at the parking stand and airplane parks
whole day at the parking stand. Firstly, settlements due to airplane that parks on the surface
of the pavement are deeper than those of no airplane. Secondly, use of longer piles is able to
reduce the settlement for two cases. Thirdly, it is very important to note that the pavement
structure behaves as a block because differential settlement on top piles and in the middle
point between piles is very small for 2 cases. Lastly, settlements on the surface and bottom of
the pavement are very small, it means that the material pavement does not deform.
50 at base of at base of
embankment, no 40 embankment, no
load of airplane load of airplane
40
at surface, with 30 at surface, with
max. load of max. load of
30 airplane
airplane
20
20 at top of piles, with at top of piles, with
max. load of max. load of
airplane 10 airplane
10
at base of at base of
embankment, with embankment, with
0 0 max. load of
max. load of
0 20 40 airplane
0 20 40 airplane
Years Years
158
60 at surface, no load 60 at surface, no
of airplane load of airplane
50 50 at top of piles, no
at top of piles, no
load of airplane load of airplane
40 40
Settlements (cm)
Settlements (cm)
at base of at base of
embankment, no embankment, no
load of airplane load of airplane
30 30
at surface, with max. at surface, with
load of airplane max. load of
airplane
20 20
at top of piles, with at top of piles,
max. load of with max. load of
airplane airplane
10 10
at base of at base of
embankment, with embankment,
0 max. load of 0 with max. load of
airplane airplane
0 20 40 0 20 40
Years Years
Fig.7-11 below presents the estimated settlements of apron pavement without using piles.
These settlements at 30th year for no loading and maximum loading are 50.4 cm and 76.6 cm
respectively. Whilst, for the 5th year, the settlements are 31.9 cm and 54.9 cm respectively. Of
course, settlements without using piles in construction give the deepest value compared to
those of using piles. All settlements in different length of piles show that settlements behave
’block’ where differential settlements at top piles and subsoil in the middle of between piles
are very small.
90
80 at surface, no load of
airplane
70
Settlements (cm)
60 at base of embankment, no
50 load of airplane
40
at surface, with max. load
30 of airplane
20
10 at base of embankment,
with max. load of airplane
0
0 10 20 30 40
Years
159
7.3.3. Effectiveness of Piles Length on Floating Piles
No doubt that use of longer piles on a construction over soft soil contributes a smaller
settlement. For construction over soft soil which creep behavior is really obvious, it is worth
to note that the effectiveness of piles length that is penetrated in soft soil must be taken into
account. Relative Settlement Reduction (RSR) can be used to describe the effectiveness of
pile length when constructing the floating pile group on soft soil. The RSR is a ratio for a
deviation of settlements between no piles and piles to settlement no piles or (So-SL)/So. Where
So is settlement without using piles and SL is a settlement with a pile at certain length of piles.
Fig. 7-12 provides the effectiveness of pile length using RSR parameter for case no load of
airplane and parking load of airplane. It shows when a ratio of pile length to soft soil
thickness is less than 0.2, the length of the piles will not quite significant reduce settlement.
Moreover, time of period and magnitude of loading are an important aspect when observing
the effectiveness of pile length. Short-term monitoring is more effective than long-term
period to know the effectiveness. Whilst, higher load implies that longer pile is needed to
reduce the excessive settlement. In other words, the use of pile length more than 0.2 of soft
soil thickness is an important consideration to overcome the creep behaviour.
1 1
0.9 0.9
Relative settlement reduction
Relative settlement reduction
0.8 0.8
0.7 0.7
30 years 30 years
0.6 0.6
20 years 20 years
0.5 0.5
10 years 10 years
0.4 0.4
5 years 5 years
0.3 0.3
3 years 3 years
0.2 0.2
1 year 1 year
0.1 0.1
0 0
0 0.5 1 0 0.5 1
Pile length / Soft soil thickness Pile length / Soft soil thickness
(a) (b)
Fig. 7-12 Effectiveness of pile length, (a) No load of airplane (b) Parking load of airplane
160
CHAPTER 8
In the framework of research on analysis of dynamic loading behaviour for pavement on soft
soil, several research topics have been presented. The research can be divided into six main
topics consisting of the effective distance of geosynthetics for reinforcement, analysis of
stress concentration ratio on high compressibility of soil, various types of load transfer
platforms on floating columns, performance of low embankment on end-bearing piles, the
influence of traffic load on arching effect, settlement on floating piles under traffic load. The
most important findings of the research are summarized in the following sections and then
followed by recommendations for further research.
8.1. Clonclusions
The research that has been presented is aimed to establish the behaviour of pavement material
under traffic load over soft soil that being supported by piles including both with
geosynthetics reinforcement and without geosynthetics reinforcement. The findings have
been discussed in the previous chapters, and they can be concluded in several topics below:
161
From laboratory work of Ping et al. (2007) and FE-calculation, h/s value around of
1.4 can be used to distinguish between low embankment and high embankment based
on the equal plane strain at the surface of the embankment.
Furthermore, for practical purpose it is important to note here that main conclusions below
are utmost important thing to keep in mind. Firstly, height of the embankment must be more
than 1.4 times of clear spacing between piles to avoid differential settlement at the surface of
the embankment. Secondly, because floating piles foundation works as a pile group, it is
important to remember that ratio pile spacing to pile diameter must be kept always less than
5. Thirdly, the use of small size for diameter of the pile is an useful way when being applied
to low embankment construction. Last thing, when we construct an embankment over soft
soil supported by floating piles using FEM and estimate settlement because of creep,
consolidation calculation type must be run without taking into consider a load compaction
during construction period. However, plastic calculation type can be used to predict
settlement with involving a static traffic load. This type of calculation is relatively faster than
consolidation calculation type.
163
The phenomenon of increasing bearing capacity of displacement piles along the
increasing time after installation is well known and these findings are based on some
empirical field experiments. However, its mechanical behaviour has not been fully
understood yet.
164
REFERENCES