THESIS DOCTOR OF PHILOSOPHY GUHA Indranil 2021 PDF
THESIS DOCTOR OF PHILOSOPHY GUHA Indranil 2021 PDF
THESIS DOCTOR OF PHILOSOPHY GUHA Indranil 2021 PDF
By
INDRANIL GUHA
BE (Mech), M. Tech (Pipeline), MIEAust, CPEng, NER
August 2020
Thesis declaration
I, Indranil Guha, certify that:
This thesis has been substantially accomplished during enrolment in this degree.
This thesis does not contain material which has been submitted for the award of any other
degree or diploma in my name, in any university or other tertiary institution.
In the future, no part of this thesis will be used in a submission in my name, for any other
degree or diploma in any university or other tertiary institution without the prior approval
of The University of Western Australia and where applicable, any partner institution
responsible for the joint-award of this degree.
This thesis does not contain any material previously published or written by another
person, except where due reference has been made in the text and, where relevant, in the
Authorship Declaration that follows.
This thesis does not violate or infringe any copyright, trademark, patent, or other rights
whatsoever of any person.
This thesis contains published work and/or work prepared for publication, some of which
has been co-authored.
Signature:
Date: 13/08/2020
i
Abstract
Pipelines are the safest, most reliable, economic and efficient means for the
transportation of petroleum fluids and water. In the last few decades, the importance of
the pipeline transportation system has increased significantly due to the hydrocarbon
industry developing resources that are further offshore. This thesis is concerned with the
structural behaviour of submarine pipelines subjected to submarine slide impact, and
thermal loading conditions.
The research has aimed to support the transition of oil and gas developments into
deeper water and more remote conditions. The principal motivations are the needs for (a)
export and tieback pipelines to negotiate regions of unstable or steeply sloping terrain,
where submarine slides may occur potentially impacting the pipeline; and (b) for high
temperature, high pressure, pipelines laid directly on the seabed in deep water to
withstand cycles of thermal and pressure loading during operation. Emphasis is placed
on the axial pipe-soil interaction and structural behaviour of the pipeline.
Analytical models are derived for the axial pipe-soil elastic stiffness and
numerical solutions using FE software ABAQUS are provided for axial, horizontal, and
vertical motions of a pipeline relative to the seabed with the aim of expressing these
parameters in terms of fundamental elastic properties of the soil. Where appropriate, the
theoretical techniques used for pile design are transferred to pipeline conditions to model
the axial slide loading condition. A parametric study of pipeline-slide interaction has
been carried out to provide insights into the dominant governing parameters and to assess
the sensitivity of the pipeline structural loading to the geotechnical (i.e. pipe-soil, slide
material) input parameters.
Thereafter, classical buckling theory has been extended to estimate the critical
lateral buckling load of on-bottom pipelines subjected to axial loading down-stream of
the submarine slide, incorporating the dominant governing geotechnical and as-laid
parameters, and to assess the sensitivity of the critical buckling load of the pipeline to
these input parameters.
Axial walking of deepwater pipelines due to thermal and pressure cycles during
their operational life is also addressed, extending the present theoretical framework by
incorporating the elastic-plastic (bi-linear) response of the soil into the analytical solution
and verifying the solution numerically. Furthermore, the walking behaviour of on-bottom
pipeline has been analysed using a velocity-dependent friction model within ABAQUS
to provide an equivalent friction factor that allows for this velocity-dependency via a
single value.
ii
Acknowledgement
I would like to thank my supervisors Prof David J White and Prof Mark F
Randolph for their guidance and support throughout the course of this research. I am
greatly indebted to both of them for their penetrating and timely criticisms. Without that,
many of the ideas of this thesis may have remained undeveloped. The research presented
here forms part of the activities of the Centre for Offshore Foundation Systems (COFS),
supported as a node of the Australian Research Council Centre of Excellence for
Geotechnical Science and Engineering (grant CE110001009) and through the Fugro
Chair in Geotechnics, the Lloyd’s Register Foundation Chair and Centre of Excellence
in Offshore Foundations and the Shell EMI Chair in Offshore Engineering.
I also owe much to my colleagues, staffs from IT and administration in the Centre
for Offshore Foundation System at the University of Western Australia for their constant
support and motivation.
I would also like to record my thanks to my wife Munmun Basak for her
continuing encouragement, love, sacrifice and support throughout the journey. She made
this journey easier for me. At the same time, I would like to thank my parents and in-
laws. Without their motivation, I could not have finished it.
Except where specific reference is made in the text to the work of others, the
contents of this dissertation are original and have not been submitted to any other
university.
iii
Authorship declaration: Co-authored publications
1. Guha, I., & Whitney, B. (2012). ‘Seismic vulnerability of Australian buried pipeline
industry’. In Proc of 11th Australia New Zealand Conference on Geomechanics (ANZ
2012), Melbourne, Australia.
This paper was based on the research proposal and chapters 1 and 2. Candidate
drafted the paper and second author Dr. Beau Whitney contributed to the final version.
2. Guha, I., Randolph, M.F., White, D.J. (2016). ‘Evaluation of elastic stiffness
parameters for pipeline-soil interaction', ASCE, Journal of Geotechnical and
Geoenvironmental Engineering, 142, 6.
DOI: 10.1061/(ASCE)GT.1943-5606.0001466
This paper was based on chapter 3. Candidate performed the analytical and numerical
modelling under the guidance of the co-authors. The candidate prepared the first draft
of the paper. Prof. Mark Randolph and Prof. David White reviewed and contributed
to the final version.
3. Guha, I., Randolph, M.F., & White, D.J., (2020). Analysis of axial response of
submarine pipeline to debris flow loading’. Accepted at Journal of Geotechnical and
Geoenvironmental Engineering, August 2020.
This paper was based on chapter 4. Candidate performed the analytical and numerical
modelling under the guidance of the co-authors. The candidate interpreted the results
reported in this paper. Prof. Mark Randolph prepared the first draft of the paper and
all authors revised and contributed to the final version.”
4. Guha, I., White, D.J., & Randolph, M.F (2020). ‘Parametric solution of lateral
buckling of submarine pipelines. Applied Ocean Research, 98.
https://doi.org/10.1016/j.apor.2020.102077
This paper was based on chapter 5. Candidate performed the analytical and numerical
modelling under the guidance of co-authors. The candidate prepared the first draft of
the paper. Prof. David White and Prof. Mark Randolph reviewed and contributed to
the final version.
5. Guha, I., White, D.J., & Randolph, M.F. (2018). ‘Subsea pipeline walking – the effect
of a bi-linear friction model’. Submitted at, ISFOG 2020.
This paper was based on chapter 6 and was accepted at OMAE 2019, but was
withdrawn (as none of the authors could attend) and resubmitted to ISFOG 2020.
Candidate performed the analytical and numerical modelling under the guidance of
co-authors. The candidate prepared the first draft of the paper. Prof. David White and
Prof. Mark Randolph reviewed and contributed to the final version.
iv
6. Guha, I., White, D.J., & Randolph, M.F. (2018). ‘Subsea pipeline walking with
velocity dependent friction’, Applied Ocean Research, 82, January 2019, Pages 296-
308.
https://doi.org/10.1016/j.apor.2018.10.028
This paper was based on chapter 7. Candidate performed the analytical and numerical
modelling under the guidance of co-authors. The candidate prepared the first draft of
the paper. Prof. David White and Prof. Mark Randolph reviewed and contributed to
the final version.
v
Structural analysis of submarine pipelines
Table of Content
INTRODUCTION ............................................................................... 1
OVERVIEW ............................................................................................................ 1
SUBMARINE PIPELINES ..................................................................................... 2
AUSTRALIAN OIL AND GAS DEVELOPMENT AND SUBMARINE
PIPELINES ............................................................................................................. 4
1.1 GEOTECHNICAL CHALLENGES .............................................................. 6
1.2 RESEARCH GOALS ...................................................................................... 7
1.3 THESIS OUTLINE ......................................................................................... 9
vii
Structural analysis of submarine pipelines
viii
Structural analysis of submarine pipelines
x
Structural analysis of submarine pipelines
INTRODUCTION
OVERVIEW
Pipelines are the most safe, reliable, economic and efficient means for the transportation
of petroleum fluids and water. In the last few decades the importance of pipeline
transportation systems has increased significantly due to their extensive use in the
developments of oil and gas infrastructure are moving towards deeper water. This
requires design and construction of long high temperature and high pressure pipelines
from deep sea to the shore. As a result, the cost of in-field flowlines and export pipelines
has increased significantly, leading to a paradigm shift in the importance of pipelines and
The first submarine pipeline was laid on the seabed in 1954. Since then offshore
production has reached water depths in excess of 2,100 m, and current exploration
has extended to almost 3,500 m depth (Kyriakides & Corona, 2007). Thousands
of kilometres of pipelines have already been laid onshore to meet the demand
supply gap for energy, and to transport water or sewage to mai ntain human
habitation throughout the land. Therefore, the onshore pipeline industry has
needed to meet challenges such as where pipelines cross earthquake prone zones,
traverse steep hill sides or in challenging regions such as the arctic. These threats
etc.
1
Structural analysis of submarine pipelines
However, the offshore industry is still evolving and with passing time new
SUBMARINE PIPELINES
investment, most of the pipelines have a life span of more than 40 years and require
limited maintenance when designed and constructed effectively. In the past three decades,
significant hydrocarbon reserves have been identified offshore, including the North Sea,
the Gulf of Mexico, the Persian Gulf, offshore Brazil, West Africa, and more recently
fields in Malaysia, Indonesia, Northwest Australia, and in the Bay of Bengal of India (Bai,
2001).
As current hydrocarbon reserves are found in much deeper water, many hundreds of
kilometres of new submarine pipelines are being proposed to transport hydrocarbon fluids
than typically encountered in shallow water. Also, hydrodynamic forces are much
reduced in deepwater, so that pipelines are generally laid directly on the seabed without
such as cyclic expansion and contraction due to fluctuating temperatures and pressure. A
thorough introduction to these issues has been provided in relevant chapters of the
offshore geotechnical engineering books of Dean (2010) and Randolph and Gourvenec
(2011), and also in the papers of Cathie et al. (2005), White and Cathie (2010), Randolph
(2013), and the recently concluded JIP paper of White et al. (2016).
2
Structural analysis of submarine pipelines
Figure 1-1: Existing and proposed pipelines at the North West Shelf of Australia
generally small and the dominant forces are from high internal temperature and pressure,
which tend to cause expansion (Bruton et al., 2008). Axial resistance between the pipe
and the seabed opposes this expansion. Excessive compressive forces lead to buckling,
but the buckling response depends critically on the lateral soil resistance. When buckling
occurs, it significantly reduces the net axial load in the pipe. On the other hand, excessive
buckling may lead to high bending strains in the pipe section. So, controlled buckling
(Figure 1-2) may be a feasible solution for relief of thermal loading. Accumulated axial
movement due to repeated thermal and pressure cycles may lead to global displacement
of pipes. This phenomenon is termed ‘walking’ (Carr et al., 2006). For design purposes,
it is very important to assess pipeline buckling and walking accurately. Recent design
approaches to control buckling and walking have necessitated predicting the available
soil resistance on pipelines undergoing movement, accounting for the associated changes
in seabed geometry and strength. The existing models are mainly derived for stability
3
Structural analysis of submarine pipelines
analyses. The challenge is to extend existing models to account for geometry changes,
Figure 1-2 Controlled lateral buckling of an on-bottom pipeline (Jayson et al., 2008)
Australia has proved to be one of world’s leading Liquefied Natural Gas (LNG) exporter
countries with estimated reserve of 170 trillion cubic feet (tcf) of gas reserve. The
offshore industry in North-West shelf of Australia ( Figure 1-1 ) is maturing rapidly due
to the huge gas reserves in this region. Much of Australia’s remote sub-sea oil and gas
resources remain undeveloped. Over the last two decades, offshore oil and gas
developments have gradually expanded infrastructure to deeper water further from shore.
This has led to a paradigm shift of focus from fixed platforms to floating production
systems, which in turn has resulted in increasing importance of offshore pipelines and
4
Structural analysis of submarine pipelines
riser systems (Randolph & White, 2008a). Australia’s gas industry relies on ultra-long
sea-bed pipelines to bring the oil and gas from remote offshore hydrocarbon fields to
The first submarine pipeline in Australian water was laid in 1983 between Burrup
Peninsula and North Rankin A platform, known as 1TL (White et al., 2013). Since then
the design and construction of submarine pipelines in Australian water has increased
under consideration offshore Western Australia (modified after Cassidy & Tian, 2012).
The gas platforms offshore Australia are connected to onshore liquefaction plants
(Liquefied Natural Gas export terminals) and also to the domestic natural gas network.
As discussed previously, constructing a pipeline network link from an offshore oil and
gas field to the onshore facility is a huge capital investment. For example, in Australia
5
Structural analysis of submarine pipelines
the construction cost of offshore pipelines is estimated to exceed $4.5 million per km
(Cassidy & Tian, 2012). That was based on the current projects at the time, Gorgon in
water depth of 1350 m, lengths 65 and 140 km; Scarborough at depth of 900 m and length
280 km; Pluto at depth of 830 m and length of 180 km; and Browse at a depth of 600 m
and lengths of 5, 24 and 400 km. With over 2000 km of offshore pipeline in design and
construction phases the estimated industry volume is expected to exceed $10 billion.
Following the other regions of the world, Australian oil and gas development is also
moving beyond the immediate continental shelf into deeper water. The scope of the
pipelines only.
From a geotechnical perspective, the transition from shallow water to deepwater pipelines
and associated infrastructure increases the prevalence of certain design challenges (see
Figure 1-3):
• Pipeline networks that are subjected to constant high pressure and high
• Export pipelines to shore, crossing steep and unstable terrain at the continental
conditions are the two major aims. By contrast in deepwater along with stability and
serviceability the design may also allow movement of the pipelines, for example
changes. Also, the response of engineered buckles is affected by the axial pipe-soil
6
Structural analysis of submarine pipelines
(a)
(b)
Figure 1-3: Some geotechnical aspects of pipeline design (a) in shallow water & (b) in deepwater (White
& Cathie, 2010)
The proposed research is concerned with the structural behaviour of submarine pipelines
subjected to impact by submarine slides, and cyclic thermal loading conditions. The
research aims to support the transition of oil and gas developments into deeper water and
more remote conditions. The principal motivation is the need for export and tieback
may occur (i.e. submarine slides), and for these pipelines to withstand other forms of
1. To improve the techniques for assessing the axial pipe-soil interaction forces
resulting from relative pipe-soil movement, including the passage of mobile slide
7
Structural analysis of submarine pipelines
for the axial pipe-soil elastic stiffness and numerical solutions will be provided.
parameters in horizontal and vertical directions. Parallels will be drawn with the
‘t-z’ techniques for assessing pile-soil interaction forces. Where appropriate, the
conditions.
interaction (and also other pipeline loading conditions), to provide insights into
the dominant governing parameters and to assess the sensitivity of the pipeline
assessment of the potential effect of slide loading (and other loading conditions)
when the pipe is subjected to axial loading down-stream of the slide, to provide
insights into the dominant governing geotechnical and as-laid parameters and to
assess the sensitivity of the critical buckling load of the pipeline to the
8
Structural analysis of submarine pipelines
velocity.
The thesis consists of eight chapters. Every chapter deals with a specific pipeline
soil interaction issue of submarine infrastructure. A brief outline of each chapter is given
below.
methods used for structural analysis of pipelines in the chosen field of the thesis. A
thorough literature review is carried out separately for each chapter as every chapter of
this thesis deals with a different pipeline structural issue. Methodologies applied in each
interaction
This chapter focuses on elastic stiffness parameters for axial, horizontal, and
vertical motions of a pipeline relative to the seabed, with the aim of expressing these
exists in the literature on the axial elastic response of on-bottom pipelines, particularly
for axial stiffness, focusing on the case of shear modulus proportional to depth. The
solution was then verified through numerical analysis. Further numerical analysis was
carried out to obtain relationships for horizontal and vertical elastic stiffnesses of on-
displacement modes were developed. Here recommendations are made for the selection
9
Structural analysis of submarine pipelines
application to the analysis of pipeline laying, buckling, walking, and on-bottom stability.
straight on-bottom pipeline. It is shown that the non-dimensional axial loads and axial
displacements depend on three non-dimensional input parameters, i.e. the driving force
in the slide zone, seabed resisting force in the passive zone, and pipe-soil stiffness. Non-
dimensional design charts are presented to show the effect of individual input parameters
on axial loads and axial displacements. The maximum axial load in the pipe is directly
proportional to the slide force, while the load at the transition from elastic to plastic soil
resistance is initially proportional to the slide force but then becomes limited. The limit
is reached for most relevant values of the slide force. Results from numerical FE analysis
to verify the analytical model are also presented, showing close agreement between
analytical and numerical solutions. Buckling was ignored in the analytical model.
However, the existing classical theory of buckling was linked to the output of the
analytical model to show the vulnerability of the pipelines towards buckling in case of
various slide loading conditions. On bottom submarine pipelines are more susceptible to
in the offshore industry due to the complexities involved in the analysis, and the potential
researchers and recent joint industry projects provide a basis for estimation of the critical
buckling load of a straight, or in some cases imperfect, pipe on either a rigid or elastic
10
Structural analysis of submarine pipelines
seabed. However, systematic solutions for the combined effects of nonlinear soil
properties and the as-laid geometry – specifically the out-of-straightness – on the buckle
initiation behaviour have not been developed previously. This chapter reports an
subjected to axial compressive loading. The seabed was modelled with lateral and axial
elastic, perfectly plastic, springs to idealise the load-displacement behaviour of the soil
and the pipe was modelled with pipe elements. Buckling was performed by a
displacement controlled finite element method with the modified RIKS algorithm that is
available in the commercial software ABAQUS. This numerical tool was used to develop
a parametric solution for the present problem in terms of the various pipe material and
geometry parameters and the lateral and axial pipe-soil interaction parameters. In
particular, the influence of the magnitude and stiffness of the lateral pipe-soil response
was investigated, highlighting the sensitivity of the pipeline response to the geotechnical
inputs. The results have been synthesised in a generic non-dimensionalised design chart
to estimate the buckling load, valid for the range of inputs covered by the parametric
study.
The objective of this chapter is to explore the gaps in the present analysis methods
proposed by various joint industry projects and others for pipeline walking behaviour.
Thereafter, existing analytical methods are extended to bridge the gap between analytical
modelling tools and numerical analysis. A literature review shows that previous
researchers have identified four major conditions that lead to pipeline walking. These are
seabed slope, Steel Catenary Risers (SCR) tension, thermal transients and the unequal
distribution of load from the product, due to separate of gas and liquid on shutdown. In
11
Structural analysis of submarine pipelines
expressed as ultimate resistance per unit length. Often this term is expressed non-
pipe weight, W. However, the elastic-plastic (i.e. bi-linear) behaviour of the soil and the
effect of this bi-linear response of the soil on the walking behaviour were poorly
additional parameter, specified as the mobilisation displacement, i.e. the amount of axial
displacement that occurs before the ultimate friction is generated, with the resistance
rising linearly with displacement up to this value. The walking behaviour is affected by
the axial friction mobilisation displacement. Numerical results were provided by previous
researchers. However, the analytical solution for the reduction of walking per cycle due
to increase in mobilisation displacement has not been attempted. This chapter will give
insight to the walking behaviour due to seabed slope, SCR tension and thermal transients.
The existing analytical solutions are extended to incorporate the elastic-plastic response
of the soil into the expression of pipeline walking and a new derivation proposed.
Accumulated axial movement due to repeated thermal cycles may lead to global
with seabed friction. In conventional analyses, seabed friction is independent of the rate
of thermal expansion but it has been recognised that the sliding resistance between a pipe
and the seabed varies with velocity, partly due to drainage effects. In this paper a
validated via single element and simple (flat seabed) pipeline cases. This model features
upper and lower friction limits, with a transition that occurs as an exponential function of
12
Structural analysis of submarine pipelines
velocity. A parametric study is performed using differing rates of heating and cool-down
in walking situations driven by seabed slope, SCR end tension and the difference between
heat up and cool down rates. The walking behaviour is compared to cases with constant
friction and solutions are proposed to express the velocity-dependent response in terms
of an equivalent constant friction. These equivalent friction values can then be applied in
existing simple solutions or more complex numerical analyses, as a short cut method to
Chapter 8: Conclusions
The findings from each chapter are summarised in this final chapter. The original
landslide and thermal loading is discussed here. In addition to the concluding remarks,
13
Structural analysis of submarine pipelines
References
during lateral buckling and pipeline walking – the Safebuck JIP. Proc. of Offshore
Carr M, Sinclair F, & Bruton D.(2006). Pipeline walking — understanding the field
layout challenges, and analytical solutions developed for the SAFEBUCK JIP.
Cathie D. N., Jaeck C, Ballard J. C. & Wintgens J-F. (2005). Pipeline geotechnics – state-
of-the-art. Proc. 2nd Int. Symp. On Frontiers in Offshore Geotechnics, Perth, 95-
114.
Cassidy, M.J., & Tian, Y. (2012). Development and application of models for the stability
78.
Jayson D, Delaporte P, Albert J-P, Prevost M. E., Bruton D & Sinclair F. (2008). Greater
Plutonio Project – Subsea Flowline Design and Performance. Proc. 31st Offshore
14
Structural analysis of submarine pipelines
London, 27-59.
White, D. J., & Cathie, D. N. (2010). Geotechnics for subsea pipelines. Proc. 2nd Int.
White, D. J., & Boylan, N. & Levy, N.H. (2013). Geotechnics offshore Australia -
White, D.J., Randolph, M.F., Gaudin, C., Boylan, N., Wang, D., Boukpeti, N., Zhu, H &
10.4043/27034-MS.
15
Structural analysis of submarine pipelines
16
Structural analysis of submarine pipelines
2.1 OVERVIEW
This thesis is concerned with the structural analysis of sub marine pipelines
stiffness
a production facility and they are relatively short in length and of moderate
gas or water for injecting into the wells, are called service lines. To transport
hydrocarbon, water etc. between offshore facilitates within a limited area infield
pipelines are used. Larger diameter and longer export pipelines, transmission
17
Structural analysis of submarine pipelines
Size wise the diameter of submarine pipelines can vary between 76 mm (3) to
the materials used for submarine pipelines are mainly carbon steel a s per
American Petroleum Industry standards Grade B to Grade X70, and even higher
Trunkline
Platform with
compression
facilities
To Domestic
market/Power Plants
Wellhead
Manifold
Figure 2-1: Schematics of different kind of pipelines in an offshore field (after Bai 2001).
Figure 2-1 shows a schematic of a typical offshore field where different pipelines
export pipelines fall under midstream, transporting hydrocarbons from the offshore
18
Structural analysis of submarine pipelines
matured significantly over the last decade, due to the extensive use of pipelines in the
challenged by extreme events, for example the geohazard of a submarine landslide (from
here on referred to as submarine slides) while crossing continental slopes from deepwater
towards the shore. Structural integrity issues also arise from the high temperature (HT)
and high pressure (HP) of the transported fluid, and the cyclic changes resulting from
operational shut down and start up cycles. Pipeline impact from a submarine slide causes
high bending and tensile stresses (and strains) in the pipe wall, leading to risk of fracture,
while HT/HP results in high compressive axial force in the pipe leading to buckling and,
This chapter summarises current methods used to analyse the effects of submarine slides,
lateral buckling and walking mechanisms, and also relevant research on pipe-soil
interaction, which is a necessary prerequisite for analysing the pipeline system as a whole.
categories - classical plasticity theory, finite element analysis and physical modelling
studies. McCarron (2011) summarised the design considerations for subsea flowlines
against lateral and upheaval buckling. Simple numerical modelling techniques of pipe-
soil interactions are also highlighted in the book. Randolph & Gourvenec (2011) also
presented the key aspects of offshore pipeline design and current practices in their chapter
interrelations have been provided by Cathie et al. (2005) and White and Cathie (2010).
19
Structural analysis of submarine pipelines
Here, a brief overview of some of the analytical methods are given that have contributed
to design approaches.
concerned primarily with the issues associated with lateral buckling, and walking. These
two issues have been the topic of a longstanding joint industry project, the SAFEBUCK
JIP (Bruton et al. 2006, 2008, Collberg et al. 2011). The axial and lateral resistances
offered by the shallow sediments are key inputs to the design and depend on the
embedment of the pipeline, and on the velocity and time of the movement relative to the
researched extensively in the past two decades by various researchers. Plasticity solutions
may be used for simple static penetration (Randolph & White 2008b, Martin & White
2012). After laying of the pipeline it undergoes hydro testing before commencing
operations, where it will be subjected to cycles of high pressure and high temperatures.
The axial expansion of a pipeline due to thermal loading mobilises axial resistance similar
to that of shaft resistance on a vertically-loaded pile (Randolph & White, 2008a). Even
though at face value the axial resistance would seem to be trivial (a sliding failure with
known vertical load), it turns out to be more complex in practice due to sensitivity to the
curved surface of the pipe and also the degree of consolidation before and during axial
movement.
under vertical and horizontal load, comparing the limiting loads with those calculated
using upper bound plasticity, and proposing yield envelopes covering different
combinations of vertical and horizontal load. Pipeline penetration and response under
20
Structural analysis of submarine pipelines
& Yu (2007) and Yu & Konuk (2007) studied the large displacement pipe-soil interaction
problem. Merifield et al. (2009) also studied the vertical penetration response of pipes
and subsequent horizontal resistance for pushed-in-place (PIP) pipes. Wang et al. (2010)
and Chatterjee et al. (2012) studied large amplitude lateral motion of a pipe using a RITSS
approach (Hu & Randolph, 1998), but implemented for the first time in finite element
software ABAQUS (see Figure 2-2). Such large deformation, two and three dimensional,
finite element analyses are relatively difficult to perform and as such are not economic
for design from a project schedule perspective. Martin et al. (2013) provided a more
connected to a 3D pipe model, but even that approach would be outside the capability of
most projects.
Figure 2-2: Large deformation finite element analysis of pipeline penetration into highly layered material
(Chatterjee et al. 2012).
Pipeline resistance during large amplitude lateral movements has been investigated
by physical modelling (White & Dingle 2011). Embedment of the pipeline occurs due to
the submerged weight of the pipe and during additional cyclic motions the lay process
21
Structural analysis of submarine pipelines
(Westgate et al., 2009, 2010). Cathie et al. (2005) summarised various models proposed
by Wagner et al. (1987), Lieng et al. (1988), Verley & Sotberg (1992) and Verley & Lund
(1995) used for assessing lateral resistance of partially embedded pipelines. Large and
small-scale tests were reported by Bruton et al. (2006) to provide the key parameters
journals, conference proceedings and technical notes, mostly linked with three major
pipeline design JIPs such as HotPipe ( Collberg et al. 2005), SMARTPIPE ® (White et al
2010) and SAFEBUCK (Collberg et al. 2011) were completed over the last two decades.
Among them the SAFEBUCK JIP (Collberg et al. 2011) was completed most recently
and aimed at pipeline-seabed interaction specifically, and some of its work has been
small and the dominant forces are from high internal temperature and pressure, which
tend to expand the pipes, increasing the net axial force (Bruton et al., 2008). Axial
resistance between the pipe and the seabed opposes this expansion. Excessive
compressive forces lead to buckling, but the buckling response depends critically on the
lateral soil resistance. When buckling occurs, it significantly reduces the axial loading.
On the other hand, excessive buckling may lead to high bending strains in the pipe section.
So, a controlled buckling may be a feasible solution for relief of thermal loading. If the
soil resistance is too high, there will be an accumulation of axial loading in the pipe.
Cumulative axial movement due to repeated thermal cycles may lead to global
displacement of pipes. This phenomenon is called ‘walking’ (Carr et al., 2006). For
design purposes, it is very important to assess pipeline buckling and walking accurately.
22
Structural analysis of submarine pipelines
Recent design approaches to control buckling and walking have necessitated predicting
the soil forces on mobile pipelines (White & Cheuk, 2008), accounting for the associated
changes in seabed geometry and strength. The existing models are mainly derived for
stability analyses. There is a need to extend existing models to account for geometry
changes, remoulding and reconsolidation effects that influence large amplitude cyclic
displacements.
It is increasingly recognised that in the axial direction, the pipe-soil ‘t-z’ interaction
response is highly sensitive to the rate of movement (which affects the drainage
conditions) and, for cyclic loading (from thermal events, for example), the pause period
between events is important (Randolph & White 2008). This means that a wide range of
axial friction factors ( typically ranging from 0.1 to unity), and corresponding maximum
load transfer forces, can apply in different conditions. There is still limited understanding
of how strongly and in what ways the adopted axial pipe-soil parameters affect the
The passive interaction of the seabed and the pipe has been studied by many groups
and there are many reports and publications available on this topic. A summary of current
research and practice in this area is given by publications related to SAFEBUCK JIP
(Collberg et al. 2011) and White et al. (2017). In conventional design approach the pipe
seabed interaction model is idealized by spring-slider systems distributed along the length
of the pipeline. The ultimate axial resistance per unit length, F, may be expressed in terms
factor, ζ, to account for wedging around the curved surface of the pipe: (Randolph and
White 2008b).
23
Structural analysis of submarine pipelines
𝐹 = 𝜇𝑁 = 𝜇𝜁𝑊 (2.3)
This is a simple formula but resolving the parameters and is very challenging. An
alternative approach used to estimate the axial resistance is called the total stress (alpha)
method, and is comparable to the equivalent technique to estimate axial pipe capacity,
Here α is the friction ratio, i.e. ratio of unit interface shear resistance τf to the undrained
strength, su, of the adjacent soil around the pipe; θD' is the half contact angle between the
pipe (of diameter D) and the seabed, so that θD'D is the contact perimeter.
The resistance of the seabed depends on the embedment or w/D ratio. When the
pipe is embedded it undergoes cyclic large horizontal movement in zones where buckling
occurs, due to thermal expansion and contraction. The SAFEBUCK JIP (Collberg, et al.
2011) highlighted that the axial breakout response can show a significant peak. The peak
in resistance that can occur falls away to a residual axial friction after breakout. A
significant peak in axial resistance can occur when the pipe moves axially for the first
time, or may be after a period of rest during which consolidation occurs. The first
movement is associated with the buckle formation. The displacement associated with the
peak is termed as the ‘mobilization displacement’. This effect has parallels in the brittle
This significant (or otherwise) of a peak in the axial response has been studied in
relation to pile capacity by Murff (1980) and Randolph (1983). They used load transfer
curves which strain-softened abruptly or progressively once the full shaft friction was
mobilized. However, the influence on the ultimate capacity of the pile of a brittle peak
and strain softening of the soil response reduces as the axial pile stiffness increases, due
24
Structural analysis of submarine pipelines
to progressive failure. Exactly the same effect is relevant to pipelines under axial loading,
The advantages and limitations of the alpha (or total stress) and beta (or effective stress)
approaches are well documented elsewhere (White and Randolph 2007; Oliphant and
Maconochie 2007; Randolph and White 2008a; Jewell and Ballard 2011; White et al.,
2011; White and Cathie, 2010). A new framework was proposed by White et al. (2012)
based on the concept of critical state of soil mechanics to incorporate the undrained and
drained conditions at the pipe-soil interface. Four elements of the frameworks are shown
in Figure 2-3.
Figure 2-3: Mechanisms affecting axial pipe-soil interaction (White et al. 2012)
Elastic stiffnesses
inconsistently in the literature, without a definitive set of rigorous coefficients to rely on.
For example, Tian et al. (2010) and Tian and Cassidy (2011) proposed taking horizontal
elastic stiffness as equal to vertical elastic stiffness in their advanced pipe–soil interaction
25
Structural analysis of submarine pipelines
model. Hodder and Cassidy (2010), in their plasticity model, considered the ratio of
horizontal elastic to vertical elastic stiffness to be 0.925, basing it on Gazetas et al. (1985)
and Gazetas and Tassoulas (1987) for strip footings on an elastic half-space. The DNV
Along with the vertical and lateral stiffnesses the pre-failure axial stiffness of the pipeline
lateral buckles or submarine slide impact. At an element level, the axial stiffness may be
estimated by assuming a simple distribution of shear stress around the perimeter of the
pile, similar to that for normal effective stress (Randolph, 2013). Where the shear
modulus of the soil increases proportionally with depth, the axial elastic pipe-soil stiffness
can be estimated by the product of the shear modulus gradient and the diameter of the
pipe.
Until the last decade very limited research was conducted on historic submarine
slide systems that lie in oil and gas development regions, such as the Storegga slide off
the coast of Norway (Bugge et al. 1998) and the margins of the Mississippi delta where
submarine mudslides have been triggered by major hurricanes (Gilbert et al. 2007). The
submarine slide hazards of Australia’s North West Shelf have been highlighted by
Hengesh et al. (2011, 2013), Zhang et al. (2015). The outcomes of a 3-year COFS-
MERIWA JIP funded by six operators on the impact of submarine slides on pipelines has
Global perspective:
As the oil and gas industry is extended into deeper water, there is an increase in the
prevalence of geohazards near submarine pipelines. The pipelines cross a wide range of
geographical terrain and must eventually cross the continental shelf break. The most
26
Structural analysis of submarine pipelines
significant geohazard in this region is the potential for landslides. Therefore, the
given a perspective on some of the challenges faced by the oil and gas operators when
siting and designing their pipelines in landslide prone area of Gulf of Mexico (see Figure
2-4).
Figure 2-4: Seabed topography and field architectures of the Mad Dog field, Atlantis field, and Mardi
Gras Export transportation system in the Southern Green Canyon area of the Gulf of Mexico
(Jeanjean et al., 2005).
There is also a plan to construct a high pressure deepwater pipeline over 1300 km (Nash
& Roberts, 2011) of highly variable geologic conditions between the Middle East and
India across the Indus River fan as shown in Figure 2-5. The performance of offshore and
onshore buried and on-bottom pipelines subjected to ground surface rupture, soil
27
Structural analysis of submarine pipelines
Figure 2-5: Morpho-tectonic map of proposed Oman-India deepwater pipelines (Nash & Roberts, 2011)
Seismic hazard analysis for the proposed Oman-India submarine pipeline was carried out
a few years ago (Campbell et al. 1996). The ground shaking hazard along the route of the
pipeline and in the Indus Canyon was found out to be relatively high. Due to this ground
shaking other potential geohazards such as liquefaction and triggering of submarine slides
may occur. Figure 2-6 shows the route of the pipeline and seismic activities in nearby
areas.
28
Structural analysis of submarine pipelines
Australian perspective:
The much discussed Greater Gorgon project of Western Australia, involves field
development of the shallow water Gorgon and deepwater Jansz fields off the north-west
coast of Australia. The pipeline between the Jansz field and the onshore LNG facility has
complex challenges including the deepwater origin (1350 m), potential debris flow off
the continental shelf and steep scarps to be crossed along the route to shallow water near
Figure 2-7: Pipeline route through continental shelf of the Jansz filed of Gorgon project (Equid,
2008).
To avoid these complex challenges, a relatively long southern route was identified.
However, the risks associated with deepwater pipelines, such as steep slopes and potential
impact from slides, were still inevitable for the new route.
There are many papers and technical reports available on the damage of shallow
water pipelines in the Gulf of Mexico from recent hurricanes, including Andrew, Lili,
Ivan, Katrina and Rita (DNV, 2006a). Gilbert, Nodine et al. (2006), Gilbert et al. (2007),
Hurricanes like Ivan in 2004 and Katrina in 2005 produced slides in this area that
29
Structural analysis of submarine pipelines
damaged many pipelines (Mineral Management Services 2005).. These mudslides are
generally localized features of several thousands of feet (1-3 km) long by hundreds of
feet (30 – 300 m) across by 50 to 100 feet (15 – 30 m) deep. Damage to the pipelines in
the vicinity of the mudslide zone are due to excessive longitudinal forces (Randolph et
al., 2010). More reliable techniques of assessing the normal and axial impact forces
generated by mudslide loading, which result in high tensile longitudinal forces in the
pipeline, would aid the prediction and mitigation of hurricane damage. However, in
contrast to the understanding of the common mudslide events that affect the dense
network of shallow Gulf of Mexico pipelines, there are very few reports in the technical
Sweeney et al. (2004) and Parker et al. (2008) presented a similar ‘string’ or ‘cable’
model of the submarine pipelines, where the pipe was assumed to resist tension but not
bending. In this model, the pipeline was subjected to uniform normal loading by the soil
inside the slide area and outside the slide area. Assuming that the normal force per unit
area is balanced in the two regions, from equilibrium, the normal force in the region
outside the slide occurs over a distance L/2 on each side, where L is the width of the slide.
It was shown from moment equilibrium that the ‘string’ deforms into a double parabola,
whose shape was defined as a function of applied load and tension at the centre of the
parabola. Both the models assume that the tension in the pipeline was uniform over a
central region, and that it decreased gradually to zero outside of this region, due to being
opposed by the axial soil resistance. Parker et al (2008) considered that the uniform
tension applies to the landslide region of width L, whereas Sweeney et al. (2004) applied
30
Structural analysis of submarine pipelines
the uniform tension to a region of length 2L (the entire region loaded laterally by the soil,
In a different mudslide loading model, Swanson & Jones et al. (1982) explicitly
accounted for the extra length (geometric slack) provided by laying the pipeline with
some curves. Their model was based on the equilibrium and compatibility equations for
an inextensible rod loaded axially and transversely allowing for the finite deflection. This
model was used to determine the maximum width of slide for different conditions by
changing the various parameters including pipe diameter, pipe thickness, pipe weight,
buried or surface laid and inclination of the slide with the pipe axis. The governing
equations were integrated in closed form between the characteristic nodes 1-4 as shown
in Figure 2-8.
Figure 2-8: Schematic of landslide impact on pipelines (Swanson &Jones et al. 1982)
The results showed that the burial of a pipeline outside the slide zone was safer because
active loading was reduced or eliminated and the passive support loading was increased.
It was also observed that with increasing inclination angle of the slide with the pipe axis,
the likelihood of failure reduced. Only downslope soil movement was considered in their
study.
31
Structural analysis of submarine pipelines
given by Summers & Nyman (1985). Their theoretical framework accounted for large
Total Potential, an approach closely related to the Principle of Virtual Work. However,
Figure 2-9: Lateral slide loading on pipeline model (Summers & Nyman, 1985)
Parker et al. (2008) verified the closed form solution with FE modelling using ANSYS
software. The pipeline was considered as a series of pipe elements supported by axial and
lateral elasto-plastic springs and the slide was modelled as a distributed load. Half of the
pipeline width was modelled, considering a fixed boundary condition outside the limits
of interest of the pipeline and a symmetric boundary condition at the centreline of the
slide. They presented their analysis of a ‘string’ model in the form of a chart showing the
maximum landslide pressure as a function of the landslide width for various pipe sizes.
They identified that the maximum landslide can be withstood, if the pipe length was
32
Structural analysis of submarine pipelines
greater than the (passive) anchor length. They also noted that, if the PLET (pipeline end
termination) was less than the minimum anchor length then the pipe would slide before
In another study Randolph et al. (2010) developed an analytical solution of lateral slide
pipe interaction. The loading on the pipe was divided into two zones, namely active
region and passive region. The active region was loaded by the slide itself and the passive
region was resisted by the surrounding soil as shown in Figure 2-10. Numerical modelling
was also carried out to verify the analysis. Firstly, normal loading was studied, and then
the analyses were expanded for the general case of a slide impacting the pipe at some
angle.
ABAQUS (2007) was used for the numerical analysis, with the pipe modelled by beam
elements and the pipe-soil response modelled using the PSI24 element available in
ABAQUS. The idealisation of the transverse slide loading on the pipe is shown in Figure
2-11. A parametric solution was presented to estimate the maximum tensile and bending
strains along with the maximum displacement of the pipe. From the parametric study it
33
Structural analysis of submarine pipelines
was revealed that for stiff pipelines or narrow slides, the structural load is dominated by
bending stresses. However, for wider slide and more flexible pipelines, as would typically
be the case in practice, the maximum stress was predominated by the tensile loading due
Figure 2-11: Finite element model used to verify the analytical solution (Randolph et al, 2010)
Liu & O’Rourke (2010) discussed the response of submarine pipelines due to slide
impact in the Mississippi Delta Region. An analytical solution was presented for pipelines
laid parallel to and perpendicular to the direction of seafloor movement (Figure 2-12). The
longitudinal effect from slide loading was based on the models discussed for buried
pipelines O’Rourke et al (1995) and the lateral loading was based on the cable model
34
Structural analysis of submarine pipelines
(a) (b)
Figure 2-12 (a) Longitudinal slide-pipe interaction (b) Lateral slide-pipe interaction(Liu &
O’Rourke, 2010)
generally small and the dominant forces are from high internal temperature and pressure,
which tend to cause expansion (Bruton et al., 2008). Axial resistance between the pipe
and the seabed opposes this expansion. Excessive compressive forces lead to buckling,
but the buckling response depends critically on the lateral soil resistance. When buckling
occurs, it significantly reduces the net axial load in the pipe. On the other hand, excessive
Hobbs (1981, 1984) was inspired by rail track buckling (Kerr, 1972). However, he
assumed the pipe as a straight beam resting on a rigid seabed with a friction factor
between the pipe and the seabed to derive the critical axial loads, maximum displacement
Taylor and Gan (1986b), Taylor and Tran (1993), Taylor and Tran (1996)
introduced initial imperfection into the analytical model proposed by Hobbs considering
deformation-dependent resistance force model (Taylor and Gan 1986a, Taylor and Gan
1987). Their results revealed that the ratio of amplitude and wavelength (/L) of the
35
Structural analysis of submarine pipelines
temperature and buckle deformation. Later Liu et al. (2014) extended the analytical
solution and provided the relationship between the critical buckling temperature,
buckling length and amplitude under different high-order global lateral buckling modes.
Hong et al ( 2015) compared results of analytical and finite element solutions for high
order lateral buckling modes with a single arch imperfection. Zhu et al. (2015)
Beam buckling in the vertical mode and in the lateral mode was studied by Chi &
Chiou (1995). Elasto-plastic beam elements were used to model a pipe with nominally
infinite length and large local deformation. Their results revealed that the buckling
behaviour of a buried submarine pipeline is influenced by the burial depth and the initial
imperfection of the pipe. The lateral buckling load is greatly reduced for shallow burial,
Torselleti et al. (1999) in the Hotpipe project carried out FE analysis with
ABAQUS of a submarine pipeline resting on a flat seabed. They employed the RIKS
method to compare the analytical buckling solution with their numerical results for initial
imperfection and pipe material properties. They also carried out analyses to select the
best element for the pipeline. The pipe was modelled using a combination of pipe and
elbow elements to incorporate the advantages of both elements. Pipe elements were used
for the straight section and elbow elements were used for the buckle section. Nonlinear
seabed pipe interaction was modelled in the lateral and axial directions. A number of
analyses were carried out to compare the available analytical buckling solution with the
numerical results. In the end they concluded that initial imperfection and lateral soil-pipe
friction forces are very important in the development of the stress-strain behaviour of the
buckle.
36
Structural analysis of submarine pipelines
influenced by lay induced imperfections. However, they failed to explore the significance
of soil properties in the buckling analyses. A buckling solution with nonlinear pipe
material was given by Peek & Yun (2004 ). They gave a scaled solution of lateral buckling
using FE analysis. They incorporated the elastic-plastic behaviour of the pipe material
and the lateral soil resistance was modelled with a series of nonlinear elastic springs
attached to each pipe node. The mobilisation displacement for this lateral soil resistance
force was taken to be small enough not to affect the results, but large enough for
numerical stability. The pipe was modelled with the elastic-perfectly plastic beam
Peek & Yun (2007) gave an analytical solution of lateral buckling triggered by flotation.
Then they used 3-D FE analysis carried out with ABAQUS to verify the analytical model.
The ABAQUS B31H element was used for the pipe and the seabed was modelled with
rigid elements. They found that from the FE analysis the bifurcation point occurred
slightly before the peak axial load. This leads to the conclusion that lateral instability
Submarine pipelines are increasingly operating at high temperature and high pressure
(HP/HT) (Carr et al. 2003). Figure 2-13 shows the trend and development of submarine
have a high tendency to Euler global buckling and walking. Recent projects have
quantifying and mitigating the more complex issues associated with walking. Pipeline
37
Structural analysis of submarine pipelines
cause a ratcheting response in the pipe lead to axial displacement (Carr et al 2003).
Figure 2-13: Trend of submarine pipelines being subjected to HP/HT (Carr et al 2003)
Palmer & Ling (1981) presented an analytical model for the pipeline expansion near
the platform, comparing the results with the two cases of North Sea pipelines. They
considered an elastic pipeline resting on a seabed that expands more towards the platform
and less away from the platform due to the fluctuating operating temperature and pressure.
The smaller expansion of the pipeline away from the platform was restrained by the
frictional resistance between the pipe and the seabed. When the operating temperature
and pressure are reduced, the movement of the pipeline reverses away from the platform.
The movement of the pipeline towards the platform was expressed separately for
exponential temperature increase with distance from the platform and also for uniform
from the platform, the expansion of the pipeline was expressed as:
38
Structural analysis of submarine pipelines
(
= 1 1 − exp(− z ) +
)
1 1 pRz
−
E 2 t
−
fz 2
4Rt
(2-1)
For a uniform temperature increase the expansion of the pipeline was given by:
2
RE (1 ) 2 t R1
p
= 1 + −
f E1 t 2 (2-2)
Poisson’s ratio, R is the radius, t is the thickness of the pipe material, f is the limiting
longitudinal friction per unit length, λ is a decay length over which the temperature
difference falls to 1/e (0.369) of its initial value, p is the internal pressure in the pipe and
Konuk (1998) developed a solution algorithm for the expansion of pipelines under
decaying ‘ratcheting’ mechanism for the pipelines subjected to cyclic temperature and
pressure conditions. A mathematical model was developed to derive the end expansion
due to uniform heating and exponential heating. A numerical algorithm was also
program EXPAND. Tørnes et al. (2000) studied the axial ‘ratcheting’ mechanism due to
thermal transient and the characterisation of ‘short’ and ‘long’ pipelines depending on
their expansion behaviours. For a ‘long’ pipeline during operating cycles the strain caused
by the axial frictional resistance due to pipe-soil interaction exactly counterbalances the
cumulative strains due to pressure and temperature. Hence, further expansion movement
is restricted and thus the pipe is fully restrained. However, for a ‘short’ pipeline the strain
caused by the frictional resistance between the pipeline and seabed is unable to balance
39
Structural analysis of submarine pipelines
the cumulative strains due to pressure and temperature even with friction has been
Perinet & Frazer (2006) illustrated the walking mechanism due to SCR tension. When
the pipes is subjected to SCR tension and undergoes a thermal cycle, the fixed point
during heat up and cool down steps differs and the pipe moves towards the SCR loading
at each cool down (resulting from a shutdown operation). The pipe expands in the
opposite direction from the fixed point with the expansion resisted by the soil resistance.
However, they did not provide an analytical solution for the behaviour of pipelines during
thermal cycles.
Carr et. al. (2006) gave expressions for the pipeline walking under various field
conditions, as part of the SAFEBUCK JIP. They noted that pipeline walking behaviour
of short pipelines occurs due to: 1) tension at each end of the pipeline, for example as
associated with a SCR; 2) global seabed slope along the length of the pipe; and 3) thermal
gradients along the pipeline during changes in operating conditions. The analytical
solution of pipeline walking per cycle, under the tension from an SCR was given by:
=
(P + S R )
− f L SR
EA f (2-3)
where ΔP is the change in fully constrained force, SR is the SCR tension, f is the axial
frictional force, L is the length of the pipeline and EA is the axial stiffness of the pipe.
Then the slope of the seabed (φ) was introduced in the above expression and the walking
=
( P + W L sin( ) − W L cos( ) ) L tan( )
EA (2-4)
40
Structural analysis of submarine pipelines
where, μ is the axial friction coefficient. The analytical model of walking due to thermal
transients was not straightforward and was solved by an incremental method. However,
a simple approximate solution for the distance walked per cycle was also given by:
f L2 f f
= 24 − − 4 if f > fθ/6
16 EA f f
(2-5)
f L2
= if f < fθ/6
8EA (2-6)
where, fθ is the force generated by a thermal transient. The above expressions were
validated with FEA analysis. The effect of the mobilisation displacement for limiting
axial friction was also studied in the paper. The results of the FEA analysis showed that
However, they mentioned that FEA analysis indicated that the mobilisation displacement
Throughout the paper the axial frictional resistance was considered to be a fixed value.
Rong et al. (2009) evaluated the axial walking behaviour of deepwater flowlines
along with the mitigation strategies. Detailed numerical analyses were carried out with
ABAQUS. In their analyses the nonlinearity of pipe-soil parameters, material of the pipe,
seabed profile, pressure and temperature profile and the connected PLET (pipeline end
termination) resistance were all considered. The pipe was modelled with 3D pipe element
(PIPE31H) and the seabed was modelled with a 3D rigid surface. Axial and lateral friction
factors were used with peak and residual values based on Carr et. al. (2006).
Quasi-static analysis was performed as this was considered reasonable for axial
walking. A detailed sensitivity study was performed to investigate the effects of friction
factor, axial mobilisation distance, PLET and any other in-line structural resistances.
41
Structural analysis of submarine pipelines
Their results showed that axial walking mainly depends on the temperature change and
the axial friction factor, with walking reduced as the friction factor is increased. Some
mitigation methods against walking were also discussed, among which anchoring was
Chudhury (2010) gave similar expressions for pipeline walking during start-up heating
and shut-down cooling cycles. The analytical tool proposed by him can be used in the
implemented at the design stage. He also validated the analytical model with FEA
analysis using ABAQUS. However, he noted that in typical shutdown and restart
operation cycles, the soil resistance during a heating up operation may be less than during
a cooling down cycle but this anomaly of soil behaviour has not been included in any of
Along with three reasons for thermal walking discussed by Carr et al. (2006), Bruton
et al. (2010) recently described other contributory reasons of walking. They identified
multiphase flow behaviour during shutdown and restart operations as one of the main
additional mechanisms of thermal walking. In multiphase pipelines that are laid on a slope,
the liquid separates out and settles at the bottom of the slope with any gas accumulating
at the top of the slope. This unequal density distribution up a slope enhances the walking
of the pipeline down the slope. During normal operating conditions the pipe moves more
in the downhill direction than uphill. Due to accumulation of the liquid down slope during
a shut down process, the increased submerged weight reduces the uphill contraction of
the pipe; the reverse situation happens for the uphill light gas filled section of the pipe.
This increases the absolute walking rate down the slope. An analytical model was
developed to estimate the effects of liquid hold up on walking for a flat seabed. Numerical
analysis was also carried out with ABAQUS to verify the analytical model. The numerical
42
Structural analysis of submarine pipelines
results matched closely with the analytical results. In addition, interaction between axial
walking and lateral buckling was also addressed, along with some mitigation methods.
The importance of combined effects of axial walking and riser load on flowline design
was investigated by Brunner et al. (2006). Their work gave an overview of the combined
flowline walking and riser dynamic loading response estimated by FE analysis using
Figure 2-14: Axial soil resistance in full heat up and cool down cycle (Brunner et. al, 2006)
Non-linear analysis was performed by applying thermal, pressure loading along with
dynamic SCR bottom tension at the touchdown point. Pipe elements were used to model
the pipe and 3D rigid elements were used to model the flat seabed. Non-linear friction
was applied between the pipe and the seabed. They found that maximum anchor loads
due to walking are generally greater than the riser dynamic loads alone, and therefore
govern the anchor pile design. They also concluded that high effective axial tension after
walking and during a cool-down cycle may also influence the minimum allowable radius
43
Structural analysis of submarine pipelines
Flowlines may be seated on undulating slopes and with cyclic thermal gradients
during operation and shut down operations. Generally global movement is associated
with short or non-cyclically constrained pipelines, which move towards the cold end or
down the slope of the seabed. Wang et al. (2010) presented an analytical model for long
and partially constrained pipelines that can move towards either end when subjected to
movement of the pipelines,. This migration of long pipelines depends on the localisation
of the transient pressure and temperature profile. The bi-linear axial friction was modelled
Figure 2-15: Typical axial friction profile considered by Wang et al. (2010)
For the analytical modelling they divided the whole length of the pipelines into a stick-
44
Structural analysis of submarine pipelines
They considered a different thermal profile for heating and cool down cycles for the stick-
slide model. The end displacement of the pipe after one cycle was plotted in Figure 2-17.
Figure 2-17: Movement of the pipe after 1st cycle with stick-slide model
The effect of friction mobilisation was studied over 40 cycles. They also noted the
earlier. In this section the available literature about pipe-soil interaction during walking
45
Structural analysis of submarine pipelines
is discussed in more detail. Advances in the study of soil-pipe interaction forces have
been described by many researchers (Cathie et al. 2005, White & Randolph 2007, and
Hill & Jacob (2008) identified knowledge gaps and resulting assumptions in the
modelling the soil behaviour in many of the previous studies. They noted that, as the
pipeline is heated from one end to the other, it moves axially relative to the soil with
variable speeds at different positions along the route of the pipeline. This range might
typically vary from 0.0002 mm/s to a maximum of approximately 0.7 mm/s. To address
the gap in the literature and available methods, they developed an in-situ testing system
hydraulic system. The axial, vertical and lateral movement of the pipe section, and
corresponding forces, were measured by transducers. The aim of the development was to
interpreting the data from a geotechnical perspective and understanding the consolidation
characteristics of the soil and the pore pressure response during pipe movement in detail.
The equipment was projected to save significant CAPEX for a project by testing in situ,
Generally axial pipe soil interaction is modelled with a simple friction coefficient
that links the available ultimate resistance to axial movement to the weight of the pipe.
However, the response may not be completely frictional and changes in pipe weight may
not completely correspond to the corresponding changes in axial resistance (Bruton et. al.
2008). A finite displacement must occur at very small loads to reach the ultimate
resistance. Due to this finite mobilisation displacement, the initiation of buckling and
46
Structural analysis of submarine pipelines
Figure 2-18: Axial soil pipe interaction schematics with ultimate resistance and mobilisation displacement
Figure 2-18 illustrates ductile, brittle and drained axial responses for pipelines, as
developed from laboratory testing. The brittle response consists of high breakout axial
resistance followed by a lower residual resistance. The breakout axial resistance can be
response is influenced by excess pore pressure generated at the pipe-soil invert. The very
small displacement at which the peak axial resistance occurs is referred to as the
Once the pipe starts moving, for example during a thermal increase, the friction
reduces to a residual value at a larger displacements. The residual friction dominates the
end expansion of the pipe and the feed-in into a lateral buckle. If the pipe displacement
occurs very slowly and no excess pore pressure is generated, the response may be termed
‘drained’ (Bruton et al 2008). Carneiro & Murphy (2011) observed the effect of soil
provided by Carr et al. (2006) is in good agreement with the FE results provided by White
et al. (2015). The difference between the results is due to the axial displacement, Xmob,
required to mobilize the full axial resistance in the FE analysis. Therefore, White et al.
47
Structural analysis of submarine pipelines
(2015) proposed an elastic correction of 2Xmob, subtracted from the predicted walk per
conventional analyses, seabed friction is independent of the rate of thermal loading and
expansion, although it has been recognised that the sliding resistance between a pipe and
the seabed varies with velocity, partly due to drainage effects. In conventional design
approaches, the rate of thermal loading during operating cycles is not considered within
the analysis. However, the velocity varies with position along the pipeline and with time
during each startup or shutdown event. Recently, it has been recognised that the sliding
resistance between a pipe and the seabed varies significantly with velocity due to drainage
effects (White et al. 2011, White & Cathie 2010), but to date much less research has been
done to incorporate this complex model (Carneiro et al. (2017) in any structural modelling
2.3 CONCLUSIONS
this chapter to lay the foundation of this thesis. A short discussion is presented here to
Numerical simulations to derive the elastic t-z stiffness of a pipe resting on the
seafloor are overly complex and a simple solution is needed. Limited prior information
exists for the axial elastic response of on-bottom pipelines in the literature, particularly
48
Structural analysis of submarine pipelines
the axial stiffness, focusing on the most relevant case for deepwater pipelines of shear
impact from submarine debris flow, has become a critical issue for the offshore pipeline
industry. The numerical modelling of the pipe structural response is also onerous –
varying slide loading, strength, impact velocity etc. which means that an analytical
solution for the structural response of a straight pipe impacted by axial slide loading is
needed.
cycles of high temperature and high pressure. This rises the possibilities of lateral
buckling and or pipeline walking issues. Classical theories to predict the critical buckling
loads are simple and easy to use without any complexity of the soil-pipe interaction.
However, in practice the geotechnical parameters play a role, along with the as laid
geometry and pipe material. Therefore, a simple correlation (including all the relevant
design chart to estimate the buckling load, valid for the range of inputs covered by the
parametric study.
soil behaviour and velocity dependent soil friction, simple analytical solutions supported
49
Structural analysis of submarine pipelines
In the concept stage of any submarine pipeline project, analytical solutions can
highlight the parameters that are most important when carrying out site investigation.
However, to validate final results of analytical solutions fine tuning of the design
In the following chapters axial pipe-soil interaction and structural analysis of submarine
pipe are presented, starting with analytical solutions followed by numerical modelling.
50
Structural analysis of submarine pipelines
References
Bai, Y. & Bai, Q. (2005). Subsea Pipelines and Risers. Elsevier Ltd.
Brunner, M. S., Qi, X., Zheng, J., & Chao, J. C. (2006) . Combined effect of flowline
walking and riser dynamic load on HP/HT flowline design. Proc. of Offshore
Bruton, D., Sinclair, F., & Carr, M. (2010) . Lessons learned from observed walking of
pipelines with lateral buckles, including new driving mechanisms and updated
Paper 20750
Bruton D. A. S., White D. J., Cheuk C. Y., Bolton M. D. & Carr M. C. (2006). Pipe-soil
displacement tests by the Safebuck JIP. Proc. Offshore Technology Conf., Houston,
during lateral buckling and pipeline walking – the Safebuck JIP. Proc. of Offshore
Bugge, T., Belderson, R.H. and Kenyon, N.H. (1998). The Storrega slide. Phil. Trans.
Campbell, K.W., Thenhaus, P.C. Mullee, J.E. & Preston, R. (1996). Seismic hazard
evaluation of the Oman India pipeline. Proc. of Offshore Technology Conference,
Houston, Paper OTC 8135.
Carneiro, D. & Murphy, D. (2011). Simple numerical models for pipeline walking
accounting for mitigation and complex soil response. Pro. of the ASME 2011 30th
Int. Conf. on Ocean, Offshore and Arctic Engineering OMAE2011.The
Netherlands. OMAE2011-49780.
51
Structural analysis of submarine pipelines
Carneiro, D., Rathbone, A., Siong, Soon, Kok., & Viecelli, G. (2017). Velocity-
10.1115/1.4034695.
Carr, M., Sinclair F, & Bruton D.(2006). Pipeline walking — understanding the field
layout challenges, and analytical solutions developed for the SAFEBUCK JIP.
Cathie D. N., Jaeck C, Ballard J. C. & Wintgens J-F. (2005). Pipeline geotechnics – state-
Chatterjee, S., White, D.J. and Randolph, M.F. (2012). Numerical simulations of pipe-
soil interaction during large lateral movements on clay. Géotechnique 62(8), 693–
705.
Chi, S. Y., & Chiou, Y. J. (1995). Buckling of offshore pipelines. Proc. 5th Int. Offshore
Collberg, Leif & Mo̸rk, Kim & Levold, Erik & Vitali, Luigino. (2005). HotPipe JIP:
Collberg, L., Carr, M. and Levold, E. (2011). Safebuck Design Guideline and DNV-RP-
Dassault Systèmes (2007). ABAQUS analysis users’ manual, Simula Corp, Providence,
RI, USA.
DNV. (2006a). Technical Report, Pipeline damage assessment from Hurricane Ivan in
52
Structural analysis of submarine pipelines
DNV. (2006b). Free spanning pipelines. Recommend practice RP-F105, Det Norske
Gazetas, G., Dobry, R., and Tassoulas, J. L. (1985). Vertical response of arbitrarily
9410(1985)111:6(750), 750–771.
9410(1987)113:5(440), 440–457.
Gerwick B.C. (2007) Construction of marine and offshore structures. CRC Press, New
York, 795 p.
Gilbert R.B., Nodine, M.C., Wright, S.G., Cheon, J.Y., Wrzyszczynski, M., Coyne, M.
Hengesh, J., Whitney, B., Rovere, A. (2011) A Tectonic Influence on Seafloor Stability
along Australia’s North West Shelf. Int. Offshore and Polar Engineering
Hengesh, J., Dirstein, J.K., Stanley, A.J. (2013) Landslide geomorphology along the
53
Structural analysis of submarine pipelines
Hill, A. J., & Jacob, H. (2008). In-situ measurement of pipe-soil interaction in deepwater.
Hobbs, R. E. (1981). Pipeline buckling caused by axial loads. J. Construct. Steel Res.
1(2), 2-10.
175-189.
Hodder,M. S., and Cassidy,M. J. (2010). A plasticity model for predicting the vertical
Hong, Z., Liu, R., Liu, W., & Yan, S. (2015). Study on lateral buckling characteristics of
61-83. https://doi.org/10.1016/S0266-352X(98)00012-3
Jeanjean, P., Hampson, K., Evans, T., Liedtke, E., and Clukey, E.C. (2005). An operator’s
perspective on offshore risk assessment and geotechnical design in geohazard-
prone areas. Proc, 2nd Int. Symp. Frontiers in Offshore Geotechnics, ISFOG 2010,
Perth, 115-143.
Jewell, R.J. & Ballard, J-C. (2011). Axial pipe-soil interaction – a suggested framework..
Proc. Offshore Tech. Conf., Houston, USA. Paper OTC 22010.
Kerr, A.D. (1972). On the stability of the railroad track in the vertical plane (No. NYU-
AA-72-35)
54
Structural analysis of submarine pipelines
OMAE98-1103
Konuk, I. & Yu, S. (2007). Continuum FE modelling of lateral buckling: study of soil
effects. Proc. of 26th Int. Conf. on Offshore Mechanics and Arctic Eng., San Diego,
OMAE2007-29376.
Liu, R., Liu, W.B., Wu, X.L., & Yan, S.W. (2014). Global lateral buckling analysis of
Martin, C.M., Kong, D. & Byrne, B.W. (2013). 3D analysis of transverse pipe-soil
Martin, C.M. and White, D.J. (2012). Limit analysis of the undrained bearing capacity of
Murff, J. D., Wagner, D. A. & Randolph, M. F. (1989). Pipe penetration in cohesive soil.
Nash, I., & Roberts, P. M. (2011). The deepwater gas route to India. Proc. 34th Annual
Offshore Pipeline Technology Conference, Amsterdam, UK.
Nodine, Mary C., Wright, Stephen G., Gilbert, Robert B. and Ward, E.G. (2006).
Mudflow and mudslide during Hurricane Ivan. Proc. Offshore Technology Conf.,
Houston, OTC 18328.
Oliphant, J. & Maconochie, A. (2007). The axial resistance of buried and unburied
pipelines. Proc. 6th Int. Conf. Offshore Site Investigation and Geotechnics.
London: Society for Underwater Technology, 125–132.
O’Rourke, M.J., Liu, X.J., and Flores-Berrones, R.,(1995). Steel pipe wrinkling due to
longitudinal permanent ground deformations. Journal of Transportation
Engineering, September/October, Vol. 121, No. 5, pp. 443-451.
Palmer, A.C., & Ling, M. T.S. (1981) Movement of submarine pipelines close to
platforms. Proc. 13th Offshore Technology Conf., Houston, Paper OTC 4067.
Parker, E., Traverso, C., Moore. R., Evans, T., Usher, N. (2008). Evaluation of landslide
impact on deepwater submarine pipelines. Proc. Offshore Technology Conf.,
Houston, Paper OTC 19459.
Peek, R. & Yun, H. (2004). Scaling of solutions for the lateral buckling of elastic plastic
pipelines, 23rd Int. Conf. on Offshore Mechanics and Arctic Engineering, Vol. 3,
Vancouver, British Columbia, Canada. OMAE2004-51054, 1-8
Peek, R. & Yun, H. (2007). Floation to trigger lateral buckles in pipelines on a flat seabed.
J. Eng. Mech. 133 (4), 442-451.
Perinet, D. & Frazer, I. (2006). Mitigation Methods for Deepwater Pipeline Instability
Induced by Pressure and Temperature Variations. Proc. Offshore Technology
Conference, Houston, Paper OTC 17815-MS.
Randolph, M.F. (1983). Design considerations for offshore piles. Proc. Conf. on Geot.
Practice in Offshore Engng, ASCE, Austin, 422-439.
56
Structural analysis of submarine pipelines
Randolph, M. F., Seo, D., & White, D. J. (2010). Parametric solution of slide impact on
pipelines. Journal of Geotech and Geoenvironmental Engineering, ASCE, July,
136(7), 940-949.
London, 27-59.
OTC 19128.
Rong, H., Inglis, R., Bell, G., Huang, Z., & Chan, R. (2009). Evaluation and mitigation
Sriskandarajah, T., Dong, S., S. Sribalachandran, S., & Wilkins, R. (1999). Effect of
initial imperfection on the lateral buckling of Subsea pipelines. Proc. 9th International
Summers, P. B., & Nyman, D. J. (1985). An approximate proceduce for assessing the
effects of mudslides on offshore pipelines. Offshore Mechanics & Arctic
Engineering Division, Transaction of ASME, Vol. 107. December, 426-432.
Sweeney, M., Gasca, A. H., Garcia Lopez, M., and Palmer, A. C. (2004). Pipelines and
landslide in rugged terrain: A database, historic risk and pipeline vulnerability.
Proc. Int. Conf. on Terrain and Geohazard Challenges Facing Oil and Gas
57
Structural analysis of submarine pipelines
Taylor, N., Gan, A.B., 1986a. Refined modelling for the lateral buckling of submarine
Taylor, N., Gan, A.B., 1986b. Submarine pipeline buckling-imperfection studies. Thin-
Taylor, N., Gan, A.B., 1987. Refined modelling for the vertical buckling of submarine
Taylor, N., Tran, V., 1993. Prop-imperfection subsea pipeline buckling. Mar. Struct. 6
(4), 325–358.
Taylor, N., Tran, V., 1996. Experimental and theoretical studies in subsea pipeline
Tian, Y., and Cassidy, M. J. (2011). Pipe-soil interaction model incorporating large lateral
Tian, Y., Cassidy, M. J., and Gaudin, C. (2010). Advancing pipe-soil interaction models
through geotechnical centrifuge testing in calcareous sand. Appl. Ocean Res., 32(3),
284–297.
Torselleti, E., Vitali, L., & Levold, E. (1999). Hotpipe Project: Snaking of submarine
pipelines resting on flat sea bottom using finite element method. Proc. 9th Int.
Verley, R. L. P. & Lund, K. M. (1995). A soil resistance model for pipelines placed on
clay soils. Proc. Offshore Mechanics and Arctic Engineering Conf, Copenhagen,
58
Structural analysis of submarine pipelines
Verley, R. L. P. & Sotberg, T. (1992). A soil resistance model for pipelines placed on
sandy soils. Proc. Offshore Mechanics and Arctic Engineering Conf, Vol V-A
Westgate, Z. W., White, D. J., Randolph, M. F. & Brunning, P. (2010). Pipeline laying
White, D. J, Bolton, M. D., Ganesan, S. A., Bruton, D., Ballard, J-C., & Langford, T.
interaction on soft natural clays. OTC 21249. Proc. Offshore Tech. Conf.,
Houston, USA.
White, D.J., Campbell, M., Boylan, N. & Bransby., M.F. (2012). A new framework for
White, D.J. Cathie, D.N. (2011). Geotechnics for subsea pipelines. Frontiers in Offshore
59
Structural analysis of submarine pipelines
White, D. J., & Cheuk, C. Y. (2008). Modelling the soil resistance on seabed pipelines
during large cycles of lateral movement. Marine Structures, 21(1), 59-79. DOI:
10.1016/j.marstruc.2007.05.001
White, D.J., Clukey, E.C., Randolph, M.F., Boylan, N.P., Bransby, M.F., Hill, A.J. and
pipeline design. Proc. Offshore Technology Conf., Houston, Paper OTC 27623.
White, D. J., & Dingle, H. R. C. (2011). The mechanism of steady friction between seabed
10.1680/geot.8.T.036
White, David, Hill, A., Westgate, Z., & Ballard, Jean-Christophe. (2010). Observations
of pipe-soil response from the first deep water deployment of the SMARTPIPE®.
10.1201/b10132-125.
White, D. J., & Randolph, M. F. (2007). Seabed characterisation and models for pipeline-
193-204.
White, D.J., Randolph, M.F., Gaudin, C., Boylan, N., Wang, D., Boukpeti, N., Zhu, H &
10.4043/27034-MS.
White, D.J., Westgate, Z., Ballard, J.-C, De Brier, C. & Bransby, M.F. (2015). Best
Yu, S. & Konuk, I. (2007). Continuum FE modelling of lateral buckling. Proc. Offshore
60
Structural analysis of submarine pipelines
Zhang J., Erbrich C.T., Finnie I., Neubecker S., White D.J., Deeks A.D., Doh G.,
Cumming G., Brown N. & Little R. (2015). Geotechnical design and construction
Zhu, J., Attard, M. M., & Kellermann, D. C. (2015). In-plane nonlinear localised lateral
61
Structural analysis of submarine pipelines
62
Structural analysis of submarine pipelines
3.1 OVERVIEW
The major focus of this research is to analyse the structural stability of submarine
pipelines due to loading from a submarine slide. To investigate the effect of axial loading
from the slide a comprehensive simple analytical model comprising an elastic straight
pipe lying on an elastic-plastic seabed was developed. To support the proposed model,
detailed numerical analysis was carried out for an elastic pipe on a non-linear seabed.
Buckling due to the compressive loading on the pipeline was neglected in the analytical
the pipe and distributed along its length. Before starting analysis of the pipeline under
loading from a submarine slide, it is necessary to assess the various input soil parameters
in terms of pre-failure (elastic) stiffness and limiting resistance for different loading
modes. The present chapter focuses on elastic stiffness parameters for axial, horizontal
(lateral) and vertical motions of the pipeline relative to the seabed, with the aim to express
Selection of elastic stiffness for pipe-soil interaction springs is covered poorly in the
literature. For example, Tian et al. (2010), Tian & Cassidy (2010) proposed taking the
horizontal elastic stiffness as equal to the vertical elastic stiffness for their advanced pipe-
soil interaction model. Hodder & Cassidy (2010) in their plasticity model considered the
ratio of horizontal elastic stiffness to the vertical elastic stiffness as 0.925, which was
63
Structural analysis of submarine pipelines
based on Gazetas et al. (1985) and Gazetas & Tassoulas (1987) for strip footings on elastic
half space. None of the above considered the axial elastic stiffness, which is typically
defined in terms of the mobilisation distance, expressed as some fraction of the pipeline
diameter, to reach the limiting axial friction, but without considering the relationship
A brief literature review was carried out to establish solutions for the elastic stiffness of
soil, under different modes of loading. (The descriptor “non-homogeneous” is used here
in the sense of soil of a single type (e.g. clay or sand), but where the elastic modulus
various in some continuous manner with depth below mudline.) Very limited information
was available for the axial elastic response of on-bottom pipelines, particularly for non-
homogeneous soil. Therefore, a new analytical method was developed for the axial
stiffness, which was verified through numerical analysis. Numerical analysis was also
carried out to obtain relationships for the horizontal and vertical elastic stiffness of on-
bottom pipelines. Thereafter, relationships amongst the elastic stiffnesses were developed
and recommendations made for the selection of proper elastic stiffnesses in all three
directions of motion.
The nonlinear behaviour of the seabed may be idealised by nonlinear springs attached to
the pipeline as shown in Error! Reference source not found.. The loads per unit length o
n the pipeline are defined as vertical load V, parallel to the z-axis, horizontal (lateral) load
H, parallel to the y-axis and axial load F parallel to the x-axis (see Figure 3-2). Likewise
the displacements along the x, y and z axes are δx, δy, and δz respectively as shown in
Figure 3-2. The stiffnesses in the three directions kv, ky, kx are defined as load (per unit
length) per unit displacement, so with typical units of kN/m2 or kPa, and are expressed as
64
Structural analysis of submarine pipelines
axial direction. The aim is to arrive at non-dimensional forms for the stiffnesses - for
soil shear modulus and Kx is a function of the pipeline embedment, normalised by the
diameter as w/D, and a factor reflecting the degree of non-homogeneity of the soil.
pipeline of nominally infinite length. Stiffness expressions for such foundations are
length (larger dimension in plan), taken as the dimension in the x-direction, and a
considered foundation’s longer side rather than by its shorter side is to approximate
Following the literature review, the axial elastic stiffness is derived for an on-bottom
that varies proportionally with depth. Then a comprehensive chart is prepared to convert
expressions for foundation stiffness from the literature into equivalent expressions for a
pipeline. These are presented in Table 3-1, Table 3-2 and Table 3-3. The dimensionless
stiffnesses of different shapes are also plotted with increasing aspect ratio for comparison,
This section reviews expressions for the vertical elastic stiffness of rigid rectangular,
circular and strip footings on the surface and embedded in soil with homogeneous and
65
Structural analysis of submarine pipelines
3.3.1 Rectangular
Solutions for rectangular embedded foundations are considered here, translating the
Pais and Kausel (1988) proposed approximate formulas to describe the dynamic
accurate numerical solutions. The dimensions considered by them are shown in Figure
3-4. The static elastic stiffness of an embedded rectangular footing was expressed as a
function of the dynamic shear modulus of the soil, G, Poisson’s ratio, , the footing semi-
K Vo (1 − ) L
0.75
(3-1)
= 3.1 + 1.6
GB B
where the curly brackets indicate that the notations used are based on the original
publication, not on those used in this thesis. Rocking and torsional stiffnesses given by
Pais and Kausel (1988) are not considered in this study. This expression is then converted
Dobry and Gazetas (1986) gave another set of equations for the vertical elastic stiffnesses
number of analytical results, which were supported by numerical analysis and interpreted
by simple physical models. Expressions for rectangular, circular and strip footings are
considered here. The vertical elastic stiffness of a rectangular footing was expressed as
66
Structural analysis of submarine pipelines
2LG (3-3)
K z = S z
1 −
For foundations of different shapes, they found that the dimensionless parameter Sz could
be explicitly related to the base-shape parameter A/4L2, where A is the base area of the
foundation, 2L and 2B are the length and width of the rectangle respectively (see Figure
3-5). The dimensionless base-shape parameter A/4L2 is equal to 1 for a square, 0.785 for
a circle and 0 for an infinitely long strip. These researchers suggested taking
(
S z = 0.73 + 1.54 A / 4 L2 )
0.75
for A/4L2 > 0.02 and Sz = 0.8 for A/4L2 ≤ 0.02. For
rectangular shapes, A/4L2 may be taken as equivalent to B/L. For long rectangular
foundations the above expression can therefore be expressed in the present convention as
kv 0.8 (3-4)
Kv = =
G 1 −
Poulos and Davis (1974) gave approximate solutions for the vertical displacement due to
symmetrical vertical loading of different shaped rigid footings on elastic soil. The
P (1 − 2 ) (3-5)
z =
z BL E
where P is the symmetrical vertical load on the rigid base, B and L are the dimensions of
the short (y-direction) and long (x-direction) sides of the rectangle respectively and βz is
a factor that depends on the aspect ratio, with 1 < βz < 1.4 for 1 ≤ L/B ≤ 10. The expression
67
Structural analysis of submarine pipelines
k (3-6)
2 b
Kv = v = x
lG 1− l
3.3.2 Circular
Poulos and Davis (1974) and also Dobry and Gazetas (1986) gave similar expression
4GB (3-7)
K z =
1 −
where B is the radius of the circular foundation. The above expression can be expressed
kv 2 (3-8)
Kv = =
G 1 −
3.3.3 Strip
Dobry and Gazetas (1986) suggested taking Sz = 0.8 in equation (3-3) for the vertical
stiffness of a strip footing, where B/L tends to zero, hence resulting in the normalised
surface of the soil. Gazetas et al. (1985) expressed the vertical elastic stiffness of an
arbitrary rigid foundation embedded in a reasonably homogeneous and deep soil deposit.
This is an extension of earlier expression by Dobry and Gazetas (1986) in equation (3-3).
The surface elastic stiffness was factored by two expressions to derive the embedded
stiffness. The two expressions represented what were termed the ‘trench’ and ‘side-wall-
68
Structural analysis of submarine pipelines
contact’ effects, and were investigated separately. The final embedded vertical stiffness
2 LGS z 1 D 4 Ab A
2/3
K emb = 1 + 1 + 1 + 0.19 s
2
(3-9)
1 − 21 B 3 4 L Ab
where Ab = is the area of the basement-soil contact surface, 2L and 2B are the length and
width of the rectangle circumscribed to the actual shape of the footing as shown in Figure
3-5. D is the embedment depth, d/2 is half the depth of side-wall contact (see Figure 3-6).
This may be approximately applied to the case of a pipeline of invert embedment w and
contact width B, with maximum w/B 0.5 for a pipeline half-embedded in the soil. In this
4 B w
(3-10)
2/3
k S 1 2w
Kv = v = z 1 + 21 B 1 + 3 L 1 + 0.19 2 B
G 1 −
Similarly Pais and Kausel (1988) also expressed the static vertical elastic stiffness for an
s o 0.25 0.8
(3-11)
K V = K V 1.0 + 0.25 + (E / B )
L/ B
where E/B is the embedment ratio. This expression can be converted into the present non-
0.25 0.8
(3-12)
K v _ emb = K V 1.0 + 0.25 + (2w / b )
l /b
may be considered to have approximately half the lateral p-y stiffness at depths sufficient
to negate the surface effects for a pile in homogeneous soil. Analytical solutions for the
latter (e.g., Baguelin et al. 1977; Osman and Randolph 2012)for undrained conditions (ν
69
Structural analysis of submarine pipelines
= 0.5) suggest a p-y stiffness of four to five times the shear modulus of the soil G, leading
K v 2.25 (3-13)
movement as
k v 1.3 (3-14)
Kv = =
G 1 −
Matsubara and Hoshiya (2000) gave spring vertical stiffness for a buried pipeline as a
kv 2 (3-15)
Kv = =
G log e
recommended practice (DNV 2019) suggests a vertical elastic stiffness derived from a
0.88G (3-16)
K v =
1 −
Converting the above expression into the present dimensionless format leads to
70
Structural analysis of submarine pipelines
k v 0.88 (3-17)
Kv = =
G 1 −
As may be seen, this is 10 % greater than the limiting case given by Dobry and Gazetas
modulus is assumed to increase linearly with depth, as shown in Figure 3-7, expressed as
G(z) = G(0)+ mz . In the limiting case where the shear modulus is zero at the surface,
G(0) = 0, the expression for shear modulus becomes G(z) = mz. The stiffness falls to zero,
apart from the specific case of incompressible soil ( = 0.5). For that case, the elastic
vertical settlement δz, due to uniform pressure q, on any arbitrary-shaped rigid footing at
z = q / 2m (3-18)
as vertical load per unit area. Therefore, for a long slender foundation with width B, q =
V / LB. Taking a reference shear modulus, G*, as the value at a depth of the pipeline
kv 2mLD' (3-19)
Kv = = =2
G* LmD'
This section gives a brief review of expressions for the horizontal or lateral (as distinguish
from axial, i.e. y-direction) elastic stiffness for rectangular, circular and strip footings,
71
Structural analysis of submarine pipelines
3.4.1 Rectangular
Pais and Kausel (1988) gave approximate formulae for the elastic stiffness in the y-
direction. However, they expressed the elastic stiffness in the lateral direction in terms of
the elastic stiffness in the axial (long) x-direction, with an additional additive term as:
K Hy
o
(2 − ) K Hxo (2 − ) L
(3-20)
= + 0.8 − 1
GB GB B
The above expression can be converted into the present convention as:
1 B
0.35
ky kx 0.4 B B (3-21)
Ky = = + 1 − = 0.4 + 3.4 + 0.8
G G (2 − ) L (2 − ) L L
In contrast to the zero asymptotic vertical stiffness reached as B/L reduces to zero, the
Dobry and Gazetas (1986) and Gazetas and Tassoulas (1987) gave an expression for the
in a reasonably deep and homogeneous soil deposit, based on boundary element analysis.
The expression for the surface elastic stiffness in the lateral direction for a rectangular
(2 − ) K y , sur A
0.85 (3-22)
S y = 2 + 2.5 2
2 LG 4 L
Similar to the vertical case the shape parameter, A/4L2 is replaced by B/L and the non-
1 B
0.85 (3-23)
ky
Ky = = 2 + 2.5
G 2 − L
A limit of 2.24 for the square-bracketed term was proposed by Gazetas and Tassoulas
72
Structural analysis of submarine pipelines
Poulos and Davis (1974) gave an approximate solution for the horizontal displacement
Q(1 − 2 )
(3-24)
h =
X BL E
where Q is the uniform horizontal load on a rectangular surface with length L and width
B and βz is a factor depending on L/B and ν, with 0.69 < βx < 1.25 for 1 < L/B <10. This
leads to a stiffness of
ky 2 B
Ky = = z (3-25)
G 1 − L
3.4.2 Circular
Poulos & Davis (1974) gave expression for the horizontal elastic stiffness of a rigid
(7 − 8 )(1 + ) p xa
x = (3-26)
16(1 − ) E
where px is the average applied pressure and expressed as horizontal load per unit area,
P/πa2, P being the load and a the radius of the circle. The above expression can be re-
ky 16(1 − )
Ky = = (3-27)
G 7 − 8
For the same condition Gazetas & Tassoulas (1987) expressed the elastic stiffness of a
8GR (3-28)
K y =
2 −
where R is the radius of the footing. The above expression can be re-written in non-
ky 4 (3-29)
Ky = =
G 2 −
3.4.3 Strip
For strip footings, Gazetas & Tassoulas (1987) gave the horizontal elastic stiffness as
K y 2.24G (3-30)
=
2L 2 −
ky 2.24
Ky = = (3-31)
G 2 −
K L = 0.76G(1 + ) (3-32)
ky
Ky = = 0.76(1 + ) (3-33)
G
investigated by Gazetas and Tassoulas (1987). The surface elastic stiffness, Ky,sur was
given by equation (3-19). Therefore, the overall expression was expressed as:
D h Aw
0.40 (3-34)
K y = K y , sur I tre I wall = K y , sur 1 + 0.15 1 + 0.52 2
B B L
where D is the embedment, h is the side wall’s centre of gravity from top of the ground
as shown in Figure 3-8 and Aw is area of the sidewall soil interface. The above expression
2 w w 2 B
0.4
(3-35)
K y _ emb
= K y 1 + 0.15
1 + 0.52 2
B B L
74
Structural analysis of submarine pipelines
where w/B will be considered as 0.5 for a pipeline embedded to half a diameter. (in-
Pais and Kausel (1988) expressed the static horizontal elastic stiffness for an embedded
rectangular footing as
s o (3-36)
K H = K H 1.0 + 0.33 +
1.34
(E / B )0.8
1+ L / B
where E/B is the embedment ratio and the same expression can be converted into present
non-dimensional form as
(2w / B )0.8
1.34 (3-37)
K y _ emb = K y 1.0 + 0.33 +
1+ L / B
which leads to a maximum increase by a factor of 1.33 compared with a surface strip
similar in geometry to a half pile in homogeneous soil, loaded laterally. Again, due to
symmetry, the pipeline horizontal (lateral) stiffness is expected to be close to half the pile
lateral “p-y” stiffness. The horizontal pipeline stiffness would therefore be identical to
Published work on shallow foundations tends to focus on vertical and lateral response
(parallel to the short side), providing only limited information on axial response (parallel
to the long side). There are, however, a few expressions in the literature for the
longitudinal stiffness of rectangular footings. Some insight may also be gained from
solutions for axial loading of piles that consider a pipeline half-embedded in the soil to
75
Structural analysis of submarine pipelines
3.5.1 Rectangular
Gazetas and Tassoulas (1987) gave an expression for the longitudinal stiffness of a
0.21 B (3-38)
K x , sur = K y , sur − GL1 −
0.75 − L
k 0.21 B 3-39)
Kx = x = Ky − 1 −
G 2(0.75 − ) L
Pais and Kausel (1988) gave an expression for the dynamic axial elastic stiffness of a
rectangular footing as
K Hx
0
(2 − ) L
0.65
(3-40)
= 6.8 + 2.4
GB B
1 B B
0.35
kx
Kx = = 3.4 + 1.2 (3-41)
G (2 − ) L L
considered by Gazetas & Tassoulas (1987), adjusting the stiffness for a surface footing
in a similar way as for the vertical and lateral stiffness. The resulting expression was
D h Aw
0.40 (3-42)
K x = K x , sur I tre I wall = K x , sur 1 + 0.15 1 + 0.52 2
B B L
76
Structural analysis of submarine pipelines
2w w 2 B
0.40
(3-43)
K x _ emb = K x 1 + 0.15 1 + 0.52 2
B B L
where w/b will be considered as 0.5 for a pipeline embedded to half a diameter. .
Pais and Kausel (1988) derived the static elastic stiffness of an embedded rectangular
s o
K H = K H 1.0 + 0.33 +
1.34
(E / B )0.8 (3-44)
1+ L / B
(2w / B )0.8
1.34 (3-45)
K x _ emb = K x 1.0 + 0.33 +
1+ L / B
pile, the axial stiffness (per unit length) may be taken as identical to that customarily
assumed in load transfer analysis of piles. Randolph and Wroth (1978) expressed the
2G a (3-46)
0 =
D
where D is the pile diameter. The parameter, , was expressed in terms of the pile length
(normalised by diameter) and Poisson’s ratio, but was typically about 4. The
D 0 (3-47)
Kx =
2 G a 4
In order to compare this expression with others for rectangular foundations, the
1999)
77
Structural analysis of submarine pipelines
L 3-48)
ln 5 + 5(1 − )
B
Kx = 2 (3-49)
Most of the solutions considered so far, with the exception of Gibson (1967), have
considered the soil as homogeneous, with constant shear modulus, G. More commonly,
for on-bottom pipelines, the soil modulus will tend to increase with depth. It is therefore
pertinent to explore further the extreme case where the shear modulus is proportional to
depth, G = mz. An analytical solution is developed below for the elastic axial stiffness for
considering concentric semi cylinders of soil, with the shear stresses on each cylinder
decreasing inversely with the radius of the cylinder. The pipeline itself is considered as
shown in Figure 3-9, with embedment, w, and contact width, D'. The contact width may
be expressed in terms of the subtended contact angle D', which is given by (White and
Randolph, 2007)
2w (3-50)
cos D = 1 −
D
where D is the diameter of the pipeline. A notional pipeline, half embedded in the soil,
78
Structural analysis of submarine pipelines
Figure 3-9 shows the notional inner pipe, of radius, R = D'/2, and a general point, P,
within the soil at a radius, r, and angle, , from the vertical. The shear stress p at the
where p,max is the maximum value at the pipe invert. This is consistent with the
Following the same approach as Randolph and Wroth (1978), equilibrium parallel to the
length of the pipe leads to the shear stress decreasing inversely with the distance from the
pipe according to
R R (3-52)
= 0 = p, max cos
r r
Just as for axial loading of a pile, the deformations in the soil can be approximated by a
du p, max R
= xr = = cos (3-53)
dr G G r
Assuming that the shear modulus varies proportionally with depth, and noting that
du p, max R
= (3-54)
dr m r2
where m is the gradient of shear modulus. This may be integrated to give the axial
du p ,max
x = dr = (3-55)
R
dr m
79
Structural analysis of submarine pipelines
Now integrating the value of τp around the pipe (ACB in Figure 3-10), leads to the
/2 /2
F 2 2R
x
=
x p Rd =
0
x
0
p , max cos d (3-56)
F p ,max
= 2R = 2mR = mD sin D (3-57)
x x
The resulting axial stiffness may be normalised by the shear modulus at a level of half
2 F D' w
Kx = = 2 sin D = 2 = 2 sin arccos1 − 2 (3-58)
mD x D D
Interestingly, this is exactly half the vertical stiffness derived for a pipeline embedded in
soil with modulus proportional to depth (Randolph, 2013). Alternatively, if the stiffness
is normalised by the shear modulus at a depth equal to the contact width, D', the non-
1 F
Kx = =1 (3-59)
mD sin D x
This compares with the non-dimensional stiffness of Kx = /4 estimated from the half-
pile solution for the case of homogeneous soil (equation (3-44)). No pipe-soil slippage
3.7 DISCUSSION
The aim of this chapter is to quantify the elastic stiffnesses of on-bottom pipelines and
the relationship between the stiffness in each direction (vertical, lateral and axial). It is
useful to compare the results of the literature review on elastic stiffnesses in each
direction, together with the approximate analytical solutions for axial stiffness of
80
Structural analysis of submarine pipelines
pipelines in non-homogeneous soil. The various values are plotted here for comparison.
In addition, relevant results from the finite element analyses presented later are also
shown. In many cases, the non-dimensional stiffnesses are a function of the aspect ratio,
Figure 3-11 shows the variation in nondimensional vertical elastic stiffness with aspect
ratio. For circular and square foundations, the Poulos and Davis (1974) and Dobry and
Gazetas (1986) solutions give similar values of 3.33. For a rectangular foundation, the
stiffnesses defined by Poulos and Davis (1974) and Pais and Kausel (1988) show similar
patterns although with diverging trends at high L/B (beyond the range of the original
studies), reducing to 0.46 and 0.83, respectively for L/B = 100 compared with the lower
limit of 1.33 proposed by Dobry and Gazetas (1986) for L/B ≥ 50. For an embedment
ratio of w/B = 0.5, the Gazetas et al. (1985) and Pais and Kausel (1988) solutions give
stiffness values of 1.5 and 1.03, respectively, for L/B = 100. For buried pipelines,
Matsubara and Hoshiya (2000) give a value of 2.7 for λ = 10 whereas Selvadurai (1985)
gives a value of 2.16 for L/B = 100. For a pipe on the seabed, the DNV (2019) guidelines
give a value of 1.46 nominally for L/B = 10, which matches closely with the Poulos and
Davis (1974) value. For nonhomogeneous soil with zero modulus at the soil surface but
a Poisson’s ratio of 0.5 rather than 0.4, the Gibson (1967) solution is equal to 2 for any
arbitrarily shaped rigid footing. The half-pile solution (Baguelin et al. 1977) gives a
slightly higher value of 2.25 for a pipe embedded by half a diameter in homogeneous soil.
The numerical solution from the finite-element analyses presented later is also plotted
vertical stiffness is found to vary from 3.0 (homogeneous soil) down to 2.62
(nonhomogeneous soil with G = mz), although both values are for a finite-layer depth of
81
Structural analysis of submarine pipelines
Figure 3-12 shows the variation in nondimensional lateral elastic stiffness as a function
of aspect ratio. For circular and square foundations, Poulos and Davis (1974) and Gazetas
and Tassoulas (1987) propose similar values of approximately 2.5 for L/B = 1. For a
rectangular foundation, the stiffnesses defined by Poulos and Davis (1974), Gazetas and
Tassoulas (1987), and Pais and Kausel (1988) shows similar patterns, reducing and
diverging when extrapolated to increasing aspect ratio. However, the stiffness values
from Gazetas and Tassoulas (1987) are higher than other values for high L/B ratios.
Gazetas and Tassoulas (1987) also propose a value of 1.4 for a strip footing between their
values for surface and embedded foundations with L/B = 100. For a pipe on the seabed,
DNV (2019) guidelines yield a value of 1.06 nominally for L/B = 10, as shown in Figure
3-12. The half-pile solution (Baguelin et al. 1977) gives a value of 2.25 for a half-pile in
homogeneous soil, plotted for comparison at L/B =100. Nondimensional lateral stiffness
for a half-embedded pipe (w/D= 0.5) is found to vary from 2.22 (homogeneous soil) down
to 2.36 (nonhomogeneous soil with G = mz) from the finite-element analyses, spanning
Figure 3-13 shows the variation in nondimensional axial elastic stiffness with aspect ratio.
For a rectangular foundation, the stiffnesses proposed by Gazetas and Tassoulas (1987)
and Pais and Kausel (1988) show similar patterns, reducing with increasing aspect ratio.
The elastic stiffness for a long-embedded foundation is higher than for a surface
foundation by about 30%, reducing with increasing L/B. The half-pile analogy from
Randolph and Wroth (1978) gives a stiffness of 0.78 for homogeneous soil, which is
slightly lower than the value of 1 from the present analytical solution for
nonhomogeneous soil. The numerical solutions for a half-embedded pipe are plotted
82
Structural analysis of submarine pipelines
forL/B = 100 in the figure and range from 1.06 for homogeneous soil down to 1.0 for G
= mz.
A series of small strain finite element (SSFE) analyses were performed in order to
validate the analytical solution for the elastic pipe-soil axial stiffness, and also derive
ratios between the vertical, lateral and axial stiffnesses. The analyses assumed the pipe
was “wished in place” at various embedment ratios (i.e. neglecting the pipe installation
process). The analytical and numerical results are compared and a parametric study
undertaken to explore the effects of soil parameters and embedment ratio on the elastic
stiffness. Simple expressions for the axial elastic stiffness for homogeneous and non-
A pipe diameter of D = 0.4 m was adopted, with both vertical and horizontal boundaries
placed at 2.5 m (10D) away from the centreline of the embedded portion of the pipe. This
was found to be sufficient to minimise boundary effects. The base of the soil model was
fixed in all three coordinate directions and the sides of the model were fixed in the
horizontal direction (Randolph et al., 2012). The portion of the pipe surface embedded in
the soil was treated as fully rough and bonded to the soil surface using the ‘tie constraint’,
which is computationally efficient and less prone to convergence problem than other
The rigid pipe section was also meshed (even though forced to displace as a rigid body)
with the number of elements around the periphery being equal to the corresponding
number of elements in the soil around the pipe (as necessary using tie constraints). A very
high Young’s modulus of 1011 MPa was adopted for the pipe to ensure essentially rigid
response. First order full integration hexahedral (C3D8) elements (as shown in Figure
83
Structural analysis of submarine pipelines
3-14) were used for both pipe and soil domains. Corresponding nodes on the two vertical
walls were tied together using the multi-point constraints available in ABAQUS, forcing
each node to have exactly the same degree of freedom. In this way, an infinitely long (but
rigid) pipe was modelled in order to evaluate the pipe-soil stiffness for any typical section
homogeneous soil (e.g. Merifield et al., 2008; 2009) has mainly been restricted to vertical
and horizontal movement of the pipe, and axial soil-pipe interaction (especially in respect
For these analyses, the embedded portion of the pipe was modelled as a rigid body, with
load and displacement related to a single load reference point (LRP) at the centre of the
pipe section. As the stresses in the pipe are not relevant for this particular study,
and embedment ratios (w/D) of 0.1, 0.2, 0.3, 0.4 and 0.5 were modelled for the analyses..
The shear modulus of the soil was assumed to increase linearly with depth, expressed as
G(z) = G(0)+ mz. In the limiting case where the shear modulus was zero at the surface,
G(0) = 0, the expression for shear modulus becomes G(z) = mz. Poisson’s ratio was taken
as = 0.4.
considered and V-H swipes were performed to establish the yield envelope in vertical-
horizontal space. Homogeneous soil with uniform shear strength of (su = 5 kPa) was
considered. Reaction forces were normalised by the shear strength of the soil and the
diameter of the pipe. Resulting envelopes were compared with those of Merifield et al.
84
Structural analysis of submarine pipelines
(2008). Exact matches were observed for all the embedment cases. The following Figure
3-15 shows the yield envelopes in V-H space for different initial embedment for the fully
bonded case.
Thereafter, numerical analyses were carried out to investigate the elastic response of
mD' (3-60)
=
G (0) + mD'
For homogeneous soil, with m = 0, = and for the limiting case of stiffness
proportional to depth, with G(0) = 0, = The value of was varied from 0 to 1 with
an increment of 0.2, by keeping m constant (at 333 kPa/m) and varying G(0) for different
values of embedment for a pipe with G/su = 167. The elastic stiffnesses were then
normalised, by the value of G at a depth D′ (=G(D′) to eliminate the effect of the arbitrary
It was found that the results of the parametric study could be fitted for both homogeneous
and non-homogeneous soils, expressing the non-dimensionless stiffness as (e.g. for the
axial direction)
b (3-61)
kx 2w
Kx = = a
G ( D ) D
The values of a and b were found using a least squares method from the numerical results
and expressed as
85
Structural analysis of submarine pipelines
a = 1 + 0.5 (1 − ) (3-62)
b = 0.15 + 0.3 (1 − )
Figure 3-16 shows the normalised values of axial stiffnesses for different values of
along with the proposed fit. The analytical solution for = 1 is also plotted on the same
graph. The comparison between analytical and proposed fit improves with increasing
embedment. This is perhaps due to the diminishing difference between the true pipe
geometry and the equivalent pipe (Figure 3-9).Also, the proposed fit shows good
Figure 3-17 shows the comparison between the actual and predicted axial stiffnesses for
In the following section finite element values of vertical and lateral stiffness, using the
same model and relationships as for the axial direction, are presented and compared with
The analytical expressions of lateral and vertical elastic stiffnesses were derived by
various researchers as discussed earlier in the chapter. However, the relationship between
axial and vertical or axial and lateral elastic stiffnesses was not addressed in the literature.
Therefore, this section aims to derive simple ratios of vertical versus axial elastic
stiffnesses and lateral versus axial elastic stiffnesses for the same group of input
86
Structural analysis of submarine pipelines
homogeneous soil ( = 0), while Figure 3-19 shows corresponding results for the extreme
Figure 3-20 illustrates the relationship between the vertical and axial and elastic
Kx (3-63)
= 0.34
Kv
Similarly, Figure 3-21 shows the relationship between the lateral versus the axial elastic
stiffness where:
Kx (3-64)
= 0.46
Ky
Just as for the axial case, the vertical and lateral elastic stiffnesses obtained numerically
may be estimated using equation (3-61) but with adjusted values of the parameter a to
reflect the above ratios. For completeness the values of a and b are summarised in Table
Figure 3-22 and Figure 3-23 show the comparison between the actual and predicted
vertical and lateral stiffnesses respectively for = to covering the full range of
agreement with the actual stiffness for both vertical and lateral loading, although with
The higher value of vertical stiffness for homogeneous soil modulus may be attributed,
at least partly, to the limited layer depth modelled numerically. However, other studies of
strip foundation vertical stiffness suggest that this limited mesh extent is likely to raise
the stiffness by only ~ 15% (Sloan & Randolph 1982, Bell et al. 1992, Yan et al. 2011).
Therefore, the discrepancy evident in Figure 3-11 at L/B = 100 between the FE results
and the other relationships is because of their inaccuracy when extrapolated to higher L/B.
87
Structural analysis of submarine pipelines
For the lateral and axial cases, the deformation is restricted to a shallower region, and any
The corresponding ratio, Ky/Kv, therefore averages 0.75. The largest outliers from each
set of ratios are those for the extreme non-homogeneous condition with G0 = 0, hence
= 1, for which the averages for the data at different embedment are about 10 % higher
than the overall averages indicated in Figure 3-24. Given the challenges of estimating
suitable ‘elastic’ shear modulus values for soil, the spread of stiffness ratios is acceptably
narrow.
3.10 CONCLUSIONS
Values of elastic stiffness for on-bottom pipelines have not been well defined in the
literature. However, based on previous studies for shallow foundations and present
analytical and numerical analyses approximate relationships for the stiffnesses and their
relative magnitudes have been proposed. An analytical model was presented to predict
the axial stiffness, validated through numerical analysis. The proposed form of
directions. These were found to be rather higher than values obtained by extrapolating
Stiffness ratios of axial to lateral (Kx/ Ky) and axial to vertical (Kx/Kv) were evaluated and
Average ratios of 0.34 and 0.46 were obtained (giving an average of 0.75 for Ky/Kv).
Generalised expressions were derived to evaluate the axial, lateral and vertical stiffnesses
directly as functions of w/D and . Overall, a robust systematic approach for selecting the
88
Structural analysis of submarine pipelines
stiffness (Kx), which is close to unity, and then deduce corresponding values for the
vertical (Kv) and lateral (Ky) stiffnesses. These results allow consistent and rigorous
This chapter lays the basic foundation for the following chapters on axial slide pipeline
interaction and then buckling analyses due to axial loading. Once the axial elastic stiffness
is known, then the ultimate axile resistance can be evaluated and axial springs can be
formulated and implemented into FE analyses of the pipeline for design calculations. The
outcomes of this chapter will be reflected in the analyses of the following two chapters.
89
Structural analysis of submarine pipelines
Rectangular K V0 (1 − ) L
0.75
𝑙 0.75 Pais and Kausel
= 3.1 + 1.6 𝑘𝑣 3.1 ( ) + 1.6 𝑏
GB B 𝐾𝑣 = = 𝑏
1988
Kv 𝑙𝐺 4(1 − 𝜈) 𝑙
90
Structural analysis of submarine pipelines
𝑘𝑣 0.8
2 LG 𝐾𝑣 = =
K z = S z 𝑙𝐺 1−𝜈
1 −
0.75
A A
S z = 0.73 + 1.54 2 for 2 0.02 Dobry and
4L 4L
Gazetas 1986
A
S z = 0.8 for 2 0.02
4L
kv 2 b
Kv = = z
lG 1 − l
P(1 − 2 )
z =
Kv z BL E Poulos and Davis
1974
where , P = total vertical load
B, L = rectangular dimensions
β z = factor dependent of L/B and
1< β z<1.4
91
Structural analysis of submarine pipelines
4GB Kv =
kv
=
2
K z = lG 1 −
1 − Dobry and
where B is equal to the radius of the circle Gazetas 1986
Circular &
2 p av a
P k 2
z = (1 − ) Kv = = v =
2 E Ga w Gb 1 − Poulos and Davis
1974
kv 0.8
Kv = =
lG 1 −
Strip K z 0.8G Dobry and
Kv =
2L 1 − Gazetas 1986
92
Structural analysis of submarine pipelines
0.88G
K v =
1 − k v 0.88
Kv = =
lG 1 −
where KV =V/Δv and V is the incremental vertical force
Pipe on seabed DNV(2019)
between pipe and soil per unit length of pipe, and Δv is the
associated incremental vertical displacement of the pipe.
1 D A
2/3
k v _ emb
K emb =
2 LGS z 4 Ab
1 + 0.19 s K v _ emb = =
1 + 1 + Gazetas , Dobry
1 − 21 B 3 4 L2 Ab lG
4 b w & Tassoulas
2/3
S 1 2w
= z +
21 b + +
3 l
where, D = embedment depth 1 1 1 0.19
1 − 2b 1985
Rectangular As = area of the sidewall soil interface = D x 2L
w is the embedment depth
with Ab = area of the base = 2L x 2B
embedment
0.25 0.8 0.25 0.8
KVs = KVo 1.0 + 0.25 + ( E / B ) K v _ emb = K V 1.0 + 0.25 + (2w / b ) Pais and Kausel
L/ B l /b
1988
where E/B is the embedment ratio w is the embedment depth
93
Structural analysis of submarine pipelines
v = q / 2m Kv =
kv
*
4
D
Kv lG D
Arbitrary shape Where q = Load / area Gibson 1967
G(z) = mz
94
Structural analysis of submarine pipelines
b
0.41 −
K (2 − ) K (2 − )
0 0
L ky kx l Pais and Kausel
y
= x
+ 0.8 − 1 Ky = = +
GB GB B lG lG (2 − ) 1988
(2 − ) K y , sur A
0..85
b
0.85
Ky Q(1 − 2 )
h =
X BL E
Poulos and
where Q = total horizontal load
B, L = rectangular dimensions
𝑘𝑦 2 𝑏 Davis 1974
β z = factor dependent of L/B and ν
𝐾𝑦 = =𝛽 √
0.692< β x<1.25
𝑙𝐺 1−𝜈 𝑙
95
Structural analysis of submarine pipelines
Gazetas and
8GR ky 4
K y = Ky = = Tassoulas 1987
2 − lG 2 −
Circular
(7 − 8 )(1 + ) p xa
x =
16(1 − ) E ky 4
Ky = = Poulos and
where lG 2 −
px = average applied pressure, P/πa2 Davis 1974
Ky
a = radius of the circle
Pipe on seabed
K L = 0.76G(1 + ) Ky =
ky
= 0.76(1 + ) DNV(2002)
lG
Where KL=FL/ΔδL, where ΔFL is the incremental
horizontal force between pipe and soil per unit length
96
Structural analysis of submarine pipelines
D h Aw
0.40
97
Structural analysis of submarine pipelines
K Hx
0
(2 − ) L
0.65
= 6.8 + 2.4 kx 1 b
0.35
b
GB B Kx = = + 1.2
3.4 Pais and Kausel 1988
lG (2 − ) l l
Rectangular
0.21 B
K x , sur = K y , sur − GL1 − kx k y 0.21 b Gazetas and Tassoulas
0.75 − L Kx = = − 1 −
lG lG 2(0.75 − ) l 1987
Kx
Rectangular K x = K x , sur I tre I wall =
0.40
with embedment D h Aw K x _ emb =
= K x , sur 1 + 0.15 B 1 + 0.52 B L2
Gazetas and Tassoulas
2w w 2 b
0.40
= K x 1 + 0.15 1 + 0.52 2
Where b b l 1987
D = embedment depth
h = side wall’s centre of gravity from top of the
ground
Aw = area of the sidewall soil interface = D x 2L
98
Structural analysis of submarine pipelines
K Hs = K Ho 1.0 + 0.33 +
1.34
(E / B )0.8
K y _ emb = K y 1.0 + 0.33 +
1.34
(2w / b )0.8 Pais and Kausel 1988
1+ L / B 1+ l / b
D 0
Kx =
2 X 2 G x 4
0 =
Pipeline as half G D To compare this with rectangular foundation
Randolph and Wroth
pile Where 4 is replaced with,
D = Pile diameter (1978)
ζ = side wall’s centre of gravity from top of the l
ground ln 5 + 5(1 − ) (Horikoshi & Randolph
Aw = area of the sidewall soil interface = D x 2L b
1999)
99
Structural analysis of submarine pipelines
F
kx = = mD' = mD sin D '
x
Where F = axial resistance per unit length 2 F D w analytical
δx= axial displacement Kx = = 2 sin D = 2 = 2 sin arcCos1 − 2
Kx Pipe on seabed mD x D D solution
Dʹ = embedded diameter
D = diameter of the pipe or
θD’ = embedment angle 1 F
Kx = =1
mD sin D x
100
Structural analysis of submarine pipelines
101
Structural analysis of submarine pipelines
kx l l w
Kx = = sin D ' = sin arcCos1 − 2
G ( z )b z z D
D w
Kx = sin arcCos1 − 2
z D
102
Structural analysis of submarine pipelines
103
Structural analysis of submarine pipelines
Reference:
Bell, R. W., Houlsby, G. T., and Burd, H. J. (1992). Finite element analysis of
Dassault Systèmes (2007) Abaqus analysis users’ manual, Simula Corp, Providence, RI,
USA.
pipelines.
135.
Gazetas, G., Dobry, R. & Tassoulas, J. L. (1985). Vertical response of arbitrarily shaped
440-457.
Hodder, M. S. & Cassidy, M. J. (2010). A plasticity model for predicting the vertical
and lateral behaviour of pipelines in clay soils. Géotechnique 60, No. 4, 247-
263.
104
Structural analysis of submarine pipelines
Horikoshi, K., & Randolph, M. F (1999). Estimation of overall settlement of piles rafts.
Soils and Foundations. Japanese Geotechnical Society. Vol. 39, No. 2, April, 59-
68.
Matsubara, K., & Hoshiya, M. (2000). Soil spring constants of buried pipelines for
Osman, A.S. and Randolph, M.F. (2012). An analytical solution for the consolidation
227.
Poulos, H.G. & Davis, E.H. (1974).Elastic solutions for soil and rock mechanics. John
Randolph, M. F., White, D. J. and Yan, Y. (2012). Modelling the axial soil resistance on
105
Structural analysis of submarine pipelines
piles. J. Geotech. Engng Div. Am. Soc. Civ. Engrs 104, GT 12, 1465-1488
Sloan, S.W. and Randolph, M. F. (1982). Numerical prediction of collapse loads using
Tian, Y., Cassidy, M. J., & Gaudin, C. (2010). Advancing pipe-soil interaction models
Yan Y. White D.J. & Randolph M.F. (2011). Penetration resistance and stiffness factors
Zhang, Y., Bienen, B., Cassidy, M. J., & Gourvenec, S. M. (2011). The undrained
106
Structural analysis of submarine pipelines
Figures:
107
Structural analysis of submarine pipelines
Figure 3-6: Dimensions of embedded footing for vertical stiffness by Gazetas, Dobry,
and Tassoulas, 1985
108
Structural analysis of submarine pipelines
Figure 3-7: Variation of shear modulus in linear and nonlinear soil, Gibson 1967
109
Structural analysis of submarine pipelines
110
Structural analysis of submarine pipelines
111
Structural analysis of submarine pipelines
Figure 3-14: Finite element model used to estimate the axial stiffness
112
Structural analysis of submarine pipelines
2.5
1.5
0.5
0
0 1 2 3 4 5 6 7
Normalised vertical resistance, V/suD
113
Structural analysis of submarine pipelines
Figure 3-17: Comparison between actual and predicted axial elastic stiffnesses
Figure 3-18: Comparison of stiffnesses in all three directions for homogeneous soil
114
Structural analysis of submarine pipelines
Figure 3-19: Comparison of stiffnesses in all three directions for non-homogeneous soil
115
Structural analysis of submarine pipelines
Figure 3-22: Comparison between actual and predicted vertical elastic stiffnesses
116
Structural analysis of submarine pipelines
Figure 3-23: Comparison between actual and predicted lateral elastic stiffnesses
0.60
=0
Ratios of axial to lateral and vertical stiffness
0.2
0.55
Kx / Ky = 0.46 0.4
0.6
0.50 0.8
1
0.45
0.40
0.35
0.30
Kx / Kv = 0.34
0.25
0.20
0 0.1 0.2 0.3 0.4 0.5 0.6
Embedment ratio, w/D
Figure 3-24: Ratios of axial to lateral and vertical stiffnesses
117
Structural analysis of submarine pipelines
118
Structural analysis of submarine pipelines
ANALYTICAL SOLUTION OF
4.1 INTRODUCTION
The offshore oil and gas industry commonly operates in deepwater, beyond the
submarine landslides, mud and volcanoes, seismicity, shallow gas and gas hydrates
(Kvalstad et al. 2001). One of the most significant geohazards on the continental slope is
the threat of submarine landslides, which typically originate from the shelf-break but may
run out several kilometres into development zones or across pipeline routes. Submarine
pipelines are the safest and most economical mode of transporting hydrocarbon products
from deepwater to the shore. Hence, such pipelines are more exposed to impact risk from
submarine slides than other seabed infrastructure. Therefore, geohazard assessment, and
in particular quantifying the response of pipelines to impact from submarine debris flow,
This chapter focuses on developing an analytical solution for axial submarine slide
pipeline interaction. The submarine slide was considered as a moving block of constant
shear strength along the longitudinal direction of the pipeline, hence developing
compressive forces in the downstream pipeline and tensile forces in the upstream of the
pipeline. The axial loading on the pipeline from the submarine slide movement can be
assessed using two different approaches: (1) geotechnical and (2) fluid dynamics. The
latter approach was not included in this thesis. Here, a geotechnical approach will be
The chapter describes the derivation of a solution for the axial load and deformation of
an elastic pipeline laid on a seabed, assuming a linear elastic – perfectly plastic t-z
response (outside the slide zone). The effect of the submarined axial slide is represented
by a uniform axial load (force per unit length) over a defined length of the pipeline. In the
analysis, buckling of the on-bottom pipeline due to the axial slide loading downstream of
the slide has been neglected. However, the peak axial load experienced by a given
pipeline due to different slide loading cases and t-z responses were compared with the
To investigate the axial slide impact on the pipeline the frictional force developed on the
pipeline due to the slide movement and seabed needs to be calculated. The following
sections start by reviewing previous work on load transfer from the slide (active loading)
and frictional resistance of the seabed (passive resistance). Thereafter, a proposed method
to estimate the net resistance on the pipe due to slide movement along the length (parallel
to the pipeline axis) is identified along with the boundary conditions. However, the case
where the pipe is inside the mudflow with a moving soil below the pipe invert has been
chapter for loading on a submarine pipeline due to slide movement parallel and normal
to the pipe. Buckling and other forms of failure of the pipe is considered later in the thesis.
The active force on any infrastructure on the seabed may be expressed in terms of the
shear strength of the sliding soil mass. A typical expression is given as (Zakeri, 2009)
120
Structural analysis of submarine pipelines
where su_slide is the undrained soil shear strength and A is the projected area in the direction
of the flow. The k-parameter has either been determined experimentally or based on field
Previous researchers (Demars, 1978; Swanson and Jones, 1982; Swanson and Jones, 1982;
Georgiadis, 1991; Bruschi et al., 2006) used a geotechnical approach to estimate the
normal force per unit length of the pipeline due to slide loading. They followed the
Georgiadis (1991) expressed the resistance factor as a power-law function of the velocity
as
n
U (4-2)
k = k r
U r
viscosity parameter.
Marti (1976) and Schapery and Dunlop (1978) followed a similar approach and expressed
n
U (4-3)
Fn = g n Ds u _ slide
D
where gn is the empirical dimensional coefficients (units of s-n). The value of gn was
derived as
g n = 11.42(125.9n ) (4-4)
n
For a buried pipeline Audibert et al. (1984), Summers and Nyman (1985) and Sweeney
et al. (2004) expressed the force on the pipe (in clays) following equation.However, for
Fn = Z C N q D (4-5)
where γʹ, is the soil effective unit weight, ZC is the embedment depth from the ground
D
Z C (1 + K 0 ) tan
(4-6)
Fa =
2
where K0, is the coefficient of earth pressure at rest and δ, is the interface friction angle.
Calvetti et al. (2004) slightly modified the previous expression (equation (4-6) ) by
studying the pipe-soil interaction for buried submarine pipeline in a sand rich slope and
Fn = ( − i w ) Z C N q ( , ) D (4-7)
For flow normal to the pipeline a hybrid approach, combining ‘geotechncial’ and ‘fluid
drag’ componets of resistance, was given by Randolph and White (2012) for undrained
conditions as
1 (4-8)
Fn = C d v n2 D + N p su ,nom D
2
Similarly for flow parallel to the pipeline (axial force) was given by
122
Structural analysis of submarine pipelines
Fa = f a su ,nomD (4-9)
In these equations Cd is the drag coefficient, ρ the density of the flowing material, vn the
velocity component, Np the bearing factor, and fa friction coefficient parallel to the pipe.
(4-10)
v/D
su ,nom = su ,ref .
ref
where su,nom is a reference shear strength at a strain rate of γref and v is the flow velocity.
The value of fa lies between 1.2 to 1.4 for typical values of β between 0.05 and 0.15.
For slide loading of onshore buried pipelines, O'Rourke et al. (1995) established five
idealised patterns for the permanent ground deformation (PGD) shown in Figure 4-1.
Their analysis showed that, for longitudinal PGD, a block pattern results in the largest
strain in an elastic pipe. As such for the rest of this chapter the slide loading will be
there are many reports and publications available on this topic. A summary of current
research and practice in this area is given by publications related to the SAFEBUCK JIP
An alternative approach used to estimate the axial resistance, referred to as the total stress
(alpha) method, is comparable to the equivalent technique to estimate the axial active pipe
force.
For buried pipes, a similar expression was proposed by American Lifeline Alliance (2001)
123
Structural analysis of submarine pipelines
where α is a ‘adhesion’ factor which is influenced by the roughness of the surface (with
α being the relative strength of pipe-soil and soil-soil shearing), and su_seabed is the
For embedded pipelines Cathie et al. (2005) and Oliphant and Maconochie, (2007)
modified the expression slightly by introducing an embedment effect, with the axial
where, θD′ is the contact angle as defined in the previous chapter (see also Figure 4-3)
The alpha method for estimating the passive frictional resistance of the seabed can be
extended by integrating the seabed shear strength, su_seabed around the embedded portion
of the pipe. For conditions of zero mudline strength intercept (i.e. seabed shear strength
/2 (4-13)
Fpassive = 2RD D cos d = DD D
0
The net axial force on a pipeline depends on the relative magnitudes of the shear strengths
of slide and seabed, and also on the pipeline embedment (expressed through D'). As the
slide material flows down the slope, the shear strength is likely to reduce significantly,
for example from 3 to 30 kPa (in situ) to 1 kPa or less as shown in Figure 4-2. By contrast
the seabed shear strength may be expected to lie in a similar range to that of the original
(in situ) submarine slide material. Modelling of the shear strength variation during the
debris flow is extremely challenging and not considered in this study. Only the effect of
the actual shear strength ratios on the pipeline structures is discussed in the next section.
124
Structural analysis of submarine pipelines
Following the alpha method the axial active force per unit length of the pipe (top), Fslide
Figure 4-3 shows the schematics of the active and passive loading on the pipe.
capture the effect of the ratio of su_slide / su_seabed on the ratio of the active and passive
The above expression is shown as a design chart in Figure 4-4 for embedment ratios 0 <
w/D ≤ 0.5 and with various values of undrained shear strength ratios within, 0.1 < su_seabed
/ su_slide ≤ 10. From the chart it is clear that, at high su_seabed / su_slide ratios (i.e. a weak slide
relative to the seabed soil), the required embedment ratio for no net load to be transferred
to the pipeline itself is small. Also, there can be no net load on the pipeline for su_seabed /
su_slide ≤ 1.
systems distributed at regular intervals along the length of the pipeline. A simplified
method of estimating axial soil resistance, Fpassive was proposed in BS 8010 (1993) as
Fpassive = N = W (4-16)
where W′, is the submerged weight of the pipe calculated as Dtpipe (pipe is density of
125
Structural analysis of submarine pipelines
enhancement factor, ζ, to account for wedging around the curved surface of the pipe as,
2 sin D (4-18)
=
D + (sin D cos D )
The factor increases to an upper limit of 1.27 for pipeline embedment of 0.5D, as shown
As discussed by White and Randolph (2007) the effective stress approach (equation
(4-17)) is preferable to the total stress method, mainly because the seabed shear strength
will change due to consolidation under the pipeline weight. Therefore, the effective stress
approach is followed in this chapter for quantifying the passive resistance, Fpassive.
Using the effective stress approach for seabed friction (equation (4-17)) and alpha method
for estimating slide loading (equation (4-14)) a dimensionless generalised expression can
be formed to capture the effect of the ratio of su_slide / su_seabed on the ratio of the active
F passive W W 2 sin D 1
(4-19)
= =
Fslide su _ slide D( − D ) su _ slide D ( D + sin D cos D ) ( − D )
Similar to equation (4-15) the above expression is shown as a design chart in Figure 4-6
for embedment ratios 0 < w/D ≤ 0.5 and with various values of pipe weight versus slide
shear strength ratios within, 0.1 < μWʹ/D su_slide ≤ 10. From the chart it is clear that, at
high μWʹ/D su_slide ratios (i.e. a weak slide relative to the seabed soil), the required
embedment ratio for no net load to be transferred to the pipeline itself is small. Also, there
126
Structural analysis of submarine pipelines
The submarine slide-pipeline-seabed interaction problem may be divided into three parts,
active slide zone, passive plastic zone and elastic zone as shown in Figure 4-7. Key loads,
displacements and frictional resistances for each zone are indicated in the schematic. The
response in each zone is solved analytically for the relevant boundary conditions in the
following sections.
The perfectly straight pipe is defined by diameter, D, wall thickness, t submerged unit
weight, Wʹ, and Young’s modulus, E, from which the axial rigidity EA can be calculated.
The slide was defined as a block pattern with length, Lslide ; to analyse the problem only
half slide length, LAB is considered here, corresponding to the ‘active slide zone’. The
length of the ‘passive plastic zone’ is LBC, beyond which point (C onwards) it is defined
as ‘elastic zone’ (Figure 4-7). The displacement at the centre of the slide, A, is uA, at the
interface of ‘active slide zone’ and ‘passive slide zone’, B, is uB and the interface of
‘passive slide zone’ and elastic zone, C is uC. The axial load generated within the pipeline
due to the slide movement along the length is defined as P. At the centre of the slide at A,
it is PA, at the interface of ‘active slide zone’ and ‘passive slide zone’, B, is PB and the
interface of ‘passive slide zone’ and elastic zone, C is PC. It is noteworthy to mention that
this idealisation assumes symmetric geometry of the slide and pipeline however, the
When a submarine landslide occurs and passes across a partially embedded pipeline with
the upper portion and at the invert of the pipe are different due to variations in the physical
properties of the slide material (debris flow material) and seabed (natural fine-grained
sediment). The net (resultant) force per unit length of the pipeline is given by
127
Structural analysis of submarine pipelines
input properties. The maximum axial load, PB, may be normalised by the axial elastic
stiffness of the pipe, EA, and presented as compressive strain, = P/EA; the axial
displacement, u, may be normalised by the slide length, Lslide, as u/Lslide. These normalised
output parameters may then be expressed in terms of normalised input parameters, i.e. the
driving force, FnetLslide/EA; resisting force per diameter in passive zone, FpassiveLslide/EA;
All the input and output parameters are defined in the previous section and tabulated in
Table 4-1.
128
Structural analysis of submarine pipelines
k x L2slide (4-26)
a3 =
EA
The following sections discuss the analytical solution of load, displacement and frictional
resistances on a pipe due to axial slide loading. Throughout the chapter it is assumed that
the passive resistance (equation (4-17)) continues to act (unchanged) within the active
The gradient kx of the initial slope of the ‘t-z’ response of a pipeline embedded within the
seabed was considered in chapter 3, focusing in particular on conditions where the seabed
closed form solution for kx was developed by considering the deformation of a vertical
(or radial) circular slice comprising a short element of pipe and the surrounding soil.
F (4-27)
kx = = mD
u
where m is the gradient of the shear modulus and D′ is the contact width, which can be
The axial load generated in the pipe due to the presence of frictional resistance of the
seabed is
dP (4-28)
= −F
dx
The compressive strain, εx, of the pipeline (assumed elastic) is written in terms of the load,
129
Structural analysis of submarine pipelines
du − P (4-29)
=− =
dx EA
Differentiating equation (4-29) and using equations (4-27) and equation (4-28) gives
d 2u F k (4-30)
= = x u
dx EA EA
dimensions m-1 . The constants C1 and C2 may be found from the boundary conditions, u
In summary, the axial displacement of the pipe can be expressed in terms of the
u( x) = u C e −x (4-32)
The total load in the pipe is expressed by substituting the value of u(x) from equation
dP (4-33)
= − F = −k xu ( x) = −k xuC e −x
dx
Integrating the above equation from x = 0 at the passive-plastic and elastic interface to
130
Structural analysis of submarine pipelines
k x uC (4-35)
PC = = k x EAuC
kx
EA
(4-36)
PC k x L2slide uC u
= = a3 C
EA EA Lslide Lslide
Since point C (Figure 4-7) represents the interface between passive plastic and elastic
zones, using equation (4-27) the elastic slip uC-slip may be expressed as
Therefore, using equation (4-35) and (4-37) the elastic load in the long pipe is given as
PC ,max a (4-38)
EA hence = 2
PC ,max = F passive
kx EA a3
In general, there will be a passive plastic zone between the active slide zone and the elastic
zone, where slip occurs between the seabed and the pipe and the resistance force per unit
length is Fpassive. The governing equations of the plastic zone are similar to those for the
dP (4-39)
= − Fpassive
dx
C C (4-40)
dP = − Fpassive dx
B B
131
Structural analysis of submarine pipelines
The axial load in the passive slide zone is expressed by integrating equation (4-40) as
PB − PC L P P 1 (4-42)
LBC = hence BC = B − C
Fpassive Lslide EA EA a2
When the passive zone LBC = 0, point B coincides with point C leading to PB = PC. In
du P( x ) (4-43)
=−
dx EA
where P is the varying load between B and C. Integrating the above equation gives the
C
1
C (4-44)
du = −
B
EA B
P( x)dx
where B is the boundary point at the active-passive intersection, C is the boundary point
at the elastic zone and P(x) is the linearly varying load between B and C.
Now integrating equation (4-44) gives the compression in the passive-plastic zone:
(PB + PC ) (4-45)
u B − uC = LBC
2 EA
can be expressed as
132
Structural analysis of submarine pipelines
𝑢𝐵 𝑢𝐶 2
1 (𝑃𝐵 −𝑃𝐶2 ) 𝐸𝐴
= +
𝐿𝑠𝑙𝑖𝑑𝑒 𝐿𝑠𝑙𝑖𝑑𝑒 2 𝐸𝐴2 𝐹𝑝𝑎𝑠𝑠𝑖𝑣𝑒 𝐿𝑠𝑙𝑖𝑑𝑒
For a small active slide load (or strong passive resistance), the elastic limit of uC-slip may
In the active zone the interaction between the pipe and the soil is assumed to be plastic.
The displacement is taken as uA at the centre (x = 0) of the slide from the symmetry. The
dP (4-47)
= − Fnet
dx
where Fnet is given by equation (4-20). Integrating equation (4-47) within the active
B B (4-48)
dP = − Fnet dx
A A
Lslide (4-49)
PA − PB = − Fnet LAB = − Fnet
2
where LAB is the half slide length, Lslide/2. Assuming the axial load at the centre of the
133
Structural analysis of submarine pipelines
Note that PB represents the largest axial force generated in the pipeline, and hence the
Similar to equation (4-43) the axial deformation can be calculated using a compatibility
du P( x ) (4-51)
= − x = −
dx EA
Integrating this, allowing for the linear variation of P between PA (= 0) and PB gives
A
1
A (4-52)
B du = − EA B P( x)dx
Therefore, using equation (4-50) the axial deformation within the slide zone is
PB Lslide F L2 uA u a (4-53)
u A − uB = = net slide hence = B + 1
4 EA 8 EA Lslide Lslide 8
For convenience the main expressions are summarized here in non-dimensional form.
PA P a P a a (4-54)
= 0; B = 1 ; C = Min 1 , 2
EA EA 2 EA 2 a 3
LBC a 1 (4-55)
= Max 0 , 1 −
Lslide 2a 2 a3
134
Structural analysis of submarine pipelines
uA a a 1 a12 1 a2
= 1 + Max 1 , −
Lslide 8 a 8 a2 2 a3
3
uB a 1 a12 1 a2
= Max 1 , − (4-56)
Lslide a 8 a2 2 a3
3
uC 1 a1 a2
= Min ,
Lslide 2 a a3
3
In the next section the effects of dimensionless parameters, i.e. driving force a1,
FnetLslide/EA ; resisting force per diameter in passive zone a2, FpassiveLslide/EA; and pipe-
displacements, u/Lslide and dimensionless passive length, LBC/Lslide of the pipeline are
design charts to show the loads and displacements for different submarine slide loading
conditions. The ranges of absolute input parameters considered relevant are tabulated in
Table 4-1, along with the corresponding ranges of the three normalised quantities a1 to
a3. The results reported in this chapter are limited to cases where the axial stress is less
than 0.003E corresponding to 630 MPa for a steel pipeline. This is just above the yield
strength of X65 grade pipe used widely for subsea pipeline projects.
For this analysis the pipeline wall thickness, D/t, has not been treated as an independent
variable, but has been kept constant at 25 (Randolph et al., 2010). To calculate the axial
rigidity of the pipe, EA, only the pipe diameter, D, was varied (between 0.1 – 1 m) keeping
Young’s modulus of the pipe, E constant at 210 GPa. The submerged weight of the pipe,
135
Structural analysis of submarine pipelines
W, was calculated as Dtpipe, with the density of the pipe material, pipe kept constant at
7900 kg/m3. A pipe embedment ratio of w/D = 0.5 and coefficient of friction, = 0.3
were adopted. The soil shear modulus gradient, m, was varied between 0.1-10 MPa/m.
Slide lengths, Lslide, ranging between 50 and 1000 m were considered along with shear
strength, su_slide , varying from 0.1 – 10 kPa. With these raw input data, relatively wide
ranges of non-dimensional input parameters (see Table 4-1) were estimated, from which
design charts for the various non-dimensional output parameters were drawn.
The non-dimensional slide force, a1, is taken as the primary influence on the non-
uC/Lslide) and non-dimensional passive length LBC/Lslide. The relationships are plotted for
different values of the secondary influence, the non-dimensional passive resistance, a2.
Then different charts are shown for different values of pipe soil elastic stiffness
parameter, a3.
0.00001 to 0.0001, for three values of pipe-soil elastic parameter, a3 = 0.01, 3 and 1000. The
force, a1, while the load at the transition from elastic to plastic soil resistance, PC/EA, is initially
proportional to a1 but then is limited (equation (4-54)). The limit is reached for most relevant
Figure 4-11, Figure 4-12, and Figure 4-13 show the corresponding variations of non-
dimensional passive length with dimensionless slide force for different values of passive
resistance, a2 = 0.00001 to 0.0001, and for pipe-soil elastic parameter, a3 = 0.01, 3 and 1000
136
Structural analysis of submarine pipelines
respectively. The length of the passive zone increases linearly with increasing slide force, with
The peak axial dimensionless axial displacement, uA/Lslide, increases with increasing slide force.
Figure 4-14, Figure 4-15, and Figure 4-16 show the variation of dimensionless displacement
with slide force for various resistances, a2 = 0.00001 to 0.0001 and for a3 = 0.01, 3 and 1000. The
axial displacement increases with increasing slide force, but at a rate that decreases with
increasing passive resistance. The elastic stiffness parameter, a3, has essentially negligible
influence.
Similarly, the normalised slide displacement, uB/Lslide, increases with increasing slide force.
Figure 4-17, Figure 4-18, and Figure 4-19 show the variation of dimensionless displacement
with slide force. Increasing passive resistance, a2, reduces the displacements in the slide and
passive zones. The values of the dimensionless peak displacements, uA/Lslide, are slightly higher (
The normalised elastic displacements increase with increasing slide loading until it reaches the
limiting value at which slip occurs, as shown in Figure 4-20, Figure 4-21and Figure 4-22. The
slip displacement increases proportionally with resistance in the passive zone, and inversely with
4.5 DISCUSSION
To compare the individual dimensionless displacements uA/Lslide, uB/Lslide and uC/Lslide, these are
plotted together on Figure 4-23 for a2 = 0.00001 and 0.0001 and a3=0.01. The peak and slide
displacements match closely for different values of a2. However, the values of elastic
displacement are much less than the peak and slide displacements. It is also useful to compare the
elastic load, Pc/EA, for different values of pipe-soil elastic stiffness, keeping a2 constant. Figure
4-24 shows the normalised elastic load as a function of the slide force for three different values
of a3 = 0.01, 3 and 1000, keeping a2 = 0.0001. For a normalised slide force of 0.006, Pc/EA =
137
Structural analysis of submarine pipelines
0.001 for a3 = 0.01 and decreases significantly for a3 = 1000 to a value of Pc/EA = 3.2x10-6. Figure
4-25 shows the effect of the pipe-soil elastic parameter on the passive length. It is evident that
the pipe-soil elastic stiffness has negligible effect on the passive length. Figure 4-26 shows a
similar observation of insignificant effect of the pipe-soil elastic stiffness on the peak
of slide load for three different values of a3 = 0.01, 3 and 1000, keeping a2 = 0.0001. For a
normalised slide load of 0.006, uc/Lslide = 0.0000001for a3 = 1000 and increases significantly to a
Normalised displacements are plotted against the normalised slide length, Lslide/D, for the whole
range of data in Figure 4-28. For an embedment ratio of w/D = 0.5 and coefficient of friction, μ
= 0.3, the normalised peak and slide displacements increase with increasing Lslide/D, whereas the
normalised elastic displacement decreases with increasing Lslide/D ratio. Normalised loads, P/EA
are also plotted against the normalised slide length, Lslide/D or the whole range of data in Figure
4-29 to show the effect of slide length. For an embedment ratio of w/D = 0.5 and coefficient of
friction, μ = 0.3, the peak loads increase with increasing Lslide/D. However, for shorter slide
lengths, elastic and peak loads merge (since no plastic zone) while for longer slides the elastic
In the next section example cases are studied to integrate all the analytical solutions together.
along the length of the pipeline. For this analysis the adopted pipeline properties were, D
= 1 m, D/t = 25, E = 210 GPa, Wʹ = 7.8 kN/m, μ = 0.3 and embedment ratio, w/D = 0.5.
Three different slide lengths with Lslide = 100, 300 and 500 m were considered, taking the
shear strength of the slide material as su_slide = 10 kPa. The corresponding non-dimensional
138
Structural analysis of submarine pipelines
The axial compressive strain, P/EA = ε (%), the dimensionless axial displacement, u/Lslide
and frictional resistance distribution along the length of the pipeline are presented here.
Figure 4-30 shows the distribution of compressive strain along the pipeline for Lslide/D =
100. In the active slide zone the axial load increases from zero (at the centre of the slide)
with a positive slope equal to the net resistance force, Fnet, until reaching a maximum
value at the edge of the slide. In the passive zone the axial load decreases linearly, thereby
decreasing the compressive strain, with a negative slope equal to the passive resistance
of seabed, Fpassive. In the elastic zone the distribution of strain decreases exponentially
The dimensionless displacement profile of the pipe along its length is shown in Figure
4-31. The displacement is maximum at the centre of the active slide zone, and decreases
in the passive slide and elastic zones. In the ‘elastic’ zone the displacements decay
The variation of normalised frictional resistance along the length of the pipe is plotted in
Figure 4-32. In the active zone applied frictional resistance force per unit length, Fnet is
very high and plotted as a positive quantity, whereas the values of passive resistance force,
Fpassive in the passive slide zone is comparatively lower and plotted as a negative quantity.
In the ‘elastic’ zone the passive resistance force of the seabed decays exponentially with
In order to show the effect of different slide lengths on the axial strain, displacement and
frictional resistance distribution along the length, a set of three cases were undertaken.
These cases examine the effects of slide loading for Lslide/D = 100, 300 and 500 (see Table
4-2). Figure 4-33, Figure 4-34 and Figure 4-35 show the variation of compressive strains,
axial displacements and frictional resistance forces along the pipe respectively. The
139
Structural analysis of submarine pipelines
forces are directly proportional to the slide length when other conditions are kept constant.
Thus, for a slide with length of 500 m the value of maximum strain is exactly five times
software ABAQUS, (Dassault Systèmes, 2007) was used to verify the analytical solution.
Figure 4-36 illustrates the schematic diagram and FE mesh used for this. The modelling
was carried out for the passive and elastic zones of the pipe, starting from B and
approximating an infinitely long straight pipe consistent with the analytical solution by
taking L∞ = 5000 m. This was the only limitation of the FEA model. The boundaries of
the pipe were left free to represent the conditions at B and at far distance from the slide.
compressive load applied to the pipe at B. The seabed was modelled with springs, with
displacement behaviour (Summers & Nyman, 1985). The pipe was modelled with
Timoshenko beam elements (ABAQUS code B21), which model linear elastic transverse
Figure 4-37 shows the distribution of normalised axial displacements for the three
different slide lengths considered in Table 4-2. The FE results are shown by solid triangles,
while the analytical solutions are represented by solid lines. The analytical solutions start
from coordinate x/D = 0, while the FE solutions start from the passive zones (point B) of
the corresponding cases. Overall analytical and FE results are in good agreement.
140
Structural analysis of submarine pipelines
compressive load is provided in the next chapter, but a brief discussion is relevant here.
In the above analytical solutions, the buckling response of the pipeline was ignored.
Analytical solutions for the critical buckling load for a straight pipe on a rigid seabed was
investigated by several researchers including Hobbs (1981), Maltby & Calladine (1995)
and Tran & Taylor (1996). Recently the SAFEBUCK JIP has introduced new analytical
solution for three modes of buckling. Also, Peek (2009) corrected the lateral buckling
solutions previously published. Therefore, in this section the critical buckling load was
estimated using the solution proposed by Peek (2009). The maximum axial load attained
in each case of slide loading was compared with the buckling initiation load to estimate
the susceptibility to buckling. The post-buckling load was estimated for a pipe of unit
weight, W, with axial and bending stiffnesses of EA and EI respectively, resting on a flat
Ppost −buckle = 8.2183076 EA EIq 2
1/ 4 (4-57)
where q = WL. The minimum axial force prior to buckling for which a buckled solution
Here only lateral buckling was considered and estimated by equation (4-58). The ratios of
peak load versus buckling load, PB/Pbuckle are plotted against the normalised slide length, Lslide/D
(for μ/μL = 1) in the same graph (Figure 4-38) for the whole range of data (Table 4-1). To see
the effect of slide loading, cases with Fnet/Fpassive = 1.07 – 40.42 are adopted here. Increasing slide
length and Fnet/Fpassive ratio increases the vulnerability of lateral buckling of the pipeline.
141
Structural analysis of submarine pipelines
4.7 CONCLUSIONS
An analytical solution was presented here for axial submarine slide loading of a straight
on-bottom pipeline. Non-dimensional expressions for axial load distribution and axial
It was shown that the non-dimensional axial loads and axial displacements depend on
three non-dimensional input parameters, i.e. the driving force in the slide zone, seabed
resisting force in passive zone, and pipe-soil stiffness. Non-dimensional design charts
were presented to show the effect of individual input parameters on axial loads and axial
displacements. The maximum axial load in the pipe is directly proportional to the slide
force, while the load at the transition from elastic to plastic soil resistance is initially
proportional to the slide force but then becomes limited. The limit is reached for most
relevant values of the slide force. The length of the passive zone was also derived non-
dimensionally and shown to increase linearly with increasing slide force, with the
gradient inversely proportional to the passive resistance. It was also shown that the axial
displacement increases with increasing slide force, but at a rate that decreases with
increasing passive resistance. The elastic stiffness parameter has essentially negligible
influence.
Numerical FE analysis was also carried out to verify the analytical model. Analytical and
numerical solutions were shown to match with close agreement. The buckling was
ignored in the analytical model. However, the existing classical theory of buckling was
linked to the output of the analytical model to show the vulnerability of the pipelines
pipelines are more susceptible to lateral buckling when impacted axially by stronger and
longer slides.
142
Structural analysis of submarine pipelines
The presented analytical framework can be used by industry practitioners to estimate the
maximum axial slide loading for a straight pipe resting on seabed. At the same time load
distribution and axial movement of the pipe along the length of the pipeline can also be
estimated.
143
Structural analysis of submarine pipelines
Lslide
parameter
pipeline, Fslide
Fnet
stiffness, kx
144
Structural analysis of submarine pipelines
parameter, a3
145
Structural analysis of submarine pipelines
References:
American Lifeline Alliance (2001). “Guideline for the Design of Buried Steel Pipe.”
ASCE, USA.
Audibert, J. M. E., & Nyman, K. J. (1977). Soil restraint against horizontal motion of
1119-1142.
Bruschi, R., Bughi, S., Spinazze,M., Torselletti, E., and Vitali, L. (2006). “Impact of
debris flows and turbidity currents on seafloor structures.” Norw. J. Geol., 86,
317-337.
BS 8010 (1993). “Code of practice for Pipelines, Part3: Pipelines subsea: design,
Dassault Systèmes Simulia. (2007). ABAQUS analysis users’ manual, Providence, RI,
USA.
Georgiadis, M. (1991). “Landslide drag forces on pipelines.” Soils Found., 31(1), 156-
161.
175-189.
Kishida, H., and Uesugi, M. (1987). “Tests of interface between sand and steel in the
simple
146
Structural analysis of submarine pipelines
Maltby, T. C., and Calladine, C. R. (1995). “An investigation into upheaval buckling of
Oliphant, J., and Macanochie, A. (2007). “The axial resistance of buried and unburied
O’Rourke, M.J., Liu, X.J., and Flores-Berrones, R. (1995). “Steel pipe wrinkling due to
451.
Peek, R. (2009). “Correction to the Infinite Mode for Lateral Buckling of Pipeline.” J.
Aeronautics, 28p.
147
Structural analysis of submarine pipelines
Randolph, M. F., Seo, D., and White, D. J. (2010). “Parametric solution of slide impact
SAFEBUCK (2008). “Safe design of pipelines with lateral buckling – design guideline,
Summers, P. B., and Nyman, D. J. (1985). “An approximate procedure for assessing the
432.
Sweeney, M., Gasca, A. H., Garcia Lopez, M., and Palmer, A. C. (2004). “Pipelines and
Proc., Int. Conf. on Terrain and Geohazard Challenges Facing Oil and Gas
Taylor, N. and Tran, V. (1996). “Experimental and theoretical studies in subsea pipeline
White, D. J., and Randolph, M. F. (2007). “Seabed characterisation and models for
Zakeri, A. (2009). “Review of state-of-art: Drag force on submarine pipelines and piles
caused by landslide or debris flow impact.” J. Offshore Mech. Arct. Eng., 131(1), ,
014001.
148
Structural analysis of submarine pipelines
Figures:
Figure 4-1: Idealized pattern for longitudinal PGD. a) Block pattern b) Ramp pattern c)
Ridge pattern d) Ramp-ridge pattern e) Asymmetric ridge pattern ( O’Rourke et al.
1995)
149
Structural analysis of submarine pipelines
Pipe position
Slope failure
Vertical scale
Figure 4-2: Schematic diagram of submarine slide stages with the shear strength
distribution along the slide length. (White et al. 2008, MERIWA project report no. M395)
Wʹ
θDʹ
θDʹ w
Seabed
Fpassive=μζWʹ
Dʹ
Figure 4-3: Schematic of frictional resistances on an embedded pipe during axial slide
loading
150
Structural analysis of submarine pipelines
2
su_seabed / su_slide = 10 - 2
1.5
Frictional resistance ratio, Fpassive / Fslide
0
0 0.1 0.2 0.3 0.4 0.5
Embedment ratio, w/D
Figure 4-4: Variation of ratio of active and passive forces with embedment ratio for
1.4
1.2
Enhancement factor, ζ
0.8
0.6
0.4
0.2
0
0 0.1 0.2 0.3 0.4 0.5
Embedment ratio, w/D
151
Structural analysis of submarine pipelines
4
μW' /Dsu_slide = 10 - 2
Frictional resistance ratio, Fpassive / Fslide
No load transfer
to the pipe
1
0
0 0.1 0.2 0.3 0.4 0.5
Embedment ratio, w/D
Figure 4-6: Variation of ratio of passive and active forces with embedment ratio for
152
Structural analysis of submarine pipelines
0.003
a3 = 0.01
0.001
0
0 0.001 0.002 0.003 0.004 0.005 0.006
Dimensionless slide force, FnetLslide/EA
Figure 4-8: Effect of dimensionless slide force on axial loads of the pipe for a3 = 0.01
0.0003
a3 = 3
Dimensionless axial load, P/EA
PB
0.0002
0.0001
Pc: a2= 0.0001 - 0.00001
0
0 0.0002 0.0004 0.0006
Dimensionless slide force, FnetLslide/EA
Figure 4-9: Effect of dimensionless slide force on axial loads of the pipe for a3 = 3
153
Structural analysis of submarine pipelines
0.0002
a3 = 1000
Dimensionless axial load, P/EA
0.00015
PB
0.0001
0
0 0.0001 0.0002 0.0003 0.0004
Dimensionless slide force, FnetLslide/EA
Figure 4-10: Effect of dimensionless slide force on axial loads of the pipe for a3 = 1000
300
a3 = 0.01
Dimensionless passive length, LBC/Lslide
200
100
0
0 0.001 0.002 0.003 0.004 0.005 0.006
Non-dimensional slide force, FnetLslide/EA
Figure 4-11: Effect of dimensionless slide force on passive zone length for a3 = 0.01
154
Structural analysis of submarine pipelines
320
a3 =3
280
200
160
120
80
40
0
0 0.001 0.002 0.003 0.004 0.005 0.006
Non-dimensional slide force, FnetLslide/EA
Figure 4-12: Effect of dimensionless slide force on passive zone length for a3 = 3
320
a3 =1000
280
Dimensionless passive length, LBC/Lslide
160
120
80
40
0
0 0.001 0.002 0.003 0.004 0.005 0.006
Non-dimensional slide force, FnetLslide/EA
Figure 4-13: Effect of dimensionless slide force on passive zone length for a3 = 1000
155
Structural analysis of submarine pipelines
0.5
a3 = 0.01
a2= 0.00001- 0.0001
Dimensionless peak displacement, uA/Lslide
0.4
0.3
0.2
0.1
0
0 0.001 0.002 0.003 0.004 0.005 0.006
Dimensionless slide force, FnetLslide/EA
Figure 4-14: Effect of dimensionless slide force on dimensionless peak displacement for
a3 = 0.01
0.5
a3 = 3
Dimensionless peak displacement, uA/Lslide
0.3
0.2
0.1
0
0 0.001 0.002 0.003 0.004 0.005 0.006
Dimensionless slide force, FnetLslide/EA
Figure 4-15: Effect of dimensionless slide force on dimensionless peak displacement for
a3 =
156
Structural analysis of submarine pipelines
0.5
a3 = 1000
0.3
0.2
0.1
0
0 0.001 0.002 0.003 0.004 0.005 0.006
Dimensionless slide force, FnetLslide/EA
Figure 4-16: Effect of dimensionless slide force on dimensionless peak displacement for
a3 = 1000
0.5
a3 = 0.01
Dimensionless slide displacement, uB/Lslide
0.3
0.2
0.1
0
0 0.001 0.002 0.003 0.004 0.005 0.006
Dimensionless slide force, FnetLslide/EA
Figure 4-17: Effect of dimensionless slide force on dimensionless slide displacement for
a3 = 0.01
157
Structural analysis of submarine pipelines
0.5
a3 = 3
Dimensionless slide displacement, uB/Lslide
0.3
0.2
0.1
0
0 0.001 0.002 0.003 0.004 0.005 0.006
Dimensionless slide force, FnetLslide/EA
Figure 4-18: Effect of dimensionless slide force on dimensionless slide displacement for
a3 = 3
0.5
a3 = 1000
Dimensionless slide displacement, uB/Lslide
0.3
0.2
0.1
0
0 0.001 0.002 0.003 0.004 0.005 0.006
Dimensionless slide force, FnetLslide/EA
Figure 4-19: Effect of dimensionless slide force on dimensionless slide displacement for
a3 = 1000
158
Structural analysis of submarine pipelines
0.012
a3 = 0.01 a2= 0.0001 - 0.00001
0.004
0
0 0.001 0.002 0.003 0.004 0.005 0.006
Dimensionless slide force, FnetLslide/EA
for a3 = 0.01
0.00004
a3 = 3
Dimensionless elastic displacement, uC/Lslide
0.00003
0.00002
0.00001
0
0 0.001 0.002 0.003 0.004 0.005 0.006
Dimensionless slide force, FnetLslide/EA
for a3 = 3
159
Structural analysis of submarine pipelines
1.2E-07
a3 = 1000
Dimensionless elastic displacement, uC/Lslide
8E-08
4E-08
0
0 0.001 0.002 0.003 0.004 0.005 0.006
Dimensionless slide force, FnetLslide/EA
for a3 = 1000
0.05
a3 = 0.01
Dimensionless displacement, u/Lslide
0.01
uC: a2 = 0.00001
0
0 0.001 0.002 0.003 0.004 0.005
Dimensionless slide force, FnetLslide/EA
160
Structural analysis of submarine pipelines
0.0015
a2 = 0.0001
PC: a3 = 0.01
0.0005
PC: a3 = 1000
PC: a3 = 3
0
0 0.001 0.002 0.003 0.004 0.005 0.006
Non-dimensional slide force, FnetLslide/EA
40
a2 = 0.0001
Dimensionless passive length, LBC/Lslide
30
10
0
0 0.001 0.002 0.003 0.004 0.005 0.006
Dimensionless slide force, FnetLslide/EA
161
Structural analysis of submarine pipelines
0.08
a2 = 0.0001
Dimensionless peak displacement, uA/Lslide
0.06
a3= 3, 1000
0.04
0.02
a3= 0.01
0
0 0.001 0.002 0.003 0.004 0.005 0.006
Dimensionless slide force, FnetLslide/EA
0.012
Dimensionless elastic axial displacement, uc/Lslide
a2 = 0.0001
a3= 0.01
0.01
0.008
0.006
0.004
0
0 0.001 0.002 0.003 0.004 0.005 0.006
Dimensionless slide force, FnetLslide/EA
162
Structural analysis of submarine pipelines
1.00E+02
μ =0.3 u_A/L_slide
u_B/L_slide
Dimensionless displacement, u/Lslide 1.00E+00 u_C/L_slide
1.00E-02
1.00E-04
1.00E-06
1.00E-08
10.00 100.00 1000.00 10000.00 100000.00
Dimensionless slide length, Lslide/D
1.00E+00
μ =0.3
P_B/EA
P_C/EA
1.00E-02
Dimensionless axial load, P/EA
1.00E-04
1.00E-06
1.00E-08
10.00 100.00 1000.00 10000.00 100000.00
Dimensionless slide length, Lslide/D
163
Structural analysis of submarine pipelines
0.003
a1 = 0.0000471
Axial compressive strain, P/EA=ε (%)
0.002
Passive slide
zone
0.001
Elastic zone
0
0 200 400 600 800 1000
Along the length, x/D
Figure 4-30: Dimensionless axial load profile for a slide of Lslide/D = 100
0.00004 a1 = 0.0000471
a2 = 0.0000149
Active slide zone
Dimensionless axial displacement,
a3 = 2.526
0.00003
0.00002
u / Lslide
Passive
slide zone
0.00001
Elastic zone
0
0 200 400 600
Along the length, x/D
164
Structural analysis of submarine pipelines
0.00006
FLslide / EA 0.00004
Passive slide
0.00002 zone
Elastic zone
0
0 200 400 600 800 1000
-0.00002
Figure 4-32: Dimensionless frictional resistance profile for a slide of Lslide/D = 100
0.015
Axial compressive strain, P/EA=ε (%)
Case A3,
0.01 Lslide/D = 500
Case A2,
0.005 Lslide/D = 300
Case A1,
Lslide/D = 100
0
0 500 1000 1500 2000
Along the length, x/D
Figure 4-33: Comparison of axial load profiles for different slide lengths
165
Structural analysis of submarine pipelines
0.00015
Dimensionless axial displacement, u / Lslide
Case A3,
Lslide/D = 500
0.0001
Case A2,
Lslide/D = 300
0
0 500 1000 1500 2000
Along the length, x/D
Figure 4-34: Comparison of axial displacement profiles for different slide lengths
0.00032
Dimensionless frictional resistance, FLslide/EA
0.00016
Case A2,
Lslide/D = 300
0.00008
0
0 500 1000 1500 2000
-0.00008
Case A1,
Lslide/D = 100
-0.00016
Figure 4-35: Comparison of frictional resistance force profiles for different slide lengths
166
Structural analysis of submarine pipelines
1m
B LBC C L∞
Axial slide loading
0.00015
Dimensionless axial displacement, u / Lslide
Case A3
0.0001
Case A2
Analytical solution ( solid lines)
0
0 200 400 600 800 1000 1200 1400 1600
Along the length, x/D
solutions
167
Structural analysis of submarine pipelines
100
Load ratio, PB/Pbuckle
Figure 4-38: Ratio of peak versus buckling load for different Fnet/Fpassive values
168
Structural analysis of submarine pipelines
169
Structural analysis of submarine pipelines
170
Structural analysis of submarine pipelines
5.1 INTRODUCTION
With the depletion of onshore and shallow offshore hydrocarbon reserves, the present oil
and gas industry is moving into the deeper sea. Deep sea developments generally include
production facilities. Deepwater pipelines are laid on the seabed and are usually left
unburied. In operation they are subjected to cycles of high temperature and high pressure.
Allan (1968) proposed a solution for the critical buckling load for vertical mode for a thin
strip with an imperfection at the centre under axial compressive load on a flat bed. Little
further work has been done to modify or include other parameters into the simple formula
Hobbs (1981, 1984) solution for a straight pipeline is still considered to be the design
guideline for buckling phenomena related to on-bottom pipelines. Both upheaval and
lateral critical buckling loads had been determined and expressed as a function of the
length of the buckle. For upheaval buckling it was assumed that the bending moment and
slope at the lift-off point is zero. For lateral buckling the lateral resistance is assumed
fully mobilised and the zero slope condition at the end of the buckle is replaced by a shear
force condition there. The buckling solutions proposed by Hobbs was based on the Kerr’s
171
Structural analysis of submarine pipelines
(1974, 1978) studies on the buckling of railway tracks, assuming a rigid foundation with
Maltby & Calladine (1995) extended Hobbs’ solution for a straight pipe on a rigid seabed
by including an initial out of straightness of the pipe. They related the critical buckling
load to the central height (or imperfection) of the pipe above the seabed. They also
proposed a solution for lateral buckling of an imperfect pipe on a seabed with a non-linear
mobilisation of lateral resistance. More details of their work has been provided in the
following sections.
Croll (1997) derived the critical vertical buckling load of a straight pipe on flat seabed,
imperfect pipe on a point protrusion support and imperfect pipe on a continuous support.
Similarly, Tran & Taylor (1996) expressed the critical upheaval buckling load of an
buckling load solution for a bar under simultaneous axial and transverse loading. In that
model the bar was considered to be straight and lying on an elastic foundation under axial
compression.
Joint industry projects (JIPs) such as the previous HOTPIPE and subsequent SAFEBUCK
pipelines ( Collberg et al. 2005, Bruton and Carr 2011). In these projects the buckling of
subsea pipelines subjected to internal pressure and temperature changes was investigated.
SAFEBUCK project outlined a design approach that addresses the inherent uncertainties
of the lateral buckling problem in a safe and quantifiable way (Bruton and Carr 2011). .
Submarine pipelines that are vulnerable to Euler buckling in the horizontal / lateral plane
are considered in this project. Hobbs (1981, 1984) solutions were adopted to provide
172
Structural analysis of submarine pipelines
analytical solution of lateral buckling and thereafter numerical solution was provided to
A key step during the design of a HP/HT pipeline is to assess whether the pipeline will
be susceptible to on-bottom buckling. If buckling is likely, then the design must ensure
that buckles will form safely and reliably at a suitable spacing to absorb the expansion,
without overstressing the pipeline or leading to fatigue. A key element of this analysis is
assessment of the critical load at which buckling will occur. This load is then compared
with the expansion forces created by the constrained expansion created by changes in
These solutions are simple and easy to use without any complexity of the soil-pipe
interaction. However, in practice the geotechnical parameters play a role, along with the
as laid geometry and pipe material, in the buckling behaviour of submarine pipelines.
This study evaluates the effects of as laid geometry of the pipe, ultimate lateral and axial
soil resistances and pipe-soil stiffness on the critical buckling load. At the end an attempt
is made to include all the relevant parameters into a single correlation, which then can be
Table 5-1 summarises the assumptions and expressions for the critical buckling load of
submarine pipelines given by various researchers over the last two decades. The aim of
this study is to provide insight into how the axial and lateral components of pipe-soil
stiffness and limiting resistance affect the lateral buckling response, and to derive new
expressions for the critical buckling load that take account of these components of pipe-
soil interaction.
173
Structural analysis of submarine pipelines
Buckling analysis of pipelines has created considerable interest in the past due to the
complexities involved in the analysis. JIP projects (e.g. Collberg et al. 2005, Burton &
Carr 2011) have provided a general overview of buckling phenomena and solutions
derived from classical theory. However, the separate effects of soil-pipeline interaction,
pipe material and geometry on the buckling initiation load were not properly identified
and discussed. This paper presents results of detailed numerical analyses carried out to
lateral buckling load of submarine pipelines. The analyses were carried out in two steps.
In the first step, the pipe response was investigated using elastic-plastic lateral springs
attached along a pipe with a lateral out-of-straightness (OOS) but without any initial
stress. The numerical solutions were found to match theory well, although it was found
that the classical theories do not allow for the effect of lateral soil stiffness. Hence, the
classical solutions were extended, based on the results of the analyses, to express the
second set of analyses, axial soil elements were introduced into the previous analyses but
the pipe was considered without any initial stress. This allowed another non-
factors representing the axial and lateral soil stiffnesses and ultimate resistances, in
addition to the initial OOS of the pipeline for a range of pipe parameters. The third set of
analyses were the extension of the second set but with the initial stress in the pipe.
The pipe was defined by diameter, D, thickness, t, Young’s modulus, E, Yield stress, σy,
with a OOS, Δ, nominal length of the OOS, L, and the bending rigidity, EI and axial
stiffness, EA can be calculated. The loads per unit length on the pipeline are defined as
horizontal limiting load per unit length H, parallel to the y-axis and axial limiting load
174
Structural analysis of submarine pipelines
per unit length F parallel to the x-axis (see Figure 5-1). Likewise the elastic slips in the x
and y axes are δx and δy, respectively as shown in the same Figure 5-1 and the
displacements in the x and y axes are δX and δY respectively. The stiffnesses in the two
directions ky, kx are defined as the limiting load (per unit length) per unit displacement,
so with typical units of kN/m2 or kPa, and are expressed as ky = H/ δy for horizontal/lateral
direction and kx = F/ δx for axial direction. The output of the analyses was the critical
buckling load and expressed as Pcr. The dimensional analysis of buckling is discussed in
The aim of this study is to provide generic design charts expressing the resulting critical
buckling load of the pipeline as a function of the various input parameters. Therefore,
dimensional analysis is a particularly suitable starting point for this study. Although this
kind of analysis alone does not usually lead to a complete solution of the problem, it
A summary of the various input parameters considered for the analyses is provided in
Table 5-2, along with the main output parameter, which is the critical buckling load , Pcr.
There are eight input variables, but only two dimensional outputs: buckling load (kN) and
nominal length of the OOS (m). As such, the output may be expressed, non-
combinations of the input variables are available, but the underlying dimensions of this
problems are length and force. Therefore, here it has been decided to use the bending
rigidity EI, as standard force and a length that approximates the pipeline diameter D. For
175
Structural analysis of submarine pipelines
pipes of moderate wall thicknesses, the nominal diameter of the pipe may be
approximated as
8𝐸𝐼 (5-1)
𝐷≈√
𝐸𝐴
For the normalisations used in the paper D refers to the external diameter and the true
The other length quantities are the out of straightness parameter , and the nominal length
L of the out of straight section of the pipeline. A more precise definition of L is provided
later, with respect to the actual pipeline profile adopted in the numerical analyses. These
Calladine (1995) for upheaval buckling of a pipe of weight (per unit length) W resting
on a rigid seabed. Here, where lateral buckling is being considered, the lateral resistance
H replaces the pipeline weight. In addition, the pipe cross-sectional rigidity EA, axial
resistance F and corresponding axial and lateral pipe-soil stiffnesses are considered. Also,
elastic slips in the lateral and axial directions were considered by dividing elastic
stiffnesses by the corresponding forces. These length quantities may be normalised by the
176
Structural analysis of submarine pipelines
Before starting the numerical analyses it is useful to discuss the theory proposed by
Maltby & Calladine (1995). They described an experimental study of some aspects of
buckling with rigid base were reviewed in their study and the experimental observations
on the growth of initial imperfection in the pipe under constant increasing axial load, and
the formation of unstable buckles were in good agreement with the findings of Tvergaad
For this analyses it is helpful to consider the Maltby & Calladine (1995) buckling solution,
as summarised in Table 5-1, replacing the pipe weight by the limiting lateral resistance
of the soil. They gave the critical load for a straight pipe resting on a rigid (smooth) seabed
as
H EI (5-3)
Pcr = 3.96
177
Structural analysis of submarine pipelines
1/4 −1/4
𝛥 1 𝐻𝛥3 𝛥 𝐻𝛥3 (5-6)
= ( ) or 4.51 ( ) =1
𝐿 4.51 𝐸𝐼 𝐿 𝐸𝐼
The boundary conditions in the Malty & Calladine (1995) solution led to a constraint on
the magnitude of the buckle length L in terms of H and EI. However, in the more general
analysis presented here, the initial length L of the out of straight section of pipeline is an
independent input variable, although also reflects closely the length of the buckle that
transpires.
An expression for the critical buckling load that reflects both H and L may be obtained
by dividing the right hand side of equation (5-5) by some power, n, of equation (5-6), to
give
1/2 𝑛/4
𝑃𝑐𝑟 𝛥2 𝐻𝛥3 𝐻𝛥3 𝛥 −𝑛 (5-7)
= 3.96 ( ) ( ) . (4.51 )
𝐸𝐼 𝐸𝐼 𝐸𝐼 𝐿
The relative contributions of H3/EI and /L may be ‘balanced’ by taking n = 2/3, which
2/3
𝑃𝑐𝑟 𝛥2 𝛥 −2/3 𝐻 ⋅ 𝛥3 (5-8)
= 1.45 ( ) ( )
𝐸𝐼 𝐿 𝐸𝐼
Justification for such balancing of the contributions from H3/EI and /L is provided later
5.4 METHODOLOGY
for the generalised relationship of equation (5-2), finite element analysis was performed
178
Structural analysis of submarine pipelines
using the commercially available software ABAQUS Version 6.9 (Dassault systèmes,
2007).
pipeline subjected to axial compressive loading. The seabed was modelled with lateral
and axial nonlinear, springs to idealise the load-displacement behaviour of the soil, and
The dimensional analysis was carried out in three stages. In the first stage only lateral
springs were used, with the axial pipe-soil response assumed frictionless. In the second
stage axial springs were introduced along with the lateral springs. In both those stages,
the out-of-straightness of the pipe was simulated without incorporating any internal
stresses in the pipe wall. A third set of analyses were then undertaken where the pipe was
‘forced’ into the out-of-straight shape, giving rise to internal stresses in the pipe. Note
that for the first, so-called “without-axial-spring” stage, the general expression for the
The fixed input parameters which were used to calculate the range of input parameters
(OOS) in the horizontal and vertical planes, either due to the natural form of the seabed
an initial stress free cosine or sine curve for the nonlinear shape of the as laid pipeline
179
Structural analysis of submarine pipelines
will generally suffice for predictions of the buckling load. The initial stress free shape
considered by Tvergaard & Needleman (1980) was modified slightly for the analyses.
The shape of the initial geometry is shown in Figure 5-2 and is given by
𝑥−𝑥0 2 (𝑥 − 𝑥 0 )𝜋
−3( ) (5-10)
𝑦= 𝛥 [𝑒𝑥𝑝 𝐿𝑠 ⋅ 𝑐𝑜𝑠 ( )]
𝐿
where Δ is the maximum out of straightness, Ls controls the value of the minor buckles to
either side of the major (central) buckle of length L ( length of the OOS section). The
quantity x0 corresponds to the distance of the centre of the buckle from the loaded end of
the pipe, with x0 taken as 2L. Note that the OOS section of the pipeline was taken
deliberately close to the axially loaded end of the pipeline rather than at the middle, to
minimise the reduction in axial load between the end of the pipe and the buckle (for
analyses that included axial springs). This helped to localise the buckle in the region of
OOS, where the axial load was high. In reporting results, actual values of Pcr are reported
from the centre of the buckle, together with actual values of Δ and L. The Δ/L ratio was
varied from 0.01 to 0.05 (Perinet & Simon, 2011). This range represents typical levels of
5.4.3 Lateral and axial soil response: stiffness and limiting resistance
The adopted ranges of lateral and axial soil resistance span typical field conditions on soft
clay. The ratio between the highest and lowest values of lateral soil resistance, H, exceeds
1000, while this ratio for axial resistance, F, exceeds 4000 (while the ratio F/H varies by
a factor of 3, because these properties are closely correlated) (Table 3). A mean estimate
of the stiffness and limiting resistance in the lateral and axial directions was made using
180
Structural analysis of submarine pipelines
the parameters given in Table 5-4, and calculation methods based on elastic and plasticity
(Guha et al., 2016, Randolph & White, 2008, Merified et al., 2008) . These typical values
were then extended to create the wider range spanned by the parametric study.
𝑘𝑥 = 𝑚𝐷𝑠𝑖𝑛𝜃𝐷 (5-11)
where kx is the axial soil stiffness, m is the constant of shear modulus (G = mz, z is the
depth below the pipe), and θD is the embedment angle of the pipe as discussed in the
previous chapter. The axial elastic stiffness, kx span a range of 0.04 to > 200 kN/m2 was
considered here.
where F is the ultimate axial frictional resistance of the soil per unit length, α is adhesion
factor, undrained shear stress varies with depth z as su_seabed = kz, D is the diameter of the
The lateral pipe-soil stiffness, ky was estimated by dividing the ultimate lateral resistance,
H, by elastic slip, δy . For this study δy in the range of 0.05D to 0.5D was considered. As
highlighted in the paper Guha et al. (2016) that lateral and vertical stiffnesses predicted
numerically are rather higher. The higher values for homogeneous soil modulus may be
attributed, at least partly, to the limited layer modelled numerically. Therefore, a range of
0-20% was proposed in the paper (Guha et al., 2016) for vertical and lateral elastic
stiffnesses. Therefore, a slightly higher value of kx/ky =0.7 was adopted in the present
chapter.
181
Structural analysis of submarine pipelines
To analyse the problem of buckling associated with material nonlinearity and geometric
nonlinearity prior to buckling an advanced method such as the Riks algorithm (Riks,
1979) is very useful. The Riks method is a static stress-based analysis method which uses
the load magnitude as an additional unknown; it solves the load and deflection
collapse behaviour, where the load-deflection response shows a negative stiffness and the
structure must release strain energy to remain in equilibrium. The Riks algorithm allows
for the complete load-deflection curve to be determined including a reversal of load path.
bifurcation and snap-through buckling and obtain a complete solution of the post-
buckling response. This method is used for cases where the loading is proportional, that
is, where a single scalar parameter governs the load magnitude. This is very important for
the analysis of structures carrying de-stabilising loads such as axial load in a pipeline,
column or pile. The advantage of the method is that the solution is viewed as the discovery
of a single equilibrium path by the loading parameters and the nodal variables. The
To reduce the complexity of the problem and make a simple solution for the same, a
loading pattern is initially defined by a load vector (Pn), this load vector represents the
relative magnitude applied to the structure. Then a proportional load factor (λ) is applied
to all of these loads. The load proportionality factor is found out as the part of the solution.
displacement (u) and load factor. The total load can be expressed by the multiplication of
the load factor with the set of unit loads over the structure at any increment (n)
182
Structural analysis of submarine pipelines
Ptotal = P n (5-13)
The response of a structural system in discrete set of nonlinear equation can be expressed
f ( , u ) = 0
(5-15)
These equations are solved in a series of steps or increments, starting from λ = 0 condition.
The nth step starts from a known solution on the equilibrium path in the solution space ζ
= (λ, u)T (Ragon et al., 2002) .This consists of two phases, Prediction phase and
Correction phase. In the prediction phase estimation is generated on the equilibrium path
𝜁 0 = (𝜆0𝑛+1 , 𝑢𝑛+1
0 )𝑇
(5-17)
With this starting point Newton-Raphson method is employed to find the next point on
the equilibrium path. The size of each step is fixed with an arc length Δl. The general set
of equations for i-th iteration executed by each algorithm at step (n+1) is expressed as
183
Structural analysis of submarine pipelines
Df ( ) = f ( )
i −1
n +1
i −1
n +1 (5-18)
( )
c ni −+11 = 0 (5-19)
where
= ni +1 − ni −+11 , i = 1, 2, … (5-20)
and
Equation (5-21) is the m x (m+1) Jacobian matrix of the system f(λ, u) = 0 and defines the
tangent stiffness of the structure in λ-u space. As mentioned above, an arc length Δl is
introduced in ABAQUS (Crisfield, 1981) which limits the length of the displacement
vector. This prevents the internal stress residual to become too large and iterating back
the solution onto the equilibrium path. In this modified algorithm a hypersphere is defined
with a locus on the λ-u plane as used by Knappett & Madabhushi (2009) for the lateral
𝑖 𝑖
(𝜁𝑛+1 − 𝜁𝑛∗ )𝑇 (𝜁𝑛+1 − 𝜁𝑛∗ ) = 𝛥𝑙 2 (5-22)
As mentioned earlier the (n+1)th increment will follow a two-stage procedure namely
prediction and correction stages. Starting from a known solution point from the nth
increment, ζ*n+1, the bounding radius or the arc-length Δl is defining the hypersphere and
184
Structural analysis of submarine pipelines
a prediction step is made to estimate ζ0n+1, on the bounding radius. At this stage
equilibrium conditions are checked for the newly defined point as the internal stresses
were generated by the external deflection. If the equilibrium condition is not satisfied,
Newton-Raphson iteration is used in the correction step for the point ζ 0n+1 to the next
point ζ1n+1 and onwards as shown in Figure 5-3. This procedure continues until the
For nonlinear buckling analysis modified Riks method is very popular. Torselleti et al.
(1999), as part of the Hotpipe project, carried out FE analysis of submarine pipeline
resting on flat seabed with ABAQUS. They employed RIKS method to compare the
analytical solution with their numerical results for initial imperfection and pipe material
lateral buckling influenced by the lay induced imperfection. In risk based design and
ABAQUS was also very popular (Maschner & Wang, 2008, Peek & Yun, 2007) with
beam or pipe elements and were mainly used for structural analysis of long pipelines for
numerical advantages (Torselleti et al.,1999; Peek & Yun ,2007; Maschner & Wang,
2008, Odina, & Tan, 2009, Chee et al., 2019, Chee et al., 2018). The details of beam and
pipe elements and the pipe-soil interaction are discussed in the following sections.
al. (1999), Zhou et al., (2019) and Wang et al. (2015), revealing that the three methods
are made so that the problem is reduced to one dimension mathematically: the primary
185
Structural analysis of submarine pipelines
solution variables are functions of position along the beam axis and the shortest distance
from the axis to any point in the continuum is smaller than typical lengths along the axis
of the element. The element library in ABAQUS contains several types of beam elements
and these are generally divided by the approach to beam theory into three parts: Euler-
Bernoulli beam, Timoshenko beam, and hybrid beam elements as shown in the Figure
5-4. Firstly, the Euler-Bernoulli beam theory is a classic and simplest in that it is based
on the assumption that cross-section does not deform in its plane or warp out of its plane
and this cross-sectional plane remains normal to the beam axis. Since this beam element
can be used to simulate beams for large axial strains as well as large rotations, this beam
theory supplies generally satisfactory results for slender beams. The ABAQUS beam
element library provides B23, B33, and so forth, as this type of beam.
The Timoshenko beam theory allows transverse shear strain that is the cross-section may
not necessarily remain normal to the beam axis. This approach can be useful for thicker
beams whose shear flexibility may be important. Also, this is effective to simulate thin
beams where Euler- Bernoulli theory is accurate. There are several element types in the
ABAQUS beam element library, such as B21, B22, B31, B32, etc. The hybrid beam
elements in ABAQUS are designed to handle very slender situations where the axial
The pipe elements in ABAQUS such as PIPE21, PIPE22, PIPE31, PIPE32, etc and their
hybrids are available with hollow, thin-walled, circular section. The hoop stress caused
by internal or external pressure loading in the pipe is included. The main purpose of these
is to provide for the subsequent difference in yield when a point in the section is subjected
to axial tension or compression stress due to its nonzero hoop stress component causing
an asymmetry in the section’s response to inelastic bending as shown in the Figure 5-5.
186
Structural analysis of submarine pipelines
The hybrid pipe elements such as, PIPE31H, PIPE32H, etc are used where it is difficult
to compute the axial and transverse shear forces in the beam by the FE displacement
method. This kind of problems of arises in geometrically nonlinear analysis when the
element undergoes large rotations and is very rigid in axial and transverse shear
deformation, such as flexible long pipe or cables. The hybrid elements use a more
generalised formulation in which the axial and transverse shear forces in the elements are
included, along with the nodal displacements and rotations, as primary variables.
Therefore, PIPE31H was selected as the appropriate element for this kind of geometric
nonlinear analysis due to the option of applicability of internal and external pressure for
these elements for future research where study can be under taken to investigate the effect
of internal and external pressure on the buckling response ( Zhou et al., 2019, Peek &
Traditionally, the buckling of submarine pipelines and the seismic fault crossing effect
on buried pipelines idealises soil behaviour using discrete springs to model the pipe-soil
interaction based on the concept of subgrade reactions proposed by Winkler (Odina &
Tan, 2009). In that model the soil is represented as series of independent springs attached
to the nodes of the pipeline. An elasto-plastic spring is generally used to idealise the force-
displacement response at the pipe surface. The spring relationships in all the three
direction are defined by the maximum force per unit length at the soil pipe interface and
the ultimate mobilization displacement. The force-displacement model considered, for all
the three directions is defined in the section below. For nonlinear structural analysis pipe
187
Structural analysis of submarine pipelines
The axial, lateral and vertical soil pipe interaction was modeled with nonlinear springs
(SPRING1) available in ABAQUS, which link each node to a fixed reference point (the
ground). The nonlinearity of the springs was defined separately in the input code. The
𝛥𝑢 = 𝑢𝑖 (5-23)
where i is defined in the spring input code and is in the local direction of the model. In
ABAQUS lateral and axial springs are directly oriented towards the global coordinates.
In this scenario if the pipe element is deformed in the global x-axis, the axial soil stiffness
will contribute to both axial and lateral stiffness component. At the same time lateral soil
stiffness may affect the axial component and lateral component. Therefore, special
computing may be employed to eliminate this issue (Sollund & Vedeld, K., 2014).
However, to predict the peak buckling load in the pipe this issue has minimum impact.
Therefore, for the simplicity purposes this has been ignored here.
To carry out FE analysis the pipe was idealised as an elastic-plastic beam resting on a
nonlinear foundation with one end on the roller and the other end was pinned. At the roller
end a displacement was given which idealises the compressive load in the axial direction.
A lateral imperfection was deliberately introduced into the pipe to investigate the effect
of OOS on the critical buckling load. The whole length of 1 kilometre of the pipeline was
divided into small elements of approximately 1 metre length and at each node of the
element bilinear springs were attached, which represents the soil restraints in the lateral
and axial directions. The effect of element size was also investigated by changing the
188
Structural analysis of submarine pipelines
element size and thereby changing the number of springs. However, the above
economy.
Based on the model described in the above sections a simple case study was performed
to check the sensitivity of the model in a series of buckling analyses. In the first case the
axial springs were not attached and the OOS was considered to be initially stress free.
Only lateral springs were attached to analyse the behaviour of lateral buckling during an
axial compressive loading. The effects of temperature, pressure and external loads were
incorporated by displacing the end of the modelled section of pipe axially in compression.
With the model discussed in the above section an example analysis was carried out to
check the effect of input parameters on the expected outcomes. For this example a pipe
with D = 0.152 m (6″) was considered with D/t ratio as 30.3 and σy / E as 0.002. The pipe
has a bending rigidity, EI of 1.24 MNm2. A bi-linear soil model was considered with
normalised limiting lateral resistance, HΔ3/EI = 4.24x10-05 and elastic stiffness, ky/H√
(8EI/EA) = 16.66 (H = 424 N/m and δy = 0.008 m). The ratio of Δ/L was varied from 0.01
to 0.05.
Axial load-displacement curves for various Δ/L ratios are compared in Figure 5-6, with
the displacement of the loading end plotted against the load at the central node of the
buckle.For smaller OOS the peak buckling load is higher and the instability is more
prominent. With increasing OOS the axial load-displacement response does not show
‘snap-through’. The summary of peak loads is plotted as shown in Figure 5-7. It depicts
the variation of peak buckling load with OOS. The peak load at Δ/L = 0.01 is 3.5 time
189
Structural analysis of submarine pipelines
higher than the peak load at Δ/L = 0.05. With increasing Δ/L ratio the load required for
buckling decreases. Therefore, it agrees with the previous findings that OOS tends to
reduce the critical buckling load as presented by others (Maltby & Calladine, 1995, Croll,
The relationships between axial load and lateral offset of the midpoint of the buckle are
plotted for various Δ/L ratios in Figure 5-8. The lateral displacement of the midpoint of
the buckle is plotted including the initial OOS values in order to show the effect of the
total lateral offset on the sustainable buckling load. Just as for the axial load-displacement
case higher OOS does not allow ‘snap-through’ buckling. After the peak load reaches its
maximum value it reduces and merges with the peak load with higher Δ/L ratios.
Figure 5-9 shows the lateral friction distribution at the buckling initiation along the length
of the pipeline for Δ/L = 0.01 to 0.05. The ultimate lateral resistance was fully mobilised
through the central lobe prior to the onset of buckling for /L > 0.01, and for /L = 0.05
the ultimate resistance was also fully mobilised in the secondary lobes when buckling
occurred. These results show the significance of initial soil condition and effect of OOS
The comparison between the initial shape and the corresponding displaced shapes of the
buckles at the peak load are shown in Figure 5-10. A smaller OOS requires a larger load
to mobilise the elastic limit. Figure 5-9 show that a larger OOS is more prone to buckle
and the corresponding lateral displacement is higher for smaller OOS. It is evident from
Figure 5-11 that the lateral offset of the pipe is very low (0.0093 m) for Δ/L = 0.01 when
buckling occurs, and has only just exceeded the elastic limit within the lateral force-
increases with increasing OOS. For Δ/L = 0.05 the offset of the midpoint of the buckle at
190
Structural analysis of submarine pipelines
To show the effect of lateral limiting friction on the peak load, FE results for three
different limiting lateral resistances (H = 100, 200 and 424 N/m) were compared. Figure
5-12 illustrates the influence of the lateral limiting resistance on the peak load of the
middle node. It is evident that the buckling initiation is significantly influenced by the
lateral limiting residence. As the limiting lateral resistance reduces from 424 N/m to 100
N/m (⁓ 4 fold decrease) the peak load reduces by a factor slightly higher than 2.
Correspondingly for a lower OOS the limiting lateral soil resistance has greater impact
on peak buckling than for a higher OOS. This is because the higher OOS does not lead to
‘snap-through’ buckling.
The analytical solution proposed by Maltby & Calladine (1995) for H = 424 N/m was
plotted on the same graph (dotted line) to compare with the FE solutions. On the same
graph (Figure 5-12) the dimensionless modified Maltby & Calladine (1995) results for H
= 424 N/m was also plotted. It is evident that the FE results are lower than the analytical
results because of the initial assumptions in the analytical modelling: (1) analytical
solution considered a fully mobilised lateral limiting friction; and (2) In the FE solution
the initial shape similar to mode 3 (of Hobbs solution) was considered predicted by
equation (5-3) whereas for the analytical solution mode 1 (of Hobbs solution) was
To show the effect of lateral soil stiffness (or reverse of elastic slip) on the peak load a
series of analyses were carried out for elastic stiffness parameter of the lateral spring, ky
ranging between 0.5 – 5 kN/m2. This was achieved by varying the elastic slip, δy between
0.05D – 0.5D for H = 424 N/m and EI = 1.24 MNm2. Figure 5-13 illustrates the effect of
elastic slip on the peak load. The elastic slip was normalised by the diameter of the pipe
191
Structural analysis of submarine pipelines
(0.152 m or 6”). With increasing elastic slip (decreasing elastic stiffness) of the soil the
Before starting the parametric study to explore the cases of ‘without axial resistance’,
‘with axial resistance and without initial stress’ and ‘with axial resistance and with initial
stress’, a simple graph was plotted to show the effect of these three cases on the
normalised peak buckling load . For this analysis similar input parameters were chosen
as above, just changing the axial soil resistance, but keeping F/H at 0.3. Also for the ‘with
axial and with initial stress’ the pipe was pre-stressed before commencing the real
buckling simulation. Figure 5-14 shows the comparison of peak buckling load varying
with OOS for three different cases. Difference between the peak loads decrease
significantly with higher OOS. It is evident from the figure that peak load in the pipe is
higher with axial resistance and highest with axial resistance and initial stress. For an out-
of-straightness, = 2.5 m (/L = 0.05) the value of peak load for ‘with axial and without
initial stress’ is 1.43 times higher than that for ‘without axial’ case due to the addition of
axial springs (resistance) along the length of the pipeline. Similarly, initial stress in the
pipe increases the peak load further 1.11 times for the same value. Based on these
results a comprehensive parametric study was carried out in the next section.
Many numerical analyses were carried out for a range of cases to establish a general
relationship for the critical buckling load, Pcr. This load was found by analysis of the
responses from >200 separate analyses using a range of pipe and pipe-soil parameters,
including different Δ/L ratios. Pipe with three different diameter sizes of 0.152 m (6″),
0.304 m (12″) and 0.912 m (36″) with D/t ratio of 30, were considered. The results were
plotted in nondimensionalised manner all cases. First a design chart was proposed based
192
Structural analysis of submarine pipelines
on a parametric study without any initial stress or an axial pipe-soil spring. Then the
individual effect of parameters such as ultimate soil resistances and OOS were
investigated through a systematic parametric study. Finally another design chart was
proposed for the case with initial stress and axial pipe-soil spring.
In this step only lateral and vertical nonlinear springs were used along the length of the
pipeline. No axial restraint was applied except for the fixed end boundary condition. The
critical buckling load for various embedment ratios and for various ultimate soil
resistances was plotted and the following nondimensionalised fit was developed:
Figure 5-15 illustrates the comparison between the actual critical buckling load and the
predicted critical buckling load. The proposed multivariate fit shows a good agreement
with the FE results. The mean absolute error |1- Pcr_predicted /Pcr_FEA| between equation
(5-24) and the FEA result was 0.14 with standard deviation value of 0.033. Whereas
Maltby & Calladine (1995) solution gives an error of 1.01 with a standard deviation value
of 0.059. Therefore, predicted critical buckling load is in good agreement with the FEA
results.
In this step along with axial and lateral nonlinear springs, axial springs were also used
along the length of the pipeline. Similar to the previous cases the critical buckling load
193
Structural analysis of submarine pipelines
for different as-laid geometry of the pipes and for various ultimate soil resistances was
0.1 −0.12
0.66 −0.1
𝑃𝑐𝑟 ∙ Δ2 Δ −0.66 H ∙ Δ3 k𝑦 8𝐸𝐼 Δ 𝐹 −.017 k𝑥 8𝐸𝐼 (5-25)
= 1.15 ( ) ( ) ( √ ) ( ) ( ) ( √ )
𝐸𝐼 𝐿 𝐸𝐼 𝐻 𝐸𝐴 √8𝐸𝐼/𝐸𝐴 𝐻 𝐹 𝐸𝐴
This solution verifies the classical buckling solution and extends it with axial and lateral
soil stiffness terms. Figure 5-16 illustrates the comparison between the calculated and
predicted critical buckling load. The proposed multivariate fit shows good agreement with
the FE results. The mean absolute error |1- Pcr_predicted /Pcr_FEA| between equation (5-25)
and FEA result was found out to be 0.63 with a standard deviation value of 0.058, which
is again an improvement on Maltby & Calladine (1995), for which the mean error was
However, the solution does not include the initial stress in the pipe due to its shape and
the axial component of the soil resistance. The following section adds these aspects.
In the field due to the pipe lay process and any uneven seabed topography, the pipe gains
an OOS and hence the pipe is pre-stressed, prior to the addition of axial stresses. Initially
a straight pipe was chosen and then it was bent into the desired shape using equation (5-
10). This is called ‘pulled in place’ case. Thereafter, the soil springs were activated i.e.
the pipe was put in contact with the soil. It was anticipated that the pipe would rebound
slightly due to the elastic slip, as well as the redistribution of elastic stresses within the
pipe. Therefore, the buckle portion, particularly the central lobe and the mid-point of the
central lobe, was pushed slightly further – so that even after rebounding the shape it would
still provide the same OOS. The net effect was that the pipe was pre-stressed before
194
Structural analysis of submarine pipelines
commencing the real buckling. At this stage the springs along the length of the pipe were
deactivated and did not impose any resistance to the lateral displacement to achieve the
deformed shape. Thereafter, all the springs were activated and buckling analyses was
carried out, via displacement control simulation. Similar to the previous cases, the ratio
Δ/L was varied between 0.01 to 0.05 to cover a wide range of OOSs.
Finally, the effect of all soil parameters and the as-laid geometry on the critical lateral
0.1 −0.18
0.66 −0.1
𝑃𝑐𝑟. Δ2 Δ −0.66 H. Δ3 k 𝑦 8𝐸𝐼 Δ F −0.26 k 𝑥 8𝐸𝐼 (5-26)
= 1.35 ( ) ( ) ( √ ) ( ) ( ) ( √ )
𝐸𝐼 𝐿 𝐸𝐼 𝐻 𝐸𝐴 √8𝐸𝐼/𝐸𝐴 𝐻 𝐹 𝐸𝐴
Figure 5-17 shows the comparison between the actual and the predicted normalised
critical buckling load of pre-stressed pipeline with similar seabed conditions and
geometry as above. The mean absolute error |1- Pcr_predicted /Pcr_FEA| between
equation
(5-26) and FEA result was 0.14 with a standard deviation value of 0.068, which is an
improvement on the Maltby & Calladine (1995) solution, for which the mean error was
0.58 and with a standard deviation value of 0.084. The multivariate fit shows good
agreement with the actual FE results. Lateral friction exists to restrain the pipe from
rebounding back for the pulled in place case, while there is no stress in the wished in
place case. During buckling the initial friction is first mobilized similar to the pulled in
place case. Then a higher axial force is needed to activate buckling. As for the wished in
place case, the axial compression directly causes the lateral friction to accumulate. As a
result, the buckling of the pulled in place case is delayed and more energy accumulates,
potentially leading to snap-through buckling (Zhou et al 2019). The increase in peak load
195
Structural analysis of submarine pipelines
due to initial stress in the pipe was also observed by the others (Zhou et al., 2019, Miles,
1998). Therefore, the results can be used by the designers to predict the buckling load of
5.7 CONCLUSIONS
multi-variate fit to the relevant dimensionless groups was developed to predict critical
buckling load. The seabed was idealised with horizontal, axial nonlinear springs to
idealise the load-displacement behaviour of the soil and the pipe was modelled with pipe
with the modified Riks algorithm that is available in the commercial software ABAQUS.
In the first step the classical theories were verified with the numerical tool and a simple
multivariate fit was proposed by carrying out parametric study. Then in the next steps
actual analysis was carried out by incorporating the axial springs and the initial stress into
the numerical model. The results have been synthesised in a generic non-dimensionalised
design chart to estimate the buckling load, valid for the range of inputs covered by the
parametric study.
196
Structural analysis of submarine pipelines
References:
Allan, T. (1968). One-way buckling of compressed strip under lateral loading. J. Mech.
Bruton, D., & Carr, M. (2011). Overview of the SAFEBUCK JIP. 10.4043/21671-MS.
Chee, J., Walker, A., & White, D. (2018). Controlling lateral buckling of subsea pipeline
DOI: 10.1016/j.oceaneng.2018.01.024
Chee, J., Walker, A., & White, D. (2019). Effects of variability in lateral pipe-soil
10.1016/j.oceaneng.2019.04.034
Collberg, L., J. Mo̸rk, K., Levold, E., & Vitali, L. (2005). HotPipe JIP: Design Guidelines
Dassault Systѐmes (2007) Abaqus analysis users’ manual, Simula Corp, Providence, RI,
USA.
Guha, I., Randolph, M. F., and White, D. J . (2016). Evaluation of elastic stiffness
5606.0001466
197
Structural analysis of submarine pipelines
Hetényi, M. (1946). Beams on elastic foundation. Univ. of Michigan Press, Ann Arbor,
MI.
Hobbs, R. E. (1981). Pipeline buckling caused by axial loads. J. Construct. Steel Res. 1(2),
2-10.
175-189.
Kerr, A. D. (1974). On the stability of the railroad track in the vertical plane. Rail Int. 5:
132–142.
Kerr, A. D. (1978). Analysis of thermal track buckling in the lateral plane. Acta Mech. 30
583-592.
Maschner, E. A., and Wang, N. Y. (2008). Risk design for stable ‘thin-walled’ pipelines
subjected to global buckling and excessive bending from seabed and pipelay
Massey, B.S. (1971). Units, Dimensional Analysis and Physical Similarity. University
Merifield, R. S., White, D. J., and Randolph, M. F. (2008). The ultimate undrained
198
Structural analysis of submarine pipelines
Odina, L., and Tan, R. (2009). Seismic fault displacement of buried pipelines using
79739.
Peek, R., & Yun, H. (2007). Flotation to triggred lateral buckles in pipelines on a flat
Perinet, D., and Simon, J. (2011). Lateral buckling and pipeline walking mitigation in
21803-MS
Ragon, S. A., Gürdal, Z., and Watson, L. T. (2002). A comparison of three algorithms
for tracing nonlinear equilibrium paths of structural systems. Int. J. Solids &
London, 27-59.
Sollund, H., and Vedeld, K. (2014). A finite element solver for modal analysis of multi-
199
Structural analysis of submarine pipelines
Sriskandarajah, T., Dong, S., S. Sribalachandran, S., and Wilkins, R. (1999). Effect of
Taylor, N., and Tran, V. (1996). Experimental and theoretical studies in subsea pipeline
Torselleti, E., Vitali, L., and Levold, E. (1999). Hotpipe Project: Snaking of submarine
pipelines resting on flat sea bottom using finite element method. Proc., of 9th
Wang, Z., Huachen, Z., Liu, H., and Bu, Y. (2015). Static and dynamic analysis on
256.
White, D. J., and Randolph, M. F. (2007). Seabed characterisation and models for
Zhou, T., Tian, Y., White, D. J., and Cassidy, M. J. (2019). Finite element modelling of
offshore pipeline lateral buckling. J. Pipeline Syst. Eng. Pract., 10 (4): DOI:
10.1061/(ASCE)PS.1949-1204.0000396.
200
Structural analysis of submarine pipelines
Tables:
Table 5-1: Summary table of estimation of critical load of buckling of submarine pipelines:
𝑤𝐸𝐼
𝑃 = 3.95√ wL/2
Δ
L/2 Allan (1968)
seabed
201
Structural analysis of submarine pipelines
L = 2
w AE
P = 4 2 EI / L2
202
Structural analysis of submarine pipelines
seabed
203
Structural analysis of submarine pipelines
𝐸𝐼𝐻
𝑃=√
Δ L Maltby & Calladine
D
where EI = bending rigidity and H = (1995)
𝐸𝐼𝐻 0.25
expressed as: 𝜆 = ( )
Δ
seabed
𝐸𝐼𝑤
𝑃 = 4.029√
Δ Croll (1997)
204
Structural analysis of submarine pipelines
and wave length was expressed as A perfectly straight pipe of unit weight w was resting on a flat, rigid
protrusion support
w
The propagation load was expressed as
D
𝐸𝐼𝑤
𝑃 = 4.029√
Δ
𝐸𝐼Δ
𝐿 = 4.427√ 𝑤
where w = weight per unit length and = A rigid protrusion was considered with amplitude of at the centre. With
displacement relative to flat seabed increasing axial load the pipe will start lifting off from the support and
205
Structural analysis of submarine pipelines
Croll (1997)
Vertical mode- imperfect pipe on
4c.
y’’
continuous support D
𝑃 = 3.007√𝐸𝐼𝑤/Δ The seabed profile was considered as dominant buckling mode. The
The maximum propagated buckle load profile of the seabed and the lift-off phenomena with propagation of
buckling load was considered here and shown in the figure. With
expressed as 3.478√𝐸𝐼𝑤/Δ and the
increasing axial load the pipe starts lifting of the seabed support. Two
corresponding length of the buckle was
design load criteria had been defined and closed form solution given for
given by 𝐿 = 4.427(𝐸𝐼Δ/𝑤)1/4
the buckling of submarine pipelines with imperfection.
206
Structural analysis of submarine pipelines
D
The critical load of buckling was
L/2 L/2
expressed as
𝑤𝐸𝐼
𝑃 = 1.5848√ Δ The empathetic model of continuous supported pipe had been analysed
= 40% P_hobbs, where = imperfection here. The critical buckling load and buckling lengths were found out
amplitude and w = weight per unit length theoretically and was supported by experiments.
207
Structural analysis of submarine pipelines
L/2 L/2
The critical load of buckling was
6 Straight bar under simultaneous axial In this model a bar was considered to be straight and on elastic foundation Hetényi, 1946
208
Structural analysis of submarine pipelines
seabed
𝑤𝜇𝐸𝐼
𝑃 = 3.86 ( )
𝐷
Modified Hobbs solution for infinite mode of buckling SAFEBUCK II 2008
where w = submerged pipe weight, μ =
pipe
𝑤𝐸𝐼
𝑃 = 3.86√
Δ
209
Structural analysis of submarine pipelines
Out-of-straightness Δ m
Buckle length/amplitude L m
Quantity Range
Out-of-straightness Δ 0.5 – 10 m
Quantity Value
k = su/z 5 kN/m3
G/su 50
211
Structural analysis of submarine pipelines
Figures:
212
Structural analysis of submarine pipelines
( ni +1 − n* ) T ( ni +1 − n* ) = l 2
Loading parameter, λ Tangential solution from ζ*n
Δl
ζ0n+1
ζ1n+1
ζ2n+1
ζ*n
Crisfield iteration
f(ζ)=f(λ,u)
Displacement, u
Figure 5-3: Modified RIKS method used in ABAQUS (Knappett & Madabhushi (2009))
213
Structural analysis of submarine pipelines
0.12
D = 0.152 m (6")
EA = 462 MN Δ/L = 0.01
L = 50 m
H = 424 N/m
0.09
Axial load, Pcr (MN)
0.02
0.06
0.03 0.04
0.03
Δ/L = 0.05
0
0 0.1 0.2 0.3 0.4 0.5
Axial displacement, δX (m)
Figure 5-6: Axial load at the central node of a buckle versus axial displacement at the
214
Structural analysis of submarine pipelines
Figure 5-8: Development of lateral displacement at buckle crown as axial load increases
215
Structural analysis of submarine pipelines
500
D = 0.152 m (6")
424 N/m EA = 462 MN
L = 50 m
Lateral resistance of the soil, H (N/m)
100
Figure 5-9: Mobilised lateral soil resistance at peak load along the length of the pipeline
4
D = 0.152 m (6") Δ/L = 0.01 ( initial shape)
EA = 462 MN
L = 50 m Δ/L = 0.01 ( Shape at peak load)
H = 424 N/m
3 Δ/L = 0.02 ( initial shape)
Δ/L = 0.05
Δ/L = 0.02 ( Shape at peak load)
Lateral position, y (m)
0
0 50 100 150 200
-1
Axial position, x (m)
Figure 5-10: Comparison between initial and displaced shapes at the peak loads for
various out-of-straightness
216
Structural analysis of submarine pipelines
0.5
0.3
0.2
0.1
Elastic limit = 0.008 m
0
0 0.01 0.02 0.03 0.04 0.05 0.06
Normalised out-of-straightness, Δ/L
Figure 5-11: Lateral offset values corresponding to peak load for various OOS
0.4
FE results for H = 424 N/m D = 0.152 m (6")
EA = 462 MN
FE results for H = 200 N/m L = 50 m
Peak buckling load, Pcr (MN)
0.3
FE results for H = 100 N/m
Modified Mantby & Calladine analytical solution for H = 424 N/m eq(8)
0.2
Analytical solution of Maltby & Calladine eq(3)
for H = 424 N/m
0.1
0
0 0.01 0.02 0.03 0.04 0.05 0.06
Normalised out-of-straightness, Δ/L
217
Structural analysis of submarine pipelines
Figure 5-14: Comparison of normalised peak buckling loads for three cases
218
Structural analysis of submarine pipelines
Figure 5-15: Multivariate fit to computed critical buckling load without axial spring
Figure 5-16: Multivariate fit to computed critical buckling load with axial spring
219
Structural analysis of submarine pipelines
Figure 5-17: Multivariate fit to computed critical buckling load with axial spring and initial stress
220
Structural analysis of submarine pipelines
221
Structural analysis of submarine pipelines
222
Structural analysis of submarine pipelines
6.1 INTRODUCTION
The increasing demand of energy is not only pushing the hydrocarbon industry towards
the deeper oceans, but also pushing the operating conditions towards the high temperature
and high pressure (HT/HP) regime. Therefore, the present deepwater submarine pipelines
are being operated at HT/HP which can require special design considerations to prevent
issues related to thermal expansion, such as buckling and walking. Theories of buckling
have been developed in last two decades and are well accepted in the industry. (e.g., DNV
2007; Collberg et al. 2011 ) The so-called pipeline walking mechanism was also well
defined in last few years by various researches and by the industry projects such as
SAFEBUCK JIP (Konuk, 1998, Tørnes et. al. 2000, Carr et. al. 2006). However,
detailed aspects of pipe-soil interaction behaviour. Therefore, this paper will give a brief
overview of the present practices to estimate the walking of pipeline along with a
Once the pipeline is laid on the seabed, it is heated or cools down during operation cycles.
Due to this heat up and cool down process the pipe tends to expand and contract
respectively. These expansion and contractions are resisted by the seabed friction. When
it cools down the pipe cannot regain the original length due to the seabed resistance. This
223
Structural analysis of submarine pipelines
phenomenon is addressed in the pipeline design guidelines. In some cases the pipeline
expansion is associated with the global displacement of the pipe because the expansion
and contraction is asymmetric between the two ends, leading to a net movement of the
pipe in one direction. This is commonly known as ‘walking’ of pipeline. Walking itself
is not a limit state, but uncontrolled walking may lead to many of the critical problems,
such as, over stressing of end connectors such as spool pieces and jumpers, loss of tension
in a SCR (Steel catenary riser), increased loading within a lateral buckle and route curve
In deepwater the pipelines are often connected to the receptor facilities by SCR. This
arrangement pulls the pipeline into tension at the touch down zone. This tension at the
end of the pipeline causes a short pipe to walk under the thermal cyclic loading conditions
towards the SCR. It is assumed that there is sufficient axial friction along the length of
the pipeline to be axially stable under the highest axial riser tension (Carr et al. 2006).
Carr et al. (2006) gave the expressions for the pipeline walking under various field
conditions in work that was part of the SAFEBUCK JIP ( Collberg et al. 2011). They
expressed that, the pipeline walking behaviour of short pipelines can occur due to (1)
tension at the end of pipelines, associated with a SCR, (2) global seabed slope along the
length of the pipe, and, (3) thermal gradients along the pipeline during changes in
operating conditions. Bruton et al. (2010) added a fourth mechanism, being (4) changes
in the distribution of weight along the pipeline due to liquid hold up on shutdown. The
following section summarises the existing theories to estimate walking triggered by these
mechanisms.
224
Structural analysis of submarine pipelines
The analytical solution of pipeline walking per cycle, under the tension of SCR was given
where, ΔP is the change in fully constrained force, TSCR is the SCR tension, f is the axial
frictional force = µW′, L is the pipeline length and EA is the axial stiffness of the pipe.
P = −( p2 − p1 ) Ai (1 − 2 ) − E AS ( 2 − 1 ) (6-2)
where p is internal pressure and subscripts 1 and 2 refer to conditions before and after the
operating change, Ai is internal cross-sectional area of the pipe, ν is Poisson’s ratio of the
pipe material, As is cross sectional area of the pipe wall, E is Youngs modulus, α is
The first term is the pressure term and assumed to be zero for these analyses and the
second term is the thermal term as defined earlier. Numerical analysis was also carried
out to verify the analytical solution. Analytical and numerical results were matched with
good agreement.
The seabed slope causes the walking of pipeline under cyclic thermal loading. Carr et al.
(2006) gave an analytical model to predict the walking of an on-bottom pipeline with
sloping seabed. Figure 6-1 shows the force diagram of a pipe element resting on a seabed
225
Structural analysis of submarine pipelines
with an angle Φ. The slope of the seabed (Φ) was introduced in the previous expression
=
P + W L sin( ) − W L cos( ) L tan( ) (6-3)
EA
In this case there is a component of pipe weight (W′sin Φ) which acts downhill and in the
same direction of expansion. When the pipe expands up the slope, this force acts against
the expansion and when the pipe expands down the slope, this force acts with the
expansion. This affects the shape of the force envelope causes the asymmetric shape of
the force profile envelope as discussed later. FE analysis was also carried out by Carr et
al. (2006) to mutually verify the two methods. Analytical and numerical results were
The heat up and cool down step is needed to be considered realistically to investigate the
to the well head and the ‘cold-end’ is considered to be closest to the receiver facility or
attached to riser. The pipe (and its contents) cools down uniformly over its length (when
there is no flow through the pipe) and heats up non-uniformly from one end to the other
with some gradient (due to the flow from one end to the other). Due to this reason,
expansion in one direction occurs during the heat up step, but the contraction during the
cool down step is more symmetric. A typical thermal transient model was shown in Figure
6 2(Carr et al. 2003 ). The important aspect of this phenomenon is the shape of the thermal
profile developed over time as the pipe is heated up from ambient. The hot fluid enters
the ‘hot-end’ at 0 and heat is lost to the surrounding and the fluid quickly cools down to
ambient temperature. With time the pipe is gradually heated up till the fluid is discharged
226
Structural analysis of submarine pipelines
at the ‘cold-end’ of the pipe. Earlier research ( Tørnes et al. 2000, Carr et al. 2003, Carr
et al. 2006) has cycle was also given by revealed that steepness of the thermal transients
play a major role in driving the walking. The analytical model of walking due to thermal
transients was not straightforward and was solved by an incremental method. However,
f L2 f f (6-4)
; T 24 − − 4 when, f > fθ/6
16 EA f f
and
𝑓⋅𝐿2 (6-5)
𝛥 𝑇 ≅ 8⋅𝐸𝐴 when, f< fθ/6
where, f is the frictional resistance of the soil, fθ (= EAαqθ ) is force associated with the
thermal gradient, qθ = θ2 – θ1. FE analysis was also carried out to verify the analytical
solution. Approximate analytical solutions were in good agreement with numerical results.
This model assumed constant thermal gradient during heat up steps until the full steady
state is reached, for simplicity. In reality the gradient reduces as the pipe heats up.
A new walking mechanism was recently been added to the existing mechanisms by
Bruton et al. (2010). In multiphase pipelines, which are laid on a slope, when the flow is
shut down, the liquid separates out and settles down quickly at the bottom of the slope
with the gas accumulating at the top of the slope. This unequal density distribution along
a slope changes the weight distribution along the pipe and enhances the walking of the
pipeline down the slope. During normal operating condition the pipe moves more in the
227
Structural analysis of submarine pipelines
downhill than uphill. Due to the settlement of the liquid down slope during shut down
process, the submerged weight reduces the uphill contraction of the pipe; the reverse
situation happens at the uphill light gas filled section of the pipe. This increases the
Figure 6-3 shows the split length, Ls and the whole length of the pipe, L. The axial
frictional force can be found out for various split lengths as:
L (6-6)
Sc = (1 + ) − (1 − )(3 + ) W0 for Ls L
L
2 2 2
L (6-7)
Sc = (1 + ) − (1 − )( − ) W0 for L Ls LB '
L
2 2 2
L L (6-8)
Sc = (1 + ) + 3 − + s (1 − ) − 1 W0 for Ls LB '
4 4
Ls LB'
L =
(Sc − P ) X ab (6-9)
EA
where, Xab is the distance between the virtual anchor points, LB′ is the unloading virtual
anchor point, ω, is the ratio of the gas filled submerged weight to the liquid filled
submerged weight (Wg /Wo), β, is the ratio of the operating submerged weight to the
liquid filled submerged weight (W/Wo). FE analysis was also carried out to verify the
analytical solution. Analytical and numerical results were matched with good agreement.
In this paper walking due to first three mechanisms were considered for the analysis and
further study. The walking due to slope with liquid hold up could be studied separately
228
Structural analysis of submarine pipelines
solution is described to estimate the end expansion of on-bottom pipeline resting on a flat
seabed with uniform resistance. This analytical expression is then verified with numerical
modelling. The pipe is considered to be elastic and resting on a perfectly flat seabed with
resistance expressed as friction coefficient. One end of the pipe is considered to be fixed
with other end remained free. Then the pipe is heated uniformly and the end expansion
was resisted by the friction between the pipe and the seabed. Then the numerical
verification of this model was carried out with ABAQUS (Dassault Systѐms, 2007). The
numerical analysis was carried out in both ways, once with a flat seabed with Coulomb
friction between the pipe and the seabed (rigid-plastic pipe-soil interaction, i.e.
mobilisation distance practically zero), and with nonlinear springs (elastic-plastic pipe-
soil interaction with a finite mobilisation distance) . Results of both the numerical analysis
Thereafter, the existing analytical solutions for various seabed slopes, SCR tensions and
thermal transients are examined with numerical solutions. The numerical results match
with very good agreement with the existing analytical solutions. However, the elastic-
plastic response of the soil was neglected at this stage. To incorporate the mobilisation
displacement of the soil into the existing methodology to predict walking of submarine
pipelines, a new analytical solution is developed. The proposed model is then verified
with numerical results. The effect of the mobilisation displacement of the soil on the
walking behaviour of the pipeline due to seabed slope, SCR tension and thermal transients
is then checked. The results of the numerical analyses are in good agreement with the
229
Structural analysis of submarine pipelines
A simple analytical solution is provided here to estimate the end expansion due to thermal
1. The pipeline remains elastic and the material properties are described by Young’s
2. The pipe can be treated as straight thin walled circular tube of thickness t and the
3. The pipeline is empty and there is no internal pressure acting inside the pipeline
4. The temperature of the pipeline is considered uniform for the present analysis.
profile was often assumed in the literature as decaying exponentially with distance
from the hot end. This effect nonlinear temperature profile will be discussed in
5. The left end of the pipe was considered to be fixed on the well and the right end
was considered to the free for expansion. This assumption was similar to analysis
The pipeline is assumed to be connected to the well head and can expand freely away
from the well head. The well head is at very high temperature. When the high temperature
oil or gas starts flowing through the pipeline it is subjected to a temperature gradient, and
due to this temperature raise Δθ, it tends to expand away towards the free end. Over the
length L, the pipeline was moving away from the well head, and so the bottom exerts on
the pipe a force f, per unit length, directed towards the well head and opposing the motion.
Due to the restriction of the movement, a compressive load P develops in the pipeline.
Figure 6-4 illustrates the schematics of thermal expansion of submarine pipeline. Hence,
230
Structural analysis of submarine pipelines
the total strain can be given by the addition of the thermal strain and the strain due to the
compressive load.
The thermal strain generated due to the uniform increase in thermal gradient is given by:
where α is the coefficient of thermal expansion of the pipeline material and Δθ is the
change in temperature.
𝜀𝑐 = − 𝑃⁄𝐸𝐴 (6-12)
dP (6-13)
=−f
dx
f = N = W (6-14)
where μ is pipe-soil of friction coefficient,W′ is the submerged weight of the pipe, and an
enhancement factor ζ, was introduced to account for wedging around the curved surface
of the pipe (White and Randolph, 2007). For simplicity the value of ζ is taken as unity.
P = − fx + C1 (6-15)
The axial load in the pipeline decays towards the free end. Therefore, at x = L , P = 0
C1 = fL (6-16)
P = f (L − x ) (6-17)
231
Structural analysis of submarine pipelines
𝐿
𝑓(𝐿 − 𝑥) (6-18)
𝛥𝑒𝑛𝑑 = ∫ (𝛼𝛥𝜃 − ) 𝑑𝑥
0 𝐸𝐴
Integrating, the maximum axial strain in the pipe due to end expansion can be expressed
as:
𝛥𝑒𝑛𝑑 𝑓𝐿
= 𝛼𝛥𝑇 −
𝐿 2𝐸𝐴
(6-19)
This is a dimensionless expression with two input groups. The first group, αΔT is the
thermal term and the second group, fL/2EA is the geotechnical term. This condition is
valid when the thermal term is greater than the geotechnical term, αΔT > fL/2EA.
Numerical finite element analysis (FEA) was carried out to support the analytical
modelling. The FEA was performed with commercial software ABAQUS (Dassault
Systѐms, 2007). The FEA method was performed in two different ways. In the first case
a straight elastic-plastic pipe resting on a rigid seabed was considered and the Coulomb
friction with rigid-plastic pipe-soil interaction, i.e. mobilisation distance practically zero
was applied between them. In the next case a different approach was followed and
nonlinear axial springs with elastic-plastic pipe-soil interaction with a finite mobilisation
distance were used to simulate the soil resistance in the axial direction.
The pipeline was modelled with pipe (PIPE31) element and the seabed was modelled with
analytical rigid element available in the ABAQUS (Dassault Systѐms, 2007) code. The
friction between the pipeline and seabed was modelled with friction force is proportional
to normal force. The general static method was used to simulate the heating of the pipeline
uniformly. Table 6-1 showing the material properties used for the pipeline. Static general
method was used to heat up the pipe uniformly from 0°C to 90°C. The left end of the pipe
232
Structural analysis of submarine pipelines
was assumed to be connected to the wellhead and fixed. The right end of the pipe was
kept free for expansion. The coefficient of friction was varied from 0.1 to 1 and the
corresponding expansion of the free end of the pipe was plotted against the coefficient of
friction.
In the second method the seabed was modelled with the only nonlinear axial springs. .
There is no change in the normal reaction force between pipe and seabed compared to the
rigid seabed case. Therefore, no vertical springs were used. The axial springs provided
the axial resistance to the expansion. The left end of the pipe was assumed to be connected
to the wellhead and fixed. The right end of the pipe was kept free for expansion. The
material property of the pipe was kept constant as before. Only the soil behaviour had
been changed in this case..Similar to the previous method the coefficient of friction in
this case can be implemented as, μ = f / W′ (Ultimate axial resistance / Weight of the pipe).
Therefore, μ was varied from 0.1 to 1 by varying the ultimate axial resistance of the soil.
The corresponding expansion of the free end of the pipe was plotted against the
coefficient of friction.
The results of end expansion of the pipe was plotted against the varying friction factor for
analytical and FE methods. Analytical solution and FE solution with rigid seabed with
Mohr-Coulomb friction model shows very close agreement. However, the difference
Figure 6-5 illustrates that end expansion with elastic plastic spring is slightly higher than
the rigid seabed case. This discrepancy can be explained by comparing the force profiles.
As very small mobilisation distance ( 0.001m) was considered, the end expansion for
233
Structural analysis of submarine pipelines
rigid seabed and elastic-plastic springs matched with good agreement with the analytical
solution.
Force profile
Figure 6-6 and Figure 6-7 delineate the force profiles of rigid seabed and elastic-plastic
springs respectively for various seabed frictions. Friction factor was varied from 0.1 to 1
keeping other parameters constant. The load in the pipe reduces along the length from the
fixed end towards the free end. With increasing friction the load in the pipe also increases.
Distributed displacement
Figure 6-8 and Figure 6-9 show the distributed expansion along the length of the pipeline
due to rigid seabed and elastic-plastic springs for various resistances. Friction factor was
varied from 0.1 to 1. It is depicted that, in both the cases expansion of the pipe reduces
with increasing friction between pipe and soil. The expansion at the free end was found
out to be high due to less friction mobilisation at the end. The elastic component means
that lower friction is mobilised near the fixed end so the overall average compression
Three cases were compared here with friction factors, μ =0.1, 0.5 & 1. Figure 6-10,Figure
6-11and Figure 6-12 show the force profiles of these three cases. As the elastic-slip of the
springs were chosen very small, the force profiles for first two values of friction factors
matched with good agreement. However, small deviation was observed for higher friction
factor near the fixed point. This is due to the presence of elastic-slip the load on the pipe
reduced from rigid case. This phenomenon was explained in detail later.
234
Structural analysis of submarine pipelines
Numerical analysis is described in this section to reproduce the existing theories and to
give an insight to the existing practiced methodologies to estimate pipeline walking rates
without recourse to numerical analysis. As mentioned earlier, walking due to SCR tension,
seabed slope and thermal transients have been studied. Force profiles with a range of
friction factors are studied and compared along with end displacements for each of the
insight to the walking behaviour of the pipeline subjected to these three triggering
mechanisms. At the end present numerical results were compared with the existing
The walking rate per cycle due to SCR tension was also verified numerically in this
section. Figure 6-13 shows the idealisation of the SCR tension (load) applied at the end
Force profile
Figure 6-14 and Figure 6-15 depicts the force profiles with the minimum and maximum
considered values of friction respectively. The force envelope with the maximum friction
is more symmetrical than the minimum friction. Also the maximum load in the pipe
End displacement
The end displacement of the loading end was also plotted for maximum and minimum
friction and for different SCR tensions. Figure 6-16 shows the comparison between end
displacements for SCR tension 100 kN and 500 kN and with the minimum friction. The
walking rate per cycle increased with increasing SCR tension. Figure 6-17 shows the
comparison between end displacements for SCR tension 100 kN and 500 kN and with the
235
Structural analysis of submarine pipelines
maximum friction. Figure 6-18 shows the comparison of end expansion over the thermal
cycle for the maximum and minimum friction with a SCR of 500 kN. The walking rate
Distributed displacement
To understand the walking mechanism of the whole length of the pipe, the distribution of
displacement at various stages in the analysis were plotted. The location of selected nodes
and the SCR tension are shown in Figure 6-19. Figure 6-20 depicts the end displacements
plotted against the thermal steps along the length. It is seen that expansion and contraction
of the nodes are associated with the heat up and cool-down steps. Five heat-up and cool-
down steps were considered here. The points P0 and P643 which are on the left side of
middle point were always moved towards the SCR tension and P1359 & P2000 moved,
which are on the right side of the middle point, moved towards the right during heat up
and towards the left during cool-down. Therefore, the cumulative walking was always
Existing analytical models have been compared with the present numerical results. Input
parameters were tabulated in Table 6-1. Figure 6-21 depicts the comparison of analytical
and numerical results for the maximum and minimum friction factors. The walking rate
is examined for a range of SCR tension. The walk per cycle is inversely proportional to
the friction factor and directly proportional to the SCR tension, keeping other parameters
remain constant. These numerical results are in good agreement with the existing
236
Structural analysis of submarine pipelines
Similar to the previous case, the walking rate per cycle due to seabed slope was also
verified numerically in this section. Figure 6-1 shows the idealisation of the problem
when a pipe is resting on a slope. A real slope was not modelled here, but the numerical
analysis was carried out by adjusting the direction of the gravity load to be equivalent to
Force profile
A seabed with 2° slope was considered here. Figure 6-22 and Figure 6-23 depict the force
profiles with the minimum and maximum friction factors respectively. The force
envelope with maximum friction was more symmetrical than with the minimum friction.
Also the maximum load in the pipe reached a higher value with the maximum friction.
End displacement
The end displacement of the loading end was also plotted for maximum and minimum
friction and for different seabed slopes. Figure 6-24 shows the comparison between end
displacements for seabed slopes of 1°, 2° and 3° and with minimum friction. The walking
Distributed displacement
Similar to the previous case of SCR tension, the distribution of pipe displacement at
various stages in the analysis were plotted. The locations of selected nodes and slope
idealisation were shown in Figure 6-25. Figure 6-26 depicts the end displacements plotted
against the thermal steps along the length. It is seen that, expansion and contraction of the
nodes are associated with the heat up and cool-down steps. Five heat-up and cool-down
steps were considered here. The points P0 and P477 which are on the left side of middle
point, p1000 were moved towards the downhill and points P1542 & P2000 moved,
237
Structural analysis of submarine pipelines
towards uphill during heat up and towards downhill during cool-down steps. Therefore,
The walk per cycle calculated analytically using the method described by Carr et. al.
(2006) has been compared with the numerical results. Input parameters were tabulated
in table 2. Figure 6-27 depicts the comparison of analytical and numerical results for the
maximum and minimum friction factors. The walking rate is examined for a range of
seabed slopes. The walk per cycle is inversely proportional to the friction factor and
directly proportional to the seabed slope, keeping other parameters constant. The present
numerical results are in good agreement with the existing analytical solution (Carr et. al.
2006).
Similar to the previous two walking triggering mechanisms, thermal transients were also
implemented in the analysis, as indicated by the examples shown in Figure 6-2. However,
simpler thermal transients can be used for numerical analysis (Carr et al. 2006). Figure
6-28 shows the type of heat up thermal transients used for this study. As discussed earlier,
cool down is a uniform heat loss process, therefore, a uniform temperature change along
the length of the pipe was considered for the cool down steps in the numerical modeling.
Force profile
The force profile of the heat up step was plotted for a thermal gradient of 20 °C/km. It
was the second and subsequent cycles that establish the steady walking behaviour of the
pipeline (Carr et al. 2006), so the force envelope of the second heat-up step is plotted in
Figure 6-29. The effects of thermal transients were shown by the intermediate force
profiles and an arrow showing the shift of force profile from cooling to a heating up step.
238
Structural analysis of submarine pipelines
Similarly intermediate steps of cool down are also plotted in Figure 6-30. The cool down
is typically at a uniform rate along the whole length of the pipe, leading to contraction
about the mid-line virtual anchor point. The pipeline unloads symmetrically about the
Thermal gradient applied to a pipeline have a significant effect on walking rate. . In the
model the left end is fixed and the right end of the pipe is free. Figure 6-31 illustrates the
walking rate as function of the axial friction force normalised by the constraint friction,
f* (f* = ∆P/L) at which cyclic constraint occurs, where ΔP is the fully constrained force,
and L is the length of the pipe. Three different thermal gradients were used to compare
the walking rate between existing analytical expression (Carr et. al. 2006) and the present
numerical solution. The results of the numerical analysis and the analytical solution match
with good agreement. The amount of walking is strongly influenced by the thermal
gradient, i.e. the walking rate of 30°C/km is 3 times that of 10°C/km for same f/f*.
Having verified that the numerical analysis being performed in this work is accurate,
through comparison with existing analytical solutions for pipeline walking, the numerical
the axial friction model. A bi-linear axial response was used in FE analysis, with the
mobilization displacement being the displacement required to mobilise the full axial
pipelines has been discussed by various authors (Tørnes et al. 2000, Carr et al. 2006,) and
239
Structural analysis of submarine pipelines
explored numerically and more recently analytically (Wang et. al. 2010). For example,
Wang et al. (2010) explored the effect of mobilisation displacement on long pipelines and
effect of mobilisation displacement on the walking of short pipelines has been developed
to date. This section gives a detail insight this behavior and thereafter proposes an
6.5.2 Approach
Firstly, numerical analysis was carried out to see the effect of mobilisation displacement
on the walking. Different mobilisation displacements with same ultimate friction were
used to compare the results. The force profiles were also compared. Thereafter, an
analytical model was developed to analyse the effect of mobilisation displacement on the
force profile and hence on the walking. The next two short sections give the detail
submarine pipelines resting on a seabed slope, and subjected to thermal cycles, a detailed
0.001 m, 0.05 m and 0.1 m were considered and are shown in the Figure 6-32. An elastic
perfectly plastic Coulomb friction model was chosen as discussed earlier. The pipe
material, size and length were kept constant as earlier cases. The seabed was considered
with a slope of, φ = 3°, and a minimum friction factor of 0.1 was considered to recover
The expansion at the down slope end of the pipe is plotted against cycle number for
different mobilisation displacements in Figure 6-33. This shows that with increasing
240
Structural analysis of submarine pipelines
The walking rate of the pipeline for various seabed slopes is also plotted for different
mobilisation displacement on the walking behaviour of the pipeline, the force profiles
during heat up and cool down steps are plotted for various mobilisation displacements,
uult = 0.001, 0.1, 0.2, 0.3, & 0.456m, with seabed slope, φ = 3°, and friction of 0.1 for the
Figure 6-35 shows the force envelopes of the pipeline for heat up and cool down steps
from very low to very high mobilisation displacement. The virtual anchor points show a
very sharp edge with uult = 0.001 m, and the crown of the force envelop became parabolic
with increasing uult. The sharper peak is due to the full mobilization of the friction along
the length of the pipeline. The phenomenon was analytically explained in the following
section.
Figure 6-36 shows the result of a force profile of a pipeline with, EA = 1.9E+10N, L = 2000
m, φ = 3°, uult = 0.456 m, ∆θ = 90°C. Only the cool down step is plotted here. Only half
of the force envelop was considered here, as this is a symmetrical case. The force profile
has been divided into four parts and the force profile for these four parts is derived to
241
Structural analysis of submarine pipelines
du P (6-20)
=
dx EA
dP W (6-21)
= Fpassive = uk = u
dx u ult
d 2P W (6-22)
2
= P = 2 P
dx EAuult
W
where, 2 = , note, λ has two values as, λup and λdown.
EAuult
and
dPdown ( x)
dx
(
= down Adown e down x − Bdown e −down x ) (6-25)
P = xu upW (6-26)
dPup ( x) (6-27)
= upW
dx
242
Structural analysis of submarine pipelines
From point x = xd :
P = (L − xd ) downW (6-28)
dPdown ( x) (6-29)
= − downW
dx
dPup ( x) (6-30)
=0
dx
dPdown ( x) (6-31)
=0
dx
xm (6-32)
= down
L − xm up
and the distance between two unloading virtual anchor points, xd and xm can be
predicted by:
u ult EA (6-33)
xd − xm =
(P − S f )
where, ∆Sf = µdownWL, the length of the unloading virtual anchor point without
mobilisation displacement is, xm, and the length of the unloading virtual anchor point with
Solution scheme:
243
Structural analysis of submarine pipelines
Bup 2 up x0
=e
Aup
(6-34)
Bdown (6-35)
= e 2down x0
Adown
upW 1 (6-36)
Aup =
up e up xu
−e
up ( 2 x0 − xu )
upW e
2 up x0
(6-37)
and, Bup = up ( 2 x0 − xu )
up e up xu
−e
downW 1 (6-38)
Adown =
down e down xu
− e down (2 x0 − xd )
downW e 2down xd (6-39)
and, Bdown =
down e down xu
− e down (2 x0 − xd )
Again, it can be shown that the new anchor length, X′ab is linked to the previous anchor
'
X ab = X ab − 2(xd − x m ) (6-40)
Therefore, the new walk per cycle can be found out by:
EA
Using equation (6-40), which is the rigid plastic walking rate solution, we get;
EA EA (P − S f )
Therefore, the ultimate walking per cycle for the elastic-plastic friction model can be
given by:
244
Structural analysis of submarine pipelines
The first term in the left hand side is the walk per cycle due to slope for a rigid plastic
seabed (i.e. without a mobilisation displacement) and the second term is the twice of the
mobilisation displacement. The same soil elastic stiffness was used for loading and
unloading cases, meaning that an increase of uult in the mobilisation displacement reduced
The above expression can be used to predict the walking of submarine pipelines, taking
into account the recoverable elastic component of the mobilisation displacement. Further
cases were analysed to explore how the mobilisation displacement affects the walking of
Numerical modelling was carried out to analyse the walking affected by the mobilisation
displacement for SCR tension and thermal transients, in addition to the case of seabed
slope that is analysed above. Analytical and numerical results were then compared to
verify the exactness of the proposed model, when applied to both the slope case (derived
above) and also the other cases. All the three walking mechanisms were studied here. The
material property of the pipe and the heating range were kept constant as given in Table
6-1.
The results in all cases showed that the correction by 2uult of the rigid-plastic solution
gives the correct walking rate for the elastic plastic friction model. Therefore, the solution
that is proven analytically above for the seabed slope case, is also applicable to the
Figure 6-39 illustrates the comparison between predicted and numerical results of walk
per cycle. The mobilisation displacements were varied from 0.001 m to 0.456 m (0.5D)
245
Structural analysis of submarine pipelines
to investigate the validity of the model. Predicted and numerical results match with good
Figure 6-40 shows the comparison between predicted and numerical results of walk per
cycle. The mobilisation displacements were varied from 0.001 m to 0.456 m (0.5D) to
investigate the validity of the model. Predicted and numerical results match with good
Figure 6-41 depicts the comparison between predicted and numerical results of walk per
cycle. The mobilisation displacements were varied from 0.001 m to 0.456 m (0.5D) to
investigate the validity of the model. Predicted and numerical results match with good
6.6 CONCLUSIONS
This chapter bridges the gap between the present analytical design tool for pipeline
walking towards a more realistic soil response. The existing analytical methods to predict
pipeline walking have been used to verify a numerical model and an analytical solution
has been developed to incorporate the elastic-plastic response into the existing solution
The solution was also shown to apply to walking due to SCR tension and thermal
transients based on full numerical analysis. A new expression was therefore proposed to
predict the walking rate on elastic-plastic soil response. This is proven analytically for the
slope case and shown numerically to apply for other cases, as a modification of the rigid-
plastic solutions.
246
Structural analysis of submarine pipelines
References:
Bruton, D., Sinclair, F., & Carr, M. (2010) . Lessons learned from observed walking of
pipelines with lateral buckles, including new driving mechanisms and updated
Paper 20750
Bruton, D., White, D. J., Carr, M., & Cheuk, J. C.Y. (2008). Pipe-soil interaction during
Carr, M., Bruton, D., & Leslie, D. (2003). Lateral buckling and pipeline walking, a
Amsterdam.
Carr, M., Sinclair, F., & Bruton, D. (2006). Pipeline Walking – Understanding the filed
layout challenges, and analytical solutions developed for the SAFEBUCK JIP.
Collberg, L., Carr, M. and Levold, E. (2011). Safebuck Design Guideline and DNV-RP-
OMAE98-1103
247
Structural analysis of submarine pipelines
Tørnes, K., Ose, B.A., Jury, J., Thomson, P. 2000. Axial Creeping of High Temperature
Pipelines caused by Soil Ratcheting. Proc. of 19th Int. Offshore Mech. and Arctic
Wang, Y. N., Maschner, E. A., & Hayes, R. (2010). A global migration risk to partially
White, D. J., & Randolph, M. F. (2007). Seabed characterisation and models for pipeline-
193-204.
Tables:
Diameter, D 0.912 m
Length, L 2000 m
Temperature rise, Δθ 90 °C
248
Structural analysis of submarine pipelines
Figures:
Figure 6-2: Example of pipeline heat up thermal transient profiles (Carr et al. 2003)
249
Structural analysis of submarine pipelines
250
Structural analysis of submarine pipelines
2.5
1.5
1
D = 0.912 m (36")
EA = 1.96 E+10 N
Analytical
0.5 Δθ = 90 °C
Rigid seabed α = 1.2 E-05 /°C
L = 2000 m
Elastic-plastic spring W′ = 6.9 kN/m
0
0 0.2 0.4 0.6 0.8 1 1.2
Friction factor, μ
Figure 6-5: Comparing results of analytical and FE analysis for end expansion of pipe
16
Rigid seabed D = 0.912 m
W' = 6.9 kN/m
μ = 1, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1 EA = 1.91 E+04 MN
12
Axial load, P (MN)
μ = 0.5
μ = 0.1
0
0 500 1000 1500 2000
Along the length, x (m)
Figure 6-6: Axial load distribution along the length of the pipeline for various friction
251
Structural analysis of submarine pipelines
16
Elastic-plastic spring D = 0.912 m
W' = 6.9 kN/m
EA = 1.91 E+04 MN
μ = 1, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1
Effective axial load, P (MN)
12
μ = 0.5
μ = 0.1
0
0 500 1000 1500 2000
Along the length, x (m)
Figure 6-7: Axial load distribution along the length of the pipeline for various friction
2.5
D = 0.912 m
W' = 6.9 kN/m
2 EA = 1.91 E+04 MN
Axial displacement, u (m)
1.5
0.5 μ = 1, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1
Rigid seabed
0
0 500 1000 1500 2000
Along the length, x (m)
Figure 6-8: Axial displacement along the length of the pipeline for various friction
252
Structural analysis of submarine pipelines
2.5
D = 0.912 m
W' = 6.9 kN/m
EA = 1.91 E+04 MN
2
Axial displacement, u (m)
1.5
0.5 μ = 1, 0.9, 0.8, 0.7, 0.6, 0.5, 0.4, 0.3, 0.2, 0.1
Elastic-plastic spring
0
0 500 1000 1500 2000
Along the length, x (m)
Figure 6-9: Axial displacement along the length of the pipeline for various friction factors
2
μ=
= 0.1 D = 0.912 m
W' = 6.9 kN/m
EA = 1.91 E+04 MN
Effective axial load, P (MN)
1.5
0.5
Rigid seabed
Elastic-plastic spring
0
0 500 1000 1500 2000
Along the length, x (m)
Figure 6-10: Comparison of force profiles between rigid seabed and elastic-plastic
253
Structural analysis of submarine pipelines
8
μ = 0.5 D = 0.912 m
W' = 6.9 kN/m
EA = 1.91 E+04 MN
Effective axial load, P (MN)
2
Rigid seabed
Elastic-plastic spring
0
0 500 1000 1500 2000
Along the length, x (m)
Figure 6-11: Comparison of force profiles between rigid seabed and elastic-plastic seabed
for, µ = 0.5
16
μ= 1 D = 0.912 m
W' = 6.9 kN/m
EA = 1.91 E+04 MN
12
Axial load, P (MN)
Rigid seabed
Elastic-plastic spring
0
0 500 1000 1500 2000
Along the length, x (m)
Figure 6-12: Comparison of force profiles between rigid seabed and elastic-plastic
seabed for, µ = 1
254
Structural analysis of submarine pipelines
Figure 6-14: Force profile with SCR Tension and for minimum friction
255
Structural analysis of submarine pipelines
Figure 6-15: Force profile with SCR Tension and for maximum friction
Figure 6-16: Comparison of end expansion for different SCR tensions with μ = 0.1
256
Structural analysis of submarine pipelines
Figure 6-17: Comparison of end expansion for different SCR tensions with μ = 0.6
Figure 6-18: Comparison of end expansion for single SCR tensions with different
friction factors
257
Structural analysis of submarine pipelines
Figure 6-19: Location of various nodes selected for distributed expansion analysis for
SCR tension
Figure 6-20: Distributed expansion at various location of the pipe for SCR tension 500
258
Structural analysis of submarine pipelines
Figure 6-21: Comparison of analytical and FE solution for walking with SCR tension
Figure 6-24: Comparison of end expansion for different seabed slope and single
260
Structural analysis of submarine pipelines
Figure 6-25: Location of various nodes selected for distributed expansion analysis for
seabed slope
Figure 6-26: Distributed expansion at various location of the pipe for seabed slope of 3°
261
Structural analysis of submarine pipelines
Figure 6-27: Comparison of analytical and FE solution for walking with sloping seabed
Figure 6-29: Force envelope due to thermal transients - second heat up step
Figure 6-30: Force envelope due to thermal transients – second cool down step
263
Structural analysis of submarine pipelines
Figure 6-31: Comparison of analytical and FE solution of pipeline walking with thermal
transients
264
Structural analysis of submarine pipelines
Figure 6-33: Effect of mobilisation displacement on a pipeline with sloping seabed and
minimum friction
265
Structural analysis of submarine pipelines
266
Structural analysis of submarine pipelines
1.2
Φ = 3ᵒ
µ = 0.1, 0.2, 0.3 L = 2000 m
1 EA = 1.91E+04 MN
W' = 6.9 kN/m
Walk per cycle, ∆ t (m) 0.8
Analytical solution
(solid lines)
0.6
0.4
0.2 FE results
0
0 0.1 0.2 0.3 0.4 0.5
Axial mobilisation displacement, uult (m)
Figure 6-37: Comparison of analytical and numerical results of pipeline walking with
mobilisation displacement
1.5
µ = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6 φ = 3ᵒ
Walk per cycle, ∆φ (m/cycle)
L = 2000 m
EA = 1.91E+04 MN
A W' = 6.9 kN/m
1
B
0.9
0.8
0.5
0
0 100 200 300 400 500
Figure 6-38: Chart to find out the effect of mobilisation displacement on pipeline
267
Structural analysis of submarine pipelines
0.7
Actual total walking per cycle, ∆ ult
0.4
0.3
0.1
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
Predicted walking per cycle, ∆ ult (m/cycle)
Figure 6-39: Verifying the effect of mobilisation displacement on walk per cycle due to
SCR tension
1.5
Actual total walking per cycle,
∆ult = ∆φ -2uult
∆ ult (m/cycle)
Figure 6-40: Verifying the effect of mobilisation displacement on walk per cycle due to
seabed slope
268
Structural analysis of submarine pipelines
0.05
10°C/km
∆ult = ∆T -2uult
0.02
Figure 6-41: Verifying the effect of mobilisation displacement on walk per cycle due to
thermal load
269
Structural analysis of submarine pipelines
270
Structural analysis of submarine pipelines
7.1 INTRODUCTION
During the operational life of a seabed pipeline, it undergoes many start-up and shut down
cycles leading to pipe-soil interaction forces that control the expansion and accumulated
movement. Accumulated axial movement due to repeated thermal cycles may lead to
global displacement, referred to as pipeline walking, and the pipeline may be designed to
buckle laterally to relieve the internal loads. The axial pipe-soil interaction forces have
significant influence on the design for walking and lateral buckling. Even though a
significant amount of experimental data on axial friction has been available (White et al.
2011), the equivalent friction coefficient for structural analysis of pipeline is often
selected without considering the significant effect of the pipeline velocity or history of
coefficient is introduced. The model originates in the study of friction between solid
materials. It has been established in tribology that the friction coefficient that opposes the
initiation of slipping from a sticking condition is different from the friction coefficient
that opposes established slipping (Rabinowicz (1965), Moore (1975) and Bhushan (1999))
The former is typically referred to as the ‘static’ friction coefficient, and the latter is
higher than the kinetic friction coefficient. One model to capture this behaviour assumes
271
Structural analysis of submarine pipelines
that the friction coefficient rates exponentially from the static to the kinetic value
where μk is the coefficient of kinetic friction, μs is the coefficient of static friction, v is the
relative velocity of the surfaces, or the slip rate, and a is the rate factor (Oden & Martins.
(1985)). When a = 0, μ = μs . With increasing a, the transition between the limiting values
coefficient on clay soils was developed Randolph et al. 2012 and matched with the
𝑇 𝑚
𝐹 𝐹 𝐹 𝐹 𝑙𝑛(2)( )
= ( ) − [( ) − ( ) ] 𝑒 𝑇50
𝑊 ′ ′
𝑊 𝑑𝑟𝑎𝑖𝑛𝑒𝑑 ′
𝑊 𝑑𝑟𝑎𝑖𝑛𝑒𝑑 ′
𝑊 𝑢𝑛𝑑𝑟𝑎𝑖𝑛𝑒𝑑
(7-2)
where F is the frictional sliding resistance, W is the weight of the pipe, T is the time
and T50 is the time at which the frictional resistance is 50% of the undrained and drained
limits. Here, the value of T50 is close to 0.05, and m = 0.5. T can be calculated as T =
cvt/D2, where cv is the consolidation coefficient for the soil and D is the diameter of the
pipe. This model captures the observed trend for low undrained friction to be mobilised
for high speed or short duration axial pipe movements, and a high drained friction to be
In the present analysis a similar model to those proposed by White et al. (2011) and
Randolph et al. (2012) has been applied to the resistance between a pipe and the seabed
during axial pipeline motion. The parameter μk is replaced by μLB, the apparent
272
Structural analysis of submarine pipelines
drained conditions. The form of the relationship has been kept constant, but the
where the term v50 is the velocity at which the friction factor lies 50% of the way between
the limiting values. The velocity parameter, v50, depends on the drainage behaviour of the
soil and can be estimated from laboratory tests and varies from soil to soil. Figure 7-2
shows the relation between friction factor and velocity given by equation (7-2); this model
is referred to throughout this paper as the friction rate model. The position of v50 is
indicated in the diagram. The constant λ is equal to ln 2 (0.693) because e-ln2 = 0.5.
7.2 OBJECTIVE
The present design practice for analysing pipeline walking behaviour was discussed in
the previous paper, and the conventional elastic-plastic pipe-soil response was analysed.
In conventional analyses, the rate of thermal loading during operating cycles is not
considered within the analysis. However, the velocity varies with position along the
pipeline and with time during the startup or shutdown event. Recently, it has been
recognised that the sliding resistance between a pipe and the seabed varies with velocity
due to drainage effects ( White et al. (2011) and Randolph et al. (2012)) but to date this
This paper develops a numerical model in which dynamic walking behaviour is studied
ABAQUS and validated via single element and simple (flat seabed) pipeline cases. This
model is then used for a parametric study exploring walking driven by seabed slope and
SCR end tension. The effect of differing rates of heating and cool-down on walking
273
Structural analysis of submarine pipelines
behaviour are also investigated. Finally, the potential to express the velocity-dependent
It has been recognised that the sliding resistance between a pipe and the seabed varies
To establish more realistic pipe-soil interaction White & Cathie (2010) and Hill et al.
(2012) performed the time dependent pipe-soil friction model tests. The tests revealed
that axial resistance is strongly influenced by the rate of movement and the pause period
between the movement events. The higher peak values of friction factors were associated
with longer waiting periods between axial sweeps and the lowest residual values were
The effect of velocity on interface shear strength is shown in Figure 7-3 (Hill et al.,2012).
The response of fine-grained soil during shearing varies between fully drained and fully
undrained conditions across a wide range of velocities ( White and Cathie (2010)). Data
from previous research by Hill et al. (2012) and Steelfelt (1993) have been summarised
in the figure, with τres the mean shear stress on the pipe surface curing steady motion at
velocity v and σ′no the initial normal effective stress. These results are all from clay
samples with an over consolidation ratio (OCR) of 1, where the OCR is the ratio of the
previous maximum effective stress to the current effective stress. A value of OCR = 1 is
typical for soft clay seabeds, particularly when disturbed and then loaded by the process
of laying and then operating a pipeline. For a high velocity, partially or fully undrained
conditions will prevail, and positive or negative excess pore pressure will be generated
274
Structural analysis of submarine pipelines
(White et al. 2012). For these normally and lightly over consolidated soils, the drained
friction represents an upper bound (UB) and the undrained friction represents and lower
bound (LB).
This paper applies a more accurate model of the response than previous studies of pipeline
walking, via a friction rate model that links these two limits.
Figure 7-4 represents the previous data fitted with the friction rate model. The actual data
is in good agreement with the friction rate model. However, in each case the v50 parameter
was varied to fit the data for particular drained and undrained limits. The data show that,
with increasing sliding rate, undrained (lower bound) friction ratio prevails.
Recently Carneiro et al. (2017) carried out numerical analyses using the quasi-static
method in the commercial software ABAQUS to investigate the effect of the rate model
studied with a flat seabed and a single pipe element, where the pipe was moved back and
forth with certain velocities. The friction factors were calculated from the output and were
compared with the predicted outcomes. This exercise confirmed that the friction rate
Thereafter, a length of pipe was considered, where the pipe was heated-up very fast and
cooled down slowly to show the operation of the rate model on a full pipeline. This was
first performed for a flat seabed case, and the variation in friction along the pipe length
275
Structural analysis of submarine pipelines
The rate model was then used to show the difference in walking behaviour for various
driving mechanisms. The effect of the rate model on walking due to SCR tension, seabed
slope and thermal transients was studied in detail. In each case an attempt has been made
to express the complex outcomes from the rate model in a simpler way by plotting the
force profiles and walking rates and thereby comparing them with the response for
conventional friction models based on the undrained and drained friction limits.
Finally, a parametric analysis was carried out to evaluate the walking rate for different
rate models. Once the walking rate was derived from the numerical analysis, a back
calculation was made to find out the equivalent velocity-independent friction factor for
each of case. The aim of this step is to illustrate to what extent the velocity-dependent
coefficient.
Firstly, the relative rate of heating and cooling was considered with a dimensionless time
tH
Thc = (7-4)
tC
where, tH is the time of heating and tC is the time of cooling. Throughout this analysis Thc,
was considered to be 0.01 unless otherwise specified. The heating of the pipe is rapid,
reflecting a sudden flooding with product, whereas the cooling is slower, reflecting heat
Another dimensionless group was formed by dividing the characteristic velocity of the
pipe expansion by v50. This group is essentially the ratio between (twice) the velocity of
the ends of the pipe if the thermal expansion is unconstrained, and the characteristic
( )L
Vp = m/s (7-5)
t H tC
where α is the coefficient of thermal expansion of the pipe material, ∆θ is the temperature
change, L is the length of the pipe. In most of these analyses, values of the parameters
These give Vp = 0.00216 m/s. Although the absolute values of tH and tC are impractically
small, the corresponding dimensionless velocities span the range between drained and
( )L / t H t C (7-6)
=
v50
Characteristic strain
The characteristic strain in the pipeline, c, is the ratio between the total pipe weight, WL,
which controls the force profile in the pipe through the mobilised friction coefficient, and
277
Structural analysis of submarine pipelines
W' L
c = (7-7)
EA
Through most of this paper the dimensionless strain was kept constant at 0.00072, based
steel pipe of outer diameter, D = 0.912 m (36″) and wall thickness, t = 0.033 m.
The characteristic strain is related to the mechanical strain that shortens the pipeline when
the friction is full mobilised. In this case, ignoring the small central region where the
limiting friction is not fully mobilised, for a constant friction coefficient , the force at
the mid-point of the pipe is WL/2 and the average force along the pipe is WL/4. The
consequent strain in the pipe is c/4 and the overall mechanical shortening is cL/4, or
cL/8 at each end relative to the centre. This mechanical shortening counteracts the ‘free
Figure 7-5 illustrates the model used for the benchmarking case with the friction rate
model in ABAQUS . A single element beam model, one metre in length, was used for the
pipe and a rigid surface element was used for the seabed. Gravity load was applied to
settle the pipe on the seabed and then dynamic steps were used to pull the pipe with
specified velocities.
Two cases A10 and A100 were studied with different decay factor a = 10 and 100 s/m
resulting in values for v50 = 0.069 m/s and 0.0069 m/s respectively, together with
uult = 0.01 m. The mobilised friction factor, μ was extracted from the results of the
analyses as μ = τ /σ′, where τ is the shear stress in the axial direction and σ′ is the normal
pressure on the seabed. The numerical results were then compared with the theoretical
278
Structural analysis of submarine pipelines
model response. A pipe with diameter, D = 0.912 m (36) with D/t = 25 and axial stiffness,
EA = 19100 MN was used for this analysis. A yield stress, σy = 410 MPa was also
specified for the pipe element, but the yield stress was not reached in this analysis and the
For the soil, throughout this study, an undrained friction factor, representing a lower
bound (LB), of μUD = 0.1 was adopted, along with a drained friction factor, representing
an upper bound (UB) of μD = 0.6. The input velocities and the rate (and rate, v50) factors
were tabulated in Table 7-1. The velocity was slowly ramped up to 10 m/s and then
reduced by a factor of 10 in six steps to 0.00001 m/s and then again increased back to 10
m/s for both cases. Every time step was of equal duration (10 s) for simplicity. Figure 7-6
Figure 7-7 shows the comparison between analytical model and the numerical results.
The velocity was varied from a very high value of 10 m/s to a low value of 0.00001 m/s
and then again increased from the low value to the high value. By doing this a wide range
of velocities were covered in order to judge the accuracy of the friction rate model. The
calculated equivalent friction factor matches closely with the theoretical expression,
At very large pipe velocities, such as those considered here, viscous effects become
important, and the friction coefficient will increase again as the shear strength of clay
For the second benchmarking case the same input parameters were chosen for a 2 km
long pipeline resting on a flat seabed. The pipe is assumed to have free ends, therefore, a
279
Structural analysis of submarine pipelines
The times of heating and cooling were varied to investigate the effect of velocity and the
resulting velocity dependent friction on the end expansion of the pipeline. For this
analysis the dimensionless time Thc was maintained at 0.01, with the time of heating taken
as tH = 10 s, and hence the time of cooling was tC = 1000 s, with 5 full cycles simulated
starting with heating by 90° and then cooling by the same amount. A waiting time of twait
= 100 s was introduced following each heating or cooling stage to stabilise the unbalanced
forces present at the end of the very fast heating. For these input parameters and taking
The dimensionless characteristic strain parameter, εc was kept constant at 0.00072. Cases
using the drained (UB) and undrained (LB) friction factors as well as the rate model were
used. The results of all of these models are compared in the following sections. Two key
outputs are the end movements and the mobilised friction at the free ends of the pipeline
which are plotted against time. In addition, the force envelopes and the mobilised friction
along the length of the pipeline at the end of each change in temperature are also compared.
Figure 7-8 shows the end displacement of the pipe at x = 0 over the five heating and
cooling cycles for different friction models, reflecting the net expansion of the pipeline.
Firstly the end displacements are plotted for both drained and undrained friction factors.
The drained friction factor leads to smaller expansion of the pipe and a lower cyclic range,
strains, and the greater mechanical extension in response to the contraction. After the first
The end expansion in response to the rate model is superimposed on the plot. The
response matches the drained solution during cool down – when the velocity is slow, and
the mobilised friction is close to the drained – and matches the undrained solution at heat
280
Structural analysis of submarine pipelines
up step – when the velocity is fast. These effects are also evident in the effective friction
Figure 7-10 shows the comparison of force profiles between the drained, undrained and
rate models and Figure 7-11 shows the corresponding variation in mobilized friction
along the pipe. The friction profile was calculated by dividing the change in force along
These results show that during the slow cool down process the mobilised friction along
the entire length of the pipeline is similar to the drained value, resulting in an almost linear
variation in force. However, the more rapid heat up process causes velocities in the range
over which the friction varies significantly. The result is a non-linear variation in force
along the pipe. The profile matches the undrained friction case near the fast-moving pipe
ends, but the mobilised friction is higher near the center of the pipe, where the velocity is
lower.
Figure 7-12 shows an enlarged view of the central portion of the pipe during cool down
steps. The drained (UB), undrained (LB) and rate model created a zig-zag portion at the
center of the pipe. These are due to the presence of mobilisation displacement in the
down and heating up the pipe experiences inward and outward displacements of the pipe
respectively. However, because of the previous heating cycle, the cooling stage does not
give sufficient inward displacement near the center of the pipe to reverse the sense of the
friction.
A comparative study is presented in this section to investigate the effect of the rate model
on steel catenary riser (SCR) tension. The force profiles, mobilised friction (expressed as
281
Structural analysis of submarine pipelines
the ratio of change in force along the length to the weight of the pipe) and walking rates
particular non-dimensionalised terms ψ = 0.0311, Thc = 0.01 and c = 0.00072. For this
case SCR tension of 500 kN at x = 0 was considered, maintaining a free end at x = 2000 m.
Force profile
Figure 7-13 shows the comparison of force envelopes for drained (UB), undrained (LB)
and rate friction models. The force profiles of the drained and rate models at the cool
down step are comparable. However, the maximum force reached with the rate model is
slightly lower than the drained value. This was due to the selection of total time for the
cool down step. The maximum value reached with the rate model in the heating step was
higher than the undrained value. However, the slopes of the undrained and rate models
Friction profile
Friction profiles are plotted in Figure 7-14. Mobilised friction was calculated from the
force profile as previously by (dP/dx)/W′. The drained (UB) and undrained (LB) values
were 0.6 and 0.1 respectively and the resulting force profiles are shown by the solid lines.
However, the results of the rate model match closely with the drained response over most
of the pipe length for the cool down step but with the undrained response towards each
Walking rate
The walking rate over the five cycles is compared for the drained (UB), undrained (LB)
and rate model in Figure 7-15. For the first cycle the walk rate was similar for all the
friction. However, the effects of friction factors were more prominent after the second
282
Structural analysis of submarine pipelines
cycles due to the stabilisation of the model from the second step onwards. The walking
rate with the undrained friction was much higher than for the drained case while results
for the rate model fall in between the drained and undrained responses.
A comparative study is presented in this section to investigate the effect of the rate model
on seabed slope with an angle Φ. Similar to previous section force profiles, mobilised
friction and walking rates were compared for drained, undrained and rate models keeping
ψ = 0.0311, Thc = 0.01 and c = 0.00072. For this case a seabed slope angle of Φ = 3° was
Force profile
Force profiles for drained (UB), undrained (LB) and rate models are compared in Figure
7-16 for the seabed slope case. The maximum force reached with the rate model during
cool down is slightly lower than for the drained case. The maximum compressive force
reached with the rate model in the heating step was higher than for the undrained case.
However, the slopes for the drained and rate models match closely towards the free end
of the pipe.
Friction profile
Figure 7-17 compares friction profiles (strictly ratios of force gradient to pipe weight) for
the drained (UB), undrained (LB) and rate models with a seabed slope. The friction values
for drained ( = 0.6) and undrained ( = 0.1) respectively are shown by the solid lines
and corresponding results for the rate model by dashed lines. Results of the rate model
match closely with the drained case over most of the pipe length, especially over the
283
Structural analysis of submarine pipelines
central part, for the cool down step but more closely with the undrained case towards each
Walking rate
Similar to the previous case of SCR tension, the walking rate over the five cycles is
compared for drained (UB), undrained (LB) and rate models in Figure 7-18. For the first
cycle the walk rate was similar for all the friction models. However, the effects of friction
factors became more prominent during subsequent cycles due to stabilisation of the model
from the second step onwards. The walking rate with the undrained (LB) friction was
much higher than for the drained case while results for the rate model falls between the
A comparative study has been carried out in this section to investigate the effect of the
rate model on walking triggered by thermal transients. For this case the pipe is assumed
to have a constant initial temperature of 10°C initially, and at the end of each cool down.
assumed. The same dimensionless parameters of ψ = 0.0311, Thc = 0.01 and c = 0.00072
were adopted.
Walking rate
Figure 7-19 shows the comparison of walking rate for drained (UB), undrained (LB) and
rate models for the thermal transient case. The walking rates over five cycles are presented
for the cold end of the pipe (x = 2000 m). The dotted line represents the displacements
with the rate model, whereas the solid lines represent the drained and undrained friction
cases. The walking rate with undrained friction was higher than for drained friction and
Force profile
The force profiles with the rate model are shown in Figure 7-20 for the thermal transient
case with the rate model. Force envelopes from the second heat up and cool down cycles
till the fifth cycle are plotted, and show that the force profile stabilized after the second
cycle.
The distributed displacement at various locations along the pipe are presented here to
illustrate the total walking behaviour of the pipeline with the rate model, for the cases of
SCR tension
Figure 7-21 shows the axial displacement at various points along the length of the
pipeline, for the case of SCR tension of 500 kN applied at x = 0. Other results for this
case have been presented previously. Towards each end of the pipeline the displacements
cycle between outward during heating and inward during cooling. However, there is a net
cumulative (negative) displacement towards the applied SCR tension. This is most
evident at the centre of the pipeline where the displacements are negative during both
Seabed slope
Figure 7-22 shows the plot of axial movement along the length of the pipe with a seabed
slope, β = 3° as indicated in the plot. The pattern of displacements is very similar to that
for the SCR tension case, but with larger net downhill (negative) cumulative
displacements. During each heating and cool down step the displacements at the pipeline
285
Structural analysis of submarine pipelines
centre are negative and even at x = 1500 m the net heating displacements become slightly
Thermal transients
Figure 7-23 shows the cyclic displacements along the length of the pipe for a thermal
transient of 10 °C/km, as indicated in the plot. The cyclic displacements follow a similar
pattern to those with no thermal gradient, being outward from the pipeline centre during
heating, and inward during cool down. The main difference relative to the previous two
cases is there is much less net walking of the pipeline. At the centre, x = 1000 m, small
cyclic displacements occur due to the additional heating for x < 1000 m compared with
x > 1000 m, but with barely discernible accumulation over successive cycles.
FRICTION
A parametric study was carried out to investigate the effect of dimensionless velocity on
the walking rate of the pipeline, to establish an approach to define an equivalent single
value of friction when using the rate model. Input parameters are tabulated in Table 7-2.
The midpoint velocity of the friction rate model was varied from v50 of 0.0069 m/s to
69.31 m/s (see Figure 7-24). This range allowed simulations to cover both extremes of
the friction model. Non-dimensional parameters and c (which was kept constant) are
shown in
286
Structural analysis of submarine pipelines
Table 7-3. All three pipeline walking mechanisms were studied in the following section
to understand how the walking rate is affected by the velocity parameter, and how a
Figure 7-25 depicts the effect of SCR tension and dimensionless velocity on the walking
rate. The SCR tension was varied from 100 to 500 kN and the dimensionless velocity was
varied from 0.00312 to 31.2. The portions of the walking profile affected by drained and
undrained friction behaviour are indicated. The walking rate became essentially constant
below the velocity limit ( ~ 0.1) controlled by drained friction. The walking rate
increased with increasing normalized velocity and reached a plateau with constant
Figure 7-26 shows the effect of seabed slope and dimensionless velocity on the walking
rate. The slope β was varied from 1° to 3° and the normalized velocity was varied from
6.83E-05 to 0.683. Arrows show the portions of the waking profile dictated by drained and
undrained friction. Similar to the previous case the walking rate was constant below the
undrained velocity limit of ~ 0.1 and then increased to a value beyond which the
walking rate was controlled by drained friction and constant for > 5. Another set of
data points could have improved the understanding of constant walking rate below the
convergence issues this was not plotted here. Section 7.5.5 below will provide details of
287
Structural analysis of submarine pipelines
Figure 7-27 depicts the effect of thermal transients and dimensionless velocity on the
walking rate. The thermal transients were varied from 10 to 30 °C/km and the normalized
velocity was varied from 6.83E-05 to 0.683. The portions of the walking profile affected
by the drained and undrained frictions are indicated, with the walking rate constant below
the drained velocity limit of ~ 1 and > 5. As noted previously with reference to Figure
7-23, the walking rate due to thermal transients is significantly lower than for typical SRC
Sets of data have been generated in the previous section which shows the relation between
the velocities of pipeline to the walking rate. However, for the design engineers a simple
static friction representing the walk per cycle affected by the dimensionless velocity
would be beneficial. Therefore, back calculations were made based on the existing
analytical models (Carr et al. 2006) to derive the equivalent friction factor for each of the
walk rate. The representative equivalent friction factor was then plotted against the
dimensionless velocity and a fit proposed to estimate the equivalent friction factor.
The equivalent friction was fitted via an exponential rate model, matching the underlying
rate model for the friction applied at each element within the simulations. Figure 7-28
depicts the representative friction factor matched closely with the back calculated friction
factors from walking rate. For SCR tension case the expression of equivalent friction can
be given by:
288
Structural analysis of submarine pipelines
Similarly, the equivalent friction factor for seabed slope (in Figure 7-29) can be evaluated
as:
The proposed models are in good agreement with the numerical results. Therefore, the
equivalent friction factors can be used to simulate the walking behaviour of pipelines
velocity of the pipe (Equation 5). This parameter can be determined a priori, without
performing any analyses, since it is related only to the free pipeline expansion, the
duration of the change in temperature and the friction model parameter, v50.
Until now the dimensionless time – being the ratio between heat-up and cool-down
duration, Thc = tH/tC – was kept constant at 0.01. However, here the effect of Thc is
investigated, varying it in the range 1000 – 0.001, by varying the heating time, tH and
keeping the cooling time tC constant at 1000 seconds. For these analyses the seabed slope,
β was kept constant at 3°, and effective characteristic strain, εc at 0.00072 as before.
Similar to previous cases equivalent friction factor was back calculated from the walking
rate for each case and compared with the proposed equivalent friction model. Figure 7-30
depicts the fitting of the proposed model for different dimensionless times, Thc for a
constant εc and β. The proposed model shows good agreement with the results of
individual analyses. A direct comparison of the proposed data and analysis results across
the range of values is shown in Figure 7-31, again showing close agreement.
The force profiles for various cases (Thc = 1000 – 0.001) are shown in Figure 7-32. The
walking rate decreases with increasing Thc, as may be inferred from the gradually
increasing offset of the maximum axial force from the pipeline centre as Thc decreases.
289
Structural analysis of submarine pipelines
The characteristic strain, εc expresses the ratio of the total pipeline weight divided by the
axial stiffness of the pipe. Hence if the weight of the pipe changes proportionally with
axial stiffness (essentially the cross-sectional area) of the pipe, εc does not change and the
walking rate is not affected. However, for a particular size of the pipe, if the weight of the
pipe is increased the walking rate is reduced. Figure 7-33 shows the effect of the
characteristic strain on the walking rate of a pipe on sloping seabeds. Two different seabed
slopes of, β = 1° and 3° were considered for three different pipe sizes of D = 0.152 m,
0.304 m and 0.912 m. The pipe length, L was kept constant at 2000 m. The dimensionless
speed, ψ was kept constant at 0.0311 and the dimensionless time, Thc was kept constant
at 0.01.
It was confirmed that, for a given value of εc, but for pipes with different diameter and
thickness, the walking rate is the same for a particular seabed slope. However, when the
weight of the pipe was doubled or halved by keeping all other parameters constant and
for a particular diameter, the walking rate was reduced and increased respectively. This
7.6 CONCLUSIONS
High temperature and high pressure seabed pipelines are often subjected to cyclic
expansion during their operating cycles. These repeated thermal cycles lead to global
the axial friction between a walking pipeline and the seabed, a velocity-dependent friction
model was introduced here. Numerical analyses were carried out to investigate the effect
velocity dependent friction model features a transition between friction limits that
290
Structural analysis of submarine pipelines
represent very slow and very fast movements, and was first compared with existing
experimental data of clay-pipe friction. Then the model was introduced into ABAQUS
and was validated with the single element and then flat seabed cases. Thereafter, walking
due to SCR tension, seabed slope and thermal transients were analysed and compared
Parametric results revealed the relationship between the velocity of the pipeline and the
walking rate. If the pipeline movements are sufficiently fast or sufficiently slow, then the
walking response matches a velocity-independent friction model tied to the lower and
upper friction limits respectively. However, there is a transitional zone in which the
pipeline movements mobilise a friction that lies between the two limits. The mobilized
friction also varies along the length of the pipeline due to the varying expansion rate at
different positions – with the ends of the pipeline always expanding or contracting the
fastest. A dimensionless velocity ratio has been defined, linking the free pipeline
For design practice, it is desirable to retain the simplicity of a simple single friction value
that can be used in structural modelling or existing analytical solutions for walking rate.
Therefore, back calculations were made based on the existing analytical models to derive
the equivalent friction factor that matches the walking rate observed in all of the velocity-
dependent friction cases. A relationship was then derived to calculate this representative
equivalent friction factor based on the dimensionless velocity ratio. This model provides
This analysis therefore assists in accounting for velocity-dependent axial pipeline friction
friction response into a single equivalent friction factor, allowing for the expected
291
Structural analysis of submarine pipelines
expansion velocity of the pipeline. This equivalent friction can be determined a priori,
without performing any numerical analyses, since it is related only to the free pipeline
expansion, the duration of the change in temperature and the friction model parameter,
v50
292
Structural analysis of submarine pipelines
References
Bhushan. B. (1999). The principles and applications of tribology. John Wiley & Sons.
Carneiro, D., Rathbone, A., Siong, Soon, Kok., & Viecelli, G. (2017). Velocity-
10.1115/1.4034695.
Carr, M., Sinclair, F., & Bruton, D. (2006). Pipeline Walking – Understanding the filed
layout challenges, and analytical solutions developed for the SAFEBUCK JIP.
Hill, A.J., White, D. J., Bruton, D.A.S., Langford, T., Meyer, V., Jewell, R., & Ballard,
J-C. (2012). New datasets and improved practice for assessment of axial pipe-soil
Geotechnics.
Oxford.
Oden, J. T., & Martins. J. A. C. (1985). Models and Computational Methods for
Randolph, M. F., White, D. J., & Yan, Y. (2012). Modelling the axial soil resistance on
Steelfelt, J. S. (1993). Sliding resistance for foundations on clay till. Proc. Wroth
293
Structural analysis of submarine pipelines
White, D. J., Ganesan, S. A., Bolton, M. D., Bruton, D. A. S., Ballard, J-C., & Longford,
OTC 21249.
White, D. J., & Cathie, D. N. (2010). Geotechnics for Subsea pipelines. Proc. of 2nd Int.
White, D. J., Campbell M., Boylan. N., & Bransby, M. F. (2012). Theoretical analysis of
axial pipe-soil interaction. Proc. SUT Conf on Offshore Site Investigation and
Geotechnics, London.
294
Structural analysis of submarine pipelines
10 10
1 10
0. 1 10
0.01 10
0.001 10
0.0001 10
0.001 10
0.01 10
0.1 10
1 10
10 10
295
Structural analysis of submarine pipelines
Diameter, D 0.912 m
pipe material, α
Change in temperature, ∆θ 90 °C
model, v50
296
Structural analysis of submarine pipelines
Parameters Value
297
Structural analysis of submarine pipelines
Figures:
µs µ = µk + (µs - µk )e-av
Friction coefficient, µ
Increasing a
µk
Velovity, v (m/s)
0.7
Drained, µUB µ = µLB + (µUB - µLB )e-λ(v/v50)
0.6
Friction coefficient, µ
0.5
0.4
0.3
0.2
0.1
v50 Undrained, µLB
0
0.0001 0.001 0.01 0.1 1 10
Velovity, v (m/s)
298
Structural analysis of submarine pipelines
0.6
0.4
Figure 7-3: Published data showing effect of velocity on interface shear strength (White
et al. 2012)
1 Soil P (SILT)
Marine clay (Hill
0.8 et al. 2012)
0.6
0.4
299
Structural analysis of submarine pipelines
Figure 7-5: Benchmarking rate model with single element test on flat seabed
10
1
Velocity, v (m/s)
0.1
0.01
0.001
0.0001
0.00001
0 20 40 60 80 100 120
Time, t (s)
300
Structural analysis of submarine pipelines
0.7
Drained, µUB µ = µLB + (µUB - µLB )e-λ(v/v50)
0.6
Friction coefficient, µ
0.5
v50 = 0.069 m/s
v50 = 0.0069 m/s
0.4
0.3
FE Results
0.2
0.1
Undrained, µLB
0
0.00001 0.0001 0.001 0.01 0.1 1 10
Velovity, v (m/s)
Figure 7-7: Comparing analytical and numerical results for single element test
Cycle number
0 1 2 3 4 5 6
0
End displacement at x = 0 (m)
-0.4
-0.6
-0.8
µUB = 0.6
-1
-1.2 µrate
ψ = 0.0311, ε c = 0.00072, T hc = 0.01
-1.4
Figure 7-8: End displacement for drained, undrained and rate models for flat seabed
301
Structural analysis of submarine pipelines
0.4
µrate
0.2
-0.2
µLB = 0.1
-0.4
-0.8
0 1 2 3 4 5 6
Cycle number
Figure 7-9: Mobilised friction at x = 0 for drained, undrained and rate models for flat
seabed
5
Cool down ψ = 0.0311
4 εc = 0.00072
UB Thc = 0.01
3
Axial force, P (MN)
2
LB
1
0
-1
-2
-3
Rate (dashed)
-4
Heat up
-5
0 500 1000 1500 2000
Position, x (m)
Figure 7-10: Force profiles for drained (UB), undrained (LB) and rate models for flat
seabed
302
Structural analysis of submarine pipelines
0.8
0.6
ψ = 0.0311
-0.2
-0.6
-0.8
0 500 1000 1500 2000
Position, x (m)
Figure 7-11: Friction profiles for drained (UB), undrained (LB) and rate models for flat
seabed
0.8
ψ = 0.0311
0.6 εc = 0.00072
Thc = 0.01
Friction factor, (dP/dx) / W'
0.4
-0.2
Figure 7-12: Enlarged view of central portion of friction profiles for flat seabed
303
Structural analysis of submarine pipelines
5
Cool down
4
UB(solid line)
3
2 LB
TSCR = 500 kN
1
0
Axial force, P (MN)
-1
-2 Rate (dashed)
-3 TSCR = 500 kN
εc = 0.00072
-4 ψ = 0.0311
Heat up
Thc = 0.01
-5
0 200 400 600 800 1000 1200 1400 1600 1800 2000
Position, x (m)
Figure 7-13: Force profiles for drained (UB), undrained (LB) and rate models with TSCR
= 500 kN
0.8
UB(solid line)
0.6
0.4
Friction factor, (dP/dx) / W'
LB
0.2
-0.2
TSCR = 500 kN Rate (dashed)
-0.4 εc = 0.00072
ψ = 0.0311
Thc = 0.01
-0.6
-0.8
0 500 1000 1500 2000
Position, x (m)
Figure 7-14: Friction profiles for drained (UB), undrained (LB) and rate models with
TSCR = 500 kN
304
Structural analysis of submarine pipelines
Cycle number
0 1 2 3 4 5 6
0
-2
-3 µrate
-4 TSCR = 500 kN
εc = 0.00072
ψ = 0.0311 µLB = 0.1
Thc = 0.01
-5
Figure 7-15: Walking responses for drained (UB), undrained (LB) and rate models with
TSCR = 500 kN
5
Cool down β=3
4 ψ = 0.0311
εc = 0.00072
3 UB (solid Thc = 0.01
line)
Axial force, P (MN)
2 β
LB
1
0
-1
-2
Rate (dashed)
-3
-4
Heat up
-5
0 500 1000 1500 2000
Position, x (m)
Figure 7-16: Force profiles for drained (UB), undrained (LB) and rate models for
305
Structural analysis of submarine pipelines
UB (solid line)
0.8
0.6
β=3
ψ = 0.0311
Friction factor, (dP/dx) / W'
0.4 εc = 0.00072
LB Rate (dashed) Thc = 0.01
0.2
-0.2
β
-0.4
-0.6
-0.8
0 500 1000 1500 2000
Position, x (m)
Figure 7-17: Friction profiles drained (UB), undrained (LB) and rate models for sloping
Cycle number
0 1 2 3 4 5 6
0
Down-hill end displacement at x = 0 (m)
-1
-2 µUB = 0.6
-3
-4 µrate
-5
β=3
β
-6 ψ = 0.0311
εc = 0.00072
µLB = 0.1
Thc = 0.01
-7
Figure 7-18: Walking response for drained (UB), undrained (LB) and rate models for
306
Structural analysis of submarine pipelines
1.2
µdecayμ(dotted line)
rate (dashed) µLB = 0.1
at cold end, x = 2000 (m)
Expansion at cold end, u (m)
1
0.8
0.6
µUB = 0.6
0.4
Displacement
0.2
εc = 0.00072
q = 10 °C/km
0
0 1 2 3 4 5 6
Cycles
Figure 7-19: Walking responses for drained (UB), undrained (LB) and rate models
for q = 10 °C/km
5
Cool down
4 Cycle 2 - 5
3
2
Axial load, P (MN)
Cycle 1
1
Rate model
0
0 500 1000 1500 2000
-1
-2
q = 10 C/km
-3
εC =ε0.00072
c = 0.00072
-4 q = 10 =
ψ 0.0311
°C/km
Heat up Cycle 2 - 5 ψ = 0.311
Thc = 0.01
-5
Along the length, x (m)
Position, x (m)
307
Structural analysis of submarine pipelines
2
TSCR = 500 kN P2000
1
Axial displacement, u (m)
-1
-2
TSCR = 500 kN
-3 εc = 0.00072
ψ = 0.0311
Thc = 0.01
-4
-500 0 500 1000 1500 2000 2500
Position, x (m)
Figure 7-21: Distributed displacements with rate model for SCR tension, TSCR = 500 kN
1
P0 P500 P1000 P1500
P2000
0
Axial displacement, u (m)
-1
-2
-3 β=3
β ψ = 0.0311
εc = 0.00072
Thc = 0.01
-4
-500 0 500 1000 1500 2000 2500
Position, x (m)
Figure 7-22: Distributed displacements with rate model for seabed slope, β = 3 °
308
Structural analysis of submarine pipelines
1.5 120
Gradient = 10°C/km
100 P2000
40 P1500
20
0.5 0
0 500 1000
Along the length, x (m)
1500 2000 P1000
-0.5
P500 10 °C/km
q == 10 C/km
-1 εCc = 0.00072
P0 ψ == 0.311
0.0311
T hc= =
0.01
0.01
-1.5
-500 0 500 1000 1500 2000 2500
Along the length, x (m)
Position, x (m)
Figure 7-23: Distributed displacements with rate model for thermal transient, q = 10
°C/km
0.7
Drained, µUB µ = µLB + (µUB - µLB )e-λ(v/v50)
0.6
Friction coefficient, µ
0.1
Undrained, µLB
0
0.0001 0.001 0.01 0.1 1 10 100 1000 10000
Velocity, v (m/s)
Figure 7-24: Range of friction models used for numerical parametric study
309
Structural analysis of submarine pipelines
10
TSCR = 100 kN Controlled by Controlled by undrained friction
Walk per cycle due to SCR tension, ∆SCR (m/cycle)
drained friction
TSCR = 200 kN
1 TSCR = 300 kN
TSCR = 400 kN
TSCR = 500 kN
0.1
0.01
ε c = 0.00072
Thc = 0.01
0.001
0.0001 0.001 0.01 0.1 1 10
Dimensionless speed, ψ
Figure 7-25: Effect of dimensionless velocity on walking rate for SCR tension
1
Controlled by
Walk per cycle due to thermal transients, ∆T
q = 10°C/km undrained
Controlled by friction
q = 20°C/km drained friction
0.1 q = 30°C/km
(m/cycle)
0.01
εc = 0.00072
T = 0.01
0.001
0.0001 0.001 0.01 0.1 1
Dimensionless speed, ψ
Drained Undrained
Figure 7-26: Effect of dimensionless velocity on walking rate for seabed slope
310
Structural analysis of submarine pipelines
10
0.01
0.001
ε c = 0.00072
Thc = 0.01
0.0001
0.0001 0.001 0.01 0.1 1 10
Dimensionless speed, ψ
Figure 7-27: Effect of dimensionless velocity on walking rate for thermal transients
Controlled by
0.2 drained friction
0.3
0.4
0.5
0.6
ε c = 0.00072
Thc = 0.01
0.7
0.0001 0.001 0.01 0.1 1 10 100
Dimensionless speed, ψ
311
Structural analysis of submarine pipelines
0.2
0.3 Controlled by
drained friction
0.4
0.5
0.6 ε c = 0.00072
Thc = 0.01
0.7
0.0001 0.001 0.01 0.1 1 10 100
Dimensionless speed, ψ
0.7
T_hc= 0.001
0.5
T_hc = 0.01
0.1 εc = 0.00072
β = 3°
tc = 1000 s
0
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
Proposed equivalent friction factor, μ_proposed
Figure 7-30: Comparison between proposed and calculated friction factors for varying Thc
312
Structural analysis of submarine pipelines
0
Controlled by undrained
friction
0.1 µeq = µLB + (µUB - µLB )e-0.95ψ
0.6 εc = 0.00072
β = 3°
tc = 1000 s
0.7
0.00001 0.001 0.1 10 1000
Dimensionless speed, ψ
Figure 7-31: Equivalent friction model for cases with varying ratios of heating and
5
Cool down T = 1000,100,10,1,0.1,0.01
4
Effective Axial Force, P (MN)
3
T=0.001
2
1
0
0 500 1000 1500 2000
-1
T=0.1
-2 T= 0.01
T=1 T= 0.001
-3 ε c = 0.00072
-4 β = 3°
Heat up T =1000,100,10
t c = 1000 s
-5
Position, x (m)
Figure 7-32: Comparison of force profiles for various ratios of heating and cooling
times, Thc
313
Structural analysis of submarine pipelines
0.3
ψ = 0.0311
β=3
(m/cycle) 0.2
0.1
β=1
0
0 0.0004 0.0008 0.0012 0.0016
Characteristic strain, εc
Figure 7-33: Effect of characteristic strain on walking rate of pipe for various weights
314
Structural analysis of submarine pipelines
315
Structural analysis of submarine pipelines
316
Structural analysis of submarine pipelines
OVERVIEW
The proposed research was concerned with the structural behaviour of submarine
pipeline subjected by submarine slides, and thermal loading conditions. The research aim
was to support the transition of oil and gas developments into deeper water and more
remote conditions. This was achieved by improving the techniques for assessing the axial
pipe-soil interaction forces resulting from relative pipe-soil movement, including the
CONCLUDING REMARKS
In this chapter the main findings, concluding remarks along with possible directions
for future work are presented. Throughout the thesis analytical models were developed to
tackle the pipe-soil interaction issues and structural problems and numerical solutions
were provided. Where appropriate, the theoretical techniques used for pile design were
Limited information exists in the literature on the axial elastic pipe-soil response of
analytical approach was developed for axial stiffness, focusing on the case of shear
modulus proportional to depth. The solution was then verified through numerical analysis.
Further numerical analysis was carried out to obtain relationships for horizontal and
317
Structural analysis of submarine pipelines
recommendation for the selection of proper elastic stiffnesses in all three directions of
motion. These recommendations allow consistent and rigorous modelling of elastic pipe–
seabed interactions with application to the analysis of pipeline laying, buckling, walking,
and on-bottom stability. They have found adoption in practice and have recently been
referenced in the new DNV code for pipe-soil interaction (DNVGL-RP-F114, 2019).
The contributions in chapter 3 lay the basic foundation for the following chapters on axial
slide-pipeline interaction and then buckling analyses due to axial loading. Once the axial
elastic stiffness is known, then the ultimate axial resistance was evaluated and axial
springs were formulated and implemented into the FE code to carry out the structural
analyses of pipeline.
In Chapter 4 an analytical solution was developed for axial submarine slide loading
of a straight on-bottom pipeline. It was shown that the non-dimensional axial loads and
axial displacements depend on three non-dimensional input parameters, i.e. the driving
force in the slide zone, seabed resisting force in passive zone, and pipe-soil stiffness. Non-
dimensional design charts were presented to show the effect of individual input
parameters on axial loads and axial displacements. The maximum axial load in the pipe
is directly proportional to the slide force, while the load at the transition from elastic to
plastic soil resistance is initially proportional to the slide force but then becomes limited.
The limit is reached for most relevant values of the slide force. Numerical FE analysis
was also carried out to verify the analytical model. Analytical and numerical solutions
were shown to match with close agreement. The buckling was ignored in the analytical
model. However, the existing classical theory of buckling was linked to the output of the
318
Structural analysis of submarine pipelines
analytical model to show the vulnerability of the pipelines towards buckling in case of
various slide loading conditions. On bottom submarine pipelines are more susceptible to
in the offshore industry due to the complexities involved in the analysis, and the potential
researchers and recent joint industry projects provide a basis for estimation of the critical
buckling load of a straight, or in some cases imperfect, pipe on either a rigid or elastic
seabed. However, systematic solutions for the combined effects of nonlinear soil
properties and the as-laid geometry – specifically the out-of-straightness – on the buckle
straight) on-bottom pipeline subjected to axial compressive loading was carried out. The
seabed was modelled with lateral and axial elastic, perfectly plastic, springs to idealise
the load-displacement behaviour of the soil and the pipe was modelled with pipe elements.
Buckling was performed by a displacement controlled finite element method with the
modified RIKS algorithm that is available in the commercial software ABAQUS. This
numerical tool was used to develop a parametric solution for the present problem in terms
of the various pipe material and geometry parameters and the lateral and axial pipe-soil
interaction parameters. In particular, the influence of the magnitude and stiffness of the
lateral pipe-soil response was investigated, highlighting the sensitivity of the pipeline
response to the geotechnical inputs. The results have been synthesised in a generic non-
dimensionalised design chart to estimate the buckling load, valid for the range of inputs
The chapter presented a new solution for the buckling force that includes both
OOS and initial stress. It was found that for the range of parameters investigated the effect
of the axial stress can be introduced simply as a new dimensionless group, without
changing the influence of the other groups in the correlation. This contribution provides
a more efficient basis for predicting buckling behaviour, and highlights in a simple way
plastic, expressed as ultimate resistance per unit length. Often this term is expressed non-
pipe weight, W. However, the elastic-plastic (i.e. bi-linear) behaviour of the soil and the
effect of this bi-linear response of the soil on the walking behaviour were poorly
is defined as the amount of axial displacement that occurs before the ultimate friction is
generated, and the resistance rises linearly with displacement up to this value. The
results were observed by previous researchers. However, the analytical solution for the
reduction of walking per cycle due to increase in mobilisation displacement has not been
insight into the walking behaviour due to seabed slope, SCR tension and thermal
transients. The existing analytical solutions were extended to incorporate the elastic-
plastic response of the soil into the expression of pipeline walking and a new derivation
was proposed. Numerical verification with ABAQUS was also done for the proposed
320
Structural analysis of submarine pipelines
expression. The resulting expression and calculation method extends existing analytical
loading and expansion but it has been recognised that the sliding resistance between a
pipe and the seabed varies with velocity, partly due to drainage effects. In Chapter 7 the
numerical model that was validated in the previous chapter is extended to explore the
implemented in commercial software ABAQUS and validated via single element and
simple (flat seabed) pipeline cases. This model features upper and lower friction limits,
was performed using differing rates of heating and cool-down in walking situations driven
by seabed slope, SCR end tension and the difference between heat up and cool down rates.
The walking behaviour is compared to cases with constant friction and solutions are
friction. These equivalent friction values can then be applied in existing simple solutions
or more complex numerical analyses, as a short cut method to account for velocity-
dependent friction.
FUTURE RESEARCH
The new developments of oil and gas infrastructure are moving into deeper water.
This requires design and construction of long high temperature and high pressure
pipelines from deep sea to shore. These pipelines are subjected to cyclic expansion during
operating cycles and geohazards such as submarine slides. Concluding remarks out of this
research have been detailed in the above section and some of the future research scopes
321
Structural analysis of submarine pipelines
1. The research in this thesis has been exclusively theoretical, building on existing
substantiate the proposed models of submarine pipeline behaviour under the slide and
thermal loading. The analytical models were limited to elastic limits of the pipe
material and the finite element analyses were confined and limited due to the SSFE
methods for elastic stiffness prediction. These solutions can be extended to multiple
directions of loading, in order to identify the coupling between the elastic and plastic
where independent springs are ‘attached’ to the pipe for each direction of motion.
This coupling is already recognised and quantified for vertical-lateral (V-H) loading
(e.g. Randolph & White 2008). To illustrate this potential future work, the FE model
developed in chapter 3 has been used to generate a V-H yield envelope for typical soil
parameters. A fully bonded pipe-soil interface was considered and V-H swipes were
soil with uniform shear strength of (su = 5 kPa) was considered. The resulting V-H
envelopes matched the existing literature for all the embedment cases. The same
methodology was used to provide the V-F and H-F yield envelopes for exactly similar
Figure 8-1and Figure 8-2, there is coupling, meaning that when a significant fraction
of either the vertical or lateral capacity is mobilised, the available axial resistance is
reduced. This may influence the response of pipeline in buckles and submarine slides,
where both the axial and lateral resistance is concurrently mobilised. For future work,
these SSFE example could be extended to large deformation FE (LDFE) analysis with
322
Structural analysis of submarine pipelines
2. In the analytical axial slide pipeline interaction problem - the slide loading was
perfectly straight and elastic in nature. However, the complexities associated with the
as-laid geometry of the pipe along with its elastic-plastic material behaviour need to
be studied in the future research. At the same time studies need to be extended for
dynamic loading conditions where the shear strength of the slide material changes
along its travel path to check the structural integrity of the pipelines along its length
3. For the lateral buckling on submarine pipeline the numerical model developed in this
thesis needs to be extended with dynamic simulation and with nonlinear user defined
solution. Furthermore the methodology can be used to provide solution in the vertical
direction, and where there is out-of-straightness present from the lay process prior to
4. All the structural analyses in this thesis involve a rigid or elastic seabed. During lateral
buckling, the pipe may undergo large lateral displacement and bending strains,
structural and geotechnical responses using a full three-dimensional (3D) model may
both structural and soil domains is a significant challenge from a project schedule and
numerical (FE) perspective. With the advent of supercomputers, and also new solver
5. For pipeline walking – the analytical solution with bi-linear soil-pipe interaction, and
the numerical solution with velocity dependent friction model need to be verified with
323
Structural analysis of submarine pipelines
the field data from various joint industry projects and also with the experimental data
through full scale laboratory test to predict walking behaviour more accurately for
carried out to provide more generic design charts for the pipeline engineers as a guide.
Regardless of limitations of the current work and scope for future research, it is
considered that the outcomes of this research present several useful contributions to the
current body of knowledge of pipe-soil interaction and structural analysis. They also
analysis of submarine pipeline and pipe-soil interaction, and have produced new analysis
324
Structural analysis of submarine pipelines
References:
Randolph M.F. & White D.J. 2008. Upper bound yield envelopes for pipelines at shallow
325
Structural analysis of submarine pipelines
Figures
1.8
w/D = 0.1 to 0.5
1.6
Normalised axial resistance, F/s uD
1.4
1.2
1
0.8
0.6
0.4
0.2
0
0 1 2 3 4 5 6 7
Normalised vertical resistance, V/suD
1.8
w/D = 0.1 to 0.5
1.6
Normalised axial resistance, F/suD
1.4
1.2
0.8
0.6
0.4
0.2
0
0 0.5 1 1.5 2 2.5 3 3.5 4
Normalised lateral resistance, H/suD
326
Structural analysis of submarine pipelines
327
Structural analysis of submarine pipelines
328