Intermediate Offshore Foundations 2021
Intermediate Offshore Foundations 2021
Foundations
The late Steve Kay was an independent geotechnical consultant with thirty-
three years’ experience as a principal engineer with Fugro, and over forty-
five years as a geotechnical specialist, mainly in the oil and gas industry, both
with contractors and consultants. His expertise was in shallow and interme-
diate (caisson, bucket, can) foundation design, with extensive worldwide
experience in offshore, nearshore and land engineering. He gave suction
foundation courses and master classes and wrote the commercially available
software package CAISSON_VHM.
Preface xiii
Notes from the co-authors xv
1 Introduction 1
1.1 Intermediate foundations 1
1.2 Matching models and data quality for good design 4
1.3 Structure of the book 6
3 Loads 25
3.1 Introduction 25
3.2 Units, sign conventions and reference point 25
3.2.1 Example – units 25
3.3 Structural to geotechnical load conversion 28
3.3.1 Example – load conversion 29
3.4 Geotechnical stresses and strains 30
3.5 Commentary 30
v
vi Contents
4 Marine geology 31
4.1 Geology, sediment types and depositional environments 31
4.2 Lateral variability top layers 31
4.3 Seafloor conditions 33
7 Investigation programs 55
7.1 Introduction 55
7.2 Desk study 55
7.3 Geophysical and geotechnical 56
7.4 Laboratory 58
8 Design basis 59
8.1 General principles 59
8.1.1 Introduction 59
8.1.2 Installation/retrieval/removal 61
8.1.3 In-place resistance 63
8.1.4 In-place resistance – non-co-planar MH loads 64
8.1.5 In-place response 68
8.2 Sign conventions, nomenclature and reference point 68
8.3 Foundation stiffness and fixity 69
8.3.1 Seafloor VHMT loads 71
8.3.2 Foundation lateral and rotational fixity 71
8.4 Load and material factors 72
8.5 Commentary 73
8.6 Closure 73
Acronyms 275
Notation 277
References 281
Index 301
Preface
xiii
xiv Preface
Susan Gourvenec
One afternoon in late 2017, I got a phone call in my office from Steve. He
said he would like to write a book on intermediate offshore foundations,
based on the materials he had provided to the ISO working group. He
wanted the information and the knowledge to be available for engineers to
use. He wanted to know if I would be keen to get involved. It seemed like a
good idea – so I said yes.
Steve called often to discuss things – always interesting, enlightening and
often amusing. Steve retained his humour during all our conversations even
in the late stages of his illness, and his desire for ‘knowledge transfer’ was
sincere.
Undoubtedly a talented geotechnical engineer, Steve also struck me as a
great enthusiast and citizen of our engineering community. Always generous
with his time and knowledge whether at conferences or committee meetings.
It is telling that he sent apologies for missing a meeting of the Technical
Panel of ISO WG10 just a few days before he passed away. A life cut too
short – but Steve certainly made an impression on offshore geotechnics and
on the people he worked with and interacted with in his professional life.
I know several geotechnical engineers of my generation, who worked
closely with and knew Steve much better than I did. They all speak highly of
his generosity in sharing knowledge and of how much they learnt from him.
Two such individuals were drafted (by Steve) into developing this book –
and I have thoroughly enjoyed the opportunity to work with both Elisabeth
and Etienne in creating ‘Steve’s book’, as it has been referred to over the last
couple of years as we worked on it.
I would like also to acknowledge the support from the publisher Taylor &
Francis, in particular Tony Moore, who enthusiastically took on the book
title and worked with Steve and me to put in a proposal. Thanks also to
Gabriella Williams, always so understanding about our various delays and
missing deadlines, and to Frazer Merritt, who took the book over the finish
line to final publication.
xv
xvi Notes from the co-authors
Elisabeth Palix
I got the chance to work with Steve at the start of my career on few projects
involving suction anchor design. We also got the opportunity to write a few
papers together. Steve was a passionate engineer who could give you a call a
few weeks after your last meeting to tell you that he had some thoughts
about your last discussion. I liked his ability to think out of the box and to
take into consideration the opinion of anyone, even a young inexperienced
female engineer. I remember that we had to put energy at first to convince
some engineers that despite the use of the word pile, anchor piles should not
be designed as a pile.
When Steve contacted me to provide him with some support for his book,
it was the first I heard about his illness. He was in a hurry to transfer the
knowledge he acquired during his career. He was afraid that everything
would be lost and disappear into the “Windows recycle bin” (his own
words).
This book has been the opportunity to strengthen the links between us at
the end of his life. I believe that working on this book gave him some energy
to fight his illness and some peace toward the end. I am very glad that I got
the chance to be part of this story with Susan and Etienne. We tried together
to finish this book as Steve would have done if he was still among us.
Etienne Alderlieste
Not too long after I started working at SPT Offshore, I met Steve. He was
contracted to provide engineering support and gave a course on suction
foundation design. For another project in 2014, Steve was asked to be
involved as reviewer for the detailed design of the suction anchors for an
FPSO in the Gulf of Guinea. Throughout the course of that project we dis-
cussed a lot of subjects, including engineering and New Zealand. He enjoyed
sharing knowledge and simultaneously challenging several design choices,
which he was good at, with the aim to end up with a sane design. The final
design accounted for potentially detrimental effects of the by-then novel
phenomenon known as trenching, without being overly cautious. Not con-
servative, that term belongs in politics.
During the years that followed we stayed in touch, discussing various
engineering matters, personal matters, but also matters concerning
CAISSON_VHM – checking functionality of a newer version, assistance
when a Windows update messed with things. With Steve’s health declining,
our contact intensified. He got me involved in this book, and asked me to
provide technical support for CAISSON_VHM to assist Richard, who
would be taking over.
Susan, Elisabeth and I happily accepted the relatively straightforward
task to “finish the book” Steve had drafted. It proved, however, significantly
more challenging for all of us than just those three words suggest. A big
thanks to Susan and Elisabeth; I really enjoyed our meetings and open
discussions!
Notes from the co-authors xvii
All co-authors
We would also like to acknowledge, Richard Kay, Steve’s son, for seeking
out and providing Steve's Mathcad files for recreating figures and analyses
presented in the book.
It has been a significant task, but also a great privilege to work on this
book. We hope you find it useful.
Susan Gourvenec,
Elisabeth Palix and Etienne Alderlieste
June 2020
Introduction
1.1 INTERMEDIATE FOUNDATIONS
1
2 Introduction
(a) (b)
(a)
(b)
Figure 1.2 Intermediate anchor foundation type – load applied beneath seafloor
level (at lug): (a) suction anchor (Dijkhuizen et al., 2003) and (b) tor-
pedo pile (De Araujo et al., 2004).
Introduction 3
Hogervorst (1980) described the first field trials of full-size suction piles,
nearshore Holland, and Senpere and Auvergne (1982) described the first
application of suction pile anchors for the Gorm Field, offshore Denmark.
Due to their simplicity, reliability and cost savings, suction foundations are
also applied to fixed platforms. Bye et al. (1995) described design and instal-
lation to support a jacket structure in very dense North Sea sands. Kolk
et al. (2001) presented studies for the first lightweight platform (Offshore
Cabinda, West Coast of Africa) founded on soft clays and subjected to ten-
sile leg loading. Andersen et al. (2005) presented a keynote paper on the
application of deepwater suction anchors, drawn from an industry spon-
sored study between the American Petroleum Institute (API) and the
Deepstar Joint Industry Project VI. The paper identified at least 485 suction
caissons had been installed at more than 50 locations in water depths up to
2000 m. These were for oil and gas applications. At the time of writing this
book over 40 complete jackets with suction buckets, and more than 500
suction foundation piles, for both oil and gas and renewable platform proj-
ects, have been installed worldwide (Tjelta, 2015). In Europe, up to the end
of 2019, 4258 monopile foundations have been installed as offshore wind
turbine (OWT) substructures, amounting to 80% of the market share by
facility (Wind Europe, 2020). Table 1.1 shows that monopiles represent
over 80% of all European OWT foundation substructures.
Intermediate foundations, which may be driven, suction installed, vibro-
or dynamically installed, have become established for a variety of offshore
applications. In many cases they are an economical alternative to piling,
with savings on materials, installation equipment and installation time.
Perhaps uniquely among foundation designs for fixed structures, intermediate
suction foundations offer a realistic option for complete removal. Apart
from drilled rock sockets for wind turbines and pylons, there are no
intermediate foundation precedents in onshore practice.
Table 1.1 S ubstructures for Offshore Wind Turbines in Europe (Wind Europe,
2020).
Cumulative
At the time of writing no (text)book exists that deals solely with interme-
diate foundations. This book attempts to fill the gap. Aspects of intermedi-
ate foundations are covered in the international standard ISO 19901-4:2016
(ISO, 2016a), but the scope is brief and limited to the petroleum and natural
gas industries. Design guidance based on existing practice for suction
installed intermediate foundations for fixed offshore wind are addressed in
the recent Offshore Wind Accelerator guidelines (Carbon Trust, 2019),
while new design methods for monopiles for fixed wind is the subject of the
PISA project (Byrne, 2017). This book aims to comprehensively address
intermediate foundations for both hydrocarbons and renewables drawing
from the historical context to current state of the art. The target audience of
this book is mainly offshore geotechnical engineers involved in installation
and design. In addition, those from other offshore disciplines may be inter-
ested in reading selected topics.
A theme throughout this book is that models used should be only as accu-
rate as the data used. In particular for offshore geotechnics, this generally
precludes using more advanced models since the required data are simply not
obtained during routine offshore investigations, and engineers have to resort
to empirical correlations. Examples of missing (or inaccurate) parameter val-
ues include (mass) coefficients of consolidation cv and permeability k, pre-
consolidation pressure p'c, and compression moduli Cc/(1 + e0), Cs/(1 + e0).
Geotechnical engineers should not be lulled into a false state of security when
(for example) a 3D finite element analysis (complete with imagery and ani-
mations) has been let loose on poor quality data.
This point – that accuracy of prediction is a function of the quality of the
method as well as the data used – was first made by Lambe in his 1973
Rankine lecture (Lambe, 1973). His original figures are reproduced in
Figure 1.3. Note that, to state his case, the accuracy contour lines in these
figures are either straight or elliptical. Other (intermediate) contour lines are
possible. For example, assuming parabolic-shaped accuracy contours
(assuming equal weights and consistent with the simple equation “accuracy
= model × data”), Figure 1.4 has been prepared. Again, this figure shows
that little is to be gained by using more accurate models with low accuracy
data. For example, using 50% accurate data and model, 25% (= 50% ×
50%) accuracy is obtained. This increases slightly to 40% accuracy (= 50%
× 80%) if the model accuracy is increased to 80%. Note that the same 40%
accuracy can be achieved if data and model accuracy are both modestly
increased to 64%.
A prime example of high-quality offshore predictions being made was for
the North Rankin “A” platform, North West Shelf, Australia. In the 1980s,
Introduction 5
Figure 1.4 A
ccuracy of predictions – equal weights and “accuracy = model ×
data.”
concerns were raised about the belled pile tip resistance in layered weak
calcarenite. High-quality data (triple tube PQ size rock coring providing
virtually 100% recovery and detailed core logging, plus stress-path-specific,
high-pressure consolidated-drained triaxial testing and in-situ plate load
tests) were matched using a double hardening constitutive soil model to
perform finite element analyses of single bells and groups of bells (Smith
et al., 1988).
6 Introduction
7
8 Offshore foundation types and mode of operation
a Resistance: ultimate under combined axial, lateral, moment and twist (VHMT) loads at seafloor
b Response: forces (including moment M and torsion T) and displacements (including rotation θr and
twist θt)
c Seafloor VHMT loads: see Figure 2.2a
d Second order effects ignored (e.g. p-delta for laterally loaded piles)
e Ratio of seafloor V load: tip V resistance can be almost zero for a pile in clay, but up to 50% for a pile
end-bearing in sand
f Flexible foundation: soil system assumed (Murff, 1975)
Offshore foundation types and mode of operation 9
(a)
(b)
Figure 2.1 In-place failure models for shallow foundations (0 ≤ L/D < 0.5):
(a) general shear, squeezing and punch-through failure modes under
pure V load (ISO 19905-01:2016b), (b) punch-through sand into clay
failure mode under VH load (Ballard et al., 2010). (Permission to
reproduce extracts from 19901-4 for a is granted by BSI.)
Offshore foundation types and mode of operation 11
(c)
(d)
Figure 2.1 Continued: In-place failure models for shallow foundations (0 ≤ L/D
< 0.5): (c) VH resistance envelopes for sand and clay (ISO 19901-
4:2016a), (d) rotational failure under VHM loading. (Permission to
reproduce extracts from 19901-4 for c is granted by BSI.)
12 Offshore foundation types and mode of operation
(e)
(f)
Figure 2.1 Continued: In-place failure models for shallow foundations (0 ≤ L/D <
0.5): (e) internal scoop mechanisms (Bransby and Yun, 2009), (f) VHM
yield envelopes for surface foundations with and without tension
(Taiebat and Carter, 2002).
there are six axisymmetric failure modes (Figure 2.1a). The modes, and
models to capture them, become more complicated for combined VH(M)
loading, especially for multi-layered soils (e.g. Figure 2.1b), and are rou-
tinely simplified (e.g. Figure 2.1c). A traditional simplification for shallow
foundations is to use the “effective area” method (ISO 19901-4: 2016a).
The effective foundation width and lengths are given by B' = B − 2 ex and
L' = L = 2 ey where ex = My/V and ey= Mx/V. The point of V load applica-
tion must lie in the middle third in order to avoid tensile stresses at foun-
dation embedment level. It is cautious, especially for embedment ratios of
around 0.5. In addition, plane strain vertical slices may be considered,
usually with shear contributions on the two out-of-plane faces to account
for 3D effects (Figure 2.1d). More possibilities occur if the foundation is
not solid but has ‘skirts’ (Figure 2.1e). A variety of VHM(T) resistance
envelopes are available for rectangular and circular foundations on und-
rained soil (Figure 2.1f). Unlike traditional bearing capacity methods,
which use linear load inclination factors, these use 3D FEA to derive
VHM(T) interaction effects more rigorously. Also, the “effective area” is
no longer needed, since the M resistance is automatically obtained for a
given VH load.
Intermediate Foundations (embedment ratio 0.5 ≤ L/D ≤ 10), Figure 2.2
(clay) and Figure 2.3 (sand):
All the modes shown in Figure 2.2 are for rigid foundations in undrained
soil and under combined VHM load. The first three modes use 2D plane
Offshore foundation types and mode of operation 13
(a)
Figure 2.2 In-place failure modes for intermediate foundations (0.5 ≤ L/D ≤ 10)
in undrained soil: (a) “shallow, intermediate and deep” caisson rota-
tional failure (Kolk et al., 2001).
14 Offshore foundation types and mode of operation
(b)
(c)
There are fewer models for sand than for clay, and these are shown in
Figure 2.3. The first (Broms, 1964a, 1964b) is for laterally loaded
caissons and the second shows the Butterfield and Gottardi (1994)
“rugby-ball”-shaped VHM envelope for shallow (not intermediate)
foundations.
Deep Foundations L/D > 10), Figure 2.4:
Unlike for shallow and intermediate foundations, the axial and lateral fail-
ure models for deep pile foundations are uncoupled. This is reasonable:
uncoupling accounts for the fact that the soil fails first during axial load-
ing whereas the pile fails first (bending stress) under lateral load. There
are no equivalent VHM(T) resistance envelopes for deep piles – only
axial resistance at pile head, Vmax versus depth profiles for various pile
diameters D.
16 Offshore foundation types and mode of operation
(d)
(e)
(f)
Figure 2.2 Continued: In-place failure models for intermediate foundations (0.5
≤ L/D ≤ 10) in undrained soil: (e) VHM envelope for clay (Kay, 2015),
and (f) anchor pile chain trenching failure model (Alderlieste et al.,
2016).
Figures 2.4b and c are for clay profiles. In sand, earth pressures are
slightly more difficult to compute and are usually assessed using pas-
sive wedge theory (Figure 2.4d).
Note that:
(a)
(b)
(c)
Figure 2.3 In-place failure models for intermediate foundations in drained soil: (a) HM load
on laterally loaded rigid pile (Randolph and Gourvenec, 2011) and (b) VHM
envelope for shallow foundations (Butterfield and Gottardi, 1994).
Offshore foundation types and mode of operation 19
(a) (b)
(c)
(d)
Figure 2.4 I n-place failure models for deep (pile) foundations (L/D > 10): (a) axial
(soil) failure, (b) and (c) lateral (pile steel) failure, (d) lateral (soil
wedge) failure (Reese et al., 1974).
20 Offshore foundation types and mode of operation
2.3.1 General
Intermediate suction foundations are closed at the top, open at the base, and
generally cylindrical. In clay, shapes other than cylindrical may be feasible
(for example, square or triangular).
Suction foundation top plates, which contain vent(s), are usually flat. Top
plates may be domed in order to reduce steel bending stresses should
extremely high underpressures be anticipated. Domed top plates have a sin-
gle central vent. A geotechnical side benefit is that domes can accommodate
more soil plug heave. If inverted T-beam stiffeners are placed below the flat
top plate, “rat-holes” are generally added to facilitate water egress. Apart
from depth markers on the outer surface (used for installation monitoring),
all steel cylinder surfaces are left unprotected in order to maximise available
friction resistance. If the surfaces are not rusty, or are painted or treated in
other ways that reduce skin friction, reduced interface friction factors should
be considered. Suction foundation embedment ratios are a function of foun-
dation type and soil conditions, and typically vary between 0.5 and 6.
Cylinder wall thicknesses are much less than driven pipe piles, typically
D/250 to D/100. Top plate thicknesses are higher than cylinder wall thick-
nesses, usually between D/100 and D/50.
With the increase of turbine size, monopiles for offshore wind farms are
experiencing an increase in diameter, but similar lengths, leading to a
decrease of their slenderness ratio (L/D). In the last decade or so, turbine
capacity increased from 4 to 12 MW and turbine suppliers have already
announced that 14 MW turbine will be on the market soon. Current proj-
ects consider monopile diameters in the range of 7 to 10 m. Due to the
conjunction of increasing turbine size and the growing maturity of this
industry (i.e. better optimisation to be linked with a significant R and D
effort in the past few years), the associated embedment ratios are decreasing
from a L/D of 5 to L/D of 3 (or even less). Limits on the size of monopiles
are likely to be driven mainly by the manufacturer’s limits in terms of
fabrication and transport (in Europe some manufacturers can fabricate a
monopile 12 m diameter and up to 120 m long and 2000t) or the offshore
installation vessels available on the market.
2.3.2 Internal stiffeners
Internal stiffeners may be added to intermediate foundations. Examples
include:
2.3.3 External stiffeners
External stiffeners (usually two vertical beams) may also be added to the
outside of the cylinder to reduce structural steel stresses for foundations
installed by skidding (instead of lifting by crane). External stiffeners
added above the top plate are not usually of interest to geotechnical
engineers.
The cylinder is usually of constant outside diameter and mainly of uni-
form wall thickness. However, in sand profiles, the wall (and stiffener) thick-
ness may increase at or near tip level; these friction breakers reduce internal
friction. Similarly, for anchor foundations, the wall thickness may increase
at or near lug level. In clays, the use of vertical stiffeners is preferred over
ring stiffeners to reduce the possibility of gapping (vertical water-filled
annuli or slots) and twist (misalignment in plan).
2.3.4 Protuberances
Protuberances include mooring and lifting padeyes/lugs, longitudinal or ring
stiffeners, changes in wall thickness, mooring chain, launching skids, water
injection tubing and others. Occasionally, fins may be added to increase
lateral resistance of intermediate anchor foundations. Fins are two steel
plates (with stiffeners to reduce bending stresses) sticking out from the cais-
son wall orthogonal to the H load direction.
2.3.5 Embedment ratio
Generally the objective is to get the foundation as deep as possible in order
to maximise foundation resistance. This is because soil generally becomes
stronger (and hence unit values of in-place resistance increase) with depth.
However, suction installation (not in-place resistance) considerations are
usually governing. If an infinitely large underpressure can be applied, and
the steel cannot yield, then installation stop criteria are usually soil plug
liquefaction (sand) and soil base failure (clay). However, intermediate
foundations installed by impact driving have the same stop criteria as deep
(pipe) piles – steel fatigue and maximum blow count.
22 Offshore foundation types and mode of operation
For monopile foundations for offshore wind farms, the diameter is gener-
ally governed by frequency and fatigue requirements. The embedment is
adjusted to limit the accumulated displacement at the top and to ensure that
a small variation of the length would not induce a significant variation of the
mudline displacement. There is no specific installation constraint that would
limit the embedment ratio apart from site specific stratigraphy (e.g. presence
of relatively shallow bed-rock). As discussed in Section 2.3.1, L/D ratio tends
to decrease with increase of turbine size, leading to more rigid foundations.
Figure 2.5 shows that (unlike suction support and anchor foundations),
offshore wind turbine monopile foundations are likely to increase in diam-
eter D (and decrease embedment ratio L/D) in the future.
(a)
(b)
Figure 2.5 E uropean monopile geometry trends. Geometry data from Negro
et al. (2017); wind farm inauguration dates from www.4coffshore.com.
Offshore foundation types and mode of operation 23
Loads
3.1 INTRODUCTION
Loads are the starting point for any foundation design and installation
assessment. This is because a lightly loaded structure (for example, a 100 kN
compressor) can be designed with a different foundation type than a heavily
loaded structure (for example, an offshore wind turbine monopile subjected
to significant seafloor lateral loads, H = 2 MN, M = 150 MNm).
In addition, offshore structural and geotechnical engineers have different
perspectives about loads. Table 3.1 shows that for support foundations: (a)
only the most critical (lowest factor of safety, FOS) load cases need
geotechnical analysis; (b) structural and geotechnical load formats are
different; (c) geotechnical engineers often assume MH loads are generally
co-planar, although less so for offshore wind, and (d) structural engineers
usually consider only in-place conditions. For anchor foundations, the
interface is simpler and perspectives are broadly the same – there are fewer
load cases, and structural loads have now only 2 degrees of freedom (DOF)
instead of 6.
3.2.1 Example – units
Note that using tonne mass eliminates the need for including gravity
(g = 9.81 m/s2) in dynamics problems. Three examples neatly illustrate this.
25
26 Loads
Perspective
Note(s): Displacement symbols φ, θ and Ψ are consistent with those commonly used for ship motions:
x = surge, y = sway, z = heave
φ (phi) = roll, θ (theta) = pitch, Ψ (psi) = yaw
Loads 27
(a) (b)
Figure 3.1 L oad and displacement sign conventions – positive shown: (a) shallow
foundations and (b) anchor mooring.
28 Loads
Figure 3.2 C
omparison of sign conventions: (a) structural analysis (b) geotechni-
cal analysis – laterally loaded foundations – HM sign convention.
Structural engineers generally supply all six structural load components (i.e.
Fx, Fy, Fz, Mx, My, Mz). Equivalent VHMT geotechnical loads need to be
derived from these. As expected, V and T loads are straightforward. However,
resultant H and M are slightly more complicated. Using Excel/MathCAD
function atan2(x,y), the equations to derive VHMT loads are given by:
V Fz (3.3)
H sqrt Fx 2 Fy 2 (3.4)
M sqrt M x 2 M y 2 ; M if H M 180 (3.5)
T = Mz (3.6)
where
αH = resultant H angle x-y plane (radians)
αM = resultant M angle x-y plane (radians)
αH = 8.8 [deg]
αM = 9.2 [deg]
|αH − αM| = 0.4 [deg]
V = −1012.78 [kN]
H = 838.76 [kN]
M = 3340.71 [kNm]
T = −135.07 [kNm]
3.5 COMMENTARY
Marine geology
31
32 Marine geology
4.3 SEAFLOOR CONDITIONS
(a)
(b)
Figure 4.1 Schematic of (a) sloping seafloor and (b) seafloor ripples.
Ad (ii) excessive seafloor slopes (see sketch in Figure 4.1a), may need
levelling.
Ad (iii) for a “sandy” seafloor, seabed ripple heights and lengths should be
determined. This is because an essentially level seafloor is required for instal-
lation in dense sands – otherwise, local piping may occur if there is limited
self-weight penetration, as sketched in Figure 4.1b.
Ad (v) current velocity is required for scour potential in sands.
Ad (vi) the survey should check for the absence of potential obstructions.
Examples and more details are given in Section 6.2 (Hazards for Intermediate
Foundations).
Chapter 5
5.1 INTRODUCTION
Effective (not total) stresses affect soil shear strength and volume change
characteristics of geomaterials.
It is misleading to think that all “sands” behave “drained” and all “clays”
behave “undrained” when loaded. These effects are important. Examples
include:
Table 5.1 lists typical soil responses for intermediate suction foundations. It
is seen that, for foundation sizing to withstand applied loads, installation
and long-term loading conditions generally do not govern. In addition, the
soil response is often “undrained” for a range of soil types and load condi-
tions. This is for both “sand” and “clay”, and also despite the fact that wave
loading has a longer time period (≈ 10 s) than boat impact or seismic load-
ing (< 1 s). Possible exceptions to the “undrained” rule are italicised in Table
5.1, and analysis methods for these special cases are given in Section 5.2.
The remarks here are solely for “intermediate” caisson foundations: “shal-
low” foundations, due to their shorter drainage path lengths, and the fact
that drainage times are generally proportional to (drainage length)2 tend to
be more “drained” than “undrained”.
35
36 Loading conditions and soil drainage
5.2 DRAINED-UNDRAINED
v D / cv (5.1)
where
v = penetration velocity
D = diameter
cv = soil mass coefficient of consolidation.
Figure 5.1 U
ndrained to drained soil response – effect of loading rate on bearing
response during spudcan foundation penetration (Finnie and Randolph, 1994).
under undrained conditions is much more dramatic than for the sand, owing
to a greater tendency for the latter to dilate at low effective stress levels.
G T k / R2 (5.2)
where
G = soil shear modulus
T = load duration
R = pile radius
γ = soil unit weight
k = soil permeability.
Figure 5.2 presents their results for closed-ended piles. Typical analyses
for open-ended pipe piles give η ≈ 1, i.e. a fully drained response (Nguyen
et al., 2012).
38 Loading conditions and soil drainage
(a)
(b)
Figure 5.2 U
ndrained to drained soil response – dynamic loading of closed-
ended piles. Centrifuge data (from Hölscher et al., 2009) effect of
loading rate on pile tip capacity (their Figure 10) and pile tip stiffness
(their Figure 8).
(a) (b)
Figure 5.3 U
ndrained to drained soil response – laterally loaded disk (from
Osman and Randolph, 2012). (With permission from ASCE.) Analytical
solutions for (a) excess pore pressure dissipation at disk/soil inter-
face T50 ≈ 0.5; T90 ≈ 5 (their Figure 7) and (b) pile consolidation
displacement (their Figure 9) T50 ≈ 10; T90 ≈ 75 versus dimension-
less time factor T = c h t/r o 2 .
results. It is noted that time factors for pore pressure dissipation are a factor
20 faster than for pile consolidation displacement. Typical analyses for large
diameter caissons (D = 10 m) in dense fine sand (“Best Estimate” parameters
k = 10−4 m/s, E' = 100 MPa, v' = 0.3, ch ≈ 1 m2/s) give dissipation times
t50 ≈ 10 s, t90 ≈ 100 s, i.e. likely to be “partially drained”–“undrained” for
wind or wave loading. For a large diameter monopile foundation (D = 5 m),
because dissipation/consolidation times are proportional to D2, dissipation
times reduce to t50 ≈ 2.5 s, t90 ≈ 25 s, i.e. “drained”–“partially drained”.
Even though ch can vary by one order of magnitude on either side of the
“Best Estimate”, these results suggest that laterally loaded intermediate
foundations may have to be considered as partially drained when either
wind or wave loaded.
5.2.4 Generic
Unlike the Finnie/Randolph and Hölscher/Huy approaches, the Zienkiewicz
et al. (1980) method can be used for both static and dynamic problems,
and it is not limited to one particular foundation type. One dimensional
40 Loading conditions and soil drainage
conditions are assumed – and this may be optimistic. The key diagram is
Zienkiewicz et al.’s Figure 3 (not reproduced here). The method first assesses
dimensionless parameters π1 and π2, using soil permeability k, drainage path
length L and loading frequency ω, and then find which Zone (I, II or III) is
associated with data point (π1, π2). Worked examples include an earth dam
response during an earthquake (“undrained” response valid only for k <
10−4 m/s at depths well below surface) and pore pressure distribution in
seafloor due to wave loading. The Zienkiewicz et al. (1980) method equa-
tions can be programmed in MathCAD or similar.
5.3 CLOSURE
No soil in the ground has a membrane around it. Hence partial drainage is
always a possibility in all soil types. Even if in-situ tests show drained condi-
tions, due to size effects, the soil may respond undrained or partially drained
in the case of large(r) foundations.
Most uncertainty is usually associated with mass cv (coefficient of consoli-
dation) values to be inserted into analytical equations. In addition, Finnie
and Randolph (1994) assume that bearing capacity modulus is proportional
to foundation width B. This may not necessarily be the case: other proposals
include B2, e.g. Hölscher et al. (2009), Huy et al. (2009) and Osman and
Randolph (2012).
Undrained shear strength may be overestimated when derived from in-
situ testing if the soil during testing behaved drained or partially drained.
This is particularly the case for CPT qc using a Nk factor.
If the foundation design is likely to include unconventional soil response,
then the geotechnical investigation/laboratory test program should make
provision for obtaining reasonably reliable parameter values. Some exam-
ples include:
If reliable triaxial CU data are not available, then possibilities include: (i)
consider a (cautious) drained analysis, (ii) use the DNV (1992) approach
using no or limited dilatancy (dilatancy parameter D = 0 corresponds to
pure shear in p − q stress space, and D = 1/3 for limited dilatancy) and (iii)
use an equivalent undrained shear strength defined by su/σ'vo = 0.25−0.35
for a NC clay.
Chapter 6
Notable case histories illustrating each of these seven hazards are given in
the following sections.
41
42 Hazards, uncertainties and risk minimisation
The 3 Girassol RTAs had been painted externally for protection against
corrosion, and their key dimensions were D = 8 m, L = 19.5 m. During the
self-weight penetration stage, actual penetrations were higher than pre-
dicted and ranged from 15 m to 16 m, compared to 12 m predicted, i.e.
only 3.5 m–4.5 m less than the full length (Dendani, 2003). Using over-
pressure, two push-out tests performed on one RTA one week and one
month post installation gave ≈ 40% and 100% friction increases. These
were in relatively good agreement with laboratory thixotropy test data. A
revised long-term axial capacity model was developed, assuming “coring”,
reduced internal friction resistance above internal ring stiffeners, external
friction resistance reductions (accounting for the painted steel surface and
an upper bound St value) and increased thixotropic factor. An additional
175 tonnes of ballast were added to each RTA to obtain the required pull-
out capacity.
All Girassol FPSO and buoy anchors had also been painted and moreover
had relatively large internal ring stiffeners (0.4 m and 0.2 m outstand respec-
tively). Higher-than-anticipated self-weight penetrations were again mea-
sured (Dendani, 2003), and one FPSO anchor was retrieved (i.e. brought
back to deck) to assist in contingency decision making. No clay was seen on
the painted steel, whereas 20 mm to 30 mm clay thickness covered the non-
painted area, confirming that unit friction of painted steel areas should be
based on the lower soil-paint (not soil-steel) adhesion. In addition to this,
another factor contributing to reduced penetration resistance was the ring
stiffener: there was probably a trapped soil wedge below, and a trapped
wedge of water above, the stiffener. As with the RTAs, satisfying in-place
resistance is more critical than understanding installation. Originally, to
meet certifying agency requirements, the FPSO and buoy anchor embedded
lengths were designed taking into account tilt and twist (both ±10o), but
ignoring reverse end-bearing. Re-analysis, incorporating various
combinations of reverse end-bearing (short term), friction “set-up” (longer
term) and without tilt (as installed) gave satisfactory anchor holding capacity
(Colliat and Dendani, 2002).
The Laminaria soil conditions comprised carbonate silty clay/clayey silt
of high plasticity (Ip ≈ 70%) with a carbonate content around 70%. Based
on in-situ T-bar and laboratory CAU triaxial compression data, a simpli-
fied “Best Estimate” undrained shear strength profile was assessed to be
su = 10 + 1.8z (Erbrich and Hefer, 2002). To moor the Laminaria FPSO, 3
clusters, each with 3 anchor piles, were used. Pile diameters D were 5.5 m,
with embedded lengths L up to 12.7 m, and 20 mm wall thickness (WT).
In addition to lug stiffeners, due to the large D/t ratio (275), 8 internal ring
stiffeners, each with 165 mm outstand, were used. Their axial spacings
varied between 2.25 m and 0.95 m (Figure 6.1). During installation, self-
weight penetrations (around 2.5 m) were as predicted. Thereafter, during
suction installation, measured underpressures were approximately 50% of
Hazards, uncertainties and risk minimisation 43
(a)
(b)
Figure 6.2 S oil flow mechanisms around internal stiffeners during installation of
suction caissons (Erbrich and Hefer, 2002).
reduced capacity under VHM loads. This is due to an unaltered lateral resis-
tance H, and load interaction effects; see Chapter 8, Design basis, Section
8.1 (General principles).
profiles were derived for various design aspects, typically su,DSS ≈ 1.25z. To
moor the FDS (Floating Distribution) semi-submersible in around 1900 m
water depth, a semi-taut mooring system with 4 clusters, each with 4 suc-
tion anchor piles, were used. To withstand factored Tlug = 10,410 kN at
θlug = 25°, final optimised anchor pile diameters D were 4.3 m, with embed-
ded lengths L = 23.8 m, i.e. a penetration ratio L/D = 5.6, as shown in
Figure 6.3a (Newlin, 2003a).
Maximum anchor pile tilt and misalignment values were both ±5°. During
installation, self-weight penetrations averaged 13 m, only slightly above the
12.2 m expected. During the subsequent suction assistance stage, all piles
remained essentially vertical (tilt < 2°), and 14 of the 16 piles had small
misalignments (< 3°). The remaining two piles (L15-P1 and L16-P3) rotated
up to 13°, requiring multiple retrieval and re-installation, sometimes with an
intentional bias (see Figure 6.3b), before they could be accepted. These two
(a) (b)
Figure 6.3 N
a Kika (a) suction pile anchor sketch (Newlin 2003a), and (b) pile twist
versus penetration (from Newlin 2003b). (With permission from ASCE.)
46 Hazards, uncertainties and risk minimisation
anchor piles were re-positioned at least 26 m away from their initial posi-
tion, along the anchor radius circumference (Newlin, 2003b).
The large rotations are unexplained and opinions differ: possible reasons
include non-vertical lug and internal stiffeners, and currents applying lateral
drag to the anchor chain (Lee et al., 2005). Another possibility is unsym-
metrical lifting points used for anchor handling. In all cases, at a given pen-
etration, all pile-soil shear has been utilised axially, so a reasonably small
torque is needed to cause rotation.
Figure 6.4 G
orm soil profiles. Anchor pile L = 8.5–9 m, D = 3.5 m (Senpere and
Auvergne, 1982).
(i.e. tipping out in the 70 kPa stiff clay) was achieved. To install the
remaining anchor piles, the top 1.5 m of the soil plug was liquefied by
changing the toe jetting system into one below the pile top. No under-
pressure data were reported.
Based on these data, it is considered that the root cause of pile refusal was
due to higher-than-expected penetration resistances, causing high suction
pressures and high upwards seepage gradients in the sand plug, leading to
48 Hazards, uncertainties and risk minimisation
Figure 6.5 E xample of local scour resulting from environmental events (copyright
Deltares, 2016).
were found in front of all nine suction piles; see Figure 6.6. The piles were
considered unfit for service and were replaced by OMNI-Max free-fall
anchors.
Detailed 3D FEA of trenches by Alderlieste et al. (2016) have shown that,
for a particular Gulf of Guinea project, the anchor pile VH capacity reduc-
tion due to trenching was approximately 7% for the design θlug angles.
More details on this topic, including possible mitigation measures, are to
be found in Section 12.12 on anchor chain trenching.
6.2 HAZARDS
Figure 6.6 A
nchor chain trench at Serpentina Field, Gulf of Guinea. Surveyed
dimensions (based on Bhattacharjee et al., 2014).
6.3 UNCERTAINTIES
6.3.1 Geotechnical data
The main causes of lack of confidence are usually:
If geotechnical data are insufficient, then there are larger than usual uncer-
tainties about:
These uncertainties are transmitted into the subsequent design phase, usu-
ally as a higher coefficient of variance (COV).
6.3.2 Geotechnical design
Selection of an inappropriate ground model and/or parameter values may
result in:
6.4 RISK MINIMISATION
6.4.1 Ground investigation
The principle of any site investigation is that it is continued until the ground
conditions are known and understood well enough for geotechnical engi-
neering design to proceed safely (Waltham, 2009).
Hazards, uncertainties and risk minimisation 53
6.4.2 Geotechnical data
Subject to economic and time constraints, intermediate foundation sites
should be investigated as comprehensively as possible. Investigation compo-
nents include bathymetric, shallow geophysical and geotechnical data, and
should be preceded by a desk study. Insufficient attention is often given to
desk studies where valuable information can be obtained at low cost.
6.4.3 Laboratory testing
Non-routine onshore laboratory tests on samples for suction foundations
normally include:
6.4.4 Geotechnical design
Judicious selection of soil layering, “Best Estimate” and “High Estimate”
parameter values and associated load/material factors account for both
installation and in-place uncertainties. Similarly, platform legs may be
designed to resist possible jack-up spudcan interaction effects and differential
settlements and designing riser–flow line connections to various types of
support foundations to resist high total settlement estimates. Large
uncertainties generally result in a more robust foundation (i.e. additional
steel weight).
6.5 CLOSURE
It has been shown that intermediate foundations are arguably more hazard
prone, especially during the installation phase, than either shallow or pile
(deep) foundations.
Notable case histories of failure (and their mitigation measures) have
been given, the majority of which occurred during installation. Engineers
54 Hazards, uncertainties and risk minimisation
have learnt from their experiences and, due to sharing these, failures have
become less frequent since suction foundations were first employed in the
early 1980s. Potential hazards for intermediate foundations, both natural
and derived from human activity, have also been listed. These may serve as
useful checklists during the necessary risk management/risk minimisation
process. Appropriate ground investigations and subsequent laboratory test-
ing also assist in reducing project uncertainty.
Chapter 7
Investigation programs
7.1 INTRODUCTION
• desk study
• geophysical (seafloor bathymetry, 3D, shallow seismic and side scan
sonar)
• geotechnical (and laboratory testing).
7.2 DESK STUDY
Desk studies are often the most cost-effective item in the site investigation
process. A desk study may reveal facts that cannot be discovered in any
other way. Despite these merits, the desk study phase tends to be skipped too
often due to lack of awareness. Good practice is to adjust the scope of work
to the intended development, and to bear in mind that more than one
foundation type may need consideration. A good example of this is a
mooring system, where a variety of anchor types are available (e.g. Senders
and Kay, 2002).
Key documents which may serve as useful checklists/templates for
geotechnical desk studies are:
55
56 Investigation programs
Figure 7.1 C
omparative geotechnical survey coverage for shallow or intermedi-
ate and deep foundations for floating and fixed structures (after Kolk
and Wegerif, 2005).
Moreover, since deep piles are no longer used, the depth investigated is
shallow, at most 30 m below seafloor.
Piled foundations are common in shallow water. In this case, geotechni-
cal engineers basically need to know where the competent layer (e.g. dense
sand or rock) is to get enough pile axial capacity. Provided that the pile
impact driving equipment delivers sufficient energy, installation is not a
problem. However, for shallow or intermediate suction installed founda-
tions, geotechnical engineers are interested in relatively subtle variations
in the soil layering within the top 30 m. For example, the weaker layer
may give in-place resistance problems. Similarly, the harder layer may
prevent installation.
In order to minimise risks associated with the hazards listed in Section
6.2, it is occasionally necessary to conduct two separate geophysical sur-
veys. The first is conventional, covering the complete project area. The
second is made at the exact foundation locations and is a high-resolution
3D seismic scan to reveal boulders, dropstones, coral rubble, etc.
58 Investigation programs
7.4 LABORATORY
Design basis
8.1 GENERAL PRINCIPLES
8.1.1 Introduction
Intermediate foundation design usually requires close cooperation between
geotechnical and other offshore disciplines, particularly structural engineers.
The following general design considerations indicate how intermediate
foundation design can be distinctive from shallow and pile foundation
design.
59
60 Design basis
(a)
(b)
8.1.2 Installation/retrieval/removal
Installation is by dead weight, generally followed by suction assistance.
Retrieval and removal are usually by overpressure, possibly with crane
assistance. Impact driving, vibratory driving or drilled-and-grouted
techniques are sometimes used to achieve high penetrations (L/D up to 10)
in competent sand and (weak) rock. Since they are well understood, these
techniques are not discussed in detail in this section.
The following general principles need consideration in assessing
installation/retrieval/removal of intermediate foundations:
Where these are critical, more complex analysis approaches are required.
Note the following:
Comparing Figure 8.2c and 8.2d, note that plug heave occurs in both cases
and, if large enough, will halt installation. However, the heave is due to the
volume of soil displaced by the foundation steel and water pocket formation
respectively. Also, water pockets can occur only in clay over sand profiles,
whereas soil displacement occurs in all soil profiles.
Design basis 63
8.1.3 In-place resistance
Following ISO 19901-4: 2016a Clause 7 for shallow foundations, the gen-
eral principles given below need to be considered in assessing in-place resis-
tance of intermediate driven and suction-installed foundations:
• Bearing failure constitutes any failure mode that could result in exces-
sive combinations of vertical displacement, lateral displacement, or
overturning rotation of the foundation.
• Foundation resistance can be analysed by various methods, including
limit equilibrium and yield surface.
• Limit equilibrium methods shall determine the shape and location of
the critical failure mechanism. These depend on the design loads, soil
stratification and foundation geometry.
• Yield surface methods ensure that factored VHMT loads lie within (or
on) a VHMT foundation resistance envelope.
• Due consideration shall be given to the possibilities of excessive dis-
placement and deformation of the foundation soil. Where these are
critical, more complex analysis approaches are required.
• Calculations using alternative methods of analysis shall include an
explanation of any possible differences due to the method adopted.
• Design loads need to be assessed with due consideration given to the
design life of the foundation.
• Seafloor gradient and/or installation tolerance has to be taken into account
in design. Tolerable foundation tilt and twist should be specified.
• Undrained calculations shall be adopted where no drainage, and hence
no dissipation of excess pore pressures, occurs during loading. This can
occur as a result of the rate of loading or the impermeable nature of the
soil. In contrast, drained calculations shall be adopted where no excess
pore pressures arise during loading. Analysis of foundations subject to
partial soil drainage during the loading event is complex, and specialist
advice shall be sought in these cases. Impact of structural openings,
stiffeners and protuberances shall be taken into account in design.
H and M loads are usually co-planar for single isolated foundations but may
be non-co-planar for unsymmetrical co-joined foundations subjected to
asymmetric VHM(T) loads (e.g. a braced tripod platform). In the latter case,
resultant H and M angles αH and αM need to be calculated, for which the
equations were given in Section 2.6. If |αH − αM|, the angle between H and
M, is significant, then it is usually cautious (i.e. overestimates in-place resis-
tance) to assume they are co-planar for overturning HM loads, and optimis-
tic if HM loads are restoring. In the latter case, options include applying a
higher factor of safety, using 3D FEA on the most critical load cases, or to
seek specialist advice. More details of the effect of non-co-planar MH loads
are given in the next section.
64 Design basis
(a)
(b)
(a)
(b)
(a)
(b)
8.1.5 In-place response
Following ISO 19901-4 (ISO 2016a), the general principles listed need to be
considered in assessing in-place response/serviceability of intermediate
foundations:
Vertical (V), horizontal (H), overturning (M) and torsional (T) loads are
centric and act at a geotechnical reference point (RP), which is at the mid-
point of the foundation at seafloor level; see Figure 8.6. This is the point of
(a) (b)
(c)
Figure 8.6 S ign conventions, nomenclature, reference point and seafloor stiffness
matrices for analysis of intermediate foundations.
Design basis 69
(d)
• Both axial V and torsion T loads are assumed uncoupled from the
lateral (HM) component. This is unlike in-place capacity (Section
10.10, Resistance under Combined VHM(T) Loads) where there is
coupling between V, T and HM. Both V and T resistances are adversely
affected by HM loads, V resistance more than T.
• The seafloor foundation stiffness matrix contains coupling between
the horizontal and rotational terms. This is because the rotation point
is located below seafloor when H and/or M is applied
• The rotation point is assumed to be on the foundation centreline – i.e.
small lateral offset is ignored.
• Torsion loads T are usually small and neglected.
(a)
(b)
Figure 8.7 S eafloor foundation lateral stiffness expressed as (a) horizontal and
moment spring stiffnesses K H and K M , and (b) non-linear stiffness
through springs with K 1 and K 2 .
Design basis 71
Rotational Stiffness k θθ
Foundation Fixity Condition Rotation θ xz M Loads
“Pinned” k θθ = 0 M=0
θ xz = maximum
“Intermediate” 0 < k θθ < ∞ 0 < M < maximum
“Fully fixed” k θθ = ∞ M = maximum
θ xz = 0
factors may be less than unity. This is useful for anchor foundations, as
sketched in Figure 8.9d. Partial load factors per guideline or regulation are
included in Tables 8.2 and 8.3. They exclude dead load factors.
8.5 COMMENTARY
8.6 CLOSURE
Design basis
Load Soil Resistance
Notes:
Refer to actual documents for detailed explanations/definitions/equations.
a API applies φ to resistance. DNV and ISO apply γ to tan(φ’) (sands, effective stress) or s (clays, total stress).γ is equivalent to 1/φ in case of clays, e.g. DNV mate-
m m u m m
rial factor φm of 1.30 is equivalent to API resistance factor γm of 0.77.
b FOS assumptions: (i) either all environmental or all static load (ii) clay (not sand) profile.
c ISO 19902: Section A.9.9.3.1 gives default load factor γL = 1.35 is for GoM. According to Section A.9.9.3.3, higher γL values may be in operation elsewhere: e.g. γL =
1.59 (NW Shelf, AUS) and γL = 1.40 (North Sea).
d ISO 19901: material factor γ = 1.25 is for combined (VHM) loading. For pure V loading γ = 1.5 (also in Section 7.3).
m m
e ISO 19902: 0.9 * γ (default load factor γ = 1.35).
L L
Design basis
75
76
Table 8.3 Partial Load and Soil Resistance Factors for Offshore Piles.
Design basis
Load Soil Resistance
Recommended
Existing Guideline or Partial Recommended Overall Factor of
Regulation (Foundation Recommended Partial Resistance Partial Material Safety b (FOS ≈ γ L .
Type) Loading Condition Load Factor (γ L or φ L ) Factor (φ m ) Factor (γ m ) a γ m)
Notes:
Refer to actual documents for detailed explanations/definitions/equations
a API applies φ to resistance. DNV and ISO apply γ to tan(φ’) (sands, effective stress) or s (clays, total stress).γ is equivalent to 1/φ in case of clays, e.g. DNV mate-
m m u m m
rial factor φm of 1.30 is equivalent to API resistance factor γm of 0.77.
b FOS assumptions: (i) either all environmental or all static load (ii) clay (not sand) profile.
c ISO 19902: Section A.9.9.3.1 gives default load factor γ = 1.35 is for GoM. According to Section A.9.9.3.3, higher γ values may be in operation elsewhere: e.g. γ =
L L L
1.59 (NW Shelf, AUS) and γ = 1.40 (North Sea).
d ISO 19902: 0.9 * γ (default load factor γ = 1.35).
L L
e API RP-2T: Bias factor B (recommended minimum 1.5) on API RP2A (WSD) FOS (i.e. FOS=2.25 for extreme/one line damaged reduced environmental loading and 3.0
Design basis
for operating).
f API RP-2SK: factor of safety for maximum anchor load determined from dynamic analysis.
77
78 Design basis
Figure 8.8 Factor of safety versus failure mode for intermediate anchor
foundations.
Figure 8.9 C
omparison of lumped factor of safety (WSD) and partial factor
design (PFD, LRFD) approaches for support and anchor foundations.
Chapter 9
9.1 INTRODUCTION
This chapter is long, emphasising the fact that installation, particularly those
methods using suction, can be more challenging than both in-place resis-
tance (Chapter 10) and in-place response (Chapter 11). Topics relating to
installation using suction assistance, driven, vibratory and drilled and
grouted installation are considered.
Suction installation topics resemble those for shallow foundations,
whereas those for intermediate offshore wind turbine monopile foundations
are essentially the same as piles.
Retrieval usually also needs to be assessed. In case requirements are not met
(e.g. target penetration depth, foundation verticality or orientation), the
installation contractor has to use overpressure to extract foundation, bump
over (possibly reposition) and reinstall.
Removal is also usually addressed. This is a consequence of site clearance
requirements. Removal is performed using overpressure, since this is less
expensive than soil removal around the foundation outside perimeter and
cutting off the foundation steel, say, 1.5 m depth below seafloor.
79
80 Installation, retrieval and removal
(a)
(b) (c)
Figure 9.2 D
rilling from a jackup using RCD (Reverse Circulation Drilling) tech-
nique (Courtesy of Fugro).
layers or upper fractured rock, a support will be required to avoid hole col-
lapse during drilling. Drilled and grouted installation has seldomly been
used for monopile installation but could become more common as more
complex sites are encountered. Two monopiles of 3.5 m diameters were
drilled and grouted on the Blyth offshore windfarm park; these two piles
were recently decommissioned. More drilled and grouted monopiles of
larger diameter (7 m diameter) should be installed on the French coast in
2022.
Table 9.1 summarises OWT monopile installation feasibility in various
soil and rock profiles. More details are given in Section 9.13.
Installation method
9.5.1 General
This section discusses possible installation measures for both under-penetra-
tion and over-penetration of suction foundations.
The objective is to install the intermediate suction foundation (almost) to
its design (target) penetration depth (TPD) below seafloor. Organisations
have differing approaches regarding target penetration depth tolerance.
These account for installation risks such as depth measurement inaccuracy,
excessive soil plug heave, encountering obstructions or exceeding the allow-
able underpressure. Examples include (a) final penetration depth cannot
exceed, say, 0.3 m above target penetration depth and (b) defining minimum
86 Installation, retrieval and removal
Figure 9.3 P
robability density functions comparing loads and (soil) resistances
for (a) in-place resistance and installation using (b) Best Estimate and
(c) High Estimate soil resistances.
Installation, retrieval and removal 87
9.5.2 Under-penetration
Under-penetration is when self-weight and underpressure alone cannot
achieve target penetration depth. It usually occurs in competent soils, such
as very dense sands or glacial till, where the soil resistance exceeds the High
Estimate value. Increasing the underpressure may cause structural damage.
None of the following measures will work if there is a genuine obstruc-
tion – e.g. a boulder causing refusal. In such cases, retrieval and re-installa-
tion is the only possible measure.
Possible measures to increase the final penetration depth include:
• impact driving
• vibration
• jetting (high pressure)
• water injection (low pressure fluidisation)a
• air lift
• dredge pumpb
• explosions
• free-fall (aka drop fall)c
• cyclic loading
• coat steelb
• friction breakers
• bevel.
Notes:
a water injection: This fluidises sand at tip level. In clays, it lubricates the
liquefied soil plug contents, and the geotechnical engineer has to design for
disturbed soil around the foundation. Similarly, paint, bitumen or other
coatings significantly reduce unit skin friction. Both measures are likely to
affect axial more than lateral in-place resistance.
c free-fall: More details, including equations and an example, are given in
9.5.3 Over-penetration
Over-penetration is when self-weight penetration exceeds target penetration
depth, i.e. when pile weight Wsub exceeds the Low Estimate soil resistance at
TPD. Over-penetration is generally associated with anchor pile foundations in
88 Installation, retrieval and removal
weak clay – the pile weight Wsub is heavy and/or the pile-soil α value is lower
than the Low Estimate (e.g. the soil is usually more sensitive than anticipated).
The soil resistance increases significantly if the top plate meets the clay
soil plug, which is above seafloor due to heave. When this happens, the fail-
ure mode changes from “coring” to “plugged”. The “plugged” : “coring”
installation resistance ratio, ηR, is a function of pile L, D, soil su profile and
pile-soil α. For a normally consolidated clay, at penetration L, ηR is typically
around 2.5 but could vary between 1.5 and 4.5. This resistance increase may
be sufficient to halt self-weight penetration (SWP). In this case (Rplugged ≥
Wsub) we have L < self-weight penetration < Ltot, and the top of the anchor
pile remains above seafloor.
If the resistance increase is not enough (i.e. Rplugged < Wsub) then the pile
would disappear below seafloor during self-weight penetration, and target
penetration depth is exceeded.
Possible measures to decrease the final penetration depth during SWP
phase include:
• closing the top plate vent(s) just before reaching target penetration depth
• adding an external protuberance (pile outstand/collar) at height above
tip equal to target penetration depth
• decrease Wsub such that penetration equals target penetration depth,
and await set-up.
Ad top plate vent closure: This can be done during offshore installation.
Ad adding external protuberance: This needs to be considered during
design. It cannot be done offshore.
Ad decreasing Wsub: This may be achieved using crane hook load tension. It
is a rather extreme measure – the time required depends on the required fric-
tional resistance increase, and the associated strength gain–time relationship.
Alternatively, buoyancy or flotation blocks can be attached to the foundation
to reduce the submerged weight.
The following paragraphs give more details of soil profiles likely to be dif-
ficult/problematic using underpressure. They are listed in easiest (feasible) to
hardest (impossible) order.
Loose sand – easy to install but in-place resistance (liquefaction)
considerations
Clay, sand and sand over clay – generally feasible
Calcareous/carbonate clays require no special considerations. This is
unlike carbonate sands (and silts), which are usually cemented and weaker
than silica/quartz sands.
Sand with occasional clay seams/layers and sand with occasional peat
seams/layers – Thin clay (or peat) seams or layers, if sufficiently thick and
present over the whole foundation footprint, may cause a flow seal and
Installation, retrieval and removal 89
Houlsby and Byrne (2005) have listed the following practical limits to the
maximum attainable suction pump underpressure:
9.8.1 Introduction
There are two main methods of assessing penetration resistance R:
CPT q c based f k f q c sand and clay ,
(9.1)
q k p q c sand and clay
where
R = penetration (or extraction) resistance
Wsteel = caisson submerged weight (and preload, if any)
Ai = caisson inner area = π Di2/4
Di = caisson inner diameter
Fortunately, since the majority of soil profiles are clay, there is no water
flow, and the resistance versus depth profile is usually simple to calculate.
Entering the penetration resistance curve with Wsteel gives the self-weight
penetration. At greater depths, the suction pressure is computed directly
from the resistance curve using equilibrium Equation (9.3).
However, if either underpressure or overpressure create fluid flow, then
penetration resistance R is a function of Δu. An iterative technique is
required since Δu appears on both sides of both load–resistance equilibrium
equations.
Values of Δu need to be compared with limiting underpressures for
base failure (clay) and liquefaction (sand), details of which are given in
Sections 9.10 and 9.11. In addition, for shallow water depths, the maxi-
mum attainable underpressure (Section 9.7) also has to be checked.
9.8.2 CPT qc method
CPT cone resistance qc data are frequently used to predict skirt penetration
resistance. This is reasonable since the prototype skirt and CPT are of simi-
lar size (as is the case here). Foundation skirt tips are usually flat, whereas
92 Installation, retrieval and removal
the CPT tool is conical. Hence coefficient kp, relating unit end-bearing qw to
CPT qc is a shape factor, and has a value of less than 1.
Penetration resistances R are assessed based on the DNV (1992) method.
For a uniform wall thickness (WT) cylindrical foundation, R is calculated
using the equation
R A pk pq c,L Asi k f q c, zdz Aso k f q c, zdz (9.5)
where
L = foundation tip depth
kf = empirical coefficient relating qc (or qt) to unit skin friction resistance
kp = empirical coefficient relating qc (or qt) to unit end-bearing resistance
Ap = foundation wall tip area = π (D2 − Di2) /4
Asi = caisson inner perimeter = π Di
Aso = caisson outer perimeter = π D
D = foundation outer diameter
Di = foundation inner diameter = D − 2 WT
qc,z = CPT cone tip resistance at depth z below seafloor
qc,L = CPT cone tip resistance at caisson embedded/tip depth L
In Equation (9.1), inner and outer skin friction components have been
separated because of differing fluid flow effects. Both terms are integrated
from seafloor to depth L. The DNV (1992) method uses average CPT cone
resistance qc profiles, plus Best Estimate and High Estimate empirical coef-
ficients kp and kf relating qc to unit end-bearing and skin friction, to derive
foundation penetration resistance. DNV RP-H103 (2014) repeats the DNV
(1992) coefficients, but a complicating factor is that assessment is done on a
single CPT. The DNV (1992) averaging approach is preferred – CPTs in
sand can vary significantly over short distances – provided that sufficient
CPTs have been made.
For sand, DNV (1992) indicates that experience has shown that coeffi-
cients kp and kf for the upper 1 m to 1.5 m should be 25% to 50% lower
than those given in DNV (1992) where local “piping” or lateral platform
movement has occurred.
Notation:
unit end-bearing resistance qw = kp qc
unit skin friction resistance f = kf qc or kf qT
qT,CPT = total cone resistance(3)
qT,T-bar= total T-bar resistance
sur,LE = low estimate soil remoulded undrained shear strength
Notes:
(1) from DNV (1992) – k values not stated by Colliard and Wallerand (2008)
p
(2) low internal f due to tip stiffener
Table 9.2 gives kp and kf coefficient values for various soil types. For
North Sea conditions (presumably competent dense sands and stiff clays),
coefficients kp and kf were first proposed by DNV (1992). Lower kf coeffi-
cients were suggested by Colliard and Wallerand (2008) for Gulf of Guinea
clay (normally consolidated, sensitive, high plasticity). Since penetration in
normally consolidated clay is governed by skin friction, the Gulf of Guinea
clay kp values are considered to be the same as for North Sea conditions.
Germanischer Lloyd (2005, 2013) requires higher High Estimate coeffi-
cients than DNV (1992). More recently, Frankenmolen et al. (2017) back-
analysed the Prelude project caisson self-weight penetration data, and their
most accurate fit values (neither Best Estimates nor High Estimates) are also
given. Carbonate soils are notoriously difficult to characterise: penetration
resistances were only slightly higher than the Low Estimate but a factor ≈ 5
lower than the High Estimate. Note that the presence of any cemented seams
within the silt would have prevented suction installation. Finally, Table 9.2
shows that there is room for improvement – especially in sand, where end-
bearing is generally the larger resistance component, and the Best Estimate
and High Estimate kp values vary by a factor 2.
Table 9.3 compares North Sea, Gulf of Guinea and South China Sea clay
characteristics. This may be useful when considering extrapolating Colliard
and Wallerand (2008) kp and kf coefficients to other sensitive, normally
94 Installation, retrieval and removal
Table 9.3 Comparison North Sea, Gulf of Guinea and South China Sea Clays.
Notation:
qc = cone resistance
qT = total cone resistance
Rf = friction ratio
Bq = pore pressure ratio
wp = plastic limit
wL = liquid limit
Ip = plasticity index
γ = unit weight
su = undrained shear strength
St = sensitivity
n/a = not available
z = depth below seafloor
Notes:
(1) Figures 2, 4 and 5 (CPT q ) and Table 1
c
(2) Probably a typographical error: R is typically 3% to 4% in (sensitive) clays
f
Installation, retrieval and removal 95
Figure 9.4 S elf-weight penetration back analyses in NC Gulf of Guinea clay: (a)
adhesion factor α and (b) k f coefficient (Colliard and Wallerand, 2008)
© IHS Markit.
where
c0 = 0.45, c1 = 0.36 and c2 = 0.48
Sand plug heave (i.e. loosening) increases the plug permeability. Assuming
head loss is reduced in inverse proportion to k, and correcting a typographi-
cal error in the denominator, the corresponding Houlsby and Byrne equa-
tion is:
u 1k f / 1 1 1k f (9.7)
where
α1 = c0 − c1 [1 – exp(− (L/D) /c2)]
kf = sand plug : soil mass permeability ratio [–]
Results for kf = 2 and 5 are also given in Figure 9.5a. Using Equation
(9.7), extrapolated αu values are 0.14 and 0.32 for kf = 1 and 3 respectively
at L/D = 1.
Figure 9.6 shows four steady-state flownets obtained for a caisson foun-
dation with L/D = 1, showing equipotentials at intervals of 10% of the
change in pore pressure. A slightly modified version of the finite element
program P72 (Smith and Griffiths, 1998) was used, which solves Laplace’s
equation over an axisymmetric region. A 15-by-15 mesh of 4-noded rectan-
gular quadrilateral elements is used, with no element integration and assem-
bly for the cut-off (foundation wall) elements.
The first flownet, Figure 9.6a, is for an infinite half space of uniform per-
meability. Note that most of the pressure drop occurs within the soil plug
between tip level and seafloor. In addition, the hydraulic gradient is reason-
ably uniform both radially and axially within the soil plug but is a maxi-
mum at seafloor. At tip level, αu ≈ 0.14 (taken by eye between the equipotential
at 0.1 and 0.2), i.e. agrees with Equation (9.7) αu = 0.14.
Installation, retrieval and removal 97
(a)
(b)
Figure 9.5 I nstallation in sand with seepage flow: (a) foundation tip pore pressure
coefficient α u versus penetration ratio L/D and (b) dimensionless flow
parameter F versus L/D. Flow rate into sand plug Q = F (kD∆u/γ sub ).
Parameter k f = soil plug: soil mass permeability ratio. Use the bottom-
most curve (k f = 1) for uniform permeability (Houlsby and Byrne,
2005).
98 Installation, retrieval and removal
(a)
(b)
(c)
(d)
Figure 9.6 C
ontinued: Flownets due to underpressure/overpressure. Axisymmetric
steady state fluid flow. Foundation penetration ratio L/D = 1.0, tip pore
pressure coefficient α u and dimensionless flow parameter F: (c) uniform
permeability k, impermeable boundary at 1.1 D bsf. αu ≈ 0.24, F ≈ 0.50. (d) con-
fined aquifer: as (c) but with impermeable suface. α u ≈ 0.32, F ≈ 0.42.
100 Installation, retrieval and removal
In the real world, the soil is usually of non-uniform permeability and the
far field boundaries are not at infinity, and sensitivity analyses of such rou-
tinely encountered cases are given on the remaining three components of
Figure 9.6. Figure 9.6b shows the flownet for the case when kf = 3 (was
kf = 1). Due to the reduced head loss in the sand plug, the tip pore pressure
coefficient αu increases, αu ≈ 0.30 (was αu ≈ 0.14). This means that the plug
seepage gradient reduces (for the same L/D and underpressure Δu). Again,
αu ≈ 0.30 agrees well with Equation (9.7) αu = 0.32.
Figure 9.6c explores the case when the foundation tip approaches an
impermeable boundary. It shows the flownet when the impermeable base is
no longer at infinity, but at 1.1D below seafloor (bsf) (i.e. 0.1D below tip
level). As expected, reduced radial transmissivity below the tip increases the
head drop in this zone, again leading to an αu increase (now αu ≈ 0.24, was
αu ≈ 0.14).
Finally, Figure 9.6d has the same boundary conditions as Figure 9.6c, plus
an impermeable top surface (outside the foundation). This case models a
confined sand aquifer overlain by clay, in which a “water pocket” (see
Section 9.12.2) may be created at the clay/sand plug interface. Since all
water flow comes from the radial far field boundary (the flownet becomes
rectangular), αu increases by a factor 1.3 (now αu ≈ 0.32, was αu ≈ 0.24 for
Figure 9.6c). Surprisingly, the decrease in dimensionless flow parameter F is
only 15% (now F ≈ 0.42, was F ≈ 0.50). This is for steady state flow. The
following commentary includes diffusion analyses of these last two cases
and shows that their excess pore pressure–time curves differ little.
These findings are also applicable to other L/D values. Hence, using infi-
nite half space solutions and uniform permeability for unconfined sand (i.e.
the first three Figure 9.6 situations) are cautious – this case under-predicts αu
and, as shown in the example that follows, also under-predicts soil plug
average seepage gradient iavg and underpressure ∆u. Since unit skin friction
and unit end-bearing decrease with increasing iavg, penetration resistance R
is over-predicted.
As an example, consider a foundation with embedment ratio L/D = 1, L =
10 m and sand γsub = 11 kN/m3 (i.e. icrit = 1.1). Using the infinite half space
solution and uniform permeability kf = 1 (i.e. the Figure 9.6a flownet), the
αu ≈ 0.15 value gives iavg = 0.85 (8.5 m head loss/10 m length) for an under-
pressure ∆u = 100 kPa. Since the hydraulic gradient iavg (0.85) does not
exceed icrit (1.1), ∆u can be increased. When iavg equals icrit, the critical under-
pressure ∆ucrit value is ≈ 130 kPa (100 kPa × 1.1/0.85). If the sand plug
permeability increases to, say, kf = 3, then the results given on Figure 9.6b
are applicable, i.e. αu ≈ 0.33 (was 0.15). The ∆ucrit value is now ≈ 165 kPa
(100 kPa × 1.1/0.67) instead of 130 kPa. This represents a 25% ∆ucrit
increase – provided that a factor 3 sand plug permeability increase can be
justified.
Regarding penetration resistance R, and defining penetration resistance
multiplication factor ηR = (1 − iavg/icrit), consider the preceding example with
Installation, retrieval and removal 101
Program P80 was used to re-analyse the four steady state cases shown on
Figure 9.6 (i.e. uniform, kr = 3, impermeable and aquifer). Figure 9.7 pres-
ents normalised caisson tip Δu pore pressure ratio at tip level (i.e. αu),
together with the corresponding t50 and t90 data points. The steady state (t =
200 s) αu results are slightly higher than the Figure 9.6 αu values. This is
probably due to program P72/Figure 9.6 averaging inner and outer nodal
Δu values, whereas program P80/Figure 9.7 used inner only.
Figure 9.7 shows that:
1. elapsed times vary between 0.03 s and 0.06 s (t50) and 3.5 s and 6.4 s
(t90) to achieve 50% and 90% of steady state Δu values
2. the fastest dissipation (lowest t90) is 3.5 s for case (b) plug (k or) cv:
mass (k or) cv ratio = 3
3. the slowest dissipation (highest t90) is 6.4 s for case (a) uniform
permeability
4. surprisingly, the Δu − t curves for cases (c) and (d) differ little.
Figure 9.7 C
aisson tip excess pore pressure – time. Axisymmetric unsteady state
fluid flow. Foundation L = D = 10 m. Soil mass c v = 5 m 2 /s for cases (a)
through (d). Soil plug c v = 15 m 2 /s for case (b).
Installation, retrieval and removal 103
In all four cases, steady state conditions are achieved rapidly (t90 within
6 seconds). This implies that (for the given geometry and parameter values)
the steady state assumption is very reasonable. Using steady state Δu values
is slightly cautious, since they are only marginally underpredicted. In addi-
tion, using case (a) is also cautious, since it has the “slowest” response.
To put this time (6 seconds) into perspective, assuming ≈ 1 mm/s penetra-
tion rate and 50 mm tip wall thickness WT, the caisson tip will have pene-
trated another 6 mm, i.e. just over 8% of WT.
Dissipation times (t90) are essentially proportional to 1/cv and L2. The
main source of error in cv (cv = k D/γwater) is the uncertainty in soil permeabil-
ity k (typically 2 or 3 orders of magnitude), followed by uncertainty in con-
strained modulus D (up to 1 order of magnitude). Note that, as the soil plug
seepage gradient increases, k increases (soil becomes more permeable), and
D decreases (soil becomes softer), which may possibly keep cv more or less
constant.
R Qd Q f Q p fAs qA p (9.8)
where
Qf = skin friction resistance along sides of caisson and any protuberances
Qp= total end-bearing resistance at caisson tip and any protuberances
f = unit skin friction resistance
As = pile inside and outside surface areas embedded in soil
q = unit end-bearing resistance
Ap = caisson tip cross-sectional area (excluding contained soil).
For undrained soil response, inner and outer unit skin friction, f, at any
point along the foundation may be calculated by the equation
f s u (9.9)
where
α = a dimensionless factor (see Tables 9.3, 9.4).
Notes: DNV (2005) recommends slightly different (usually higher) α values, e.g. API RP2GEO (2011)
Main Text for pull-out – but they are for driven (not pushed) piles
Pull-out excludes twist considerations, for which α should be decreased slightly
Table 9.6 P
enetration Resistance in Silica Sand. Skin Friction and End-Bearing
Parameters (ISO 19902:2007).
Notes:
n/a : not applicable
Unit skin friction f = min(β σ'vo, flim)
Unit wall end-bearing qw = min(Nq σ’vo, qlim)
9.9.1 Landing on seafloor
When landing on a sandy seabed, local scouring at the edge of the founda-
tion and piping around the skirts are almost inevitable unless the lowering
velocity is very low and the areas for water evacuation through the founda-
tion base are large.
Section 7.3.7.6 of DNVGL-RP-C212 (DNV-GL, 2017, 2019) discusses
piping when landing on sandy seabeds:
Velocities higher than 0.2 m/s occur during uncontrolled self-weight pene-
tration (both caissons and open-ended pipe piles), and their free-fall pene-
tration is discussed in Section 9.9.3.
Top plate vent design and suction pump design are discussed in Sections
12.6 and 12.7 respectively.
Occasionally, when competent soils (e.g. dense sand and very stiff clay) are
present at seafloor, this penetration value may not be achieved.
For a given foundation geometry, possible measures to increase self-
weight penetration are
9.9.3 Free-fall penetration
Free-fall penetration occurs without crane assistance and a heavy founda-
tion, and can also occur with a decreasing soil resistance (“hard over soft”)
profile. Since penetration is uncontrolled, risks need to be minimised. High
free-fall penetration values are associated with heavy intermediate and pile
foundations in clay. There is no free-fall risk during the subsequent suction
assistance stage. This is because suction pumps operate at a semi-fixed flow
rate; hence penetration is displacement, not load, controlled.
Uncontrolled penetration generally exceeds controlled self-weight pene-
tration – this effect is used to install torpedo piles well below seafloor. Free-
fall is best modelled using Newton’s Second Law, forward Euler integration
with time and Archimedes’ Principle for buoyancy effects. Assuming “cor-
ing” in clay and ignoring viscous drag effects on the wall tip, the corre-
sponding pseudo code to model free-fall in soil is:
program free_fall
initialise pile velocity and tip penetration depth;
subscripts 0 and 1 refer to previous and current timesteps
v1 = max(v0, 1e−3 m/s)
z1 = 0 m/s
∆t = 0.01 s
loop while pile velocity v1 > 0:
v0 = v1
z1 = z1 + v0∆t
SF: s
kin friction resistance by integrating su from 0 to
z1
EB: end-bearing resistance using su at z1
SF1 = ∫ Perim α su,z dz
EB1 = Awall su,z1 Nc,strip,deep
RF: vent (nozzle) resistance using v0
108 Installation, retrieval and removal
The above pseudo code is easily modified to include free-fall in the water
column above seafloor, and viscous effects. Viscosity is velocity dependent,
and various torpedo pile researchers have proposed dimensionless drag
coefficients, both above and below seafloor velocity dependent and shaft/tip
viscosity models (e.g. O’Loughlin et al., 2004). For a given vent area, the
vent resistance Rn is a function of the nozzle velocity and inverse discharge
coefficient squared, i.e. Rn α (vn/Cd)2. Section 12.6 (on top plate vent design)
derives the Rn equation and also gives a worked example.
Both the free-fall model and 1-D drivability (wave equation using lumped
masses and springs) analyse pile dynamics using integration with time. It is
noted that pile buoyancy effects are missing from the original Smith (1960)
drivability model, but it is possible in GRLWEAP (Pile Dynamics, 2010) to
include pile buoyancy effects by changing the pile gravity. Office drivability
practice generally excludes buoyancy – all piles are submerged, whether on
land or offshore. Including buoyancy would make piles slightly harder to
drive; this is because the additional upwards resistance increases soil resis-
tance to driving (SRD).
Soil: su = 2 + 1.5 z [kPa, m], ρsoil = 1.6 tonne/m3, ρwater = 1.0 tonne/m3
Pile – Soil – Fluid:
Figure 9.8 compares the resulting Case A (coloured red/thick) and Case B
(blue/thin) acceleration and velocity data versus depth. The following para-
graphs discuss the free-fall and self-weight penetration results.
Figure 9.8 F ree-fall penetration example: (a) foundation acceleration – depth and
(b) velocity – depth. Foundation L = 20 m, D = 5 m, M = 67.7 tonne.
Seafloor velocity = 0.2 m/s. Free-fall penetrations: 7.41 m (suction
foundation with top plate vent) and 13.12 m (open top OWT mono-
pile). Self-weight penetration = 7.37 m. Suction foundation accelera-
tions have been increased by a factor 100.
110 Installation, retrieval and removal
have that skin friction (SF) + end-bearing resistance (EB) = pile buoyant
weight. This can be written in pseudo-code notation as:
SF + EB + BF = M g.
The three terms on the left-hand side – SF, EB and BF (buoyancy force) –
are all penetration depth dependent, and BF decreases slightly with increas-
ing penetration into a higher-density material than water. By trial and error
it is found that the SWP is 7.37 m. This value can be checked by verifying
that the forces balance. At 7.37 m depth, we have SF = 520.6 kN, EB = 38.3
kN, BF = 105.3 kN, giving SF + EB + BF = 664.2 kN. As expected, this bal-
ances the 67.733 tonne pile, whose weight M g = 664.2 kN.
9.9.3.5 Example – commentary
For this foundation, the Case B (OWT monopile) free-fall value (13.1 m) is
just over 175% of the self-weight penetration (7.51 m). Incorporating vis-
cous drag effects are likely to reduce the Case B free-fall value. The Case A
(suction foundation) free-fall value (7.41 m) is only 40 mm deeper than the
self-weight penetration (7.37 m).
9.10 INSTALLATION IN CLAY
The key design challenge for installation analyses in clay is usually calcula-
tion of installation and allowable underpressure. The effect on allowable
underpressure when encountering firm to stiff clay layers overlying sand and
weaker clay layers can also be a challenge.
The allowable underpressure may be limited by maximum suction pump
pressure (usually atmospheric or possibly vacuum, see Section 12.7) and
structural considerations (steel cylinder buckling, etc.).
where
Fi = inner skin friction resistance (= α su,av Aint,surface)
Ai = plan view inside area where underpressure is applied
Nc,circle = circular foundation bearing capacity factor ≈ min[6 (1 + 0.2
L/D), 9]
su,av,tip = average su,z between caisson tip (L) and depth (L + D αD,su)
αD,su = su averaging parameter (αD,su ≈ 0.25)
γm = material factor on end-bearing = 1.5
Suction anchor embedment ratios L/D are usually 6 or less for normally
consolidated clays (e.g. Andersen et al., 2005). Equation (9.13), which is
sensitive to inner friction adhesion parameter αi, usually provides L/D values
well in excess of 6.
ISO 19901-4 (ISO 2016, equation A.77) also uses static equilibrium but
recommends a minimum safety factor of 1.5 on Δuall, i.e. on both REB and
Fi. Consequently, ISO 19901-4 (ISO 2016) gives lower Δuall values (and
lower limit embedment ratios) than DNV RP-E303 (DNV, 2018).
Note that using either the Skempton (1951) or the DNV (2005) equation
Nc,circle is likely to be cautious (i.e. they underestimate Δuall) and provide an
additional reserve on material factor γm. Section 10.5.2 gives reasons for
Nc,circle underestimation.
Randolph and Gourvenec (2011) considered the soil plug stability ratio
(the ratio of underpressure to cause plug base failure to the underpressure
required to penetrate the caisson to that depth). They present design curves
of plug stability ratio and embedment ratio, along with an equation for lim-
iting embedment ratio in normally consolidated soil conditions, as shown in
Figure 9.9.
Section 9.10.3 gives an installation example for clay using the DNV
RP-E303 (2005) equation. The pseudo code may be easily modified to
include other base failure options.
At the other end of the scale, low embedment ratios (L/D ≈ 0.5) have not
yet been investigated in detail. This is a typical potential problem in compe-
tent clays, where installation (rather than in-place capacity) governs.
Conventional reverse end-bearing (REB) theory using Nc,circle = 6 may not be
sufficiently cautious. Numerical 3D FEA may be required to verify that the
REB failure model is indeed applicable. The problem may be in fact non-
axisymmetric. A possibility is an essentially circular/spheroidal failure arc,
beginning inside the caisson and subjected to underpressure (i.e. negative
Installation, retrieval and removal 113
Figure 9.9 D
esign curves for suction foundation installation in clay (Randolph
and Gourvenec, 2011).
surcharge), intersecting the caisson at two points on the caisson tip, and
exiting at seafloor outside. Research may be required.
The former is because a mixture of water and clay can be trapped between
internal ring stiffeners, and that trapped water can give a significant contri-
bution to soil plug heave (Andersen et al., 2005).
The latter is because the amount of soil entering the intermediate founda-
tion is likely to be strongly affected by the proximity to internal base failure.
During the early stages of suction assisted penetration, the factor of safety
against base failure is high, and soil flow around the foundation tip is con-
sidered to be similar to that during the self-weight penetration stage.
However, near the end of penetration, the soil plug is closer to failure requir-
ing mobilisation of reverse end-bearing resistance for stability, and hence a
greater proportion of soil met at the foundation tip is likely to flow inside.
Since installation will fail if the soil plug meets the top plate before the
intermediate foundation has penetrated to target depth, the foundation
overall height should be increased by at least the height of the plug heave.
The additional height is equal to (part of) the foundation steel volume
divided by its footprint area, plus the full volume of any internal stiffeners.
For a 50%–50% split during both self-weight penetration and suction-
assisted penetration stages, and no internal stiffeners, the additional height
is given by ∆Hheave/L ≈ 2 WT/D (Romp, 2013). Hence, for wall thickness/
diameter ratios between 1/300 and 1/100, height increases could be up to
2% of the embedded length. For large intermediate foundations (D = 5 m, L
= 30 m), ∆Hheave could be around 0.75 m. More soil plug heave can be
accommodated by domed (than flat) top plates.
R suction Fi Fo Q w (9.14)
Installation, retrieval and removal 115
where
Fi Asi i su, zdz (9.20)
and
αi = clay-steel inner adhesion factor (unit friction fi = αisu)
αo = clay-steel outer adhesion factor (unit friction fo = αosu)
Aplug = caisson internal (suction) and clay plug area
= π Di2/4
Asi = caisson inner perimeter
Aso = caisson outer perimeter
Awall = caisson wall tip area
L = caisson tip penetration into clay
Lcan = caisson total length
Nc,strip = strip foundation bearing capacity factor (clay)
≈ min[5 (1 + 0.2 L/D), 7.5]
Nc,strip,deep = deep strip foundation bearing capacity factor (clay) ≈ 7.5
R = soil resistance
su,av,tip = average su,z between caisson tip (L) and depth (L + Di αD,su)
(note: Di, not D)
αD,su = su averaging parameter (αD,su ≈ 0.25)
Wsub,steel = caisson submerged weight (including ballast)
The corresponding pseudo code for the above installation in clay is:
program install_clay
initialise and evaluate constants:
2 WT, Asi = π Di, Aso = π D, Ap = π (D2 - Di2) /4 and
Di = –
Aplug = π Di2/4
116 Installation, retrieval and removal
∆L = 0.1 m
nL = Lcan/∆L + 1
zero matrix ANS (size nL rows by 12 columns)
k = 0
k = k + 1
L = (k – 1) ∆L
Fi: inner skin friction resistance by integrating su
from 0 to L
Fo: outer skin friction resistance by integrating su
from 0 to L
Fi = Asi αi ∫ su,z dz and Fo,clay = Aso αo ∫ su,z dz
Qw: end-bearing resistance using su at depth L
Qw = Awall Nc,strip,deep su,L
Rsuction = Fi + Fo + Qw
Rretrieval = Rsuction
∆usuction = (Rsuction − Wsub,steel)/Aplug
∆uretrieval = (Rretrieval + Wsub,steel)/Aplug
su,av,tip = average su,z between caisson tip (L) and
depth (L + Di αD,su)
∆uall,install = Fi/Aplug + Nc,circle su,av,tip/γm,install
∆uall,rerievall = Fi/Aplug + Nc,circle su,av,tip/γm,retrieval
store in kth row of matrix ANS:
k, L, and L/Lcan in columns 1 through 3
Fi, Fo, Qw, Rsuction and Rretrieval in columns 4 through 8
∆uall,install and ∆usuction in columns 9 and 10
∆uall,retrieval and ∆uretrieval in columns 11 and 12
Installation plots:
col# 07 (Rsuction) versus col# 02 (L) installation resistance – depth
col# 10 (∆usuction) and col# 09 (∆uall,install) versus col# 02 (L) underpressure
– depth
Retrieval plots:
Using the data just presented, at the end of penetration (L = 17 m), we have
k = 171
L = 17 m
L/D = 4.25
and
and
118 Installation, retrieval and removal
Figure 9.10 shows two resistance and pressure profiles. Figure 9.10a is for
penetration. It is seen that, at final tip penetration (17 m bsf), the required
underpressure ∆usuction (blue/thin line) is ≈ 175 kPa. The corresponding
allowable suction pressure Δuall (red/thick line) is just over 350 kPa. Since
Δuall exceeds ∆usuction, there is no major risk of soil plug base failure during
suction assisted penetration. The self-weight penetration value is 7.6 m, and
is more than sufficient to start the suction assisted stage. The high self-weight
penetration (nearly half the embedded pile length) is not uncommon for
anchor piles in weak clay, and self-weight penetration would decrease for αi
and αo values in excess of 0.3.
Figure 9.10b shows the corresponding retrieval profiles. These were
obtained using the same αi = αo = 0.3 value and the same soil su profile, but
with a γm = 1.25 (was 1.5). Because there is (reverse) end-bearing, the
retrieval and penetration resistance profiles are identical. Unlike penetra-
tion, retrieval resistance (and overpressure) drops once the pile starts to
move upwards. Hence the overpressure profiles (both required and allow-
able) are valid only at/near the penetration depth at which retrieval com-
mences. Since pile submerged weight opposes retrieval (instead of assisting
penetration), the overpressure magnitude is higher than the corresponding
underpressure: at 17 m bsf, ∆uretrieval is now 295 kPa (compared to ∆usuction
= 173 kPa). Despite the higher pressure, there is still no base failure risk at
all depths, except apparently above 3 m bsf. This is unlikely to occur – this
zone is well above the SWP depth (7.6 m bsf), and a vessel crane is used to
keep the suction anchor vertical and lift the it from the seafloor.
9.10.3.2 Commentary
Interested readers are encouraged to perform sensitivity analyses. For a
given su profile, key input variables include αi, αo and Wsub,steel. The model
assumes a uniform pile geometry. In practice, anchor piles are rarely
Installation, retrieval and removal 119
Figure 9.10 I nstallation and retrieval resistance in clay. Example using f = αs u method.
Anchor pile D = 4 m, L = 17 m in LOC clay (s u = 5 + 2z [kPa, m]).
BE) α values for retrieval and lower (LE or BE) α values for installation.
Another reason may be that, during the preliminary design phase, when
both installation and retrieval assessments are made, installation data are
available from previous projects nearby, thereby reducing the uncertainty
bandwidth. Other organisations may have other γm,retrieval opinions. In any
case, γm,retrieval values of less than 1.5 should be justified.
Since the tip resistance contribution is small, and unit friction f = αsu, it is
obvious that various α and su combinations can be used to compute the
same f (and hence assess resistances and pressures). In the example, su has
been kept constant.
Even though there is limited risk of plug base failure during suction-
assisted penetration, base failure could be possible during the preceding self-
weight penetration stage if high overpressures occur in the “water plug” due
to undersized vents. Section 12.6 (Top Plate Vent Design) below includes an
example calculation.
For suction foundations in clay, two resistance estimates are usually pro-
vided at various time periods after installation:
9.10.5 Boulders in clay
Foundation tip integrity (buckling initiation) may be adversely affected by
large-sized material. Such materials include boulders in glacial till, drop-
stones in weak clays, coral rubble in carbonate clay, and flintstones in chalk.
Like sand and gravel, boulder angularity can vary from angular (sharp edges
with relatively plane sides with unpolished surfaces) to rounded (smoothly
curved sides and no edges); see Figure 9.11.
This section outlines a reasonable model and design procedure should
there be a possibility of encountering boulders (or similar) during suction
foundation installation in clay. Assumptions made include (a) the boulder is
spherical with diameter Dboulder and (b) a constant coefficient of friction μ
models the boulder–steel interface and (c) diameter Dboulder < D should the
boulder be within the pile. The design procedure consists of three stages: (i)
check for refusal, then, if there is no refusal, (ii) check foundation tip integ-
rity and modify foundation tip geometry if necessary. Then (iii), repeat
stages (i) and (ii) if the geometry has changed. Since the forces acting on the
pile tip are needed for all stages, these are first derived in the following
paragraphs.
Figure 9.11 B
oulders, Whiterose Development, Grand Banks, Offshore
Newfoundland. PanGeo (2010). This area is unsuitable for suction
foundations - ground conditions include hard clays and very dense
granular glacial till with boulders.
122 Installation, retrieval and removal
Fx = Fnfac1 (9.23)
Fz = Fnfac2 (9.24)
Figure 9.12 B
oulder in clay – left, boulder and pile wall contact point, and right,
forces acting on the boulder.
Installation, retrieval and removal 123
Figure 9.13 B
oulder in clay – boulder polygon of forces and equivalent pile tip
forces.
If 90 o ≤ θxz ≤ 180°
Fx = Fnfac3 (9.25)
Fz = Fnfac4 (9.26)
where
and
Rp F x
2
Fz 2 (9.32)
and
2
Fn R s 2 1 2 – Wsub,boulder fac1 Wsub,boulder fac2 / 1 2 (9.33)
Tn Fn (9.34)
Vboulder 4 / 3 Dboulder / 2
3
(9.39)
Figure 9.14 shows how pile tip forces Fx, Fz and Rp vary with θxz (Example
1 following, L = 17 m, μ = 0.6), and Table 9.7 lists force variations for
Table 9.7 Boulder and Pile Tip Forces – Sensitivity to Friction Coefficient μ.
Notation:
μ = boulder-pile friction coefficient
δ = boulder-pile friction angle, atan(μ)
Fx,max = maximum pile tip lateral force
Fz,max = maximum pile tip axial force
Rp,min= minimum pile tip resultant force, Rp =√(Fx2 + Fz2)
Rp,max= minimum pile tip resultant force, Rp =√(Fx2 + Fz2)
Rs = soil resistance on boulder, Nc Aboulder su
Wsub,boulder = boulder submerged weight,Vboulder γsub,boulder
∆Rboulder = net soil resistance on boulder, Rs - Wsub,boulder
Notes:
(1) Refusal Example 1 – weak clay with 2 m diameter dropstone
(2) R = 995kN
s
(3) W
sub,boulder = 57 kN
(4) ∆R
boulder = Rs - Wsub,boulder = 938 kN
friction coefficient μ varying between 0.2 and 1.0. This figure and table are
key items. Note that:
(a) Fx, Fz and Rp kinks at θxz = 90° are due to the sign change in Tn
(b) maximum lateral Fx has almost the same magnitude as maximum
axial Fz. (Example 1: Fx,max = 827 kN, Fz,max = 938 kN)
(c) axial Fz = Fn when θxz = 90°
(d) lateral Fx = Tn when θxz = 0° and 180°
(e) maximum axial Fz occurs when θxz = atan(1/μ) (Example 1: 59° for
μ = 0.6)
(f) resultant Rp magnitude is essentially constant (Example 1: Rp,min = 938
kN, Rp,max = 964 kN)
(g) maximum resultant Rp.max ≈ ∆Rboulder (Example 1: Rp,max = 964 kN,
Rs = 995 kN, Wsub,boulder = 57 kN, ∆Rboulder = Rs – Wsub,boulder = 938 kN,
Rp,max/∆Rboulder = 1.028)
Findings (f) and (g) make possible a reasonably accurate simplification for
both refusal and integrity assessments. That is, ∆Rboulder (= Rs − Wsub,boulder)
the net soil resistance, can be applied as a force Rp.max on the pile tip at any
angle θxz. – anywhere from vertical to horizontal. Since the force initiating
lateral buckling is 50% of that for axial buckling (Aldridge et al., 2005), this
means that lateral is more important than axial.
126 Installation, retrieval and removal
9.10.5.2 Refusal
Encountering boulders causes “spikes” in soil resistance–depth profiles. This
is similar to CPTu encountering gravel in sand, where the qc value temporar-
ily increases as the obstruction is pushed aside. However, underpressure ∆u
increases whereas CPT u2 decreases. Refusal occurs if ∆u exceeds allowable
Δuall,install.
It is seen that, if a single boulder is met at final tip penetration (17 m bsf),
the required underpressure ∆usuction,boulder rises by 46 kPa to 219 kPa (was
173 kPa). The corresponding allowable suction pressure Δuall was 353 kPa.
Hence, the design is satisfactory – there is no refusal (∆usuction,boulder ≤ Δuall).
Had the boulder been encountered at a shallower depth, or the initial con-
tact point was not directly on top of the boulder, then both Rboulder, and
∆usuction,boulder would be smaller. Note that ∆Rboulder, the net soil resistance,
subtracts Wsub,boulder from Rboulder.
Foundation tip integrity assessments will subsequently demonstrate that
the tip WT has to be increased from 20 mm to at least 35 mm to preclude
local buckling.
Because both pile and soil data are different, a separate calculation using the
Section 9.10.4 program is needed. The fine details are not presented. At tar-
get penetration depth (L = 10 m) and no boulder, we have:
It is seen that, if two boulders are met at final tip penetration (10 m bsf),
the required underpressure ∆usuction,boulder is 1.36 MPa (was 0.92 MPa). The
corresponding Δuall value is 1.46 MPa. There is again no refusal (∆usuction,boulder
≤ Δuall), but is marginal (∆usuction,boulder/Δuall = 0.93). This implies that γm,install
is closer than usual to, but still exceeds 1.5. Major challenges are likely had
both boulders been larger than 2 m diameter. Foundation tip integrity assess-
ments will subsequently demonstrate that, for the current pile and soil data,
WT has to be increased from 50 mm to 100 mm to preclude local buckling
and buckle propagation. Underpressures are extremely high, and structural
steel checks are required. A domed top is necessary.
Fz 2.8 y WT 2 (9.41)
where
y pile yield stress 345 MPa for steel (9.42)
D / WT
3
5 1 soil 2 E pile / E soil (9.43)
where
Since both examples do not refuse, pile tip integrity assessments should be
made. Again, take forces Fx and Fz equal to net boulder resistance ∆Rboulder.
Example 1 and 2 ∆Rboulder values were 0.566 MN and 6.14 MN respectively.
Installation, retrieval and removal 129
Notation:
n/a not applicable - no local buckling
Notes:
(1) point load from a single boulder
(2) pile: σ = 345 MPa, E
y pile = 210e3 MPa
soil = 400 su,tip, νsoil = 0.5
(3) soil: E
The friction coefficient μ value (whether static or sliding), and where the
point of contact is on the boulder, are both irrelevant: it is reasonable to
assume that the Rp magnitude remains essentially the same. It is also consid-
ered sufficiently accurate to use the net force (∆Rboulder = Rs − Wsub,boulder) for
both axial and lateral integrity assessments.
The model described here has assumed that the bearing capacity factor of
a sphere equals Nc,circle. Ball penetrometer research indicates higher Nc values
for deep (full flow) conditions, namely 11 < Nc,ball < 15 depending on the
interface α (0 − 1) value (Randolph et al., 2000). No solutions are available
for “wished-in-place” near-surface spheres at shallow depth, but Nc values
are also expected to be higher than those for shallow circular foundations.
As usual, opinions differ, but 6 < Nc,circle < 9 used herein is considered reason-
able: major model inaccuracies include assumptions of shape (e.g. an ellipse
instead of sphere) and ULS (the boulder is pushed aside before achieving the
peak soil resistance).
In practice, since the geometry has changed for both the examples, Stage
(iii) (= Stage (i) re-check for refusal followed by Stage (ii) reassess pile tip
integrity) should be carried out. Stage (i) implies additional ”no boulder”
installation analyses with (a) increased Wsub,steel and (b) increased WT values.
However, increasing Wsub,steel further increases Δuall, which is beneficial. The
”no boulder” soil resistance Rsuction value will not increase markedly, since
the Qwall increase will be offset by a smaller Fi value. It is also advisable to
perform analyses with other alpha values (αLE, αBE, αHE), and to check that
there is no risk of over-penetration (particularly for the weak clay Example 1)
due to increased foundation weight, Wsub,steel.
where
Achain cross sectional area
1 / 2 z lug / tan volume z lug Dboulder / 2 (9.48)
Vchain may be compared with soil volumes swept out by an anchor pile.
Plugged implies that the soil plug contains internal stiffeners/protuberances,
Installation, retrieval and removal 131
Typically, Vchain ≈ 0.5 Vpile,coring, but, as usual, swept volumes are sensitive
to pile, chain and boulder geometry. Nevertheless, assuming boulders are
randomly distributed in 3D space, then, even though the probability of
boulder contact is greater for the pile than the chain, due diligence implies
that the chain–boulder contact should also be analysed. The major issue is
that θlug will increase, decreasing the FOS on axial pull-out failure. Whatever
analysis option is selected, the crux is to select parameter values giving the
required p (soil normal/bearing resistance per unit chain length) and f (soil
frictional resistance per unit chain length) (Vivatrat et al., 1982) and to uni-
formly distribute p and f over an embedded chain length, ∆Lchain modelling
the contact zone. A reasonable ∆Lchain value is Dboulder /4. Using the pile tip
model, p = (Rs − Wsub,boulder) /∆Lchain and f = μp. Analysis options include (a)
introducing a competent soil layer at the boulder depth and (b) modifying
the chain analysis program. These options are left to interested engineers to
investigate.
9.11 INSTALLATION IN SAND
(a) to sufficiently reduce inner friction and foundation tip resistance such
that penetration can occur
132 Installation, retrieval and removal
without incurring
(e) the maximum underpressure should not exceed the lowest astronomi-
cal tide (LAT) water depth (Section 9.7).
Ad (a) reduced inner friction and foundation tip resistance. Unlike clay,
most resistance is from the foundation tip, and the primary effect of under-
pressure is to reduce this tip resistance rather than to increase the load
(which is a secondary effect).
Ad (b) sand plug liquefaction failure. Underpressure is traditionally lim-
ited to that which gives zero vertical effective stress, i.e. seafloor icrit (critical
hydraulic gradient) = γsub /γwater.
The foregoing implies that underpressure cannot exceed this critical value.
However, there is some evidence from carefully executed laboratory perme-
ability tests (e.g. Fleshman and Rice, 2014; Panagoulias, 2016; Panagoulias
et al., 2017), and from foundation installation, both in the field and centri-
fuge (summarised by Panagoulias, 2016; Panagoulias et al., 2017) that i >
icrit in medium-dense to very dense sands. Figure 9.15 shows suction
Figure 9.15 N
ormalised suction penetration curves (Panagoulias et al., 2017)
(With permission of SUT, London)
Installation, retrieval and removal 133
penetration data for six North Sea locations, together with the theoretical
(Houlsby and Byrne, 2005) icrit curves. No sand plug liquefaction occurred.
In addition, no excessive sand plug heave was reported in any of the six
cases considered. Hence, provided that due care is taken, gross failure is
unlikely at underpressures (slightly) exceeding this critical hydraulic gradi-
ent. However, at the time of writing, a credible model has yet to be pub-
lished, and installation procedures have not yet been generally accepted.
Upward seepage, even if at pressures less than critical, also creates some
loosening/density reduction within the soil plug (Tran, 2005). This leads to
internal sand heave, which, if not allowed for by stick-up/over height, may
prevent the foundation penetrating to its design depth. In addition, Tran’s
small-scale model tests showed that no major adverse hydraulic conditions
(only 10% heave) occurred, even with L/D ratios up to 2.0 and critical gra-
dient > 1. Another factor mitigating severe sand plug liquefaction is proba-
bly passive arching in aged sands – arching theory (for the no-flow case) has
been given by Paikowsky (1989).
The loosening/density reduction effect needs to be accounted for in in-
place resistance assessments – the sand plug strength is reduced from origi-
nal in-situ values.
Ad (c) piping failure. More important than heave is piping. If sand “boil-
ing” occurs, piping channels may be locally formed, especially along the
wall/soil plug interface. These will create a hydraulic leak – the water plug
beneath the top plate is in direct contact with the seawater via a continuous
channel along the foundation inner and outer wall. In this case, the suction
pump is merely removing seawater. There is no water flow within the soil to
reduce inner friction and end-bearing resistance. Also, a mixture of sand and
water may come into the suction pump. In this case, installation will fail.
From the very limited amount of data available, it appears that piping
failure is more likely when the foundation is restrained from penetration.
The basic reason is that the foundation tip can no longer penetrate into
weakened soil. Hence, obstructions in sand, such as cobbles or boulders, can
quickly initiate local piping.
Ad (d) excessive soil plug heave. This is greater in sand than in clay (for a
given foundation geometry). This is because there is an additional compo-
nent of sand expansion due to stress decreases from upwards flow and other
effects. If high installation pressures are anticipated, and domed top plates
are required, these can take more soil plug heave than flat top plates. In the
Gorm Field caisson installation (Senpere and Auvergne, 1982) excessive
sand heave and piping failure was observed. Excess sand had to be removed
by water jetting to enable the caissons to reach design penetration depth. No
installation details were released – e.g. underpressure magnitudes.
Ad (e) maximum underpressure should not exceed the LAT water depth.
This is a potential problem only in shallow water. European examples
include the southern North Sea, Irish Sea and the Baltic.
134 Installation, retrieval and removal
• instrumentation
• careful installation
• retrieval/reinstallation
• pump flowrate (as high as possible, e.g. Tran, 2005)
• increased preload (gives smaller underpressure)
• reverse flow (i.e. overpressure) at end of installation. This is similar to
hydraulic sand fill placement causing compaction. Reverse flow
pushes back the sand plug and increases plug relative density/strength,
but a potential disadvantage is that it increases/creates a void under
baseplate.
Installation data for fewer than 20 projects are in the public domain (e.g.
Panagoulias, 2016; Alderlieste and Van Blaaderen, 2015). Hence, the models
described here (and associated parameter values) are likely to evolve/
improve as more data become available and additional back-analysis is per-
formed. In addition, attempts are made to numerically investigate this phe-
nomenon; see e.g. Martinelli et al. (2020).
Figure 9.16 compares critical underpressure as a function of embedment
ratio based on published equations. Only two (Houlsby and Byrne, 2005;
Sturm et al., 2015) account for a lower soil plug permeability. All equations
are for a water flow source at seafloor (not for a source at a near-infinite
radius), i.e. they cannot be used for confined sand aquifers.
Installation, retrieval and removal 135
Figure 9.16 C
omparison critical suction pressures (based on Senders and
Randolph, 2009). Critical underpressure, ∆u crit /(γ'D) versus embed-
ment ratio L/D, k r : inner soil plug/outer soil permeability ratio, [–].
where
Fi, suction Fi,no _ flow 1 isand (9.55)
Fo,suction = Fo,no _ flow (9.56)
Q w,suction Q w,no _ flow 1 isand (9.57)
Fi ,no _ flow Asik f ,sand q c, zdz (9.58)
where, noting that subscripts “suction”, “liq” and “no_flow” mean “with
suction assistance and fluid flow”, “required for full sand plug liquefaction”
and “with suction assistance but no fluid flow” respectively:
αu = excess pore pressure ratio at tip level
= Δusuction (tip)/Δusuction
Δuno_flow = water (under) pressure – no fluid flow
Δuliq = water (under) pressure – full liquefaction
Δusuction = water (under) pressure – fluid flow
γsub = sand plug submerged unit weight
γwater = (sea) water unit weight
Aplug = soil plug area
Ap = caisson wall tip area
Asi = caisson inner perimeter
Aso = caisson outer perimeter
D = caisson outer diameter
Di = caisson inner diameter
Installation, retrieval and removal 137
The corresponding pseudo code for the installation in the sand model is:
program install_sand
initialise and evaluate constants
– 2 WT, Asi = π Di, Aso = π D, Ap = π (D2 - Di2) /4 and
Di = D
Aplug = π Di2 /4
∆L = 0.1 m
nL = Lcan /∆L + 1
zero matrix ANS (size nL rows by 21 columns)
k = 0
k = k + 1
L = (k − 1) ∆L
Fi,no_flow: inner skin friction resistance by integrating qc
from 0 to L
Fo,no_flow: outer skin friction resistance by integrating qc
from 0 to L
Fi,no_flow = Asi kf,sand ∫ qc,z dz and Fo,clay = Aso kf,sand ∫ qc,z dz
Qw,no_flow: end-bearing resistance using qc at depth L
Qw,no_flow = Awall kp,sand qc,L
Rno_flow = Fi,no_flow + Fo,no_flow + Qw,no_flow
∆uno_flow = (Rno_flow − Wsub,steel) /Aplug
Fi,liq = 0
Fo,liq = Fo,no_flow
Qw,liq = 0
Rliq = Fi,liq + Fo,liq + Qw,liq
∆uliq = (Rliq − Wsub,steel) /Aplug
αu = αu_fun (L, D)
iterate using ∆usuction to balance resistance Rsuction and
load Lsuction
isand = isand_fun (L, D, ∆usuction, γsub)
Fi,suction = Fi,no_flow (1 − isand)
Fo,suction = Fo,no_flow
Qw,suction = Qw,no_flow (1 − isand)
Rsuction = Fi,suction + Fo,suction + Qw,suction
Lsuction = Aplug ∆usuction + Wsub_steel
exit iteration with ∆usuction, etc., when Rsuction = Lsuction
store in kth row of matrix ANS:
k, L, and L/Lcan in columns 1 through 3
Fi,no_flow, Fo,no_flow, Qw,no_flow, Rno_flow and ∆uno_flow in columns 4
through 8
αu and isand in columns 9 and 10
Fi,liq, Fo,liq, Qw,liq, Rliq and ∆uliq in columns 11 through 15
Fi,suction, Fo,suction, Qw,suction, Rsuction, Lsuction and ∆usuction in
columns 16 through 21
Installation, retrieval and removal 139
Plot col# 07 (Rno_flow), col# 14 (Rliq), and col# 19 (Rsuction) versus col# 02 (L)
resistance – depth
Plot col# 08 (∆uno_flow), col# 15 (∆uliq), and col# 21 (∆usuction) versus col#
02 (L) underpressure – depth
Pile – Soil:
kp,sand, kf,sand = 0.3, 0.001 (BE) and kp,sand, kf,sand = 0.6, 0.003 (HE)
k = 61
L=6m
L/D = 0.6
and
and
Fi,liq = 0 kN
Fo,liq= Fo,no_flow = 6840 kN
Qw,liq = 0 kN
Rliq = Fi,liq + Fo,liq+ Qw,liq = 0 kN + 6840 kN + 0 kN = 6.84 MN
∆uliq = (Rliq − Wsub,steel) /Aplug = (6.84 MN − 6 MN) /77.0 m2 = 10.9 kPa
and
Figure 9.17 shows results for 10 m diameter caissons tipping out in dense
sand (Dr ≈ 80%). Figure 9.17a is for Best Estimate DNV kp and kf values. It
is seen that, if there is no flow, the Best Estimate underpressure ∆u (solid
blue/thin line) at 12 m penetration is moderate, just under 300 kPa. Since
there is no flow, there is no risk of the soil plug liquefying (the major risks
are cavitation if the WD is less than 30 m and steel cylinder implosion if the
Installation, retrieval and removal 141
Figure 9.17 I nstallation in sand. Example using DNV CPT approach. Caisson 10 m diam-
eter. 1.5 m of loose sand (qc = 2z [MPa,m], Dr ≈ 30%) underlain by dense
sand (qc = 10 + 1.5z [MPa,m], Dr ≈ 80%). (a) BE kp and kf, (b) HE kp and kf.
Figure 9.17b shows the corresponding High Estimate resistance and pres-
sure depth diagrams. The situation changes dramatically: full plug liquefac-
tion first occurs at around 11 m bsf (i.e. L/D ≈ 1.1), and the final ∆u ≈ 175
kPa (was ≈ 100 kPa) at 12 m penetration. In such cases, where liquefaction
occurs only for the last metre of penetration, and only for High Estimate kp
and kf values, it is tempting to accept the risk and adopt suitable due dili-
gence procedures (Section 9.11.1) during installation.
9.11.3.2 Commentary
Large diameter caissons (D ≈ 10 m) are necessary to achieve embedment
ratios around L/D ≈ 1 to 1.5 in competent sands. This is largely due to the
fact that load is proportional to D2, whereas resistance is proportional only
to D.
All profiles are straight lines because CPT qc varied linearly in each layer.
In the real world, greater variability is expected.
The model can be modified for non-uniform caissons, other αu and isand
relationships. Unlike the clay installation example (Section 9.10.3), retrieval/
removal has been omitted, and interested engineers can easily include this.
In this case, note that (unlike clay), all Qwall terms practically become zero.
Again, Wsub,steel counteracts overpressure. Hopefully, friction “set-up” is
allowed for by judicious kf selection.
SWP was 1.5 m for both BE and HE assessments. This coincides with the
top surface of the underlying dense sand. Since this depth exceeds (say) 0.5
m, this is satisfactory and is sufficient to start the suction-assisted stage. Had
the surface loose sand layer been absent, then SWP ≈ 0, and there is a high
risk of local piping occurring around the can tip. Special measures (listed in
Section 9.9.2) need to be considered.
With the exception of Lehane et al. (2014), all these studies were for
driven pile (not intermediate) foundations.
9.11.8 General
Installation studies are difficult and opinions/models/results differ widely.
Key challenges usually include:
(a) unit inner skin friction and end-bearing reduction in the sand aquifer
due to upwards seepage in soil plug
(b) water pocket in soil plug at the clay/sand interface. This is because
installation stops if separation occurs and the clay plug meets the
foundation top plate (see Figure 8.2d)
(c) hydraulic fracture of the clay plug (see Figure 8.2f). This is perceived
to be generally associated with competent clay and moderate
underpressures.
146 Installation, retrieval and removal
(a) The clay plug is impermeable. Hence, there is no fluid flow (i.e. no
head drop) within the clay plug
(b) The sand response is fully drained. Hence there are no excess pore
pressures (effective stress = total stress, pore pressures are steady
state).
The objective is to find the sand aquifer flowrate Qflow at any penetration
Lsand. From this, and knowing the time to penetrate the sand aquifer, ΔHpocket,
the required water pocket height, i.e. clay plug heave, ∆Hheave, can be
computed.
The two steady state flow approximate equations are:
iplug 1 – 1 / HLD / Lsand water ubase (9.66)
There are also two equilibrium equations and another two linking sand
resistances to the hydraulic gradient in the sand plug:
(b) forces on caisson steel
utop A plug Wsub,steel Fi,clay Fo,clay Fi,sand Fo,sand Q w,sand (9.67)
Figure 9.19 W
ater pocket. Steady state radial fluid flow – sand plug and aquifer
annulus components are connected with no change in hydraulic head.
Figure 9.20 W
ater pocket. Comparisons Q flow versus Δu base and i plug versus Δu base
approximate relationships with axisymmetric FEA at H sand /D = 0.4
and 1.0. Diameter D = 18.5 m, k sand = 1e-4 m/s.
ubase Q flow / 8 / 13 water ksandHsand / HLD (9.72)
iplug
min 1 1 / HLD / Lsand water ubase ,1 (9.73)
150 Installation, retrieval and removal
If the aquifer flow rate Qflow is negative, then there is no water pocket, and
the no-seepage values should be used, namely:
Q flow = 0 (9.77)
ubase 0 (9.78)
iplug = 0 (9.79)
utop Fi,clay Fo,clay Fi,sand,no _ seep Fo,sand,no _ seep Q w,sand,no _ seep
Wsub,steel / A plug (9.80)
where, noting that subscript “no_seep” means “with suction assistance but
no seepage/fluid flow”:
The corresponding pseudo code for the water pocket model is:
program water_pocket
initialise and evaluate constants; assume DNV CPT qc model
for resistances
F*,clay: skin friction resistances by integrating qc from 0 to
Hclay
Fi,clay = π Di kf,clay ∫ qc,z dz and Fo,clay = π D kf,clay ∫ qc,z dz
penrate,sand,avg: caisson average penetration rate
Aplug = π Di2/4
penrate,sand,avg = Qsuction /Aplug
Wsub,clay = Aplug Hclay γsub,clay
NLsand = 21
zero matrix ANS (size NLsand rows by 10 + 4 = 14 columns)
k = 0
loop tip penetration Lsand from 0 to Hsand in k equal steps:
k = k + 1
Lsand = max(Lsand, 0.001 m)
FI,sand,no_seep and FI,sand,no_seep: skin friction resistances by
integrating qc from Hclay to Hclay + Lsand
Qw,sand,no_seep: end-bearing resistance using qc at depth (Hclay
+ Lsand)
Fi,sand,no_seep = π Di kf,sand ∫ qc,z dz and Fo,sand,no_seep = π D kf,sand
∫ qc,z dz
Qw,sand,no_seep = Awall kp,sand qc,z
αHLD = (1 + 2.5 Hsand Lsand/D2)
Qflow, Δubase, iplug, Fi,sand, Qw,sand and Δutop using Equations (9.71)
through (9.76)
Rno_seep = Fi,clay + Fo,clay + Fi,sand,no_seep+ Fo,sand,no_seep+
Qw,sand,no_seep
Fo,sand = Fo,sand,no_seep
Rsuction = Fi,clay + Fo,clay + Fi,sand+ Fo,sand+ Qw,sand
store in kth row of matrix ANS:
Lsand, Lsand /Hsand, Rno_seep and Rsuction in columns 1 through 4
Qflow, iplug, Δutop and Δubase in columns 5, 6, 9 and 10
Fi,sand and Qw,sand in columns 7 and 8
end loop on tip penetration Lsand
k = 0
loop tip penetration Lsand from 0 to Hsand in k equal steps:
k = k + 1
Qflow,avg = ∑ Qflow (i = 1 to k)
Qflow,avg = Qflow,avg /k
Lsand = ANS(k,1)
tpen,Lsand = Lsand/penrate,sand,avg
Installation, retrieval and removal 153
col# 01: L sand caisson tip penetration into sand aquifer [m]
col# 02: L sand / H sand penetration/aquifer height ratio [–]
col# 03: R no_seep soil resistance without suction assistance [MN]
col# 04: R suction soil resistance with suction assistance [MN]
col# 05: Q flow sand aquifer flowrate (into base of sand plug) [m 3 /s]
col# 06: i plug sand plug (vertical) average seepage gradient [–]
col# 07: F i,sand caisson inner wall friction – sand, suction [MN]
col# 08: Q w,sand caisson wall tip resistance – sand, suction [MN]
col# 09: Δ u,top water (under) pressure at clay plug top [kPa]
col# 10: Δ u,base water (under) pressure in water pocket [kPa]
col# 11: Q flow,avg average flowrate Q whilst in aquifer at L sand [m 3 /s]
col# 12: t pen,Lsand time taken to penetrate sand aquifer to L sand [s]
col# 13: V water water pocket volume at time t pen,Lsand [m 3 ]
col# 14: ΔH pocket water pocket height (= clay plug heave [m]
∆H heave )
Plot col#14 (x, ΔHpocket) versus col# 01 (y, Lsand) water pocket height (=
clay plug heave) – penetration.
Sand permeability can range between ksand = 1e-4 m/s and ksand = 1e-6 m/s;
HE taken.
The sand layer has been analysed at nLsand = 21 equally spaced depth
intervals.
and
Fo,sand = 5.87 MN
Rno_seep = Fi,clay + Fo,clay + Fi,sand,no_seep + Fo,sand,no_seep + Qw,sand,no_seep = 46.8 MN
Rsuction = Fi,clay + Fo,clay + Fi,sand+ Fo,sand + Qw,sand = 23.1 MN
Qflow,avg = 6.43e-3 m3/s
tpen,Lsand = 2270 s
Vwater = 17.95 m3
ΔHpocket = 0.23 m
Figure 9.21 graphs the results. Because of the model complexity, the basic
components are also plotted to gain confidence before discussing the key
plot (ΔHpocket versus Lsand) on Figure 9.21b.
The two right hand underpressure versus penetration plots on Figure 9.21a
show that:
around 100 kPa to just under 200 kPa when the tip penetrates the sand,
increasing to 222 kPa halfway through (Lsand = 3 m) and around 260
kPa at the base of the sand aquifer (Lsand = 6 m). These magnitudes are
normally reasonable for structural engineers. The red/thick lines denote
results had no water pocket been formed. In this case, underpressures
are significantly higher (around 350, 500 and 700 kPa respectively), for
which shell buckling and/or top plate design would be of concern.
• aquifer flow rates, just over 20 m3/h for the first 3 m of sand penetra-
tion, are creating the water pocket. They are less than 10% of Qsuction
(300 m3/h).
Figure 9.21 W
ater pocket example. Foundation D = 10 m, WT = 40 mm. Soil
profile: 4 m of 100 kPa clay over 20 MPa sand, k = 1e-4m/s, Q suction =
250 m 3 /hour: (a) soil resistances and underpressures versus
penetration.
156 Installation, retrieval and removal
Figure 9.21 C
ontinued: Water pocket example. Foundation D = 10 m, WT = 40
mm. Soil profile: 4 m of 100 kPa clay over 20 MPa sand, k = 1e-4m/s,
Q suction = 250 m 3 /hour: (b) water pocket flowrate and height versus
penetration and time.
• because there is a water pocket, iplug, the sand plug seepage gradient is
non-zero and increases with increasing penetration. This is caused by
increasing underpressure. A critical state “liquefaction” (iplug = 1) is
present from just over 5 m penetration into the sand.
• both inner friction and wall tip components decrease due to seepage,
and are zero below 5 m.
• the differential water pressure (Δutop - Δubase) is ≈ 135 kPa and con-
stant. This is as expected – when (Δutop − Δubase) Aplug is the force mov-
ing the clay plug upwards.
Installation, retrieval and removal 157
Figure 9.21 C
ontinued: Water pocket example. Foundation D = 10 m, WT = 40
mm. Soil profile: 4 m of 100 kPa clay over 20 MPa sand, k = 1e-4m/s,
Q suction = 250 m 3 /hour:(c) water pocket flowrate, sand plug gradient,
sand resistances and underpressures versus penetration.
Revised/
Revised original
Original and revised ΔH pocket ΔH pocket
Parameter(s) parameter value(s) [m] [–] Remarks
Figure 9.22 R
everse punch-through failure model: Forces acting on the clay plug.
No reverse-end-bearing (Cotter, 2009).
The model assumes zero reverse end-bearing (REB), and no water pocket
development, at the clay/sand interface. There are two equilibrium equa-
tions: forces on the caisson steel (during installation) and forces on the com-
plete clay plug (reverse punch-through). Note that Wsub,steel assists penetration
but Wsub,clay,RPT opposes reverse punch-through. From these equations, the
corresponding underpressures Δu can be computed. Reverse punch-through
failure is possible when ΔuRPT > Δuinstall. The four equations are:
uinstall A plug Wsub,steel Fi,clay,adhesion Fo,clay,adhesion Q w,clay (9.83)
where
σvo = total soil vertical stress at depth z = (γsub,clay + γwater) z
γsub,clay = clay submerged unit weight
γwater = (sea) water unit weight
Aplug = soil plug area
Ap = caisson wall tip area
As,i = caisson inner perimeter
As,o = caisson outer perimeter
Fi,clay,adhesion = caisson inner wall friction – clay adhesion above tip
Fo,clay,adhesion = caisson outer wall friction – clay adhesion above tip
Fi,clay,shear = caisson inner wall friction – clay shear below tip
Hclay = clay layer and clay plug thickness
kf,clay = DNV empirical coefficient – unit skin friction resistance in clay
kp,clay = DNV empirical coefficient – unit end-bearing resistance in clay
160 Installation, retrieval and removal
where
grad1 = α π (Hclay/D) + π Nc,strip (WT/D)
intercept1 = (Hclay/D) (π WT/D - π/4)
X = Wsub,steel/(γsub,clay D3)
Y = su,av/(γsub,clay D)
su,av,Hclay = average su in clay layer
Nc,strip = bearing capacity factor at caisson tip
Due to assuming Δuinstall acts over π D2/4, not π Di2/4, equation 9.12.19 is
slightly approximate. Reverse punch-through will not occur if yfRPT < 0.
The calculations to find and compare underpressures Δuinstall and ΔuRPT are
straightforward. At Lclay = 2 m (half way through the clay), we have:
(a)
(b)
Figure 9.23 R
everse punch-through example. Foundation D = 10 m, WT = 50 mm,
W sub = 6 MN. Soil profile: 4 m of ≈ 70 kPa clay over sand. Soil resistances
and underpressures versus penetration: (a) High Estimate k p, k f = 0.6,
0.05. Reverse punch-through from 3.98 m and SWP = 1.47 m, (b) Low
Estimate k p, k f = 0.4, 0.03. No reverse punch-through. SWP = 2.48 m.
9.11.10.2 Commentary
Based on these results, the probability of reverse punch-through occurring is
extremely small. The model excluded water pocket type fluid flow effects in
the sand. If they had been included, then increases would occur in the down-
wards-acting force on the punched-through plug base (currently assumed to
Installation, retrieval and removal 163
be zero) and hence also increase the underpressure ΔuRPT. This would make
reverse punch-through even more unlikely.
The reverse punch-through example presented here was for competent
clay. Had weak NC clay been studied, then (as with using lower kp, kf coef-
ficients in the example) the increased SWP values would decrease Δuinstall and
the likelihood of reverse punch-through. This exercise is left to the interested
reader.
9.12.1 Impact driving
Driving is not the most common way to install a foundation in a full massive
of weak rock. Traditionally, weak rock concerns have been encountered
mainly for jacket pile or anchor pile installation in cemented carbonate sand
(mainly Persian Gulf and Australian projects) or in chalk (CIRIA, 2002).
Case studies related to pile driving in other types of weak rock refer mainly
to driving through limited thicknesses of few meters (e.g. Puech et al., 1990);
see Figure 9.24.
The effect of impact driving on the rock mass will depend on rock char-
acteristics (e.g. strength, porosity, mineralogy and in-situ fractures), but also
on the foundation geometry. Stevens et al. (1982) considered that driving in
rock will severely fracture the rock layers and reduce the rock to granular
material. The shaft friction used in soil resistance to driving predictions
(SRD) is therefore defined using sand parameters and end-bearing propor-
tional to the rock unconfined compressive strength (UCS).
In chalk, the above approach proposed by Stevens et al. (1982) proved to
be inappropriate (Dührkop et al., 2017 and Wood et al., 2015) for low to
medium density chalk. A zone of remoulded material is observed around the
pile annulus when driving open-ended piles in these types of chalk. SRD is
linked to the grade and density of the chalk and needs to account for “fric-
tion fatigue” (unit skin friction decreases with increasing penetration) and
set-up effects. Figure 9.25 shows evidence of remoulding and fracturing
around a pile after driving in chalk.
In practice, the ground will need to accommodate the volume of steel pen-
etrating. At the start of driving, the easiest way will be to use the available
void above the ground, leading to ground uplift around the foundation and
Figure 9.25 C
halk-pile interface model – St. Nicolas pile tests, D = 0.76 m,WT = 44.5 mm
(Muir Wood et al., 2015).
Installation, retrieval and removal 165
creating wedges of fractured material. When the pile starts to penetrate, the
ground will accommodate the steel penetration through a combination of:
• Rock crushing around the pile wall. Depending on the rock type, the
crushed material will be either in a pasty or remoulded state (e.g. for
grade C and D chalk) or as a granular material.
• Radial fractures up to a few foundation diameters from the founda-
tion edge. In the presence of natural fractures, blocks will be pushed
radially and fractures will tend to close.
9.12.1.2 Risk of buckling
Open-ended piles are vulnerable to tip distortion. An initial imperfection
from the theoretical pure cylindrical geometry can progressively grow with
increasing pile penetration when the stiffness of the surrounding soil exceeds
the elastic pile stiffness (Figure 9.27). Care is required during fabrication,
transportation and installation to avoid creation of an initial imperfection
166 Installation, retrieval and removal
Figure 9.26 I chthys pile tip stiffening arrangement, D/WT ratio over 110 (Erbrich
et al., 2017).
that can lead to buckling. Aldridge et al. (2005) present equations, and both
Erbrich et al. (2010) and Aldridge et al. (2005) describe conditions for prop-
agation of an initial imperfection during pile penetration (i.e. extrusion
buckling). The propagation of a small initial deformation requires that the
soil be both strong enough and stiff enough to deform the pile further.
For monopile or large diameter anchor piles, the ratio between the diam-
eter and the pile tip thickness (i.e. D/WT ratio) is becoming very high com-
pared to classical flexible piles (see Figure 9.28). Figure 9.13 shows that for
D/t greater than 40, buckling was observed. Risk of insufficient hoop stiff-
ness or heterogeneous conditions across the diameter are significantly
greater for stubby large diameter foundation. The ring stiffness is inversely
proportional to cube of D/t ratio.
Installation, retrieval and removal 167
Figure 9.27 P
lastic hinge mechanism assumed to cause local tip buckling (Aldridge
et al., 2005).
Figure 9.28 P
ile tip geometries, oil and gas full blue dots and monopiles empty
dots, red circles correspond to recorded tip damage (Randolph, 2018).
168 Installation, retrieval and removal
9.12.2 Drilling
9.12.2.1 Drive drill drive
Drilling can be used as a contingency measure when driving is expected to
be difficult (as the drilling tool is generally already mobilised together with
the hammer) or in case of unexpected pile refusal (e.g. presence of boulder).
After the end of driving, the drilling head is inserted inside the pile and drill-
ing starts from the seafloor. Drilling can either stop after the particularly
difficult layer (or boulder) or be continued to the pile target depth. In the
latter case, hole stability needs to be guaranteed. The driving process is then
continued until reaching the target depth or the next difficult layer.
Drive Drill Drive (DDD) operations often require under-reaming capabil-
ity, i.e. the possibility for the drilling head to enlarge its diameter once it
reaches the pile tip. Drilling more than the pile diameter removes resistance
beneath the tip in a particularly hard layer (or boulder).
The switch of tool between the hammer and drilling tool is time consum-
ing. Drive Drill Drive (DDD) installation is often scheduled only if one pass
Installation, retrieval and removal 169
Figure 9.29 D
rilled and grouted installation sequence with sacrificial casing
(Courtesy of Fugro).
170 Installation, retrieval and removal
In all cases, the time during which the hole is kept unfilled needs to be
minimised.
Unfortunately, there are not many references applicable for shallow hole
stability for foundations (i.e. low drilling fluid pressure due to relatively low
depth relative to e.g. well-boring). Most of the published work covers well-
bore drilling, tunnelling or mining activities. For cylindrical geometry, zero
drainage boundary conditions and homogenous rock mass, hole stress con-
ditions can be derived based on elastic solutions.
For weak rock where discontinuities are present, the complexity of the
stress distribution and geometry of the discontinuities requires FE (Finite
Element) or DE (Discrete Element) methods. It should, however, be noted
that FE and DE methods will not capture the chemical interactions and time
dependency.
The main difficulty remains to get a detailed enough idea of the geometry
of nearby discontinuities. Offshore, the available information is often lim-
ited to the logging of joints and fractures on 100 mm samples from a single
borehole. It might be necessary, where faults are prevalent, to perform
Installation, retrieval and removal 171
several boreholes per foundation, enlarge the depth of boreholes (to com-
pare geological markers) and/or consider downhole logging such as natural
gamma-ray, televiewer and calliper. These tools can provide valuable infor-
mation on the nature, frequency and inclination of discontinuities (crossing
the borehole).
During grouting operations, it will be necessary to avoid a grout flow
inside the foundation. This can be done by either creating a grout plug at the
bottom of the foundation before filling the annulus, through a physical
boundary (e.g. internal membrane) or through a mechanical closure (e.g.
using external packers).
In the presence of fractured rock or karstic voids, the volume of grout to
be injected can largely exceed the annulus volume (volume between the
foundation and the rock). Grouting operations offshore remain a delicate
operation that is difficult to control.
9.12.3 Vibratory
Pile installation using vibratory means would not be considered if the stra-
tigraphy includes cemented material or hard (even small) layers.
Figure 9.10 has given an example for an anchor pile in LOC clay. Figure 9.17
presents a typical example for caisson installation in sand. Note the
following:
9.14.1 Suction foundations
Procedures for retrieval and removal are similar to those for installation
(Section 9.2 on penetration resistance assessment), except that:
9.14.2 OWT monopiles
Driven monopiles cannot be removed in their entirety – their pull-out capac-
ity (resistance plus submerged weight) is significant. For example, a 5 m
diameter 30 m embedded length monopile pull-out capacity could easily
exceed 60 MN, depending on the exact soil profile and other factors. There
are, however, on-going developments to investigate removal of such piles,
see e.g. Balder et al. (2020). Marine legislation generally stipulates that piles
be cut off at a certain depth (typically between 3 m and 5 m) below seafloor.
This is generally achieved by explosives, mechanical or abrasive cutters.
Jetting may be necessary to remove the soil beforehand. Pile cutting may be
either internal or external; usually the former approach is selected, in view
of the soil volumes to be excavated. Grout overflow or other seafloor
obstructions may hinder access.
Drilled-and-grouted monopiles cannot be removed entirely either. For the
decommissioning phase, the requirement is generally to cut the monopile 1
m below seabed and to fill it with rock. Only a few offshore windparks have
already been decommissioned. Few relatively small drilled and grouted
Installation, retrieval and removal 173
9.15 CLOSURE
In-place resistance
10.1 INTRODUCTION
This chapter covers in-place resistance for both support type and anchor
pile type intermediate foundations.
Figure 2.2 (Section 2.2, on in-place resistance modes) showed six failure
modes for intermediate foundations under combined VHM loads in und-
rained soil. All assume co-planar HM (horizontal load and moment) and
negligible T (torsion). Figure 2.3 showed the corresponding failure modes
for drained soil response. All modes are suitable for preliminary/routine
design calculations, but possibly need to be complemented by FEA for final
design.
Figures 2.2 and 2.3 are both for combined VHM loading. Figure 10.1
(Section 10.5 on maximum axial resistance) gives the corresponding failure
models for pure axial (pull-out) V load.
Resistance envelope models for VHM(T) loading are discussed further in
Section 10.10.
175
176 In-place resistance
Tension cracks (also known as gaps) can occur behind the trailing edge
(active side) of foundations under lateral load. This has been realised as a
potential problem for sheet pile walls for many years. The crack depth
depends upon many factors, including the lateral load magnitude and the
in-situ stress profile.
The occurrence of either tension cracking or gapping reduces VHMT
resistance. For offshore foundations, if a tension crack develops to skirt tip
level, then in-place undrained resistance capacity at tip level (including
reverse end-bearing) may be severely affected by loss of suction in fine-
grained soil. In addition, in the zone above the tension crack, a small reduc-
tion in bulkhead net resistance (passive minus active) occurs.
Tension cracking is not generally considered for normally consolidated
clay subjected to short-term loads (Randolph and Gourvenec, 2011).
Jeanjean et al. (2006) found no evidence of gapping in more competent OC
clay (su ≈ 25 kPa at 9 m depth). Gapping is, however, expected to be non-
reversible in cemented formations or rock. In the case of drilled and grouted
piles, the gap will most likely occur between the steel pile and the grout
annulus. This gap should normally be limited to the upper 1–2 meters due
to the limited displacements.
For intermediate anchor foundations in normally consolidated clay, ten-
sion cracking is not usually considered. This is because locating the lug
slightly below optimum depth causes restoring loads and decreases the risk
of forward rotation. However, tension cracking and gapping may need con-
sideration for intermediate support foundations under overturning loads
and/or in stiff clays.
Mitigation measures for foundations include a flexible mat (a “gap
arrester” or “mud liner”) around the perimeter (e.g. Keaveny et al., 1994;
Mana et al., 2013). This inhibits a tension crack by preventing water supply,
and hence suction loss, at the soil–foundation interface.
Figure 10.1 F ailure models for maximum axial tensile resistance, V max , in und-
rained soil. Reverse end-bearing is included only in the “plugged”
model (Senders and Kay, 2002).
Figure 10.1 shows the three possible failure modes – “plugged”, “coring”
and “leaking” – under tensile V loads. Identification of these modes have
been based on previous work. Fuglsang and Steensen-Bach (1991) reported
results of 1 g model and 40 g centrifuge pull-out tests on 65 mm and 80 mm
diameter L/D = 2 caissons in kaolin clay. To derive end-bearing and inner/
outer skin friction components, they used three model caisson types:
“plugged” (Fo + Qbase), “coring” (Fo + Fi) and “coring without the clay plug”
(Fo). They verified that there was indeed reverse end-bearing failure. In sev-
eral 1 g “plugged” tests, but none of the centrifuge tests, brittle tensile failure
was observed – the clay at the caisson tip abruptly failed in tension. This was
possibly the “leaking” mode – the tests were carried out underwater.
Randolph and House (2002) considered the same three “plugged”, “cor-
ing” and “leaking” failure modes as Fuglsang and Steensen-Bach (1991),
using instead the term “base-vented” when a hydraulic short circuit (i.e.
“leaking”) prevents underpressure at the caisson base. For the “coring”
mode, reverse end-bearing resistance (Qwall) under the caisson wall is
included. Although admittedly small, it is optimistic to include this compo-
nent. Senders and Kay (2002) presented the same three models, but with
component Qwall omitted from the “coring” mode. The following paragraphs
give fuller details and update the aforementioned work.
In the “coring” failure model, there is no reverse end-bearing under the tip
wall, and axial resistance (Vmax) equals outer friction (Fo) and inner friction
(Fi). The corresponding Vload (axial pull-out load) equals Fo + Fi + Wsteel. In
the “leaking” failure model, Vmax simply equals outer friction (Fo), and Vload
(axial pull-out load) equals Fo + Wsteel + Wplug.
Hence, “coring” Vmax and Vload are given by the equations:
Vmax Fo Fi (10.3)
Vmax = Fo (10.5)
Note that for a normally consolidated clay (su/σ'vo≈ 0.25) with clay-steel
inner adhesion factor αi ≈ 0.6, “leaking” is more critical than “coring”
failure when embedment ratio L/D ≈ 3. That is, if “leaking” failure can
occur, then “coring” resistance is probably governing for short, stubby
foundations, whereas “leaking” is more likely to be critical for longer
foundations.
Finally, for routine FEA, it is very difficult to model pull-out using either
a “plugged” failure model with a user-assigned Nc value (i.e. reduced REB;
Clukey and Morrison, 1993; Clukey et al., 1995), or a “leaking” failure
model. In the remaining “plugged” and “coring” cases, design practice
should be to obtain and interpret the sum of soil forces along the inner and
outer shaft areas (i.e. mobilised friction Fo and Fi) and over the base area (i.e.
mobilised bearing resistance Qbase). If “plugged” failure is occurring, then
one would expect Vload ≈ Fo + Qbase and Vload ≈ Fo + Fi if “coring”. For the
“plugged” condition, the back-figured Nc,circle value should be compared
with the Skempton (1951) value for circular foundations. In particular,
Nc,circle is expected to be in excess of 9 for L/D = 2.5, for the reasons given in
Section 10.5.2.
180 In-place resistance
Failure Mode
Note(s): H = M = 0.
182 In-place resistance
Note that these L values are solely for pure V load. Placing the lug level at
optimum depth below seafloor for resisting combined load gives signifi-
cantly smaller L values – between 12 m and 14 m (see the Example in
Section 10.7).
N c,strip 7.5 for wall end bearing strip (Skempton, 1951) (10.7)
N c,circle min 6 1 0.2 L / D , 9 (Skempton, 1951)) (10.8)
N c,circle min 6.2 1 0.34 arctan L / D , 9 (DNV, 2018) (10.9)
• For (reverse) end-bearing over the full foundation base area, a reason-
able approach is to take Nc = Nc,circle and su averaged over (say) 0.25 D
below tip level. This approach is consistent with that used for jack-up
spudcan leg penetration assessments, see. ISO 19905-1:2016 (ISO
2016b), Section A.9.3.2.2. Other organisations may have different
opinions regarding averaging depth etc.
• For (reverse) end-bearing over the foundation wall tip area, take Nc =
Nc,strip and su at tip level.
• Other methods are available for shallow foundations; see Section
3.2.1 on in-place resistance modes.
In-place resistance 183
Note that it is cautious to assume that bearing capacity factor Nc,circle has a
maximum value of 9 for closed-ended and plugged open-ended piles. See
Section 10.5.5 on axial myths.
Commentary:
Table 9.6 gives values of ISO 19902 pile bearing capacity factor Nq.
10.5.4.2 Tensile loads
Assuming a “coring” failure mode, and ignoring reverse end-bearing below
the wall tip, then Vmax is given by:
Note that friction resistance for intermediate foundations using ISO 19901-4
(2016) for driven piles is considered to be optimistic. The reasons for
decreased friction include:
10.5.5 Axial myths
10.5.5.1 Bearing capacity factor Nc
It is a geotechnical myth that Nc,circle has a maximum value of 9 (Skempton,
1951) for closed-ended and plugged open-ended piles. This value is cau-
tious. The DNV (2005) equation is merely a “better” (non-linear) fit to
the original Skempton (1951) Figure 2 curve. Skempton’s work was based
on theoretical considerations, small-scale laboratory, and two field tests
on screw cylinders/screw piles. All these boundary conditions are not
In-place resistance 185
N c , circle Reference
VHM loads are at seafloor (see Figure 3.1a). Care should be taken in distin-
guishing between free and fixed head resistances H: “Fixed head” Hmax
occurs when the foundation cannot rotate (rotation θxz = 0) and V = 0. “Free
head” Ho is when V and M are both zero. Figure 10.3a shows their positions
in the MH plane. Normalised Hmax and Ho values versus foundation embed-
ment ratio for various undrained shear strength profiles (Figure 10.3b) are
given in Figure 10.3c and 10.3d. For embedment ratios L/D > 3, reasonable
Np values are 10.5 and around 3 for most clay profiles, and hence the ratio
Hmax/Ho ≈ 3.5.
(a)
Figure 10.3 M
aximum lateral resistance in clay. Normalised resistance factor N p
= H/(LDs u,av,L ) (a) MH plane, lateral resistances “fixed head” H max and
“free head” H o (Kay and Palix, 2010).
188 In-place resistance
(b)
Figure 10.3 C
ontinued: Maximum lateral resistance in clay. Normalised resis-
tance factor N p = H/(LDs u,av,L ) (b) Soil profiles.
(c)
(d)
Figure 10.3 C
ontinued: Maximum lateral resistance in clay. Normalised resis-
tance factor N p = H/(LDs u,av,L ) (c) Horizontal capacity of suction
caissons (Supachawarote, 2006), N p versus L/D (α o = 0.65) (Kay,
2015).
190 In-place resistance
These paragraphs are about suction anchor pile lateral capacity. These are
essentially rigid caissons. Analyses of laterally loaded piles are not discussed.
The terms anchor chain and anchor line are synonymous.
Lateral performance, particularly moment failure of the steel at design
loads, may be governing (rather than soil failure). Providing additional steel
stiffeners at or near the lug level should prevent this. Note that additional
steel increases installation resistance.
Unless load and chain inclination is provided at lug level, a separate anal-
ysis is required of the anchor chain below seafloor. This is to assess the
decreased load and increased uplift angle, θlug, at lug level (see Figure 3.1b).
Note that competent soils, such as sand and stiff clay, cause significant
changes to the chain profile, particularly a further θlug increase. In such cases,
the lug level is usually restricted to relatively shallow depths below seafloor
in order to avoid axial pull-out. Equilibrium equations for an embedded
anchor chain have been given by Degenkamp and Dutta (1989) and
Neubecker and Randolph (1995).
If there is no restriction about lug depth, then maximum horizontal resis-
tance (Hmax) is obtained such that, in the ultimate limit state, the caisson
behaves as a “fixed head” body, translating laterally with zero rotation as
shown on Figure 8.1b. In this case, optimum lug level (on the caisson shaft) is
determined by the intersection of the line of action of the lug chain load and the
caisson centreline. In practice, the lug level is just below optimal (say, around
0.5 m for a 3 m diameter caisson). This ensures backwards rotation at failure,
reducing potential for a tensile crack to open on the caisson trailing edge.
The optimum lug level is close to the weighted average of the shear
strength profile, and the general expression for any shear strength profile is:
L L
ez
0
su, z z dz / su, z dz
0
(10.13)
where
L = caisson embedded length
D = caisson outer diameter
su,av,L = average su,z between caisson head (seafloor) and caisson tip (L)
and
Example – Fixed Head Pile. Use the same data as used previously in Vmax
example (Section 10.5.1), basically a 5 m diameter anchor pile to withstand
7500 kN factored load at lug level (Tlug) acting at 30° (θlug). Additional data
included:
Figure 10.4 Effect of anchor lug level on lateral resistance (Randolph and
Gourvenec, 2011).
Note that these L values are for combined VHM load. Putting the lug level
at seafloor, under pure V load, gives significantly higher L values – between
15.6 m and 22.3 m (see the Example in Section 10.5).
Figure 10.5 shows the corresponding three Vlug-Hlug envelopes, together
with the actual Tlug load components, namely Vlug = Tlug sin(θlug) and
Hlug = Tlug cos(θlug). The three shapes are dissimilar to V-Hmax diagrams – see
In-place resistance 193
Notes:
zlug,opt,fix: caisson optimum lug level depth (below seafloor)
Lopt,fix: caisson optimum embedded length (“fixed head” condition)
Figures 8.9 (ISO 19901-4) and 10.12 (Supachawarote et al., 2004). For
example, taking the “plugged” failure model, (a) Hlug first increases, and
then decreases, as θlug decreases from 30° to 0° and (b) Vlug first increases,
and then decreases, as θlug increases from 30° to 90°. The latter is beneficial:
Tlug exceeds 7500 kN should θlug exceed 30°..
The Vlug- Hlug envelope shape is the locus of points where Vlug and Hlug
loads (at constant rlug, zlug, varying θlug) intersect the “tongue”-shaped VHM
resistance envelope in 3D VHM space shown in Figure 2.2e. This does not
have the same shape as ellipsoidal V-Hmax diagrams, which are vertical slices
through the same envelope at Hmax.
where
Nc,prot = bearing capacity factor
In-place resistance 195
For “undrained” clay profiles, usually the “plugged” (not “coring” mode)
is critical, and hence Tmax = Touter + Tbase
These paragraphs are about intermediate anchor pile design for tilt (non-
vertical or out-of-plumb) and twist (misalignment or misorientation). Both
reduce holding capacity compared to a perfectly installed and orientated
anchor. For resistance calculations, the most unfavourable situation (i.e.
maximum tilt and maximum twist occurring at the same time) is often stud-
ied in detailed design. Typical specifications for installation tolerance are
±5.0° for pile tilt and ±7.5° for pile twist. Generally, maximum tilt (+5.0°)
reduces capacity more than maximum twist (+7.5°).
Tilt is accounted for by increasing or decreasing the chain load angle at
lug level (θlug) by the tilt angle. This approach is similar to the common
assumption for inclined piles, namely that the axial and lateral responses are
largely unaffected by small angles of inclination from the vertical (e.g.
Poulos and Davis, 1980). Adding (e.g.) 5° tilt to θlug increases the V load
component on the anchor, and subtracting 5° increases the H load compo-
nent. Usually, anchors are more sensitive to V (than H) loading, and hence
+5° tilt is more unfavourable than −5°. A non-horizontal seafloor, especially
if it is sloping down in the same direction as the mooring line(s), should be
considered. A reasonable procedure is to take the design tilt to be the greater
of slope angle and tilt tolerance angle. For example, if the seafloor slope
angle = 10°, and tilt tolerance = ±5.0°, then take design tilt = +10°.
Twist influences suction anchor foundation resistance less than tilt. This is
seen in typical tolerance values: namely < ±5.0° for tilt and < ±7.5° for twist.
Due to twist, part of the maximum available friction resistance will be used
in torsional shear when the chain pulls on the suction anchor. The basic
principles are given in Figure 10.6. Assuming the torsion load is fully applied
first, a simple calculation can be done: knowing the applied torsion load and
the maximum available resistance (on the suction anchor perimeter and
base), the amount of shear utilised can be calculated. Vector analysis then
gives the amount of shear available, and hence the reduction percentage, on
the maximum frictional resistance for axial loading. To simplify design cal-
culations, the pile-soil outer adhesion factor (αo) for axial loading can there-
fore be decreased by the reduction percentage. Typically, this αo decrease is
small – around 5% (Senders and Kay, 2002). In the worked example that
follows, the αo decrease is around 2.5% (was 0.65, now 0.635). Interested
196 In-place resistance
Figure 10.6 E ffects of twist on suction anchor holding capacity (Taiebat and
Carter, 2005).
engineers may care to repeat the example calculation for 15° (was 7.5°)
twist. This occurred at Na Kika, Gulf of Mexico, details of which were given
in Section 6.1.2. Large diameter foundations (D > 15 m, say) will probably
give αo decreases > 5%, and, if Tload/Tmax is large, Tload – VHM resistance
interaction effects (Section 10.10) may need to be considered.
10.9.1 Example – twist
To illustrate the difference twist (misalignment) makes to clay-steel outer
adhesion factor αo, consider a typical anchor pile in a NC soil profile. The
relevant data are:
Pile: D = 5 m, L = 30 m
Soil: su = 2 + 1.5 z [kPa, m], αo= 0.65
(su,av,L = 2 kPa + 1.5 kPa/m 15 m = 24.5 kPa, su,tip = 2 kPa + 1.5 kPa/m
30 m = 47 kPa)
Lug: Tlug = 3 MN, θlug = 30°
Twist: ψlug = 7.5°
The calculations to find αo* (i.e. adhesion αo allowing for torsion load Tload)
are straightforward. They are as follows:
For routine design, these two factors are considered using numerical analy-
sis models, either based on limit equilibrium and plastic limit analysis
(Section 10.11) or a resistance envelope approach (Section 10.12). The 3D
FEA is generally not used for preliminary (and the majority of routine) geo-
technical design of intermediate foundations.
For routine design, it is considered reasonable to assume that Tload
affects only axial V (not lateral HM) resistance. Justification for this is
given by Taiebat and Carter (2005) – as shown in Figure 10.7 – Saviano
and Pisanò (2017) – shown in Figure 10.8 – and also Suroor and Hossain
(2015), not shown.
198 In-place resistance
(a)
(b)
Figure 10.7 E ffect of torsion (T) load on vertical (V) and horizontal (H) resis-
tance (Taiebat and Carter, 2005).
10.11.1 General
Besides VHM(T) coupling, base shear and possible tension cracking need
to be considered – this is as for shallow foundation design, but unlike pile
foundation design. Unlike shallow foundations in clay, there is not gener-
ally a need to consider gap formation below the foundation base. Analysis
methods include Limit Equilibrium, Plastic Limit and Resistance Envelope
(aka Yield Function). The latter is a comparatively recent development in
geotechnical engineering.
(d)
Figure 10.8 E ffect of torsion (T) load on vertical (V) and horizontal (H) resis-
tance (Saviano and Pisanò, 2017): (a, b, c) effects of misalignment on
VH capacity envelope, (d) HT capacity diagram.
same design soil profile – and (b) differing vertical and horizontal Factors of
Safety on lug load Tlug.
Table 10.5 compares intermediate foundation resistance software func-
tionality. Only CAISSON_VHM uses a yield function (details given in
Section 10.12), and additional functionality includes:
• caisson tilt inclusion by adding tilt angle to θlug – see Section 10.9
• axial resistance Vmax using either (a) “plugged”, “coring” or “leaking”’
failure models (Section 10.5.2) or (b) user-defined
• caisson buoyant weight Bw using either WT/D ratios or user-defined
• embedded chain profile using Degenkamp and Dutta (1989); these
chain element equations were first stated by Vivatrat et al. (1982), and
are identical to those given in Appendix A3 of DNV EP-R301 (2000)
• separate λT load factors on V and H components of lug load Tlug
• resistance modifications – see Section 10.12.7
• caisson twist (misalignment): reductions are usually small, less impor-
tant than tilt (Senders and Kay, 2002), and can be simply done by
decreasing αo – see Section 10.9.
Table 10.5 C
aisson In-place Resistance Software – Comparison Functionality
(Kay, 2015).
Software
Functionality A B C D E F
Solution method PL ULS ULS YF FE 2D FE 3D
Support foundations ⨯ ✓ ✓ ✓ ✓ ✓
Anchor foundations ✓ ⨯ ✓ ✓ ✓ ✓
Chain profile ⨯ ⨯ ✓ ✓ ⨯ ⨯
Rotational failure ✓ ✓ ✓ ✓ ✓ ✓
Zero rotation failure ✓ ⨯ ✓ ✓ ✓ ✓
Internal scoop failure ⨯ ✓ ✓ ⨯ ⨯ ✓
Tension crack ✓ ✓ ✓ ⨯ ✓ ✓
High axial load ⨯ ⨯ ✓ ✓ ✓ ✓
Reduced REB ✓ ⨯ ⨯ ✓ ⨯ ⨯
L, D optimisation ⨯ ⨯ ⨯ ✓ ⨯ ⨯
Multiple L,D ⨯ ⨯ ⨯ ✓ ⨯ ⨯
Multiple load cases ⨯ ⨯ ⨯ ✓ ⨯ ⨯
Multiple soil layers ⨯ ✓ ✓ ✓ ✓ ✓
Tilt ✓ ⨯ ✓ ✓ ✓ ✓
Twist ⨯ ⨯ ✓ ✓ ✓ ✓
Soil reactions ⨯ ⨯ ⨯ ⨯ ✓ ✓
Notes:
A: AGSPANC (AG, 2003); FALL16 (OTRC, 2008); B: CANCAP2 (Fugro, 2009); C: CAISSON (Kennedy
et al., 2013); D: CAISSON_VHM (Kay, 2015); E: BIFURC (Jostad, 1997); SPCalc (XG Geotools, 2014),
F: ABAQUS, PLAXIS 3D etc.
ULS: limit equilibrium (neither LB nor UB); PL: plastic limit analysis (usually UB);YF: yield function (based
on FE); FE: finite element analysis (slight UB), either 2D or 3D; LB: lower bound; UB: upper bound.
In-place resistance 201
10.12.1 General
This section can be skipped by engineers uninterested in foundation opti-
misation. The majority of this section presents methods for undrained
soils. This is because, as noted, a generic envelope for intermediate (and
also for deeply embedded surface) foundations in drained sand under
VHM loads is not yet available.
Table 10.6 R
esistance Envelopes for Intermediate Foundations in Undrained
Soil (Clay).
Note: NC, uniform, stepped and crust soil su profiles given on Figure 10.3b.
based on around 5500 quasi 3D (2D Fourier) finite element analyses for a
large number of caisson L/D values (0.5 ≤ L/D ≤ 6) and soil su profiles
(Figure 10.3b).
The resistance envelope is “tongue”-shaped (see Figure 10.9) in VHM
space, and consists of a rotated ellipse in M-H space (Figure 10.9b), and
recovers the familiar ellipsoidal shape of VH diagrams for shallow founda-
tions in V-Hmax space (Figure 10.9c). The complete envelope is described by
three simple equations, and all parameter values are smooth functions of
foundation embedment ratio (L/D) and non-dimensional soil shear strength
profile (ez,su/L). For a multi-layered soil profile, interpolation is used between
the su profiles shown in Figure 10.3b. Major assumptions included no ten-
sion crack and no internal soil scoop failure. Torsion loading, usually small,
can be accounted for by decreasing su (Senders and Kay, 2002; DNV, 2005).
Benchmarks (including Andersen et al., 2005; Randolph and Houlsby, 1984;
Taiebat and Carter, 2010), plus design optimisation examples for support
and anchor foundations, are given by Kay and Palix (2010, 2011). Van Dijk
(2015) gave resistance envelope equations for open-ended caissons, consid-
ering internal scooping.
The results can be applied to any type of rigid intermediate foundation –
either support or anchor pile. For anchor foundations, since the VHM load
reference point is on the centreline at seafloor (see Figure 3.1), and the foun-
dation is essentially rigid, chain loads and inclination angle at lug level (Tlug
and θlug) need to be transformed to seafloor VHM load using the three
equations
In-place resistance 203
(a)
(b)
(c)
V = Vlug (10.17)
H = H lug (10.18)
where
Vlug = Tlug sin(θlug)
Hlug = Tlug cos(θlug)
rlug = lug radial offset from foundation centreline
zlug = lug level depth below seafloor
For θtilt = +5°, the seafloor V and H loads remain unchanged, cos(θtilt) = 0.996
and sin(θtilt) = 0.087, and the transformed lug coordinates are given by
10.12.3 MH ellipses
For intermediate foundations, it is convenient to work in terms of dimen-
sionless M and H values given by the equations (Figure 10.10)
H H / L Dsu,av,L (10.22)
M M / L2Dsu, av, L (10.23)
Ad (ii): since Hmax lies close to the semi-major axis (see Figure 10.11)
Hence, one can obtain ΦMH, aHM and bHM values from
10.12.4 V-Hmax ellipsoids
Senders and Kay (2002), Supachawarote et al. (2004) and Kay and Palix
(2010) all found that V-Hmax resistance envelopes for caissons in clay (L/D
in the range 0.5–6.0) were ellipsoidal and could be fitted with equations of
the type
206 In-place resistance
Figure 10.10 M
H ellipses at V = 0. Summary final data points, caisson D = 5 m,
L/D = 1.5, 2, 3, 4, 5 and 6. Soil s u profiles: (a) and (b) Constant, e z,
s u /L = 1/2, (c) and (d) NC, e z,su /L = 2/3 and (e) and (f) Stepped, e z,su
/L = 3/4. (Kay and Palix, 2010)
In-place resistance 207
Figure 10.11 M
H ellipse resistance envelope. Key data points H 0 , M 0 and H max
(Kay and Palix, 2010).
aVH bVH
H max,V V
1 (10.27)
H max Vmax
Figure 10.12 shows dimensionless V-Hmax envelopes for different L/D val-
ues. As L/D increases, the less V-H interaction occurs, i.e. circles become
more square. Hence, ellipsoid curvature parameters aVH and bVH are again
functions of L/D and the su profile.
For example, for normally consolidated clay, Supachawarote et al. (2004)
gives
The high bVH value implies that vertical resistance is less affected by hori-
zontal load than vice versa.
208 In-place resistance
Figure 10.12 V
-H max resistance. Dimensionless envelopes for differing caisson
embedment ratios, L/D. Normally consolidated clay (Supachawarote
et al., 2004).
Mload L H
2
max,V aMH cos t sin MH bMH sin t cos MH (10.31)
L Dsu, av,L H max
2
cos MH H sin MH M
fVHM
aMH
2
1 (10.34)
bMH Vmax
The first two terms of the yield function, Equation (10.34), are independent
of V. They represent the geometrical form of a rotated ellipse in non-dimen-
sional (M*, H*) space when V = 0. The third term, which is independent of
H and M, has a non-dimensional axial load term, |V/Vmax|, and double expo-
nents. This term “shrinks” the yield surface (i.e. the M*, H*ellipses) for
non-zero axial load. For (almost) pure V load, especially when |V/Vmax| > 1,
numerical problems (imaginary numbers) with exponents are avoided by
using a simpler form of Equation (10.34), namely
(i) “gapping” in the active soil zone behind caisson for both support and
anchor foundations (e.g. AG, 2003; DNV, 2005)
(ii) “anchor chain trenching” – a tapered wedge-shaped trench occurring
in front of caisson due to anchor chain motions (e.g. Alderlieste et al.,
2016; Colliat et al., 2018); this cannot occur for support foundations,
as they are top loaded and have no chain
(iii) “vertical slot” above pad-eye for anchor foundations, which may
occur during installation
(iv) “reduced end-bearing” for tension loaded anchor foundations (Clukey
and Morrison, 1993; Clukey et al., 1995)
(v) “internal scoop” for support foundations at shallow L/D (Bransby
and Yun, 2009; van Dijk, 2015)
(vi) “wings” for anchor foundations (Dührkop and Grabe, 2009)
(vii) non-co-planar HM load.
210 In-place resistance
Apart from (vi), and possibly (vii), all decrease in-place resistance.
Since all geotechnical foundation models are inaccurate to a certain
degree, and the CAISSON_VHM resistance model is no exception, it is rea-
sonable to consider it to be also valid for these exceptions without a signifi-
cant decrease in model accuracy.
Hence, four reduction factors (η) are available. Figure 10.13 shows how
ηaMH, ηbMH and ηφMH modify the MH ellipse parameters aMH, bMH and φMH,
and any ηFo and ηREB combination modify Vmax. Table 10.7 gives η examples/
senses – values come from design specifications, literature or numerical
analyses. Seafloor trenching reduces Ho more than Hmax resistance, with a
smaller Vmax reduction which is mainly related to Fo (Alderlieste et al., 2016).
Because of this, and the MH ellipse shape, η order is ηaMH < ηbMH < ηφMH
< ηFo < 1. For “internal scoops” (possible only when L/D < 1.5), moment
reduction is more severe than Hmax (van Dijk, 2015). Hence η order is ηbMH
< ηaMH < 1 = ηFo = ηREB. Adding “wings” increases axial friction less than
lateral resistance; hence η order is 1 = ηREB > ηFo > ηaMH ≈ ηbMH. Note that,
since a seafloor trench increases optimum lug level (OLL) depth, whereas
wings have the opposite effect (OLL decreases), ηφMH > 1 and ηφMH < 1
respectively.
Lateral and axial resistances are modified by applying reduction factors to
auto-computed aMH, bMH and Vmax values. The shape of the yield surface (i.e.
MH ellipses and V-Hmax ellipsoids) is assumed unchanged, and the size of
the yield surface is defined by parameters Vmax and aMH, bMH.
Reduction factors ηaMH and ηbMH values < 1 may be used to decrease lat-
eral resistance due to the possibility of a “tapered wedge”-shaped soil gap
occurring in front of the caisson anchor foundation caused by anchor chain
motions. Because optimum lug level (OLL) depth also increases, an ηφMH
value > 1 is also required. To assess project-specific ηaMH, ηbMH and ηφMH
values, 3D FEA are needed (Alderlieste et al., 2016).
The same principle is also applied to “gapping”. Preliminary reduction
factor ηaMH sources possibly include COFS (2003) API/Deepstar project and
pairs of cases C1 and C3 and C2 and C4 from Andersen et al. (2005). From
these, indicative reductions in Vmax (giving ηFo and/or ηREB) and Hmax (ηaMH)
may be assessed. A simpler (very cautious) approach is to simply reduce su
by 50% over the crack depth and use η values = 1. Again, 3D FEA is neces-
sary for detailed/final design.
Lateral resistance. MH ellipses: (a) ηaMH < 1, (b) ηaMH < 1, (c) ηaMH and
ηaMH < 1, (d) ηφMH > 1, (e) ηφMH < 1
Ad (a) through (c):
• Use ηaMH and ηbMH to shrink major and/or minor axes, and hence lat-
eral resistance.
• Hmax resistance is proportional (directly related) to ηaMH.
• Ho resistance is a function of both ηaMH and ηbMH..
In-place resistance 211
Figure 10.13 In-place resistance reduction. Top two rows (a) through (e): MH
ellipses. Bottom row (f) through (h): V-H max ellipsoids.
Ellipse rotation (i.e. ηφMH /= 1) alters both Hmax and Ho resistances; ellipse
shape is unchanged.
Axial resistance. V-Hmax ellipsoids (f) ηFo < 1, (g) ηaMH < 1, (h) ηFo and ηaMH
<1
Use either ηFo and/or ηREB to reduce axial resistance Vmax
212 In-place resistance
Notes:
ηaMH MH ellipse resistance semi-major axis length reduction factor
ηbMH MH ellipse resistance semi-minor axis length reduction factor
ηφMH MH ellipse resistance rotation angle increase factor
ηFo caisson outer skin friction reduction factor
ηREB caisson tip reverse end-bearing resistance reduction factor
REB caisson tip reverse end-bearing resistance
HM non co-planar HM load
+ overturning M load
− restoring M load
10.12.8 Resistance comparisons
CAISSON_VHM was verified internally (Kay and Palix, 2010–2015). The
following paragraphs give details of external (third-party) comparisons
made with other programs.
10.12.8.1 Support foundations
Organisation A requested CAISSON_VHM analyses for three projects,
comprising two support foundations and one anchor foundation, during
their free trial period.
Organisation A made spot checks at L/D ≈ 0.5, 3 and 1, see L and zlug
italicised data in Table 10.8, and showed reasonable agreement. In addition,
less than 1 s computer time was needed for all CAISSON_VHM analyses (3
runs, 1 or 2 load cases, 24 caisson diameters).
ORGANISATION A C AISSON_VHM
Parameter Organisation CASSION Remarks
A Analysis VHM
Analysis
Project 1 – Spreadsheet Support foundation, free head.
SUBSEA Single soil layer. Two VHM
TEMPLATE load cases. Optimum
Diameter D [m] 5.5 5.5 caisson embedded lengths
Embedded length 3 2.4 L opt found for seven caisson
L [m] diameters (2.5 m to 5.5 m).
Project 2 – Internal Support foundation, free head.
SUBSEA Software Two soil layers. Factored
STRUCTURE loads VHM(T) = 2.27, 2.0,
(MONOPILE) 20.6 and 15.9 (MN,m).
Diameter D [m] 8.5 8.5 Reduced clay-steel outer
Embedded length 27.5 27.5 friction adhesion factor
L [m] α o,twist used to account for T
load. Optimum caisson
embedded lengths L opt found
for ten caisson diameters
(4.0 m to 8.5 m).
Project 3 – PLAXIS 3D Anchor foundation, fixed
SUCTION head. Two soil layers.
ANCHOR Factored chain load 17 MN
(MOORING applied at seafloor with an
SYSTEM) (1) angle of 24°. Optimum
Diameter D [m] 12 12 caisson geometry
Embedded length 15 14 (embedded lengths L opt,fix
L [m] and optimum lug level
Lug level depth, 27.5 27.5 depths z lug,opt,fix ) found for
z lug [m] seven caisson diameters
(9 m to 15 m).
Note: PLAXIS 3D analysis was not optimised
Figure 10.14 Comparison anchor pile Tlug – BIFURC, CAISSON_VHM and PLAXIS 3D.
In-place resistance 215
extremely low embedment ratio L/D implies that cross plates may be neces-
sary in order to prevent internal scoop failure. These calculations assume
that the caisson is not heavily loaded axially – Figure 10.12 shows that
V-Hmax ellipsoids are almost vertical until around V/Vmax = 0.5.
Figure 10.15 shows the example embedded length L calculations. 5 m
diameter caisson in NC clay (su(kPa) = 2 + 1.5z) with factored lateral load
H = 500 kN.
• there are two equations to be solved – one for the resistance envelope and
another for the chain load at lug level – as summarised in Figure 10.17
• expressions for maximum caisson axial and lateral resistance, plus 3
MH ellipse constants (aMH, bMH and ΦMH) and two VH ellipsoid con-
stants (aVH and bVH) all have to be programmed in terms of caisson
geometry (and shear strength profile).
Figure 10.16 Vmax-Hmax ellipsoid example, optimum lug level (after Kay and Palix, 2010).
In-place resistance 217
(a)
Figure 10.18 C
aisson support foundation design example using CAISSON_VHM
and CANCAP2: (a) soil s u profile, NC clay, s u ≈ 2z; (b) MH plot-load
cases and envelopes for given caisson geometries (V = 3000 kN);
(c) required optimum caisson embedded lengths, L opt , and embed-
ment ratios, L opt /D, for given caisson diameters D (factored VHM
loads = 2795, 1079, 7267 [kN, m]) (Kay and Palix, 2011).
to give an insufficient FOS against base failure when L/D > 6). Based on the
resulting Lopt-D curve, and the requirement to limit 30 year settlements to
less than 0.25 m, the final caisson dimensions were D = 5.0 m, embedded
length L = 11.4 m, total length = 12.4 m, side WT = 25.4 mm. The caisson
weight was 73 tonnes.
PLAXIS 3D (PLAXIS, 2013) was used to verify the final caisson geometry
for two of the most critical VHM load cases. The caisson was modelled as
totally rigid with elastic material (steel). The clay-steel outer adhesion factor,
αo, was nominally 0.65, but with an appropriate reduction for torsion
moment. The failure points (solid blue diamonds) obtained by PLAXIS 3D
for the two load cases were added to the CANCAP2 MH resistance enve-
lope and load case MH data (Figure 10.19). It can be seen that the resistance
from the PLAXIS 3D FEA model is very slightly higher than the CANCAP2
resistance envelope. The PLAXIS 3D analysis needed 24 hours to run (2010)
for each load case, whereas all the CAISSON_VHM resistance envelope
analyses took less than 10 seconds.
In-place resistance 219
(b)
(c)
Figure 10.19 C
aisson support foundation design example using CAISSON_VHM:
PLAXIS 3D verification. Caisson D = 5 m, L = 11.4 m.
caisson: D = 3 m, 3.5 m and 4 m. Wall and top plate thicknesses are D/200
and D/100
lug: on caisson outer face, zlug at optimum depth
soil: su = 2 + z (0 ≤ z ≤ 5 m), su= 7 + 1.73 (z - 5) (5 m ≤ z ≤ 20 m) kPa, γsub
= 4 kN/m3
soil-pile: αo = 0.6 (including twist allowance), αi = 0.65.
pull-out Vmax: lesser of plugged and unplugged conditions, ηREB = 0.8
chain: 81 mm diameter, submerged weight 1.12 kN/m.
loads and FOS: see Table 10.9.
Table 10.10 compares detailed results for both tilt angles. Figure 10.20a
summarises the caisson geometry optimisation mode, which finds the
optimum “fixed head” zlug,opt,fix and Lopt,fix for given caisson diameter(s) D,
In-place resistance 221
T z0 θ z0 λT
Load Case [kN] [deg] [–]
1: All Lines Intact 2045 7.5 2.00
2: One Line Damaged 3377 13.7 1.50
3: Transient 4211 15.5 1.25
seafloor chain loads Tz0 and inclinations θz0 (i.e. Tlug and θlug at caisson lug
level zlug). A 5° tilt angle was used. Figure 10.20b shows that LC3 always
has the largest embedment length L and is therefore critical. The associ-
ated maximum embedment ratio L/D is just over 6, suggesting that base
heave installation may be an issue had D been less than 3 m. Figure 10.20c
and Figure 10.20d provide optimum buoyant weight Bw and geotechnical
efficiency (defined as η = Tlug/Bw) data used to finalise foundation geome-
try. LC3 efficiency η values decrease from ≈ 22 (D = 3 m) to 17 (D = 4 m),
showing that “deeper is better”. Final design was based on 3.5 m OD
caissons (5° tilt, LC3), see Table 10.12. Additional checks subsequently
made included capacity, lug level, tilt and twist. The sensitivity mode used
(vary θlug and/or zlug, find Tlug) is shown on Figure 10.21a. Figure 10.21b
and 10.21c show Tlug sensitivity to θlug and zlug respectively.
Less than 5 s computer time was needed for the two runs for Figures 10.20
and 10.21.
Alderlieste et al. (2016) give details of a more complicated anchor chain
trenching example. This used PLAXIS 3D to first assess lateral resistance
reduction factors, followed by CAISSON_VHM analyses for in-place resis-
tance. Details are given in Section 12.12.
222
In-place resistance
Figure 10.20 F PSO anchor design example using CAISSON_VHM, geometry optimisation mode (zero rotation). 5° tilt: (a) optimisa-
tion solution – (a) find L opt,fix and z lug,opt,fix for given D, T lug and θ lug (b) L opt,fix , and L/D versus D, (c) B w and L opt,fix versus D,
(d) efficiency (T lug /B w ) and L/D – D (Kay, 2015).
In-place resistance 223
(a)
(b) (c)
Figure 10.21 F PSO anchor design example, T lug sensitivity mode (non-zero rota-
tion). D = 3.5 m, LC3, λ T = 1.25; (a) optimisation solution – find T lug
for given L, D, z lug and θ lug , (b) T lug λ T and T lug versus θ lug at z lug,opt,fix ,
(c) T lug versus z lug at θ lug (Kay, 2015).
embedment. Govoni et al. (2016) provide h0, m0 and a values for strips in
30° sand at L/D = 0 and 0.5.
10.14.1 General
The axial resistance and initial stiffness of weak rock will depend greatly on
the mode of installation. Driving in weak rock will lead to the creation of a
remoulded or crushed zone around the wall of the structure and potentially
will lead to radial cracks around the structure (Section 9.13). The axial and
lateral resistance of piles in weak rock are discussed in the following section
for both installation cases.
Qs .UCS. A s fl . As (10.37)
Where:
α = adhesion factor
UCS = Unconfined Compressive Strength.
0.11. UCS
0.5
(10.38)
Figure 10.22 A
dhesion factor (α) against rock material UCS (Terente et al.,
2017).
conducted in the UK (St Nicolas site), comprising axial and lateral load tests
(Jardine et al., 2019).
Hole stability is one of the first factors to consider for a drilled and grouted
foundation. It is noted that when sub-layers of granular material are encoun-
tered in the rock matrix, the drilling fluid circulation will most probably
wash out these materials. Wash out will produce large voids that could lead
to stress relief in the rock mass, and also require a significant amount of
additional grout to be injected to fill these voids. If pressurised, the injection
of grout can also lead to amplification or creation of fractures.
In-place resistance 227
• through a shearing bond between the grout and the steel pile
• sliding friction between the grout shaft and the rock (when the drilled
hole is smooth and grout infiltration limited)
• dilation of an unbounded rock-grout interface until asperities shear off.
For stubby piles submitted to a large number of cycles of loading, the axial
resistance should not rely on the passive shear capacity between the steel
and cast grout (c.f. DNVGL RP 0419 recommendations for grouted connec-
tions). Pile movement during grout curing could jeopardise the steel/grout
bonding. In addition, cyclic lateral loading will tend to create a gap between
the grout and the steel. The shear capacity of the interface can be increased
Figure 10.23 P
luto – drilling head (3m diameter) with underreamer. (Courtesy of
Large Diameter Drilling (LDD))
228 In-place resistance
by using shear keys, i.e. protuberances on the steel surface, on the interface
that when subsequently cast into the grout provides a mechanical resistance
to relative sliding between the steel and grout. API RP 2A WSD (2000) and
ISO 19902 (2011) provide guidance on the number of shear keys to be
considered.
In non-carbonate weak to moderately strong rock, high friction values
have been recorded. Most published methodologies define the unit shaft
friction of drilled and grouted piles as a function of the Unconfined
Compressive Strength (UCS), independent of the type of friction developed
(i.e. sliding or dilation). Table 10.11 lists the main methodologies proposed
for shaft friction, fs, where
fs . UCS (10.39)
m
where:
UCS: Unconfined Compressive Strength
m = constant usually taken as 0.5
α: adhesion factor comprised between 0.2 and 0.45.
Table 10.11 U
ltimate Friction Formulation for Drilled and Grouted Piles in
Weak to Moderately Strong Rock.
n r (10.40)
SRC c
1 ds
where:
∆r is the mean roughness height
ds is the socket diameter
ηc is the construction method reduction factor
n is the ratio Em/UCS of rock mass modulus over unconfined compressive
strength.
The effect of the SRC on the adhesion factor is provided in Figure 10.25.
For slightly cemented carbonate rock (UCS < 5MPa), ARGEMA (1994) sug-
gest considering the ultimate shaft friction at the carbonate rock–grout
interface proposed by Abbs and Needham (1985; see Figure 10.26). For
chalk, CIRIA (2002) differentiate the axial resistance depending on the den-
sity. For low to medium density, the friction is a function of the average
effective stress, fs = 0.8 σv'. For high-density chalk, the axial resistance is
proportional to the UCS: fs = 0.1 UCS.
230 In-place resistance
Figure 10.25 E ffect of shaft resistance coefficient (SRC) on socket adhesion fac-
tor (Seidel and Collingwood, 2001).
Figure 10.26 U
ltimate shaft friction at the carbonate rock–grout interface (from
Abbs and Needham, 1985).
In-place resistance 231
evolution of normal stresses applied on the interface during rock socket pile
displacement (Figure 10.27 and CFMS, 2019).
The imposed stiffness should be representative of the actual rigidity of
the rock mass. The rock–grout interface is represented by a regular saw-
tooth with a chord length and an asperity height corresponding to the
roughness of the in situ hole. Guidance for characterising rock masses for
designing drilled and grouted offshore pile foundations can be found in
Puech and Quiterio-Mendoza (2019), and testing conditions are described
by Stavropolou et al. (2019).
In the case of weak rock (i.e. grout being much stronger than the rock),
the interface response will be governed by the rock behaviour. The failure
surface will pass through the weak rock. For stronger rocks (i.e. strength
similar to the grout), it is more the “geometrical” dilatancy that will govern
the interface response.
The main challenge of this experimental methodology (based on CNS
tests) is to get access to proper design parameters and installation conditions
at an early stage in the design process. The measured shear force is very
sensitive to the shape of the asperities. It remains, however, the best way to
understand the degradation of asperities under cyclic loading without
requiring full-scale tests.
(a)
(b)
Figure 10.28 C
haracteristic shape of p-y curves for weak rock according to (a)
Fragio et al. (1985) and (b) Reese (1997) methods. (With permis-
sion from ASCE)
234 In-place resistance
Figure 10.29 C
oncept of “chipping” and characteristic shape of p-y curves (after
Erbrich, 2004).
10.15 GROUP RESISTANCE
10.15.1 General
Groups of suction foundations may be installed. Examples are given in Figure
10.30. They include two or more suction anchors for FPSO mooring points,
and braced foundation clusters to ensure a near-vertical subsea template.
10.15.2 Considerations
Considerations include:
Like pile groups, the primary consideration is the spacing between individ-
ual foundations. At close spacings, when foundations interact, bracing (if
any) imposes kinematic restraints – an unbraced self-installing platform
behaves differently than a braced tripod.
Group load (capacity and response) is normally considered when the centre-
to-centre spacing is < 8 D. Generally, for foundations in undrained “clay”,
the group capacity may be less than the single suction foundation capacity
multiplied by the number of foundations in the group. Axial and lateral
group deflections are normally larger than that of a single suction founda-
tion subjected to the average load of the group.
10.15.3 Design procedures
There are no established procedures for designing groups of shallow or
intermediate foundation systems.
Therefore, engineers resort to ad hoc (case-by-case) approaches.
10.15.4 Design tools
Design tools (listed in order simplest first) include:
(a)
(b)
(c)
(d)
Figure 10.30 C
ontinued: Examples of suction foundation groups: (c) cluster pile
for Ceiba project (Courtesy SPT Offshore), (d) L6-B platform dur-
ing sail-out (from Alderlieste and Van Blaaderen, 2015).
238 In-place resistance
Figure 10.31 Laterally loaded anchor. Braced pile groups (Cox et al., 1984).
In-place resistance 239
10.16 CLOSURE
This chapter has covered in-place resistance, both for support and anchor
pile type intermediate foundations, and the following sections subdivide
comments into those for “clay” and “sand”.
10.16.1 Clay
As for installation, models for undrained “clay” are more developed than
for drained “sand”. Reasonably accurate models have been given to assess
maximum available Vmax, Hmax and Ho (“fixed” and “free head”) and Tmax
resistances. For anchor foundations, pull-out Vmax has three possible failure
models – “plugged”, “coring” and “leaking” (see Figure 10.1), for which
both applicability criteria and equations are provided.
VHM(T) interaction effects are important: HM load affects V resistance,
and, to a lesser extent, T load affects V resistance. Like shallow foundations,
VHM(T) yield envelopes for undrained soil are a promising development
for preliminary/routine design. Equations and software are available. PLET
and FPSO anchor design examples are also given.
Both FEA and centrifuge tests suggest that the original (Skempton, 1951)
bearing capacity factor Nc,circle is cautious. Hence Qtip is probably underesti-
mated. This has implications for both reverse end-bearing (in-place, this
chapter) and plug heave (installation, Chapter 9) failure.
For anchor piles, Vmax-Hmax and Vlug-Hlug, envelope shapes differ; compare
Figure 10.16 with Figure 10.5. This should be considered during detailed
design when making sensitivity analyses.
10.16.2 Sand
Models for “sand” are less mature than for “clay”. This is as expected: effec-
tive stress governs strength, and pore water build-up/dissipation is an addi-
tional complicating factor.
For sustained load, sands behave drained, and calculations can be easily
made to assess both the compressive resistance Vmax (which is typically very
large) and the tensile pull-out resistance (which is very small). There is also
a good understanding of near-surface foundations to combined VHM loads,
largely based on experimental data. Yield surfaces defining possible VHM
load combinations have been published, but these are for either shallow
foundations or jack-up spudcans. Hence, as a start, a nice research project
would be to derive a reasonably cautious model for assessing intermediate
foundation resistance in drained sand under static VHM loading. Besides
foundation geometry, which must take account of embedment, key soil
parameters include friction angle φ' and dilatancy parameter ψ'. This
research would hopefully lead to a drained resistance envelope and yield
240 In-place resistance
function similar to, but far more complex than, that developed for interme-
diate foundations in undrained soil.
For short-term loading (e.g. a single wave load or boat impact), resistance
estimates are challenging because the equivalent undrained sand strength is
very difficult to assess. Undrained strengths are principally determined by
dilatancy parameter D (DNV, 1992) and the onset of pore fluid cavitation.
Hence both vertical compressive and tensile resistances are typically very
large, but difficult to assess. The resistance to combined moment and hori-
zontal loads is also considered to be large. Like clays, VHM interaction
effects must be considered. Resistance assessments become even more com-
plex for cyclic load. An even more challenging research project would be to
initially assess resistance in undrained sand under static VHM loading.
Chapter 11
In-place response
11.1 INTRODUCTION
11.2.1 General
Calculation of foundation displacements (usually axial settlements) should
consider the following components:
• immediate displacements
• primary consolidation settlement
• secondary compression (creep) settlement
• settlements induced by cyclic loads
• regional settlements.
The first three components are the same as for shallow foundations (ISO
19901-4:2016, section 7.5.2.1). Differential settlements (induced by spatial
soil variability, moments, torque and eccentricity) are not assessed for indi-
vidual intermediate support foundations, as they are essentially rigid.
241
242 In-place response
(a) Sands (and stiff clays and rock): settlement magnitudes are small and not
problematic. Hence, settlements are not considered for routine design.
(b) Clays (fine-grained lithogenous sediments) are the most abundant of
all marine sediments (about 70% by volume). These clays are usually
normally consolidated (or at best lightly overconsolidated) and give
the largest settlements.
(c) Organic soils (peats) are rarely encountered offshore.
Table 11.1 summarises the relative importance of each of these three settle-
ment components for different soil types.
Settlement considerations include:
Table 11.1 R
elative Importance of Settlement Components for Different Soil Types.
Intermediate Support Foundations – Routine Design
Soil Type
11.2.2 Immediate displacement
References for methods to assess immediate, elastic displacements that
account for non-uniform soil profiles (e.g. linearly increasing soil strength)
and foundation embedment are provided in this section. The results are
generally used to provide a range of foundation stiffnesses for use by
structural engineers or to check the initial portion of “load-settlement”
response from geotechnical finite difference/element analyses.
The complete 4 × 4 seafloor foundation stiffness matrix (Section 8.3 and
Figure 8.6c) should be assessed. Note that the lateral stiffness sub-matrix
should not be derived using either p-y curves for deep (pile) foundations or
elastic solutions for flexible piles. Reasons include the fact that intermediate
foundations are essentially rigid and generally larger in diameter. Instead,
either elastic solutions or numerical analysis should be used (DNV OS-J101,
2014a).
In order to evaluate the terms appearing in the stiffness matrix, sources
for useful elastic solutions (all continuum based) for statically loaded
circular intermediate foundations are listed in Table 11.2). It is seen that
only Doherty and Deeks (2003) provide all six stiffnesses – and in non-
homogeneous soil – but only for embedment ratios L/D ≤ 2. Of these
solutions, only Baguelin et al. (1977) include a weaker zone. Their equation
12 can be used to demonstrate that a disturbed soil annulus has only a
limited effect on lateral stiffness kxx. For example, a disk with diameter D =
5 m, outside boundary radius R = 15D, soil Young’s Modulus E' = 100 MPa,
and Poisson’s ratio ν' = 0.3 has an intact stiffness kxx = 121 MN/m/m length.
For a 0.5 m wide disturbed zone with E' = 40 MPa, kxx reduces to 113
MN/m/m length, i.e. a small 7% decrease. This has implications for laterally
loaded OWT monopiles installed by vibratory methods, where both the
annulus width and soil modulus are likely to be greater than for driven
foundations. That is, since lateral stiffness is largely due to the far field
displacement regime, and most lateral displacements at depth remain in the
elastic zone, driven and vibrated monopile lateral responses under HM load
will not differ significantly.
Until Gupta and Basu (2016) and Arany et al. (2017), there were no elastic
continuum solutions for laterally loaded OWT type rigid caissons
(embedment ratio L/D ≤ 10) in non-homogeneous soil. Carter and Kulhawy
(1992) were valid for L/D ≤ 10 but are only for uniform soil and Doherty
et al. (2005) included non-homogeneous soil but for L/D ≤ 2.
There are no elastic solutions for deeply embedded rigid circular founda-
tions under MH loading – this situation occurs at the base of a 1-D Winkler
support model for OWT monopiles: Poulos and Davis (1991) give solutions
for V load only and for MH load on an embedded vertical plate (their sec-
tion 15.4), and Doherty and Deeks’ (2003) caisson geometry (b) is valid
only for L/D = 2.
244 In-place response
(a) (b)
(c) (d)
(e)
Figure 11.1 P
rimary Consolidation Settlement (support foundation): (a) simple
two-spring model, (b) through (e) corresponding (SLS, permanent)
loads and resistances; Soil Plug submerged weight not shown for
clarity.
In-place response 247
mobilised inner skin friction Fi, which acts upwards on the steel and down-
wards on the soil plug, has a value equal to Qbase. Considering loads on the
clay soil plug (Figure 11.1d), the top spring is compressed upwards from tip
level by a vertical effective stress increase ∆σv' = Qbase/Abase. Due to inner
friction resistance Fi, ∆σv' decreases with distance above tip level, and the
corresponding uplug (soil plug spring compression) final primary consolida-
tion settlement value can be assessed using textbook equations (e.g. Holtz,
1991) for 1-D oedometer compression and soil parameters Cc, Cs (compres-
sion indices) and p'c (preconsolidation pressure).
Additional settlement will occur due to increased vertical soil stresses in
the vicinity of the foundation. Therefore, settlement assessment also requires
calculation of the stress distribution in the soil below the foundation. Hence,
utip (soil settlement below tip level) can be computed using loads at tip level
and conventional load spread theory or similar. Distributions of ∆σv' (which
may be assumed to be uniform) due to two loads which need to be superim-
posed. Figure 11.1e shows that they are (i) base load Qbase acting over a
circle diameter D and (ii) mobilised outer skin friction Fo spread out with a
load spread factor ns to act over an annulus with inner diameter D and outer
diameter D + 2 L/ns.
For a fixed geometry, consolidation settlement magnitudes are sensitive
to Fo, Cc, p'c and Cs (in that order). If uncertainties exist, then analyses
should consider a high estimate of uplug. This is achieved by using a low
estimate of Fo. The model presented here is useful for foundation geome-
tries where plug final settlement may be significant, generally embedment
ratios L/D less than around 3. Longer intermediate foundations behave
more like pipe piles – tip settlement component utip dominates but is small
at working loads.
For axially loaded support foundations in clay, settlement rates are a
product of 1D vertical drainage within soil plug plus 2D drainage (com-
bined vertical + radial outwards) below foundation tip level. The major
unknown is the mass coefficient of consolidation value.
11.2.6 Regional settlement
Causes of regional settlement include (most common first):
Do not use either pile cyclic p-y data (ISO 19901-4:2016a Section 8.5, soil
reaction for piles under lateral loads) or elastic solutions for flexible piles for
intermediate foundation design.
For OWT intermediate foundations (aka monopiles) the main geotechni-
cal design challenge is that tilt should not exceed the allowable threshold.
Typical SLS criterion is that accumulated tilt at seafloor < 0.5°. A secondary
challenge is eigen frequency assessment of complete OWT system (i.e. foun-
dation, substructure, mast and rotor nacelle).
Monopile design requires a good knowledge of the soil structure interaction
at small displacements. Offshore wind turbines are commonly designed so
that their first natural frequency avoids frequencies of operational and envi-
ronmental loads, see Figure 11.2 (CFMS, 2019). Design also needs to ensure
that the permanent inclination of the foundation does not exceed a fixed limit
at the turbine end of life (the inclination typically should remain < 0.5°).
From back analysis of monitored data, Kallehave et al. (2012) found cur-
rent design models fail to reproduce the natural frequency although it is a
key parameter for the design. Most soil-pile stiffness models derived for
flexible piles tend to underestimate the natural frequency.
The method normally used to design piles under lateral and overturning
loads is based on a Winkler modelling approach, commonly termed the p-y
method. P-y curves proposed in standards API RP2A (API, 2000) or
DNV-GL RP C212, (DNV-GL, 2017, 2019) are generally regarded as being
satisfactory for jacket pile design but fail to provide reliable response for
large diameter piles at small displacements. Several semi-empirical
approaches have been proposed to stiffen the initial lateral response for
“large diameter” piles, notably:
• Kallehave et al. (2012), Sorensen et al. (2010) and Kirsch et al. (2014)
for sand
• Stevens and Audibert (1979) and Kirsch et al. (2014) for clay.
The initial response could also be calibrated on the basis of in situ measure-
ments such as High Pressure Dilatometer (HPD) tests, or using finite
In-place response 249
Figure 11.3 P
ISA model: (a) soil reaction components incorporated in the PISA
design model; (b) 1D FE model employed in the PISA analysis.
Figure 11.4 Illustration of chalk-pile load deflection concept (Muir Wood et al.,
2015).
et al., 2020), highlighted that the so-called p-y method was originally devel-
oped for relatively slender piles and should not be sufficient to model the
soil structure interaction of rigid piles with L/D < 6. However, in the case of
weak rock, field experiences tend to show that properly calibrated p-y curves
can predict reasonable pile rotation.
In the case of weak carbonate rock, pile installation (see Section 9.13) can
lead to the creation of a crushed or remoulded zone around the pile. In that
case, the subsequent lateral response will be softer than if the pile was sur-
rounded by intact rock. To account for the presence of this remoulded zone,
Muir Wood et al. (2015) proposed to consider two springs in series
(Figure 11.4). The overall deflection of the pile (y) under a given load is
assumed to be the sum of the deflection of the remoulded chalk (yr) and the
deflection of the intact chalk (ys).
A similar concept was applied for driven and drilled and grouted piles in
calcarenite (Lovera, 2019) where the grout annulus, crushed zone and intact
rock were modelled as three springs in series (Figure 11.5).
The equivalent initial local stiffness can be deduced using
p Eoed
Eequiv (11.1)
y a grout Eoed crushed Eoed
Dsteel Egrout Dcrushed Grock
where:
For the grout,
252 In-place response
Figure 11.5 M
odel for pile response in weak rock with grout anulus, crushed
zone and intact rock modelled as three springs in series (Lovera
2019).
a parameter equals zero for a driven pile and 1 for drilled and grouted piles;
λ parameter is smaller than one, taking into account that the pressure does
not apply on the all circumference.
Dcrushed = Dsteel for driven pile and Dgrout drilled hole diameter for drilled and
grouted piles;
δcrushed is the crushed zone thickness;
Eoed is the oedometer modulus of the crushed zone material at low stresses.
Figure 11.6 S ketch of the rheological model considered for piles driven in soft
rock (Lovera et al., 2020).
11.5 CLOSURE
This chapter has discussed in-place response and has shown that different
intermediate foundation types have their own challenges. Anchor founda-
tion movements are usually non-critical and are rarely assessed.
Support foundations may be conveniently subdivided into those for oil
and gas and renewable projects. The challenges are listed below in frequency
order (most commonly encountered first).
Support platform projects require 4 × 4 foundation stiffness matrices to
pass to the structural engineer. Simple (elastic) solutions have been listed to
evaluate the stiffness terms kzz (axial), kxx, kxθ, kθx, kθθ (lateral, coupled) and
kΨΨ (torsional). ISO/API type “p-y” curves are unsuitable for intermediate
foundation lateral response assessments due to the rigid body rotational
behaviour of intermediate foundations.
254 In-place response
12.1 INTRODUCTION
255
256 Miscellaneous design considerations
12.2 SCOUR PROTECTION
12.3.1 Introduction
Intermediate foundations may interact with the following adjacent infra-
structure (the most important are listed first):
• jack-up spudcans
• adjacent intermediate foundations
• adjacent shallow foundations and pipelines
• conductor wells.
12.3.2 Jack-up spudcans
Minimal-facilities hydrocarbons platforms normally require subsequent loca-
tion of a jack-up rig to install well conductors. Similarly, offshore wind tur-
bine monopiles or tripod substructures may be installed using a jack-up rig
fitted with spudcans. Figure 12.2 shows examples. Hence there may be a risk
of overstressing platform foundations. Such situations require consideration
of interaction between the jack-up spudcan(s) and platform foundations, par-
ticularly in weak soil profiles. Geotechnical analysis involves assessment of
soil movements induced by spudcans and the effects of these movements on
(a)
(b)
the foundation, both lateral and axial (Mirza et al., 1988; Siciliano et al.,
1990). Normal practice for pipe pile foundations in normally consolidated
clay profiles is to ignore interaction effects provided that the closest soil gap
distance is at least one spudcan diameter (ISO 19905-1:2016b).
However, for intermediate foundations, this gap distance (1 spudcan
diameter D) should be critically reviewed. This is because the foundations
are essentially rigid: any load due to lateral soil displacement is transferred
directly into the platform leg, and not partially absorbed by pile bending.
A key issue is predicting lateral soil movements induced by the spudcan at
the intermediate foundation. These movements have both horizontal and
rotational components. Due to smaller foundation embedment depths (than
piles), rotational effects are usually more critical than lateral effects –
foundation rotations could induce additional BM/overstress the platform
leg. This situation is similar to kicking one leg of a three-legged stool.
Jack-up spudcan–foundation interaction analysis is difficult. This is
because the accuracy of any prediction is a function of the accuracy of both
the data and the model (Lambe, 1973), and larger-than-usual inaccuracies
generally exist both for data (e.g. soil shear strength profile, spudcan geom-
etry and preload, foundation geometry and gap distance) and the model (i.e.
spudcan penetration tip depth and resulting lateral soil movements around
the spudcan).
Spudcan installation induces vertical foundation movements. However,
industry perception (for pile foundations) is that axial (downdrag/upthrust) is
not so critical as lateral. The majority of axial pile resistance is derived from
soil well below the spudcan and is rarely (possibly never) studied in detail.
Hence, instead of complex 3D large strain Finite Element Analyses, a rea-
sonable design approach may be to:
Note that care has to be taken with orientating foundation p-y data. This is
because the offsets need to be applied along a line connecting the spudcan
and the jacket leg(s), whereas the platform structural model automatically
rotates these to coincide with the storm/wind direction. However, since
spudcan penetration is unlikely to occur during a “small” (10-year summer
storm) event, then orientation may not be critical.
Further information on this approach is given by Xie (2009) and Xie et al.
(2017) on lateral soil movement data, by Tho et al. (2013 and 2015) on
Miscellaneous design considerations 259
Figure 12.3 O
phir well head platform, soil properties (Alderlieste and Dekker,
2018).
Figure 12.4 T
op view (left) and side view (right) of work-over jack-up rig at
Ophir well head platform. Suction pile diameter D = 6 m, target
penetration L = 13.7 m. Spudcan diameter D = 17 m, spudcan design
penetration h = 12.75 m, minimum spudcan-suction pile clear soil
gap S = 19.4 m (Alderlieste and Dekker, 2018).
Figure 12.5 S uction pile displacements due to spudcan penetration (not to scale)
at Ophir well head platform. Left: first spudcan installation. Right:
due to second spudcan installation (Alderlieste and Dekker, 2018).
Miscellaneous design considerations 261
12.3.5 Conductor wells
Conductor installation techniques should consider the possibility of soil
removal (wash-out) occurring in sand and appropriate measures (e.g. using
driven isolation casings) taken. This is to prevent loss of resistance at nearby
intermediate foundations.
Intermediate foundations also interact with the conductor wells. Generally,
the only potential problem is that, when laterally loaded, intermediate foun-
dations may induce additional displacements in the conductors. However,
due to their flexibility, they will probably accommodate the extra lateral soil
movement without significant bending moment increases.
12.4 SHALLOW GAS
• Site investigation: gassy soils are difficult to sample and test in the
laboratory. Therefore, site investigations should include in-situ tests.
262 Miscellaneous design considerations
• Installation: since such soils are “softer”, and the vertical stress reduc-
tion may cause gas to come out of solution, more heave than usual
may be anticipated. However, this may be offset by the fact that instal-
lation underpressures (i.e. tensile forces on the soil plug) are likely to
be low in normally consolidated soft clays.
• In-place response: in time, gas migration could create a gas “bell” under
the foundation top plate, causing an upwards force. This potential
problem is removed by providing open vents in the top plate. Assuming
vents are present, foundation compressive resistance can still accom-
modate bearing pressures under the foundation top plate. Maximum
axial tensile pull-out resistance, Vmax, cannot rely upon “passive suc-
tion” (since the top plate vent is permanently open). Instead, pull-out
should be based on outer friction plus the minimum of (inner friction
or soil plug weight) plus foundation submerged weight only, i.e. the
“coring” and “leaking” failure models shown in Figure 10.1.
pump system can provide around 80 litres/s. Both ∆V and Q are considered
reasonable. Better (and hopefully lower) Q estimates may be possible by
careful examination of suction pump instrumentation data (penetration rate
and pump outflow rate) during installation, especially at/near final penetra-
tion, where steady state conditions are most likely to have been achieved,
and the water plug volume decrease is small: subtracting the water plug
volume decrease from the pump outflow gives Q, the amount of water enter-
ing the soil plug.
Laboratory model tests on “vacuum anchors” have been reported by sev-
eral researchers during the 1980s, including Goodman et al. (1961), Brown
and Nacci (1971), Wang et al. (1975, 1977, 1978) and Helfrich et al. (1976).
Wang et al. (1977) give an equation for (vertical) breakout capacity.
Permanent underpressure is currently not a relied-on solution.
Vent(s) in the top plate are required to expel the “water plug” (the volume
of water trapped between the top plate and seafloor) during the self-weight
penetration phase. The vent outlet pipes should be of sufficient diameter to
prevent high overpressures occurring in the “water plug” for the specified
range of penetration velocities. If water pressures are high, they may cause
piping along the foundation perimeter (sand profile) or base failure (clay soil
profile). If the foundation design relies on in-place “passive suction”, or
extraction using overpressure is envisaged, then vent design should incorpo-
rate a valve. The valve has to be open during installation, and subsequently
closed. This is to ensure a permanent seal.
For intermediate foundations installed by crane, an additional vent design
consideration may be while the foundation is being lowered through the
wave splash zone: vents should also be sufficiently large to allow essentially
free air egress/ingress and hence avoid “snatching” (large variations in lift-
ing wire loads).
The basic equation for incompressible fluid flow through a single vent
(nozzle) is given by the equation:
v n Cd 2 g H (12.1)
where
Hence, for a single vent, the overpressure increase, ∆u, and corresponding
vent (nozzle) resistance Rn, are given by:
R n u Ai – An kN (12.3)
where
∆H = vent (nozzle) differential head loss [m]
= [1/(2 g)] (vn2/Cd2)
γwater = fluid (sea water) density (= 10.05) [kN/m3]
Ai = pile internal area (π Di2/4) [m2]
An = vent (nozzle) area [m2]
g = acceleration due to gravity [m/s2]
vn = vent (nozzle) velocity [m/s]
= vpile Ai/An
vpile = foundation velocity [m/s]
Cd = vent (nozzle) discharge coefficient [–]
Equation (12.2) can also be used to assess the vent (nozzle) area An
required to prevent the overpressure ∆u causing base end-bearing (“push-
out”) failure. In addition to the foundation velocity Vpile, top plate vent area
An results are sensitive to the coefficient of discharge, Cd. The minimum Cd
value should be taken to maximise vent (nozzle) resistance Rn. If the vent is
sharp-edged and circular, the “ideal” Cd value is 0.611 (= π/(π + 2)). Hence a
cautious value is usually taken to be 0.61. Since entry into the seafloor is
usually the most critical stage, then it is usually sufficient to analyse this situ-
ation only. For controlled lowering, a reasonable value for foundation
touch-down velocity, vpile, is 0.2 m/s (see Section 9.9).
Based on these results, the vent design is satisfactory (∆u < ∆uall).
As stated earlier, the seafloor situation is assumed to be the most oner-
ous. More complicated calculations are necessary should velocity increase
with depth below seafloor, for example due to free-fall caused by uncon-
trolled lowering (Section 9.9.3, free-fall penetration, and Figure 9.8, free
fall velocity – depth curve). In such cases, the Rn, the corresponding vent
(nozzle) resistance, has to be computed, and interested engineers may care
to check that Rn is 143 kN when vpile = 0.2 m/s.
During installation, the suction pump maximum achievable flow rate should
be sufficient to discharge both the “water plug” volume plus the anticipated
outflow rate at seafloor. This is especially important in sands and gravels of
high permeability and large diameter foundations (e.g. > 10 m).
Furthermore, the suction pump used during installation should be capa-
ble of generating adequate underpressure/overpressure.
To calculate the “water plug” flow rate component, a penetration velocity
of 1 mm/s is generally used. This is a factor 20 slower than a CPT. Houlsby
and Byrne (2005) give an equation for assessing seafloor outflow rate due to
underpressure. A high estimate of soil permeability should be used.
12.8 FOUNDATION INSTRUMENTATION
• penetration
• ambient water pressures acting externally on the top plate and inter-
nally (either under the top plate or at the suction outlet pipe, rather
than at the pump, in order to avoid errors due to pipe friction and
Venturi effects)
• soil plug heave
• inclination and north-seeking gyroscope (intermediate anchor
foundations)
• position (x, y, z co-ordinates)
• preload.
Preload measurements should be made both in air and under water using a
calibrated load cell mounted between the crane wire and the hook. The
foundation steel surface condition (both inside and outside) should also be
documented.
Long-term foundation instrumentation is a growing area for offshore
wind farms, where extensive monitoring systems, both short and long term,
are being used to optimize foundation design (e.g. Devriendt et al., 2013;
Byrne et al., 2015a, 2015b; Shonberg et al., 2017). In general, instrumenta-
tion should also be considered to improve geotechnical knowledge about
foundation response during storm loading (Svanø et al., 1997).
Miscellaneous design considerations 267
12.9 STEEL DESIGN
Structural steel designers should also check the foundation geometry. Critical
items normally include stress concentrations and fatigue in the vicinity of
the load application point, and buckling of the steel cylinder and top plate
under applied underpressure/overpressure.
In addition, foundation steel tip integrity assessments (e.g. Aldridge et al.,
2005; Erbrich et al., 2010) should be considered if there is a risk of the foun-
dation tip locally encountering an isolated gravel (or cobble or boulder-
sized) fragment (see Section 9.10.5), and this fragment cannot be fractured
(i.e. split longitudinally) or pushed to one side (i.e. bearing capacity failure
of the soil surrounding the fragment).
If required, like soil reaction stresses (Section 12.10), steel stresses may be
assessed by 3D non-linear FEA. However, the steel response is non-linear,
whereas it is generally sufficient to assume that the soil remains elastic. Usually
structural analysis software (SACS, Nastran, SESAM, ANSYS, etc.) is used
instead of geotechnical analysis software (e.g. PLAXIS or OptumG2/G3) to
determine steel stresses. Soil continuum elements (not a bed of Winkler springs)
should be placed both inside and outside the foundation. If the complete inter-
mediate foundation is being modelled, then a finite element mesh extending 3L
radially and 2L vertically should be sufficiently accurate (Kuhlemeyer, 1979). If
infinite elements are used, then the finite mesh radial extent can be reduced to
2L (Kay and Palix, 2010). Table 12.1 provides reasonable soil E and ν values.
For model verification, mesh nodes have 6 DOF (loads Fx, Fy, Fz, Mx, My
and Mz and displacements δx, δy, δz, φx, θy and Ψz), Figure 3.1a gives their
sign convention, and the RP is at seafloor on the foundation centreline.
Assuming lateral loading is dominant, a simple way to verify the model is to
run two load cases. The first applies a “fixed head” unit lateral displacement
(δx = 1 mm, remaining 5 DOF = 0) to the RP to find pile head lateral seafloor
stiffness matrix terms kxx (= Fx/ δx) and kxθ (= My/δx). The second load case
finds the missing two terms kθx (= Fx/θy) and kθθ (= My/θy) by applying “fixed
head” unit rotation (θy = 1 milliradian, remaining 5 DOF = 0). Off-diagonal
terms kxθ and kθx should be essentially identical. All four lateral stiffness
matrix terms can then be compared with solutions in Table 11.2 (Elastic
Continuum Solutions – Static Load).
12.10 SOIL REACTIONS
Notation:
E = Young’s Modulus
ν = Poisson’s Ratio
s u = soil undrained shear strength
q c = CPT cone tip resistance
σ c = unconfined compressive strength
Note:
a For laterally loaded piles.
12.11 UNDERBASE GROUTING
At the end of installation, a water-filled void may exist beneath the top plate.
This void generally occurs in sand profiles. The presence of a void is gener-
ally due to using a high estimate stickup value, primarily to ensure that the
target depth could be successfully reached. It is usually not practical to
embed the foundation beyond its target depth. Reasons include (a) the void
may be non-uniform and the (sand) plug may be in contact with the top
plate underside and/or internal stiffeners. In such cases, additional penetra-
tion may induce high underpressures and steel stresses, and (b) the jacket
will be (slightly) lower than intended, which may cause issues with pipelines
and other connections.
Because of the limited foundation embedment ratio (usually less than
one), there is a perceived risk of additional axial settlement occurring
during cyclic storm loads. Hence, underbase grouting is generally applied.
The objective is to replace the water in the void with cement grout.
Miscellaneous design considerations 269
Figure 12.6 A
nchor chain trench 1D wide. V lug -H lug envelopes (Hernandez-
Martinez et al., 2015).
Figure 12.7 A
nchor chain trench D/3 wide. V lug -H lug envelopes and governing ALS
loads (Alderlieste et al., 2016).
TOTAL projects, each with an FPSO and off-loading terminal (OLT) buoy
taut line mooring system. Instead of vertical sided trenches modelled previ-
ously, due to erosion and slope instability, the trenches are wider at the top
than the base, with slope angles increasing to near-vertical for ≈ 3 m above
the base. Figure 12.9a shows a cross-section 20 m away from an anchor pile
(D ≈ 4.5 m, L ≈ 17 m).
Sassi et al. (2018) performed preliminary small scale (75 g–80 g) centri-
fuge tests on D = 4.5 m, L = 15 m, zlug = 9.5 m suction piles jacked into either
blocks extracted from undisturbed core samples (su = 1 + 1.4 z), or reconsti-
tuted clay (slightly higher su profile, gradients between 1.4 kPa/m and 2.0
kPa/m. A spatula, base width = D = 4.5 m (Figure 12.9b) excavated large
trenches reaching the suction anchor and lug level depth. From the first series
of tests (pure vertical pull-out; to model no reverse end-bearing (REB) the
anchor tip was in contact with a sand layer), the presence of a trench had
only a minor reduction (as expected). The “with REB” Qpl (= Fo + Qbase) and
“without REB” Qun (= Fo + Fi) resistances were ≈ 3.0 MN and ≈1.6 MN
respectively. Back-figured average α and Nc,circle values were low (0.3, 8).
272 Miscellaneous design considerations
(a)
(b)
Figure 12.8 A
nchor chain trench – “Fixed Head” incremental plastic strains and
failure mechanism (Alderlieste et al., 2016).
For the remaining two test series, soil pore pressure data were obtained at ≈
0.2 D, 1.2 D and 1.4 D below the tip, and these data confirmed that passive
suction (i.e. REB) was developed even when trenches were present. Anchor
holding capacity reductions of 11% were obtained for θlug between 30o and
40o. Surprisingly, this value is less than FEA by Hernandez-Martinez et al.
(2015) and Arslan et al. (2015) (20% and 45%) who had 1D wide trenches
and included full REB, but similar to that of Alderlieste et al. (2016) (≈ 7%)
who had a smaller trench width (D/3) and excluded REB (i.e. they had a
”plugged” soil response). 3D FEA predictions type C1 (Lambe, 1973) were
also made.
Miscellaneous design considerations 273
(a)
(b)
Figure 12.9 O
bserved anchor chain trench geometry in Gulf of Guinea clay
(Colliat et al., 2018; Sassi et al., 2018).
12.12.1 Commentary
Reverse end-bearing capacity accounts for ≈ 50% of anchor axial pull-out
resistance. Excluding reverse end-bearing capacity adds something of the
order of 2 m length for a 5 m diameter anchor pile (see the Fixed Head Pile
Example in Section 10.7).
To date, most of the aforementioned work has been confined to the Gulf
of Guinea (GoG). It is hoped that this will eventually extend to other areas
of the world where anchor chain trenching has been observed.
One possible explanation for vertical trenches is that the Gulf of Guinea
clays may be similar to those at Bothkennar UK. Both have strong in-situ
structure, higher than usual sensitivity, and are lightly overconsolidated.
Bothkennar soils have been subject to bioturbation by worms and exhibit
rapid strain-softening due to destructuring (Atkinson et al., 1992; ICE,
1992). It is possible that vane shear test, T-bar and CPT testing all underes-
timate Gulf of Guinea clay peak undrained shear strength, su,peak. This may
explain why vertical trench heights have been observed exceeding those
given by classical wedge theory. (Mini) pressuremeter testing, plus a strain-
softening material model (Kay and Van Woensel, 2002) may be useful.
Mitigation measures for anchor chain trenching include:
12.13 CLOSURE
275
276 Acronyms
The terms caisson and foundation, lug and padeye are synonymous.
277
278 Notation
ez,su ∫ s
= analytical ez based on su,z = u,z
0
z dz
∫
/
0
su,z dz
Abadie, C.N., Byrne, B.W. and Houlsby, G.T. (2017) Modelling of monopile response
to cyclic lateral loading in sand. Proc. of the 8th Int. Conf. Offshore Site
Investigation and Geotechnics. OSIG 2017, London, Society of Underwater
Technology, pp. 1046–1053.
Abbs, A.F. (1983) Lateral pile analysis in weak carbonate rocks. Proc. Conference on
Geotechnical Practice in Offshore Engineering: University of Texas at Austin,
Austin, Texas, USA, 27–29 April, pp. 546–556.
Abbs, A.F. and Needham, A.D. (1985) Grouted piles in weak carbonate rocks. Proc.
Offshore Technology Conference, Houston, Texas, USA, Paper OTC-4852-MS.
Achmus, M., Kuo, Y-S. and Abdel-Rahman, K. (2009) Behavior of monopile
foundations under cyclic lateral load. Computers and Geotechnics, Vol. 36,
pp. 725–735.
Advanced Geomechanics (2002) Suction Pile Analysis Code, AGSpanc Version 4.2,
User’s Manual. Doc. No. AGR-1176, Rev 5.
AG (2003) AGSPANC: Suction Pile Analysis Code, User Manual Version 4.2
Advanced Geomechanics, Australia.
AGS (2000) Guidelines for Combined Geoenvironmental and Geotechnical
Investigations, Association of Geotechnical and Geoenvironmental Specialists,
Beckenham, UK.
Ahn, J., Lee, H. and Kim, Y-T. (2014) Holding capacity of suction caisson anchors
embedded in cohesive soils based on finite element analysis. International
Journal for Numerical and Analytical Methods in Geomechanics, Vol. 38, no.
15, pp. 1541–1555.
Alderlieste, E.A. and Dekker, M.J. (2018) Suction pile design and installation chal-
lenges for the Ophir WHP. Proc. 1st Vietnam Symposium on Advances in Offshore
Engineering (VSOE 2018), Hanoi, Vietnam, Springer, pp. 329–335.
Alderlieste, E.A., Romp, R.H., Kay, S. and Lofterød, A. (2016) Assessment of sea-
floor trench for suction pile moorings: a field case. Proc. Offshore Technology
Conference, Houston, Texas, USA, Paper OTC-27035-MS.
Alderlieste, E.A. and van Blaaderen, E.A. (2015) Installation of suction caissons for
an asymmetrical support structure in sandy soil. Proc. 3rd International
Symposium on Frontiers in Offshore Geotechnics, Taylor & Francis, Oslo,
Norway, pp. 215–220.
281
282 References
Aldridge, T.R., Carrington, T.M. and Kee, N.R. (2005) Propagation of pile tip dam-
age during installation. Proc. 1st International Symposium on Frontiers in
Offshore Geotechnics, Taylor & Francis, Perth, Australia, pp. 823–827.
Andersen, K.H. (2015) Cyclic soil parameters for offshore foundation design. Proc.
3rd International Symposium Frontiers in Offshore Geotechnics, June, Oslo,
Norway, pp. 5–82.
Andersen, K.H. and Jostad, H.P. (1999) Foundation design of skirted foundations
and anchors in clay. Proc. Offshore Technology Conference, Houston, Texas,
USA, Paper OTC-10824-MS.
Andersen, K.H. and Jostad, H.P. (2004) Shear strength along inside of suction anchor
skirt wall in clay. Proc. Offshore Technology Conference, Houston, Texas, USA,
Paper OTC-16844-MS.
Andersen, K.H., Jostad, H.P. and Dyvik, R. (2008) Penetration resistance of offshore
skirted foundations and anchors in dense sand. Journal of Geotechical and
Geoenvironmental Engineering, Vol. 134, no. 1, pp. 106–116.
Andersen, K.H., Murff, J.D., Randolph, M.F., Clukey, E.C., Erbrich, C.T., Jostad,
H.P., Jansen, B., Aubeny, C., Sharma, P. and Supachawarote, C. (2005) Suction
anchors for deepwater applications. Proc. 1st International Symposium on
Frontiers in Offshore Geotechnics, Taylor & Francis, Perth, Australia, pp. 3–30.
API (1997) Recommended practice for planning, designing and constructing Tension
Leg Platforms RP 2T, Ed 2 American Petroleum Institute, Washington, DC. (Ed 3
2015).
API (2000) Recommended practice for planning, designing and constructing fixed
offshore platforms - working stress design RP-2A, American Petroleum Institute,
Washington, USA.
API (2005) Recommended practice for design and analysis of station keeping sys-
tems for floating structures, API RP-2SK, American Petroleum Institute,
Washington, USA.
API (2008) Recommended practice for design and analysis of station keeping sys-
tems for floating structures, API RP-2SK, American Petroleum Institute,
Washington, USA.
API (2011) Recommended practice for geotechnical and foundation design consider-
ations RP 2GEO, Ed 1 American Petroleum Institute, Washington, DC. Errata (2014).
Arany, L., Bhattacharya, S., Macdonald, L. and Hogan, S.J. (2017) Design of mono-
piles for offshore wind turbines in 10 steps. Soil Dynamics and Earthquake
Engineering, Vol. 92, pp. 126–152.
ARGEMA (1994) Foundations in Carbonate Soils, English edn, Technip, Paris.
Arslan, H., Peterman, B.R., Wong, P.C. and Bhattacharjee, S. (2015) Remaining
capacity of the suction pile due to seabed trenching. Proc. 25th International
Symposium on Ocean and Polar Engineering Conference, Hawaii, pp. 924–931.
ASTM (2009) Standard practice for description and identification of soils (Visual-
Manual Procedure), ASTM D2488-09a. USA.
ASTM (2011) Standard practice for classification of soils for engineering purposes
(Unified Soil Classification System), ASTM D2487-11. USA.
Atkinson, J.H., Allman, M.A. and Böese, R.J. (1992) Influence of laboratory sample
preparation procedures on the strength and stiffness of intact Bothkennar soil
recovered using the Laval sampler. Géotechnique, Vol. 42, no. 2, pp. 349–354.
Attewell, P.B. and Farmer, I.W. (1976) Principles of engineering geology, Chapman
and Hall, London.
References 283
Baguelin, F., Frank, R. and Saïd, Y.H. (1977) Theoretical study of lateral reaction
mechanism of piles. Géotechnique, Vol. 27, no. 3, pp. 405–434.
Balder, T., de Lange, D.A., Elkadi, A.S.K., Egberts, P.J.P., Beuckelaers, W.J.A.P.,
Coronel, M., van Dijk, J., Atkinson, R. and Luger, H.J. (2020) Hydraulic pile
extraction scale tests (HyPE-ST): experimental design and preliminary results.
Proc. 4th International Symposium on Frontiers in Offshore Geotechnics, article
#3484, publication #1069 (IC-ISFOG21), 8–11 November 2021 (original dates
16–19 August 2020), Austin, Texas, USA, ASCE Geo-Institute and DFI (Proc pub-
lished in August 2020).
Ballard, J-C., Delvosal, O. and Yonatan, P.H. (2010) Simplified VH equations for
foundation punch-through sand into clay. Proc. 2nd International Symposium on
Frontiers in Offshore Geotechnics, Perth, Australia, pp. 650–655.
Barbosa, P. Geduhn, M. Jardine, R.J. Schroeder, F.C. and Horn, M. (2015) Offshore
pile load tests in chalk. Proc. 16th Eur. Conf. Soil Mech. Geotech. Eng, Edinburgh,
Scotland, pp. 2885–2890.
Bhattacharjee, S., Majhi, S.M., Smith, D. and Garrity, R. (2014) Serpentina FPSO
mooring integrity issues and system replacement: unique fast track approach.
Proc. Offshore Technology Conference, Houston, Texas, USA, Paper
OTC-25449-MS.
Bransby, M.F. and Yun, G-J. (2009) Undrained capacity of skirted strip foundations
under combined loading. Géotechnique, Vol. 59, no. 2, pp. 115–125.
Broms, B.B. (1964a) Lateral resistance of piles in cohesive soils. Journal of the Soil
Mechanics and Foundations Division, Proc. American Society of Civil Engineers,
Vol. 90, no. SM2, pp. 27–63.
Broms, B.B. (1964b) Lateral resistance of piles in cohesionless soils. Journal of the
Soil Mechanics and Foundations Division, Proc. American Society of Civil
Engineers, Vol. 90, no. SM3, pp. 123–156.
Broughton, P., Davies, R.L., Aldridge, T. and Carrington, T. (2002) Foundation
design for the refloat of the Maureen steel gravity platform. Proc. ICE, Geotechnical
Engineering Journal, Vol. 155, no. 2, pp. 111–118.
Brown, G.A. and Nacci, V.A. (1971) Performance of hydrostatic anchors in granular
soils. Proc. Offshore Technology Conference, Houston, Texas, USA, Paper
OTC-1472-MS.
BSI (1999) Code of practice for site investigations, British Standard BS 5930:1999.
British Standards Institution, UK.
Burd, H., Byrne, B., McAdam, R., Housby, G., Martin, C., Beuckelaers, W.,
Zdravkovic, L., Taborda, D., Potts, D., Jardine, J., Gavin, K., Doherty, P., Igoe, D.,
Skov Gretlund, J., Pacheco Andrade, M. and Muir Wood, A. (2017) Design aspects
for monopile foundations. Proc. Workshop TC209 19th ICSMGE, Seoul, 2017,
pp. 35–44.
Burd, H., David, M.G., Taborda, D., Zdravkovic, L., Abadie, C., Byrne, B., Housby,
G., Gavins, K.G., Igoe, D.J.P., Jardine, R.J., Martin, C., McAdam, R., Pedro, A.M.G.
and Potts, D.M. (2020) PISA design model for monopiles for offshore wind tur-
bines: application to a marine sand. Géotechnique. doi:10.1680/jgeot.18.P.277.
Burt, N.J. and Harris, R.P. (1980) Design, installation and testing of belled pile foun-
dations. Proc. Offshore Technology Conference, Houston, Texas, USA, Paper
OTC-3872-MS.
Butterfield, R. and Gottardi, G. (1994) A complete three-dimensional failure enve-
lope for shallow footings on sand. Géotechnique, Vol. 44, no. 1, pp. 181–184.
284 References
Bye, A., Erbrich, C., Rognlien, B. and Tjelta, T.I. (1995) Geotechnical design of
bucket foundations. Proc. Offshore Technology Conference, Houston, Texas,
USA, Paper OTC-7793-MS.
Byrne, B., Housby, G., Burd, H., Gavin, K.G., Igoe, D.J.P., Jardine, R.J., Martin, C.,
McAdam, R., Potts, D.M., David, M.G., Taborda, D., Zdravkovic, L. (2020) PISA
design model for monopiles for offshore wind turbines: application to a stiff gla-
cial clay till. Géotechnique. Published Online doi:10.1680/jgeot.18.P.255.
Byrne, B.W. (2017) New design methods for offshore wind turbine monopiles. Proc
8th International Conference on Offshore Site Investigation and Geotechnics,
held 12–14 September 2017, Royal Geographical Society, London.
Byrne, B.W., McAdam, R., Burd, H.J., Houlsby, G.T., Martin, C.M., Zdravkovic, L.,
Taborda, D.M.G., Potts, D.M., Jardine, R.J., Sideri, M., Schroeder, F.C., Gavin, K.,
Doherty, P., Igoe, D., Muir Wood, A., Kallehave, D. and Skov Gretlund, J. (2015a)
New design methods for large diameter piles under lateral loading for offshore
wind applications. Proc. 3rd International Symposium Frontiers in Offshore
Geotechnics, June, Oslo, Norway, pp. 705–710.
Byrne, B.W., McAdam, R., Burd, H.J., Houlsby, G.T., Martin, C.M., Gavin, K.,
Doherty, P., Igoe, D., Zdravkovic, L., Taborda, D.M.G., Potts, D.M., Jardine, R.J.,
Sideri, M., Schroeder, F.C., Muir Wood, A., Kallehave, D. and Skov Gretlund, J.
(2015b) Field testing of large diameter piles under lateral loading for offshore
wind applications. Proc. XVI ECSMGE Geotechnical Engineering for
Infrastructure and Development, Edinburgh, Scotland, pp. 1255–1260.
Carbon Trust (2019) Suction installed caisson foundations for offshore wind: design
guidelines. Carbon Trust. https://prod-drupal-files.storage.googleapis.com/docu-
ments/resource/public/owa-suction-caisson-design-guidelines-report.pdf.
Carter, J.P. and Kulhawy, F.H. (1992) Analysis of laterally loaded shafts in rock.
ASCE Journal of Geotechnical Engineering, Vol. 118, no. 6, pp. 839–855.
Cathie Associates (2014) OPILE instruction manual, Cathie Associates SA/NV,
Belgium.
CEN (2004) Eurocode 7: geotechnical design - part 1: general rules, European
Standard EN 1997-1:2004. (With Corrigendum EN 1997-1:2004/AC, February
2009, European Committee for Standardization.
Cerfontaine, B., Collin, F. and Charlier, R. (2016) Numerical modelling of transient
cyclic vertical loading of suction caissons in sand, Géotechnique, Vol. 66, no. 2,
pp. 121–136.
CFMS (2019) Recommendations for planning and designing foundations of offshore
wind turbines, French workgroup ‘Foundations of Offshore Wind Turbines,
Version 2.1. http://www.cfms-sols.org/sites/default/files/Rapport-cfms-eoliennes-
offshore-EN.pdf.
Chan, S.F. and Hannah, T.H. (1980) Repeated loading on single piles in sand. ASCE
Journal of Geotechnical and Geoenvironmental Engineering, Vol. 106, no. GT2,
pp. 171–188.
Chatzivasileiou, I. (2014) Installation of suction caissons in layered sand, MSc.
Thesis, Delft University of Technology, The Netherlands.
Chen, W. and Randolph, M.F. (2007a) External radial stress changes and axial capac-
ity for suction caisson in soft clay. Géotechnique, Vol. 57, no. 6, pp. 499–511.
Chen, W. and Randolph, M.F. (2007b) Uplift capacity of suction caissons under
sustained and cyclic loading in soft clay. ASCE Journal of Geotechnical and
Geoenvironmental Engineering, Vol. 133, no. 11, pp. 1352–1363.
References 285
Devriendt, C., Van Ingelgem, Y., De Sitter, G., Jordaens, P.J. and Guillaume, P. (2013)
Monitoring of resonant frequencies and damping values of an offshore wind tur-
bine on a monopile foundation, Offshore Wind Infrastructure Application Lab,
Belgium, presentation to EWEA 2013.
Dijkhuizen, C., Coppens, T, and van der Graaf, P. (2003) Installation of the Horn
Mountain Spar using the enhanced DCV Balder. Proc. Offshore Technology
Conference, Houston, Texas, USA, Paper OTC-16367-MS.
DNV (1992) Foundations, Classification Notes No. 30.4. Det Norske Veritas.
DNV (2008) Design of Offshore Steel Structures, General (LRFD Method), Offshore
Standard DNV-OS-C101. Det Norske Veritas.
DNV (2014a) Design of Offshore Wind Turbine Structures, Offshore Standard
DNV-OS-J101. Det Norske Veritas.
DNV (2014b) Modelling and Analysis of Marine Operations, Recommended
Practice DNV-RP-H103. Det Norske Veritas.
DNV (2018) Geotechnical Design and Installation of Suction Anchors in Clay,
Recommended Practice DNV-RP-E303. Det Norske Veritas.
DNV-GL (2016) Support Structures for Wind Turbines, Standard DNVGL-ST-0126.
DNV-GL (2017, 2019) Offshore Soil Mechanics and Geotechnical Engineering,
Recommended Practice DNVGL-RP-C212.
Doherty, J.P. and Deeks, A.J. (2003) Elastic response of circular footings embedded
in a nonhomogeneous half-space. Géotechnique, Vol. 53, no. 8, pp. 703–714.
Doherty, J.P., Houlsby, G.T. and Deeks, A.J. (2005) Stiffness of flexible caisson foun-
dations embedded in nonhomogeneous elastic soil. ASCE Journal of Geotechnical
and Geoenvironmental Engineering, Vol. 131, no. 12, pp. 1498–1508.
Dührkop, J. and Grabe, J. (2009) Design of laterally loaded piles with bulge. Proc.
28th International Conference on Ocean, Offshore and Arctic Engineering,
ASME, Honolulu, Hawaii, USA. OMAE2009‑79087.
Dührkop, J., Maretzki, S. and Rieser, J. (2017) Re-evaluation of pile driveability in
chalk. Proc. of the 8th Int. Conf. Offshore Site Investigation and Geotechnics.
OSIG 2017, Society of Underwater Technology, London, Vol. 2, pp. 666–673.
Dyson, G.J. and Randolph, M. (2001) Monotonic lateral loading of piles in calcare-
ous sediments. ASCE Journal of Geotechnical and Geoenvironmental Engineering,
Vol. 127, no. 24, pp. 346–352.
Edwards, D.H., Zdravkovic, L. and Potts, D.M. (2005) Depth factors for undrained
bearing capacity. Géotechnique, Vol. 55, no. 10, pp. 755–758.
Engin, H.K., Khoa, H.D.V., Jostad, H.P. and Alm, T. (2015) Finite element analyses
of spudcan – subsea template interaction during jack-up rig installation. Proc. 3rd
International Symposium Frontiers in Offshore Geotechnics, June, Oslo, Norway,
705, pp. 1275–1280.
Erbrich, C.T. (2004) A new method for the design of laterally loaded anchor piles in
soft rock. Proc. Offshore Technology Conference, Houston, Texas, USA, Paper
OTC-16441-MS.
Erbrich, C.T., Barbosa-Cruz, E. and Barbour, R. (2010) Soil-pile interaction during
extrusion of an initially Deformed Pile. Proc. 2nd International Symposium on
Frontiers in Offshore Geotechnics, Perth, Australia, pp. 489–494.
Erbrich, C.T. and Hefer, P. (2002) Installation of the Laminaria suction piles – A case
history. Proc. Offshore Technology Conference, Houston, Texas, USA, Paper
OTC-14220-MS.
References 287
Erbrich, C.T., Lam, S.Y., Zhu, H., Derache, A., Sato, A. and Al-Showaiter, A. (2017)
Geotechnical design of anchor piles for Ichthys CPF and FPSO. Proc. of the 8th
Int. Conf. Offshore Site Investigation and Geotechnics. OSIG 2017, Society of
Underwater Technology, London, pp. 1186–1197.
Erbrich, C.T. and Tjelta, T.I. (1999) Installation of bucket foundations and suction
caissons in sand – geotechnical performance. Proc. Offshore Technology
Conference, Houston, Texas, USA, Paper OTC-10990-MS.
Feld, T. (2001) Suction buckets: a new innovation foundation concept, applied to
offshore wind turbines, PhD. Thesis, Aalborg University, Denmark.
Feng, X., Gourvenec, S. and White, D.J. (2019) Load capacity of caisson anchors
exposed to seabed trenching. Ocean Engineering, Vol. 171, no. 1, pp. 181–192.
Feng, X., Randolph, M.F., Gourvenec, S. and Wallerand, R. (2014) Design approach
for rectangular mudmats under fully three dimensional loading. Géotechnique,
Vol. 64, no. 1, pp. 51–63.
Finnie, I.M.S. and Randolph, M.F. (1994) Punch-through and liquefaction induced
failure of shallow foundations on calcareous sediments. Proc. 7th International
Conference on the Behaviour of Offshore Structures, BOSS ’94, Vol. 1, Pergamon,
Oxford, pp. 217–230.
Fisher, R. and Cathie, D. (2003) Optimisation of gravity based design for subsea appli-
cations. Proc. 2nd British Geotechnical Association International Conference on
Foundations, 2–5 September, University of Dundee, Scotland, UK, IHS BRE Press.
Fleshman, M.S. and Rice, J.D. (2014) Laboratory modeling of the mechanisms of
piping erosion initiation. ASCE Journal of Geotechnical and Geoenvironmental
Engineering, Vol. 140, no. 6. doi:10.1061/(ASCE)GT.1943-5606.0001106.
Fragio, A.G., Santiago, J.L. and Sutton, V.J.R. (1985) Load tests on grouted piles in
rock. Proc. Annual Offshore Technology Conference, Houston, Texas, USAs,
Paper OTC-4851-MS.
Frankenmolen, S.F., Erbrich, C.T. and Fearon, R.E. (2017) Successful installation of
large suction caissons and driven piles in carbonate soils. Proc. of the 8th Int.
Conf. Offshore Site Investigation and Geotechnics. OSIG 2017, London, Society
of Underwater Technology, pp. 539–548.
Fuglsang, L.D. and Steensen-Bach, J.O. (1991) Breakout resistance of suction piles in
clay. Centrifuge 91: Proc. International conference centrifuge, Boulder, Colorado,
13–14 June 1991, A.A. Balkema, Rotterdam, pp. 153–159.
Fugro (2009) CANCAP2 User’s Manual, Fugro Engineers BV Internal Document
No. FEBV/CDE/MAN/083.
Gazetas, G. (1991) Foundation vibrations. In: Foundation engineering handbook,
Fang, H-Y. (Ed.), 2nd edn, Van Nostrand Reinhold Company, New York, pp.
553–593.
Germanischer Lloyd (2013) Rules for the certification and construction, IV Industrial
services, Part 7 Offshore substations, Chapter 2 Structural Design, Germanischer
Lloyd, Hamburg.
Germanischer Lloyd Wind Energie GmbH (2005) Rules and guidelines: IV –
Industrial services, Part 2 – Guideline for the certification of offshore wind tur-
bines, Germanischer Lloyd, Hamburg.
Gerolymos, N., Zafeirakos, A. and Karapiperis, K. (2015) Generalized failure enve-
lope for caisson foundations in cohesive soil: static and dynamic loading. Soil
Dynamics and Earthquake Engineering, Vol. 78, pp. 154–174.
288 References
Gilbert, R.B. and Murff, J.D. (2001) Identifying uncertainties in the design of suction
caisson foundations. Proc. International Conference on Geotechnical, Geological
and Geophysical Properties of Deepwater Sediments honoring Wayne A. Dunlap,
OTRC, Houston, Texas, pp. 231–242.
Goodman, L.J., Lee, C.N. and Walker, F.J. (1961) The feasibility of vacuum anchor-
age in soil. Géotechnique, Vol. 1, no. 3, pp. 356–359.
Gottardi, G., Houlsby, G.T. and Butterfield, R. (1999) Plastic response of circular
footings on sand under general planar loading. Géotechnique, Vol. 49, no. 4, pp.
453–469.
Gourvenec, S. and Jensen, K. (2009) Effect of embedment and spacing on co-joined
skirted foundation systems on undrained limit states under general loading.
International Journal of Geomechanics, Vol. 9, no. 6, pp. 267–279.
Gourvenec, S. and Steinepreis, M. (2007) Undrained limit states of shallow founda-
tions acting in consort. ASCE, International Journal of Geomechanics, Vol. 7, no.
3, pp. 194–205.
Gourvenec, S.M., Feng, X., Randolph, M.F. and White, D.J. (2017) A toolbox for
optimizing geotechnical design of subsea foundations. Proc. Annual Offshore
Technology Conference, Houston, Texas, USA Paper OTC-27703-MS.
Govoni, L, Gottardi, G., Antoninia, A., Archettia, R. and Schweizer, J. (2016) Caisson
foundations for competitive offshore wind farms in Italy. Procedia Engineering,
Vol. 158, pp. 392–397.
Griffiths, D.V. (1982) Elasto-plastic analyses of deep foundations in cohesive soil.
International Journal for Numerical and Analytical Methods in Geomechanics,
Vol. 6, pp. 211–218.
Griffiths, D.V. (1985) HARMONY – A program for predicting the elasto-plastic
response of axisymmetric Bodies subjected to non-axisymmetric loading, Report
to Fugro-McClelland Engineers B.V., University of Manchester, UK.
Guo, W., Chu, J. and Kou, H. (2016) Model tests of soil heave plug formation in suc-
tion caisson. Proc. Institution of Civil Engineers, Geotechnical Engineering, Vol.
169, no. GE2, pp. 214–223.
Gupta, B.K. and Basu, D. (2016) Design charts for laterally loaded rigid piles in elas-
tic soil. Proc. ASCE Geo-Chicago 2016 GSP 270, Chicago, Illinois, USA, pp.
393–406.
Hansteen, O.E., Jostad, J.P. and Tjelta, T.I. (2003) Observed platform response to a
“monster” wave. Proc. 6th Int. Symposium on Field Measurements in Geomechanics,
Swets and Zeitlinger, Lisse, The Netherlands.
Heins, E. and Grabe, J. (2017) Class-A-prediction of lateral pile deformation with
respect to vibratory and impact pile driving. Computers and Geotechnics, Vol. 86,
pp. 108–119.
Helfrich, S.C., Brazill, R.L. and Richards, A.F. (1976) Pullout characteristics of a suc-
tion anchor in sand. Proc. of Offshore Technology Conference, Houston, Texas,
USA, Paper OTC-2469-MS.
Hernandez-Martinez, F.G., Rahim, A., Strandvik, S., Jostad, H.P. and Andersen, K.H.
(2009) Consolidation settlement of skirted foundations for subsea structures in
soft clay. Proc. 17th International Conference on Soil Mechanics and Geotechnical
Engineering, pp. 1167–1172.
Hernandez-Martinez, F.G., Saue, M., Schroder, K. and Jostad, H.P. (2015) Trenching
effects on holding capacity for In-service suction anchors in high plasticity clays.
References 289
Proc. 20th Offshore Symposium, Texas Section of the Society of Naval Architects
and Marine Engineers, Houston, Texas, USA.
Hogervorst, J.R. (1980) Field trials with large diameter suction piles. Proc. of
Offshore Technology Conference, Houston, Texas, USA, Paper OTC-3817-MS.
Holeyman, A., Peralta, P. and Charue, N. (2015) Boulder-soil-pile dynamic interac-
tion. Proc. 3rd International Symposium on Frontiers in Offshore Geotechnics,
Taylor & Francis, Oslo, Norway, pp. 563–568.
Hölscher, P., van Tol, A.F. and Huy, N.Q. (2009) Influence of rate effect and pore
water pressure during rapid load test of piles in sand. Proc. Rapid Load Testing on
Piles, van Tol, F. and Hölscher, P. (Eds.), CRC Press, pp. 59–72.
Holtz, R.D. (1991) Stress distribution and settlement of shallow foundations. In:
Foundation Engineering Handbook, Winterkorn, H.F. and Fang, H.Y. (Eds.), 2nd
edn, Van Nostrand Reinhold Company, New York, pp. 166–222.
Horvath, R.G. (1978) Field load test data on concrete to rock bond strength,
University of Toronto, Canada, Publication No. 78-07. 23.
Horvath, R.G. and Kenney, T.C. (1979) Shaft resistance of rock-socketed drilled
piers. Symposium on Deep Foundations, ASCE National Convention, Atlanta,
Georgia, USA, pp. 182–214.
Hossain, M., Lehane, B., Hu, Y. and Gao, Y. (2012) Soil flow mechanisms around
and between stiffeners of caissons during installation in clay. Canadian
Geotechnical Journal, Vol. 49, no. 4, pp. 442–459.
Houlsby, G.T. and Byrne, B.W. (2005) Design procedures for installation of suction
caissons in sand. Proc. Institution of Civil Engineers, Geotechnical Engineering,
Vol. 158, no. GE3, Paper 13818, pp. 135–144.
Houlsby, G.T., Kelly, R.B. and Byrne, B.W. (2005) The tensile capacity of suction
caissons under rapid loading. Proc. 1st International Symposium on Frontiers in
Offshore Geotechnics, Taylor & Francis, Perth, Australia, pp. 405–410.
Houlsby, G.T., Kelly, R.B., Huxtable, J. and Byrne, B.W. (2006) Field trials of suc-
tion caissons in sand for offshore wind turbine. Géotechnique, Vol. 56, no. 1,
pp. 3–10.
Houlsby, G.T. and Martin, C.T. (2003) Undrained bearing capacity factors for coni-
cal footings on clay. Géotechnique, Vol. 53, no. 5, pp. 513–520.
Houlston, P., Manceau, S., Benson, A. Alobaidi, I., Evans, R., Pacheco Andrade, M.
and Meincke, J.H. (2017) Laboratory testing based design method for assessing
cyclic response of driven piles in chalk. Proc. of the 8th Int. Conf. Offshore Site
Investigation and Geotechnics. OSIG 2017, Society of Underwater Technology,
pp. 449–456.
Huy, N.Q., van Tol, A.F. and Hölscher, P. (2009) Rapid model pile load tests in the
geotechnical centrifuge. Proc. Rapid Load Testing on Piles, van Tol, F. and
Hölscher, P. (Eds.), CRC Press, pp. 103–127.
Ibsen, C.B. and Thilsted, C.L. (2010) Numerical study of piping limits for suction
installation of offshore skirted foundations and anchors in layered sand. Proc. 2nd
International Symposium on Frontiers in Offshore Geotechnics, Perth, Australia,
pp. 421–426.
ICE (1992) Bothkennar soft clay test site: characterisation and lessons learned, 8th
Géotechnique Symposium in Print. Géotechnique, Vol. 42, no. 2, pp. 163–375.
IHC (2010) IHC Hydrohammer. Pile Driving Equipment, IHC Merwede, The
Netherlands, 13 pp.
290 References
Irvine, J., Terente, V., Lee, L.T. and Comrie, R. (2015) Driven pile design in weak
rock. Proc. 3rd International Symposium on Frontiers in Offshore Geotechnics,
Taylor & Francis, Oslo, Norway, pp. 569–574.
ISO (2006) Petroleum and natural gas industries – fixed concrete offshore structures,
International Organization for Standardization, ISO 19903-1:2019.
ISO (2007) Petroleum and natural gas industries – fixed steel offshore structures,
International Organization for Standardization, ISO 19902:2007.
ISO (2011) Petroleum and natural gas industries – specific requirements for offshore
structures – part 7: stationkeeping systems for floating offshore structures and
mobile offshore units, International Organization for Standardization, ISO/DIS
19901-7.
ISO (2014a) Petroleum and natural gas industries – specific requirements for off-
shore structures – part 7: stationkeeping systems for floating offshore structures
and mobile offshore units, ISO 19901-7:2013.
ISO (2014b) Petroleum and natural gas industries – specific requirements for off-
shore structures – part 8: marine soil investigations International Organization for
Standardization, ISO 19901-8:2014.
ISO (2016a) Petroleum and natural gas industries – specific requirements for off-
shore structures – part 4: geotechnical and foundation design considerations,
International Organization for Standardization, ISO 19901-4:2016.
ISO (2016b) Petroleum and natural gas industries – site specific assessment of mobile
offshore units - part 1: jack-ups, International Organization for Standardization,
ISO 19905-1:2016(E).
Janbu, N. (1985) Soil models in offshore engineering. Géotechnique, Vol. 35, no. 3,
pp. 241–281.
Jardine, R., Chow, F.C., Overy, R.F. and Standing, J.R. (2005) ICP design methods
for driven piles in sands and clays, Thomas Telford Ltd., London.
Jardine, R.J., Buckley, R.M., Byrne, B.W., Kontoe, S., McAdam, R.A., Andolfsson, T.,
Liu, T., Schranz, F. and Vinck, K. (2019) Improving the design of piles driven in
chalk through the ALPACA research project. Revue Française de Géotechnique,
Vol. 158, no. 2. doi:10.1051/geotech/2019008.
Jardine, R.J., Standing, J.R. and Chow, F.C. (2006) Some observations of the effects
of time on the capacity of piles driven in sand. Géotechnique, Vol. 56, no. 4, pp.
227–244.
Jeanjean, P. (2006) Set-up characteristics of suction anchors for soft Gulf of Mexico
clays: experience from field installation and retrieval. Proc. Offshore Technology
Conference, Houston, Texas, USA, Paper OTC-18005-MS.
Jeanjean, P., Znidarcic, D., Phillips, R., Ko, H-Y., Pfister, S., Cinicioglu, O. and
Schroeder, K. (2006) Centrifuge testing on suction anchors: Double-wall, over-
consolidated clay, and layered soil profile. Proc. Offshore Technology Conference,
Houston, Texas, USA, Paper OTC-18007-MS.
Jostad, H.P. (1997) BIFURC - Version 3 Undrained Capacity Analyses of Clay,
Technical Report, Norwegian Geotechnical Institute.
Kallehave, D., LeBlanc, C., Liingaard, M.A. (2012) Modification of the API P-Y
formulation of initial stiffness of sand, Integrated Geotechnologies - Present and
Future. Proc. of the 7th Int. Conf. Offshore Site Investigation and Geotechnics.
OSIG 2017, Society of Underwater Technology, London, pp. 465–472.
Karlsrud, K., Jensen, T.G., Lied, E.K.W., Nowacki, F. and Simonsen, A.S. (2014)
Significant ageing effects for axially loaded piles in sand and clay verified by new
References 291
field load tests. Proc. Offshore Technology Conference, Houston, Texas, USA,
Paper OTC-25197-MS.
Kay, S. (2013) Torpedo piles – VH capacity in clay using resistance envelope equa-
tions. Proc. 32nd International Conference on Ocean, Offshore and Arctic
Engineering, OMAE2013, 9–14 June, Nantes, France.
Kay, S. (2015) CAISSON_VHM: a program for caisson capacity under VHM Load
in undrained soils. Proc. 3rd International Symposium on Frontiers in Offshore
Geotechnics, Taylor & Francis, Oslo, Norway, pp. 277–282.
Kay, S. (2018) User’s Manual CAISSON_VHM Version 00.12, Computer Program
for Caisson Foundation Resistance in Undrained Soils during VHM Load, S.Kay
Consultant Report SKA/MAN/003, March, 850 pp.
Kay, S. and Palix, E. (2010) Caisson capacity in clay: VHM resistance envelope Part
2: VHM envelope equation and design procedures. Proc. 2nd International
Symposium on Frontiers in Offshore Geotechnics, Perth, Australia.
Kay, S. and Palix, E. (2011) Caisson capacity in clay: VHM Resistance Envelope –
Part 3: extension to shallow foundations. Proc. 30th International Conference on
Ocean, Offshore and Arctic Engineering, OMAE2011, 19–24 June, Rotterdam,
the Netherlands.
Kay, S. and Van Woensel, F. (2002) Cone pressuremeter tests in clay. Proc. 5th
European Conference Numerical Methods in Geotechnical Engineering, Paris,
4–6 September 2002, Laboratoire Central des Ponts et Chaussées (LCPC) and
Presses de l’École Nationale des Ponts et Chaussées, Paris, pp. 99–105.
Keaveny, J.V., Hansen, S.B., Madshus, C. and Dyvik, R. (1994) Horizontal capacity
of large-scale model anchors. Proc. 13th Int. Conf. Soil Mechanics and Foundation
Engineering, New Delhi, India, pp. 677–680.
Kellezi, L., Kudsk, G. and Hansen, P.B. (2005) FE modelling of spudcan-pipeline
interaction. Proc. 1st International Symposium on Frontiers in Offshore
Geotechnics, Taylor & Francis, Perth, Australia, pp. 551–557.
Kelly, R.B., Byrne, B.W., Houlsby, G.T. and Martin, C.M. (2003) Pressure chamber
testing of model caisson foundations in sand. Proc. 2nd British Geotechnical
Association International Conference on Foundations, 2–5 September, University
of Dundee, Scotland, UK, IHS BRE Press, pp. 421–432.
Kelly, R.B., Byrne, B.W., Houlsby, G.T. and Martin, C.M. (2004) Tensile loading of
model caisson foundations for structures on sand. Proc. 14th International
Offshore and Polar Engineering Conference, Toulon, Vol. 2, pp. 638–641.
Kelly, R.B., Houlsby, G.T. and Byrne, B.W. (2006a) A comparison of field and labo-
ratory tests of caisson foundations in sand and clay. Géotechnique, Vol. 56, no. 9,
pp. 617–626.
Kelly, R.B., Houlsby, G.T. and Byrne, B.W. (2006b) Transient vertical loading of model
suction caissons in a pressure chamber. Géotechnique, Vol. 56, no. 10, pp. 665–675.
Kennedy, J., Oliphant, J., Maconochie, A., Stuyts, B. and Cathie, D. (2013) CAISSON:
a suction pile design tool. Proc. 32nd International Conference on Ocean,
Offshore and Arctic Engineering, OMAE2013, 9–14 June, Nantes, France.
Kirsch, F., Richter, T., Coronel, M. (2014) Geotechnische Aspekte bei der
Gründungsbemessung von Offshore-Windenergieanlagen auf Monopfählen mit
sehr großen Durchmessern, Stahlbau Spez. 2014 – Erneuerbare Energien 83.
Kolk, H.J., Kay, S., Kirstein, A. and Troestler, H. (2001) North Nemba flare bucket
foundations. Proc. Offshore Technology Conference, Houston, Texas, USA, Paper
OTC-13057-MS.
292 References
Kolk, H.J. and Wegerif, J. (2005) Offshore site investigations: new frontiers. Proc.
1st International Symposium on Frontiers in Offshore Geotechnics, Taylor &
Francis, Perth, Australia, pp. 145–161.
Kuhlemeyer, R. (1979) Static and dynamic laterally loaded floating piles. ASCE
Journal of Geotechnical and Geoenvironmental Engineering, Vol. 105, 289–304.
Lambe, T.W. (1973) Predictions in soil engineering. Géotechnique, Vol. 23, no. 2, pp.
151–201.
Langford, T.E., Solhjell, E., Hampson, K. and Hondebrink, L. (2012) Geotechnical
design and installation of suction anchors for the Skarv FPSO, Offshore Norway.
Proc. 7th Int. Conf. Offshore Site Investigation and Geotechnics. OSIG 2012,
Society of Underwater Technology, London, pp. 613–620.
LeBlanc, C., Houlsby, G.T., Byrne, B.W. (2010) Response of stiff piles in sand to
long-term cyclic lateral loading. Géotechnique, Vol. 60, no. 2, pp. 79–90.
Lee, Y.C., Audibert, J.M.E. and Tjok, K.-M. (2005) Lessons learned from several
suction caisson installation projects in clay. Proc. 1st International Symposium on
Frontiers in Offshore Geotechnics, Taylor & Francis, Perth, Australia, pp.
235–241.
Lehane, B.M., Elkhatib, S. and Terzaghi, S. (2014) Extraction of suction caissons in
sand, Géotechnique, Vol. 64, no. 9, pp. 735–739.
Lehane, B.M., Schneider, J.A. and Xu, X. (2005) A review of design methods for
offshore driven piles in siliceous sand, University of Western Australia, Perth,
UWA Report, GEO:05358.
Lin, S.-S. and Liao, J.-C. (1999) Permanent strains of piles in sand due to cyclic lat-
eral loads. ASCE Journal of Geotechnical and Geoenvironmental Engineering,
Vol. 125, No. 9, pp. 798–802.
Lovera, A. (2019) Cyclic lateral design for offshore monopiles in weak rocks. PhD
Thesis, Université Paris-Est, Environmental Engineering, France. NNT:
2019PESC1016.
Lovera, A., Ghabezloo, S., Sulem, J., Randolph, M.F. and Palix, E. (2020) Extension
of the p-y curves framework for cyclic loading of offshore wind turbines mono-
piles. Proc. 4th International Symposium on Frontiers in Offshore Geotechnics,
Austin, Texas, USA.
Lunne, T., Robertson, P.K. and Powell, J.J.M. (1997) Cone penetration testing in
geotechnical practice, Blackie Academic & Professional, London.
Lunne, T. and St John, H.D. (1979) The use of cone penetrometer tests to compute
penetration resistance of steel skirts underneath North Sea gravity platforms.
Proc. 7th European Conference on Soil Mechanics and Foundation Engineering,
Brighton, England, September, Vol. 2, British Geotechnical Society, London, pp.
233–238.
Lunne, T.A., Myrovll, F. and Kleven, A. (1982) Observed settlements of five North
Sea gravity platforms, pp. 25, NGI Publication 139. ISBN: 0078-1193. A
Norwegian Geotechnical Institute (NGI) publication, Norway.
Madsen, S., Andersen, L.V. and Ibsen, L.B. (2012) Instability during installation of
foundations for offshore structures. Proc. 16th Nordic Geotechnical Meeting,
Copenhagen, Denmark, Vol. 2, Dansk Geoteknisk Forening, pp. 499–505.
Mana, D.S.K., Gourvenec, S. and Randolph, M.F. (2013) A novel technique to miti-
gate the effect of gapping on the uplift capacity of offshore shallow foundations.
Géotechnique, Vol. 63, no. 14, pp. 1245–1252.
References 293
Martinelli, M., Alderlieste, E.A., Galavi, V. and Luger, H.J. (2020) Numerical simula-
tion of the installation of suction bucket foundations using MPM. Proc. 4th
International Symposium on Frontiers in Offshore Geotechnics, Austin, Texas,
USA.
McManus, K.J. and Davis, R.O. (1997) Dilation-induced pore fluid cavitation in
sands. Géotechnique, Vol. 47, no. 1, pp. 173–177.
McNamara, A.P. (2000) Behaviour of soil around the internal stiffeners of suction
caissons and the effect on installation and pullout resistance, Geomechanics
Group Report No. G1512, University of Western Australia, Australia, June.
Meigh, A.C. and Wolski, W. (1979) Design parameters for weak rock. Proc. 7th
European Conference on Soil Mechanics and Foundation Engineering, British
Geotechnical Society, London, UK.
Metrikine, A.V., Tsouvalas, A, Segeren, M.L.A., Elkadi, A.S.K., Tehrani, F.S., Gómez,
S.S., Atkinson, R., Pisanò, F., Kementzetzidis, E., Tsetas, A., Molenkamp, T., van
Beek, K. and de Vries, P. (2020) GDP: a new technology for gentle driving of
(mono)piles. Proc. 4th International Symposium on Frontiers in Offshore
Geotechnics, article #3483; publication #1069 (IC-ISFOG21), 8–11 November
2021, Austin, Texas, USA, ASCE Geo-Institute and DFI (Proc. published in August
2020).
Mirza, U.A., Sweeney, M. and Dean, A.R. (1988) Potential effects of jack-up spudcan
penetration on jacket piles. Proc. of Offshore Technology Conference, Houston,
Texas, USA, OTC-5762-MS.
Muir Wood, A., Mackenzie, B., Burbury, D., Rattley, M., Clayton, C.R.I., Mygind,
M., Wessel Andersen, K., LeBlanc Thilsted, C. and Alberg Liingaard, M. (2015)
Design of large diameter monopiles in chalk at Westermost Rough offshore wind
farm. Proc. 3rd International Symposium on Frontiers in Offshore Geotechnics,
Taylor & Francis, Oslo, Norway, pp. 723–728.
Murff, J.D. (1975) Response of axially loaded piles. Journal of the Geotechnical
Engineering Division, Proceedings of the American Society of Civil Engineers, Vol.
101, no. GT3, pp. 356–360.
Murff, J.D. (1994) Limit analysis of multi-footing foundation systems. Proc. Int.
Conf. on Computer Methods and Advances in Geomechanics, Morgantown, 1,
pp. 223–244.
Murff, J.D. and Hamilton, J.M. (1993) P-ultimate for undrained analysis of laterally
loaded pile’. ASCE Journal of Geotnical Engineering, Vol. 119, no. 1, pp.
91–107.
Negro, V., López-Gutiérrez, J-S., Dolores Esteban, M., Alberdi, P., Imaz, M. and
Serraclara, J-M. (2017) Monopiles in offshore wind: preliminary estimate of main
dimensions. Ocean Engineering, Vol. 133, pp. 253–261.
Neubecker, S.R. and Randolph, M.F. (1995) Profile and frictional capacity of embed-
ded anchor chain. ASCE Journal of Geotnical Engineering, Vol. 121, no. 11, pp.
787–803.
Newlin, J.A. (2003a) Suction anchor piles for the Na Kika FDS mooring system. Part
1: site characterization and design. Proc. International Symposium on Deepwater
Mooring Systems: Concepts, Design, Analysis and Materials, American Society of
Civil Engineers, Houston, Texas, USA, pp. 28–54.
Newlin, J.A. (2003b) Suction anchor piles for the Na Kika FDS mooring system. Part
2: site installation performance. Proc. International Symposium on Deepwater
294 References
the 8th Int. Conf. Offshore Site Investigation and Geotechnics. OSIG 2017,
Society of Underwater Technology, London, pp. 570–577.
PanGeo (2010) Mitigating risk for offshore installations, presentation to The
Hydographic Society, UK, 21 October 2010, 54 pp. https://www.ths.org.uk/event_
details.asp?v0=234.
Peck, R.B. (1969) Advantages and limitations of the observational method in applied
soil mechanics. Géotechnique, Vol. 19, no. 2, pp. 171–187.
Pells, P.J.N., Rowe, R.K. and Turner, R.M. (1980) An experimental investigation into
side shear for socketed piles in sandstone. Proc. International Conference on
Structural Foundations on Rock, Sydney, 7–9 May, Balkema, Vol. 1, pp. 291–302.
Pile Dynamics, Inc. (2010) GRLWEAP, wave equation analysis of pile driving,
Computer Program Package.
Poulos, H.G. and Davis, E.H. (1980) Pile foundation analysis and design, John Wiley
and Sons, New York, Series in Geotechnical Engineering, 410 pp.
Poulos, H.G. and Davis, E.H. (1991) Elastic solutions for soil and rock mechanics,
Series in Soil Engineering, Originally published 1974 by John Wiley & Sons (ISBN
10: 0471695653 ISBN 13: 9780471695653), republished in 1991 by Centre for
Geotechnical Research, Sydney University. Available online http://www.usucger.
org/PandD/PandD.htm.
Poulos, H.G. and Randolph, M.F. (1983) Pile group analysis: a study of two
methods. ASCE, Journal of Geotechnical Engineering, Vol. 109, no. 3, March,
pp. 355–372.
Puech, A., Becue, J.P. and Colliat, J.-L. (1988) Advances in the design of piles driven
into non-cemented to weakly cemented carbonate formations. Proc. 1st
International Conference on Calcareous Sediments, 15–18 March 1988, Perth,
A.A. Balkema, Rotterdam, pp. 305–312.
Puech, A., Colliat, J-L., Nauroy, J.F. and Meunier, J. (2005) Some geotechnical speci-
ficities of Gulf of Guinea deepwater sediments. Proc. 1st International Symposium
on Frontiers in Offshore Geotechnics, Taylor & Francis, Perth, Australia, pp.
1047–1053.
Puech, A., Poulet, D. and Boisard, P. (1990) A procedure to evaluate pile drivability
in the difficult soil conditions of the southern part of the Gulf of Guinea. Proc.
Offshore Technology Conference, Houston, Texas, USA, Paper OTC-6237-MS.
Puech, A. and Quiterio-Mendoza, B. (2019) Caractérisation des massifs rocheux
pour le dimensionnement de pieux forés en mer (Characterization of rock masses
for designing drilled and grouted offshore pile foundations). Revue Française de
Géotechnique, Vol. 158, no. 5. doi:10.1051/geotech/2019011.
Randolph, M.F. (1981) Piles subjected to torsion. Journal of the Geotechnical
Engineering Division, ASCE, Vol. 107, no. GT8, pp. 1095–1111.
Randolph, M.F. (2018) Potential damage to steel pipe piles during installation.
International Press-in Association (IPA), Vol. 3, no. 1, March 2018, pp. 8, Special
Contribution. https://www.press-in.org/en/page/specialcontributions.
Randolph, M.F. and Gourvenec, S. (2011) Offshore geotechnical engineering, Spon
Press, London.
Randolph, M.F. and Houlsby, G.T. (1984) The limiting pressure on a circular pile
loaded laterally in cohesive soil. Géotechnique, Vol. 34, no. 4, pp. 613–623.
Randolph, M.F. and House, A.R. (2002) Analysis of suction caisson capacity in
clay. Proc. Offshore Technology Conference, Houston, Texas, USA, Paper
OTC-14236-MS.
296 References
Randolph, M.F., Martin, C.M. and Hu, Y. (2000) Limiting resistance of a spherical
penetrometer in cohesive material. Géotechnique, Vol. 50, no. 5, pp. 573–582.
Randolph, M.F., O’Neill, M.P. and Stewart, D.P. (1998) Performance of suction
anchors in fine-grained calcareous soils. Proc. Offshore Technology Conference,
Houston, Texas, USA, Paper OTC-8831-MS.
Randolph, M.F. and Wroth, C.P. (1978) Analysis of deformation of vertically loaded
piles. Journal of the Geotechnical Engineering Division, American Society of Civil
Engineers, Vol. 104, no. GT12, pp. 1465–1488.
Reese, L.C. (1997) Analysis of laterally loaded piles in weak rock. ASCE Journal
of Geotechnical and Geoenvironmental Engineering, Vol. 123, no. 11, pp.
1010–1017.
Reese, L.C., Cox, W.R. and Koop, F.D. (1974) Analysis of laterally loaded piles
in sand. Proc. Offshore Technology Conference, Houston, Texas, USA, Paper
OTC-2080-MS.
Romp, R.H. (2013) Installation-effects of suction caissons in non-standard soil con-
ditions, MSc. Thesis, Delft University of Technology, The Netherlands.
Rosenberg, P. and Journeaux, N.L. (1976) Friction and end bearing tests on bed-
rock for high capacity socket design. Canadian Geotechnical Journal, Vol. 13,
pp. 324–333.
Rowe, R.K. and Armitage, H.H. (1987) A design method for drilled piers in soft
rock. Canadian Geotechnical Journal, Vol. 24. pp. 126–142.
Salgado, R., Lyamin, A.V., Sloan, S.W. and Yu, H.S. (2004) Two and three-dimen-
sional bearing capacity of foundations in clay. Géotechnique, Vol. 54, no. 5, pp.
297–306.
Sassi, K., Kuo, M.Y.-H., Versteele, H., Cathie, D. and Zehzouh, S. (2017) Insights
into the mechanisms of anchor chain trench formation. Proc. 8th Int. Conf.
Offshore Site Investigation and Geotechnics, OSIG 2017, Society of Underwater
Technology, London, pp. 963–970.
Sassi, K., Zehzouh, S., Blanc, M., Thorel, L., Cathie, D., Puech, A. and Colliat-
Dangus, J-L. (2018) Effect of seabed trenching on the holding capacity of suction
anchors in soft deepwater Gulf of Guinea clays – Centrifuge testing and numerical
modelling. Proc. Offshore Technology Conference, Houston, Texas, USA, Paper
OTC-28756-MS.
Saviano, A. (2015) Effects of misalignment on the undrained HV capacity of suction
anchors in clay, MSc. Thesis, TU Delft, The Netherlands.
Saviano, A. and Pisanò, F. (2017) Effects of misalignment on the undrained HV
capacity of suction anchors in clay. Ocean Engineering, Vol. 133, pp. 89–106.
Schneider, J.A., White, D.J. and Lehane, B.M. (2007) Shaft friction of piles in sili-
ceous, calcareous and micaceous sands. Proc. 6th Int. Conf. Offshore Site
Investigation and Geotechnics. OSIG 2007, Society of Underwater Technology,
London, pp. 367–382.
Seibold, S. and Berger, W.H. (2010) The seafloor: an introduction to marine geology,
3rd edn, Springer-Verlag, Germany, 272 pp.
Seidel, J.P. and Collingwood, B. (2001) A new socket roughness factor for predic-
tion of rock socket shaft resistance. Canadian Geotechnical Journal, Vol. 38, no.
1, pp. 138–153.
Senders, M. (2005) Tripods with suction caissons as foundations for offshore wind
turbines on sand. Proc. 1st International Symposium on Frontiers in Offshore
Geotechnics, Taylor & Francis, Perth, Australia, pp. 397–405.
References 297
Senders, M. (2008) Suction caissons in sand as tripod foundations for offshore wind
turbines, PhD thesis, School of Civil and Resource Engineering, UWA, Australia.
Senders, M. and Kay, S. (2002) Geotechnical suction pile anchor design in deep
water soft clays. Proc. 7th Annual Conference, Deepwater Risers, Moorings and
Anchorings, IBC, London, UK, 50 pp.
Senders, M. and Randolph, M.F. (2009) CPT-based method for the installation of
suction caissons in sand. ASCE Journal of Geotechnical and Geoenvironmental
Engineering, Vol. 135, no. 1, January, pp. 14–25.
Senders, M., Randolph, M.F. and Gaudin, C. (2007) Theory for the installation of
suction caissons in sand overlain by clay. Proc. 6th Int. Conf. Offshore Site
Investigation and Geotechnics. OSIG 2007, Society of Underwater Technology,
London, pp. 429–438.
Senpere, D. and Auvergne, G.A. (1982) Suction anchor piles – a proven alternative to
driving or drilling. Proc. Offshore Technology Conference, Houston, Texas, USA,
Paper OTC-4206-MS.
Shonberg, A., Harte, I., Aghakouchak, A., Brown, C.S.D, Andrade, M.P. and
Liingaard, M.A. (2017) Suction bucket jackets for offshore wind turbines: applica-
tions from in situ observations. Proc. TC 209 Workshop - 19th ICSMGE, 20
September 2017, Seoul, pp. 65–77.
Siciliano, R.J., Hamilton, J.M., Murff, J.D. and Phillips, R. (1990) Effects of jackup
spudcans on piles. Proc. Offshore Technology Conference, Houston, Texas, USA,
Paper OTC-6467-MS.
Simons, N., Menzies, B. and Matthews, M. (2002) A short course in geotechnical site
investigation, Thomas Telford Publishing, London, ISBN 9780727729484.
Skempton, A.W. (1951) The bearing capacity of clays. Proc. Building Research
Congress, Vol. 1, pp. 180–189.
Smith, E.A.L. (1960) Pile-driving analysis by the wave equation. Journal of the Soil
Mechanics and Foundations Division ASCE, Vol. 86, no. SM4, pp. 35–61.
Smith, I.M. and Chow, Y.K. (1982) Three-dimensional analysis of pile drivability.
Proc. First Int. Conf. Num. Methods in Offshore Piling, ICE, London, UK, 20 pp.
Smith, I.M. and Griffiths, D.V. (1998) Programming the finite element method, 3 rd
edn, Wiley, USA, 534 pp.
Smith, I.M., Hicks, M.A., Kay, S. and Cuckson, J. (1988) Undrained and partially
drained behaviour of end bearing piles and bells founded in untreated calcarenite.
Proc. 1st International Conference on Calcareous Sediments, 15–18 March 1988,
Perth, Vol. 2, A.A. Balkema, Rotterdam, pp. 663–679.
SOLCYP (2017) Design of piles under cyclic loading – SOLCYP recommendations,
Puech, A. and Garnier, J. (Ed.), ISTE – WILEY, France.
Sorensen, S.P.H., Ibsen, L.B., Augustesen, A.H. (2010) Effects of diameter on initial
stiffness of p-y curves for large-diameter piles in sand, in: numerical Methods in
Geotechnical Engineering. CRC Press, pp. 907–912.
Sørlie, E. (2013) Bearing capacity failure envelope of suction caissons subjected to
combined loading, MSc. Thesis, Norwegian University of Science and Technology,
Trondheim.
Stavropolou, E., Dano, C. and Boulon, M. (2019) Experimental investigation of
the mechanical behaviour of soft carbonate rock/grout interfaces for the design
of offshorewind turbines. Proc. E3S Web of Conferences 7th Int. Symp.
Deformation Characteristics of Geomaterials, Vol. 92, IS-Glasgow. doi:10.1051/
e3sconf/20199213006.
298 References
Tjelta, T.I. (1995) Geotechnical experience from the installation of the europipe
jacket with bucket foundations. Proc. Offshore Technology Conference, Houston,
Texas, USA, Paper OTC-7795-MS.
Tjelta, T.I. (2015) The suction foundation technology. Proc. 3rd International
Symposium on Frontiers in Offshore Geotechnics, Taylor & Francis, Oslo,
Norway, pp. 85–93.
Tran, M.N. (2005) Installation of suction caissons in dense sand and the influence of
silt and cemented layers, PhD. Thesis, Department of Civil Engineering, University
of Sydney, Australia.
Tran, M.N., Randolph, M.F. and Airey, D.W. (2004) Experimental study of suction
installation of caissons in dense sand. Proc. 23rd International Conference on
Offshore Mechanics and Arctic Engineering, 20–25 June 2004, Vancouver,
Canada, Paper OMAE2004-51076.
Tran, M.N., Randolph, M.F. and Airey, D.W. (2007) Installation of suction caissons
in sand with silt layesr. ASCE Journal of Geotechnical and Geoenvironmental
Engineering, Vol. 133, no. 10, October 1, pp. 1183–1191.
Van den Berg, P. (1994) Analysis of soil penetration, PhD. Thesis, TU Delft, The
Netherlands.
Van Dijk, B.F.J. (2015) Caisson capacity in undrained soil: failure envelopes with
internal scooping. Proc. 3rd International Symposium on Frontiers in Offshore
Geotechnics, Taylor & Francis, Oslo, Norway, pp. 337–342.
Versteele, H., Jaeck, C., Silvano, R. and Bignold, D. (2015) Sinkage analysis of
A-shaped jack-up mudmats using the coupled Eulerian Lagrangian approach.
Proc. 3rd International Symposium on Frontiers in Offshore Geotechnics, Taylor
& Francis, Oslo, Norway, pp. 1323–1328.
Versteele, H., Kuo, M.Y.-H., Cathie, D., Sassi, K. and Zehzouh, S. (2017) Anchor
chain trenching – numerical simulation of progressive erosion. Proc. 8th Int. Conf.
Offshore Site Investigation and Geotechnics, OSIG 2017, Society of Underwater
Technology, London, pp. 971–977.
Vicente, V., López-Gutiérrez, J-S., Dolores Esteban, M., Alberdi, P., Imaz, M. and
Serraclara, J-M. (2017) Monopiles in offshore wind: preliminary estimate of main
dimensions. Ocean Engineering, Vol. 133, pp. 253–261.
Villalobos, F.A. (2007) Installation of suction caissons in sand. Proceedings 6th
Chilean Conf. on Geotechnics (Congreso Chileno de Geotecnia), Chile.
Vivatrat, V., Valent, P.J. and Ponterio, A.A. (1982) The influence of chain friction on
anchor pile design. Proc. Offshore Technology Conference, Houston, Texas, USA,
Paper OTC-4178-MS.
Vulpe, C., Gourvenec, S. & Power, M. (2014) A generalized failure envelope for
undrained capacity of circular shallow foundations under general loading.
Géotechnique Letters, Vol. 4, no. July–September, pp. 187–196.
Waltham, T. (2009) Foundations of Engineering Geology, 3rd edn, CRC Press,
London. doi:10.1201/9781315273488.
Wang, M.C., Demars, K.R. and Nacci, V.A. (1977) Breakout capacity of model suc-
tion anchors in soil. Canadian Geotechnical Journal, Vol. 14, no. 2, pp. 246–257.
Wang, M.C., Demars, K.R. and Nacci, V.A. (1978) Applications of suction anchors
in offshore technology. Proc. Offshore Technology Conference, Houston, Texas,
USA, Paper OTC-3203-MS.
Wang, M.C., Nacci, V.A. and Demars, K.R. (1975) Behavior of underwater suction
anchor in soil. Ocean Engineering, Vol. 3, no. 1, pp. 47–62.
300 References
Watson, P.G., Gaudin, C., Senders, M. and Randolph, M.F. (2006) Installation of
suction caissons in layered soil. Proc. 6th International Conference on Physical
Modeling in Geotechnics, Hong Kong, pp. 685–692.
Whitehouse, R.J.S., Harris, J.M., Sutherland, J. and Rees, J. (2011) The nature of
scour development and scour protection at offshore windfarm foundations.
Marine Pollution Bulletin, Vol. 62, no. 1, pp. 73–88.
Williams, A.F. and Pells, P.J.N. (1981) Side resistance rock sockets in sandstone,
mudstone, and shale. Canadian Geotechnical Journal, Vol. 18, no. 4, pp.
502–513.
Wind Europe (2018) Offshore wind in Europe - Key trends and statistics 2017,
February, 37 pp.
Wind Europe (2020) Offshore wind in Europe - Key trends and statistics 2019,
https://windeurope.org/data-and-analysis/product/offshore-wind-in-europe-
key-trends-and-statistics-2019/.
XG Geotools (2014) SPCalc Manual.
Xie, Y. (2009) Centrifuge model study on spudcan-pile interaction, PhD. Thesis,
Department of Civil and Environmental Engineering, National University of
Singapore, Singapore.
Xie, Y, Leung, C.F. and Chow, Y.K. (2017) Centrifuge modelling of spudcan-pile
interaction in soft clay overlying sand. Géotechnique, Vol. 67, no. 1, pp. 69–77.
Zienkiewicz, O.C., Chang, C.T. and Bettess, P. (1980) Drained, undrained, consoli-
dating and dynamic behaviour assumptions in soils. Géotechnique, Vol. 30, no. 4,
pp. 385–395.
Index
Page numbers in italics refer to figures and those in bold refer to tables.
A in-place response, 7, 8
internal stiffeners, 21, 42–44, 119
accuracy of prediction, 4–5, 210, 258 loads and load transfer, 25, 26, 27,
active suction, 23 28
adhesion factors, 95, 104 optimum lug level (OLL), 60, 216,
outer for torsional loading, 196–197 217
pile-rock, 225, 226, 228, 229, 229 partial factor design approach, 72,
pile-soil, 182, 185 76–77, 78
adjacent infrastructure, 256–261 riser tower anchors (RTA), 41–42;
alpha method, 120 see also suction foundations
anchor chain, 190, 221 anchor holding capacity, 41, 269–270,
boulder contact, 130–131 271
trenching, 48–49, 50, 178, 209, 210, API RP 2A, 74, 76, 165, 228, 248
212, 269–274 API RP 2GEO, 9, 56, 74, 76, 103, 106
failure model, 17, 272 API RP 2SK, 9, 56, 76
anchor foundations, 1, 2, 41–46 API RP 2T, 76
anchor piles, 42–44, 239 axial failure models, 15–17, 19,
geometry, 118–119, 166 176–182, 184, 192–194
groups, 238 axial friction, 61
in-place resistance, 212–213 axial pull-out load, 178, 179
installation hazards, 44–48 axial resistance, 187–194, 209–212
misalignment, 45, 45 drilled and grouted piles in rock,
over-penetration, 87–88 226–232
pull-out capacity, 59–60 driven piles in weak rock, 224–226
tilt, 45, 220–223 tensile, 262
buoy anchors, 42 axial shear stress, 249
design considerations, 7, 59, 60 axial soil resistance, 91
external stiffeners, 21
factors of safety, 72, 78
FPSO anchors, 42, 220–223 B
in-place resistance, 7, 8, 176–179,
180 back analysis, 95, 143, 248, 266
analysis, 198, 202, 209, 212–213 base failure, 23, 111–112, 117–120
caisson anchor, 65–66, 67 base moment, 249, 250
FPSO anchor design example, bearing capacity factor, 182, 184–186,
220–223 191
resistance reduction, 209–211, bearing capacity method, 90–91,
212 103–105, 114–120
301
302 Index
T
V
target penetration depth (TPD), 85,
87–88, 114, 268 vacuum anchors, 263
tensile V loads, 176–179, 181–182, 184 vent design, 263–265
Index 315