Omae2017 61282
Omae2017 61282
net/publication/317616074
Advanced 3-D FEA Modelling for a Modern and Multidisciplinary Pipeline Design
Approach
CITATIONS READS
2 352
5 authors, including:
Some of the authors of this publication are also working on these related projects:
All content following this page was uploaded by Lorenzo Marchionni on 19 October 2018.
OMAE2017-61282
ADVANCED 3-D FEA MODELLING FOR A MODERN AND
MULTIDISCIPLINARY PIPELINE DESIGN APPROACH
Gravel
Berms
3.2 Pre-Processing
Using the pre-processing module into the GUI interface (red
section in Figure 1) it is possible to select all the relevant
Figure 4 – Sleepers Used as Triggers for Lateral Buckling. parameters needed for the ISB, OBR or PW analysis. In
particular, the following data are considered:
Mitigation Measures Against Pipeline Walking Pipeline data (size, material and mesh characteristics);
If the pipeline is susceptible to pipeline walking, artificial Residual Lay-Pull as a function of the KP;
anchoring points shall be placed along the route to limit Design Conditions in terms of pressure and differential
pipeline ends displacements within acceptance criteria for temperature (steady or transient conditions);
spools or termination structures (PLET/FLET). Non-linear pipe-soil relationships considering the pipeline
Non-linear springs and/or local masses can be used to in empty, flooded or in operation condition;
simulate the effects of ancillary devices such as suction piles Material physical and mechanical properties (even
and chains, mattresses or in-line structures. temperature dependent);
Laying direction;
Lateral Counteracts Artificial lateral OOS along the pipeline route;
Lateral counteracts (LaCs) are modelled as analytical rigid Mitigation measures type and location (berms, mattresses,
surfaces. Any specific force-displacement pattern simulating sleepers, counteracts, suction piles or in-line structures).
the interaction between the pipeline and LaCs can be
managed in the model by use of contact pair (Figure 5). 3.3 Post-Processing
The post-processing module (blue section in Figure 1) permits
to plot all the relevant parameters along the pipeline route at
one or more time increments/steps representing different
operating scenarios (typically as-laid, hydrotest and operating
conditions, but also intermediate instants or predefined cycles
for PW analysis). In particular, the following outcomes are
available:
- Horizontal and vertical pipeline configuration;
- Free-span height/length/clearance;
- Steel and effective axial force;
- Axial (feed-in), lateral and vertical pipeline movements
with respect to the as-laid configuration;
- Longitudinal, hoop and Von Mises stresses and strains;
Figure 5 – 3-D Seabed Model With ABAQUS – Lateral - Bending and torsional moments;
Counteracts Along the Pipeline Route. - Local Buckling Unity Check as per DNV-OS-F101 [14] or
pipeline integrity check according to SAFEBUCK [10]
Loading Sequence and Solution Algorithm design criteria (modified DCC check).
The FE model is subject to the following loading steps: - Vertical and lateral reactions of berms (used for ISB
1. Step 1: Lay the interested pipeline stretch on the sea analyses) or forces applied to chains (suction piles) or
bottom considering: in-line structures (used for PW analyses).
o Application of the relevant residual lay tension to the
straight pipeline section under investigation; The use of such tools significantly reduces the time to prepare
and process a set of pipeline global analyses for the expert
engineer, thus providing the basis for decision making and
Vertical
Residual Axial Friction Factor 0.5 -- -0.5
0.732 – 0.924
KP (km)
P/L Submerged Flooded Weight t/m 1.4
1.2
Unity Check
Table 1 – Case 1: Main Input Parameters. 0.8
0.6
0.4
During flooding and hydrotest conditions, the pipeline sags into 0.2
0
5.6 5.65 5.7 5.75 5.8 5.85 5.9 5.95 6
the Deep Depression (Figure 6) due to the increased weight and KP (km)
Flooded KP5752
Flooded KP5732
Flooded KP5700
-515
Flooded KP5650
Vertical Seabed Profile
particular, during hydrotest, the pipeline flooding operation is
-520
controlled by a pig from shallow to deep. The pipeline
-525
configuration along the free-span generated through the Deep
-530
-535
Depression is strongly affected by sequential increase of weight
-540
(Figure 9). It is fundamental in these circumstances to simulate
5.6 5.65 5.7 5.75 5.8 5.85 5.9 5.95 6
KP (km)
the actual load sequence since an incremental uniform
Figure 9 - Sequential Water Flooding along the Deep application of loads along the overall pipeline route leads to
Depression. unrealistic results.
In order to prevent pipeline overstress along the Deep
This scenario was expected to not overstress the pipeline since Depression (Figure 10b), alternative solutions have been
the rock berms were accurately positioned to control the investigated considering:
sagbend. However, the vertical displacement is not fully Additional gravel supports to limit the pipe sagging after
recovered once the line is depressurized and dewatered due to flooding (Figure 11),
the axial resistance provided by the soil at the shoulders The use of batches of thicker pipe to improve the pipeline
(Figure 10a). Therefore, during the following load step when capacity,
the internal pressure is set to the incidental one and the design A combination of the above solutions.
temperature is applied (named mechanical design step), the
pipeline continues sagging, leading to an excessive bending The selected solution was to include additional gravel supports
concentration at the bottom of the Deep Depression to limit the pipe sagging. In particular, the position and height
(Figure 10b). of two additional post-lay gravel berms has been optimized.
-490
Figure 11a shows the vertical seabed and pipeline profiles
-495
resulting from the analysis run including the additional
-500 intervention works, while Figure 11b shows the applied vertical
-505 bending moments and local buckling unity checks (now within
-510 acceptance criteria).
Z-coordinate (m)
-515 -490
-520 -495
-525 -500
As Laid -505
-530
Flooded
Hydrotest
Z-coordinate (m)
Depressurization -510
-535
Mechanical Design
Vertical Seabed Profile
-515
-540
5.6 5.65 5.7 5.75 5.8 5.85 5.9 5.95 6
KP (km)
-520
0.5
0
distribution and hence the pipe-LaCs interaction. The most
Vertical
0.6
0.4
0.2
0
5.6 5.65 5.7 5.75 5.8 5.85 5.9 5.95 6
KP (km)
Lateral
stabilization pursued with the installation of Lateral -10000
-5000
Pre-IW. 0.8
Unity Check
0.6
areas where the pipeline was partially free spanning and the 0.2
0
2 2.2 2.4 2.6 2.8 3 3.2 3.4 3.6 3.8 4 4.2 4.4 4.6 4.8 5 5.2 5.4 5.6 5.8
contact with the soil occurred only in correspondence of KP (km)
transferred to the counteracts and to the GB. 2 2.2 2.4 2.6 2.8 3 3.2 3.4 3.6 3.8 4 4.2 4.4 4.6 4.8 5 5.2 5.4 5.6 5.8
Lateral Displacements (m)
-2
The position of lateral counteracts has been defined based on a
-4
simplified analytical approach (2-D model). This model is 2 2.2 2.4 2.6 2.8 3 3.2 3.4 3.6 3.8 4 4.2 4.4 4.6 4.8 5 5.2 5.4 5.6 5.8
Vertical Displacements (m)
0
based on the equilibrium between the horizontal forces along -2
the pipe, which tend to straight the curve, and the lateral -4
-6
resistance provided by lateral counteracts. The horizontal forces 2 2.2 2.4 2.6 2.8 3 3.2 3.4 3.6 3.8 4
KP (km)
4.2 4.4 4.6 4.8 5 5.2 5.4 5.6 5.8
have been evaluated during the installation phase, neglecting Figure 12c - Pipeline Axial, Lateral and Vertical
the beneficial contribution of soil friction, and the flooded Displacements.
condition. 300
As Laid
Maximum Counteract Reaction (kN)
perform the on-bottom roughness analysis has been re-used to: 200 Mechanical Design
100
0.2
uneven seabed morphology and the local forces from pipe- As Laid
Flooded
Hydrotest
0.15
0.05
The analysis has been carried out for all the design conditions
from empty to normal operating / mechanical design. Attention 0
3285 3508 3635 3689 4400 4460 4530 4593 4976
KP Counteracts (m)
5135 5190 5454 5490 5560 5603
has been paid to the temporary conditions and, in particular, the Figure 13 - Maximum LaCs Displacements and Reactions.
water filling process and the depressurization after the
1
In general the P/T (pressure and temperature) transient profiles
0
can greatly vary depending on the pipeline service (e.g. oil, gas, First End
Second End
water, etc…) and the types of cycle the flowline is subjected to -1
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
Walking Rate (m/cycle)
0.1
during field life. In this case, two types of start-up sequences 0.08
0.06
0.04
were provided (i.e. cold restart, warm restart preceded by hot 0.02
0
-0.02
oil circulation in both directions). As an example, Figure 14 -0.04
-0.06
First End - Walking Rate
Second End - Walking Rate
shows the cold restart temperature profiles along the flowline at -0.08
-0.1
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24
different time instants, which shall be applied cyclically in the Cycle Number (#)
FEM simulation along with shut-down profiles. The in-house Figure 15 - Pipe Ends Axial Displacements and Walking Rate
Matlab tool allows the engineer to produce the input file for Without Mitigation Measures.
ABAQUS which combines any start-up or shut-down sequence As-Laid
8
in a very flexible manner, drastically reducing the time spent Hydrotest
Pre-Commissioning
for input file generation. 7 Half Cycle - Cycle #1
End Cycle - Cycle #1
Lateral Displacements (m)
10
3.75
4.00
Figure 16 - Lateral Displacement Profiles Without Mitigation End Cycle - Cycle #12
Half Cycle - Cycle #13
0 4.25
4.50
Measures. End Cycle - Cycle #13
Half Cycle - Cycle #14
End Cycle - Cycle #14
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 4.75 Half Cycle - Cycle #15
5.00 End Cycle - Cycle #15
KP [km] 10.00 To minimize the load required to stop PW care must be taken Half Cycle - Cycle #16
50.00 End Cycle - Cycle #16
of the effective axial force and of the potential for lateral Half Cycle - Cycle #17
End Cycle - Cycle #17
Figure 14 - Cold Restart Temperature Profiles at Different buckling occurrence along the pipeline. If no lateral buckle Half Cycle - Cycle #18
End Cycle - Cycle #18
Time Instants (in hours) along the Flowline. occurs, the best theoretical position for a line which does not Half Cycle - Cycle #19
End Cycle - Cycle #19
Half Cycle - Cycle #20
reach fully restrained conditions is mid-line. If lateral buckles End Cycle - Cycle #20
Half Cycle - Cycle #21
End Cycle - Cycle #21
Half Cycle - Cycle #22
End Cycle - Cycle #22
Half Cycle - Cycle #23
End Cycle - Cycle #23
9 Copyright © 2017 by ASME
Half Cycle - Cycle #24
End Cycle - Cycle #24
occur, the best theoretical position is at virtual anchor point Due to the evidenced complexity of the assessment, fully
(uphill with respect to cumulated displacement direction). non-linear FE analyses are required to verify all the possible
However, it is worth noting that: interactions between PW and lateral buckling.
For a mitigation measure which reacts in both directions In order to optimize the mitigation measures, most promising
(e.g. In-Line Structures), the best theoretical position is strategy were investigated including mattresses, suction piles,
between the virtual anchor point and the lateral buckle dead man anchors and post-installed In-Line Structures (ILS).
apex; however since the mitigation measure itself might Mattresses aim at increasing the vertical force exerted by the
change the effective axial force profile and since the pipe on soil and therefore increasing the pipe axial reaction.
transient profiles change the load on the mitigation Dead man anchors and suction piles are linked to the pipeline
measures as the P/T cycles develop, the real optimal via forged piece/collar and chains and represent direct
position cannot be determined in advance. constraints for pipeline displacements. ILS is installed along the
Robustness checks must be carried out to verify that the pipeline and assembled with an in-line anchor (ILA)
mitigation measure strategy for PW does not induce lateral incorporated in a quad-joint with the functionality to transfer
buckle activation between pipeline end and the mitigation the axial load from the pipe to the mitigation device. The ILS
device itself, in which case an uncontrolled feed-in might shall sustain the axial reaction assessed by FE walking global
still trigger PW. analyses.
When PW and lateral buckling occur at the same time, a During the optimization activities, mattresses and dead man
number of related phenomena can take place, all of which have anchors were discarded since they were shown to be not really
been observed during the optimization analyses, e.g.: effective on soft soils and not suitable for the high required
While fixing first end of the pipeline as shown below, loads to suppress walking propensity, while suction piles were
second end might migrate feeding-in into the lateral buckle discarded due to tight installation schedule, layout constraints
in the opposite direction; and presence of hard ground soils in depth. Post-installed
Walking at second end might still occur in case a route in-line structures have been proven to be the most flexible and
curve positioned e.g. close to the second virtual anchor reliable solution for quite to very severe applications. In
point shown in Figure 17 slightly moves laterally; addition, this kind of mitigation allows to take the advantages
Fixing walking at first end might “transfer” the propensity of the “wait and see” approach. Actually, they could be
to accumulate axial displacements to the second end. installed at a later time in case walking develops during field
life with a minimum preliminary set-up, (ILA installation at
predetermined KP along the line) and they could be adapted to
provide a range of axial reactions as needed.
Effective Axial Force
Slope/
T-transients effects and reacting in the longitudinal direction, being the other end of
Lateral buckle the spring fixed.
A number of FE analyses have been carried out, including
Mitigation measure sensitivities on pipe-soil lateral friction factors, P/T transient
(uphill) profiles, ILS position along the route, ILS gap
Reaction
Virtual anchor points (100 mm, -100 mm, +/-50 mm), with the aim to minimize the
KP number and dimension of the post-installed ILS, optimize their
positions to limit PW and verify the robustness of the chosen
mitigation measures.
Effective Axial Force
Weak slope/
T-transients effects life. Thanks to this ILS, the walking rate at first end was
negligible at the end of the simulation. The maximum
Lateral buckle
registered reaction was lower than 200 kN (see effective axial
Mitigation measure force discontinuities at KP 0.714 in Figure 19), well below the
(uphill) structure capacity (~600 kN).
Reaction
Virtual anchor points
KP
Figure 17 - Effective Axial Force With Lateral Buckle and
Seabed Profile (Schematic).
Displacement (m)
End Cycle - Cycle #1
Half Cycle - Cycle #2 0.5
600 End Cycle - Cycle #2
Pipe Axial
Half Cycle - Cycle #3
End Cycle - Cycle #3
Half Cycle - Cycle #4 0
400 End Cycle - Cycle #4
Half Cycle - Cycle #5
End Cycle - Cycle #5
Half Cycle - Cycle #6
End Cycle - Cycle #6
-0.5
200 Half Cycle - Cycle #7
End Cycle - Cycle #7
Half Cycle - Cycle #8
End Cycle - Cycle #8 -1
0 Half Cycle - Cycle #9 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
End Cycle - Cycle #9 0.02
Half Cycle - Cycle #10 First End - Walking Rate
Figure 18 - Effective Axial Force Profiles With One ILS Half Cycle - Cycle #20
End Cycle - Cycle #20 Figure 20 - Pipe Ends Axial Displacements and Walking Rate
(KP 0.714). With One ILS (KP 0.714).
As-Laid
6 Hydrotest
Pre-Commissioning
The position of the ILS has been optimized taking into account Half Cycle - Cycle #1
End Cycle - Cycle #1
Half Cycle - Cycle #2
5 End Cycle - Cycle #2
all the constraints given by installation and geotechnical Half Cycle - Cycle #3
End Cycle - Cycle #3
Half Cycle - Cycle #4
disciplines. Figure 19 shows the effective axial force profile
KP 0.714). Figure 20 and Figure 21 show the pipe ends axial End Cycle - Cycle #10
Half Cycle - Cycle #11
End Cycle - Cycle #11
displacements, walking rate and pipeline lateral displacements 1 Half Cycle - Cycle #12
End Cycle - Cycle #12
Half Cycle - Cycle #13
along the route, respectively. At the second end a small walking End Cycle - Cycle #13
Half Cycle - Cycle #14
End Cycle - Cycle #14
0
rate (< 5mm/cycle) was registered towards mid-line see 0 0.2 0.4 0.6 0.8 1
KP (km)
1.2 1.4 1.6 1.8 - Cycle #152
Half Cycle
End Cycle - Cycle #15
Half Cycle - Cycle #16
Figure 20. Actually from Figure 21 it is evident that the pipe is Figure 21 - Lateral Displacement Profiles With One ILS End Cycle - Cycle #16
Half Cycle - Cycle #17
End Cycle - Cycle #17
moving into the buckle amplifying the lateral displacement as (KP 0.714). Half Cycle - Cycle #18
End Cycle - Cycle #18
Half Cycle - Cycle #19
600
Half Cycle - Cycle #1
End Cycle - Cycle #1
Half Cycle - Cycle #2
only in case during field life an excessive cumulated
End Cycle - Cycle #2
Half Cycle - Cycle #3
displacement at second end occurs. It has been recommended to
End Cycle - Cycle #3
400
Half Cycle - Cycle #4
End Cycle - Cycle #4
the Company to check this condition with regular surveys.
200
Half Cycle - Cycle #5
End Cycle - Cycle #5
Half Cycle - Cycle #6
With two ILS, it has been proven that PW susceptibility is
0
End Cycle - Cycle #6
Half Cycle - Cycle #7 successfully and completely suppressed as shown in Figure 22,
End Cycle - Cycle #7
Half Cycle - Cycle #8
End Cycle - Cycle #8
with walking rate at the end of the simulation almost negligible
Half Cycle - Cycle #9
-200
End Cycle - Cycle #9
Half Cycle - Cycle #10
and while maintaining a reaction on the structure lower than the
-400
End Cycle - Cycle #10
Half Cycle - Cycle #11 allowable value. Including the sensitivities, the overall
End Cycle - Cycle #11
Half Cycle - Cycle #12
End Cycle - Cycle #12
maximum ILS reaction was lower than 350 kN.
-600 Half Cycle - Cycle #13
End Cycle - Cycle #13
Half Cycle - Cycle #14
-800 End Cycle - Cycle #14 0.5
Pipe Axial Displacement (m)
0 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 Half2Cycle - Cycle #15
KP (km) End Cycle - Cycle #15
Half Cycle - Cycle #16
Figure 19 - Effective Axial Force Profiles With One ILS End Cycle - Cycle #16
Half Cycle - Cycle #17 0
End Cycle - Cycle #17
(KP 0.226, Non-Optimized Case). Half Cycle - Cycle #18
End Cycle - Cycle #18
Half Cycle - Cycle #19
End Cycle - Cycle #19 -0.5
Half Cycle - Cycle #20
End Cycle - Cycle #20
First End
Second End
-1
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
0.02
First End - Walking Rate
Walking Rate (m/cycle)
-0.01
-0.02
0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
Cycle Number (#)