Next Article in Journal
Bioinspired Central Pattern Generator and T-S Fuzzy Neural Network-Based Control of a Robotic Manta for Depth and Heading Tracking
Next Article in Special Issue
CFD Investigation on Secondary Flow Characteristics in Double-Curved Subsea Pipelines with Different Spatial Structures
Previous Article in Journal
Laboratory Study of Turbulent Mass Exchange in a Stratified Fluid
Previous Article in Special Issue
Opportunities and Challenges to Develop Digital Twins for Subsea Pipelines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lateral Buckling of Subsea Pipelines Triggered by Sleeper with a Nonlinear Pipe–Soil Interaction Model

1
Ocean College, Zhejiang University, Zhoushan 316021, China
2
Centre for Marine Technology and Ocean Engineering (CENTEC), Instituto Superior Técnico, Universidade de Lisboa, 1049-001 Lisbon, Portugal
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2022, 10(6), 757; https://doi.org/10.3390/jmse10060757
Submission received: 22 April 2022 / Revised: 20 May 2022 / Accepted: 22 May 2022 / Published: 30 May 2022
(This article belongs to the Special Issue Subsea Pipelines)

Abstract

:
Buckle-initiation techniques, such as sleepers, are usually installed to trigger lateral buckling at pre-designated locations to release the axial compressive forces induced by thermal loading. Taking the nonlinear pipe–soil interaction model into account, a mathematical model is proposed to investigate the lateral buckling of subsea pipelines triggered by a sleeper. The numerical solution is validated by comparing the model with solutions in the literature, and the model shows good agreement. The discrepancy between them is analysed by presenting the effect of mobilisation distance during buckling. The influence of the breakout resistance, sleeper height, and sleeper friction coefficient on the buckled configuration, post-buckling behaviour, and minimum critical temperature difference is discussed parametrically. The results show that the deformation of the buckled pipeline shrinks, and both the minimum critical temperature difference and the maximum stress along the buckled pipeline enlarge when the nonlinear pipe–soil interaction model is incorporated. However, the influence of the nonlinear pipe–soil interaction reduces with increasing sleeper height.

1. Introduction

Subsea pipelines may buckle laterally due to the excessive axial compressive force due to high-temperature and high-pressure conditions. Lateral buckling occurs when the axial compressive force reaches critical levels. Lateral buckling, if not controlled, can lead to serious accidents involving local buckling, fracture, and fatigue [1]. To control this phenomenon, buckle initiation techniques, such as sleepers, are employed along pipelines to trigger buckles at predesigned locations. A sleeper is a pipe segment that is installed underneath and perpendicular to the pipeline, which typically has a low friction surface to reduce the lateral friction force. Thus, the pipeline is uplifted vertically. A combination of the vertical out-of-straightness and low lateral resistance results in reduced critical buckling force. When the sleeper is used as the buckle initiation facility, part of the pipeline is suspended. The pipeline segment at the end of the suspended section has a larger embedment into the seabed, since a vertical concentrated force exists. This embedment affects the lateral breakout resistance, which is a key design parameter governing the initiation of the lateral buckle. Thus, a nonlinear pipe–-soil interaction model is considered in the mathematical model to investigate the influence of breakout resistance on post-buckling behaviour.
The global buckling of subsea pipelines was investigated by numerous researchers. Hobbs’ solutions for a straight pipeline were derived by assuming specific buckling mode shapes and constant lateral soil resistance [2]. Based on this, an analytical model was proposed by Taylor and Gan [3] with a consideration of initial imperfection. A simplified analytical model was proposed by Croll [4] for upheaval buckling. The interaction between propagation buckling and global buckling in subsea pipelines was investigated by Karampour et al. [5]. This interaction leads to a significant reduction in buckle design capacity. The lateral buckling of imperfect pipelines was studied by Liu et al. [6], using FEM. The analytical solutions for the high-order lateral buckling of a pipeline with symmetric and anti-symmetric initial imperfection were derived by Hong et al. [7] and Liu et al. [8], respectively. Lateral buckling was investigated by Konuk [9,10] with coupled lateral and axial pipe–soil interactions. Zhang et al. [11,12] derived unified formulas for the critical buckling forces of the upheaval and lateral buckling of subsea pipelines with different types of initial imperfection. The influence of the pipe length on the lateral buckling behaviour of imperfect pipelines was investigated through FEM [13].
More recently, researchers investigated the influence of the nonlinear pipe–soil interaction model on lateral buckling. Zeng and Duan [14] used a quintic polynomial formula to simulate nonlinear pipe–soil interactions. Incorporating the tri-linear pipe–soil interaction model, Chee et al. [15] investigated the effect of imperfections on the buckling response through FEM. Considering both the initial imperfection and the nonlinear lateral soil resistance model, the critical force of the lateral buckling was analysed by assuming that the length of the buckled region equals the wavelength of initial imperfection [16].
To increase the reliability of buckle formation predictions, buckle initiation facilities were incorporated into the mathematical models. A single buoyancy or distributed buoyancy with a specific length installed along the pipeline was considered to derive some simple analytical solutions [17]. The critical load of lateral buckling triggered by a single buoyancy was investigated by Shi and Wang [18]. Moreover, dual distributed buoyancy sections with a gap between them were employed to initiate lateral buckling [19]. A new way to trigger lateral buckling is to introduce a pre-deformed section along the pipeline before installation [20]. Lateral buckling triggered by a sleeper was investigated experimentally by Silva-Junior et al. [21] and de Oliveira Cardoso and Solano [22]. Bai et al. [23] studied the lateral buckling triggered by dual sleepers through FEM. Analytical solutions for antisymmetric buckling modes triggered by a sleeper were obtained by Wang and Tang [24]. They found that the symmetric buckling mode was more likely to occur with lower sleeper friction or smaller sleeper height. Hong and Liu [25] investigated the vertical deflection of a pipeline on a sleeper by FEM.
By assuming constant lateral soil resistance, analytical solutions were derived for the lateral thermal buckling triggered by a sleeper in [26]. In practice, the pipeline always has an initial embedment into the soil and the lateral soil resistance is not constant. However, there are no studies about lateral thermal buckling triggered by sleepers that consider nonlinear lateral soil resistance.
The innovative aspect of this study is the nonlinear pipe–soil interaction model that is incorporated into the governing equations. In previous published studies about lateral buckling triggered by sleepers, the lateral soil resistance, f ( w 2 ) , is assumed to be constant. However, in this study, this function is nonlinear and it includes the effect of breakout resistance. This is due to the fact that in practice, pipe–soil interactions are nonlinear.

2. Mathematical Modelling

To avoid rogue buckles along subsea pipelines, buckle-initiation techniques, such as installing sleepers along the pipeline, are usually employed to trigger the pipeline to buckle in a controlled way at the predesignated location. For a pipeline laid on a sleeper and subjected to a temperature difference T 0 , the axial compressive force is accumulated. The axial compressive force, P 0 , is expressed as
P 0 = E A α T 0
where E is the elastic modulus, A is the cross-sectional area of the pipeline, and α is the coefficient of linear thermal expansion.
When P 0 is larger than the critical value, lateral buckling can be triggered at sleeper. The configuration and load distribution of lateral buckling are illustrated in Figure 1. From Figure 1a, it is clear that part of the pipe segment within l 1 x l 1 is uplifted by the sleeper. In l 1 x l 1 , the soil resistances are zero. However, there are concentrated contact forces F s between pipeline and sleeper at the sleeper and F t between the pipe and seabed at the end of the suspended section, respectively. The vertical configuration of the pipeline laid on a sleeper was solved by Wang et al. [26]. From their derivation, F s and F t can be expressed as
F s = 4 3 W p i p e l 1   and   F t = 1 3 W p i p e l 1
where W p i p e is the submerged weight per unit length and l 1 is the half-length of the free span, solved by
l 1 = 72 E I v o m W p i p e 4
where I is the moment of inertia and v o m is the sleeper height. Therefore, the value of l 1 can be obtained by Equation (3) when v o m is specified. Furthermore, F s and F t can be solved.
After the pipeline buckles, additional pipe coming from the thermal expansion is fed into the buckled section. Therefore, the axial force reduces partially due to the release of axial strain (see Figure 2). The axial force within the suspended region l 1 < x < l 1 , denoted by P , is constant. At x = ± l 1 , there is a jump in axial force with an amplitude of f A t induced by F t . Within the region where the pipeline makes contact with the seabed, the axial force increases because of the restraint of axial soil resistance. The axial force will reach P 0 at x = ± l s . From Figure 2, the axial force distribution P ¯ ( x ) is
P ¯ ( x ) = { P   ( 0 x < l 1 ) P + f A t + f A ( x l 1 )   ( l 1 x l s )
where f A = μ A W p i p e is the axial soil resistance per unit length and μ A is the axial friction coefficient. The force f A t = μ A F t is induced by F t .
The axial force at x = l s is
P ¯ ( l s ) = P 0 = P + f A t + f A ( l s l 1 )
Linear beam theory is used to simulate pipeline buckling. Thus, the equilibrium equations governing lateral deformation are [27]:
{ E I d 4 w 1 d x 4 + P d 2 w 1 d x 2 = 0       ( 0 x < l 1 ) E I d 4 w 2 d x 4 + P ¯ ( x ) d 2 w 2 d x 2 = f ( w 2 )       ( l 1 x l 2 )
where w 1 and w 2 are lateral deflections, E I is bending stiffness, and f ( w 2 ) is the nonlinear lateral soil resistance determined by the nonlinear pipe–soil interaction model, as shown in Figure 3. Here, the variation in the axial force within the buckled region l 1 x l 2 is ignored when solving lateral deformations. This assumption is acceptable [26,28].
Here, the nonlinear pipe–soil interaction model proposed by Chatterjee et al. [29] is employed. It can simulate breakout resistance and is given by
μ = w | w | ( μ b r k ( 1 e a 1 ( | w | D ) a 2 ) + ( μ r e s μ b r k ) ( 1 e a 3 ( | w | D ) a 4 ) )
where μ , μ b r k and μ r e s are the equivalent friction coefficients, and D is the external diameter of the pipeline. The quantities f ( w ) = μ W p i p e , F b r k = μ b r k W p i p e and F r e s = μ r e s W p i p e are, therefore, respectively, the nonlinear lateral soil resistance, the breakout resistance, and the residual resistance. The value of the coefficient a 3 in [29] is given as
a 3 = a 5 ( W p i p e V m a x ) + a 6
where W p i p e is the weight of the pipe and V m a x is the vertical bearing capacity. The values of a 5 and a 6 are calculated by
a 5 = 8.2 v i n i t D 4.9 ,   a 6 = 5.8 v i n i t D + 4.5
where v i n i t is the initial embedment of the pipe into the soil. In [29], a 1 = 25 and a 4 = 1.5 are employed, but a 2 = 1 is used here in order to have a finite linear resistance, which is physically realistic. V m a x = 5 W p i p e and v i n i t = 0.3 D are adopted so that a 3 = 2.272 , and set μ r e s = 0.5.
Due to symmetry, half a pipeline is considered. The slope of the deflection at x = 0 is zero, while the shear force f o w = μ s F s / 2 at x = 0 is induced by the friction force μ s F s . Here, μ s is the friction coefficient between pipeline and sleeper. The displacement, slope, and moment at x = l 2 are also zero. The boundary conditions at x = 0 and x = l 2 are
{ d w 1 d x ( 0 ) = 0 d 3 w 1 d x 3 ( 0 ) + f o w E I = 0 w 2 ( l 2 ) = 0 d w 2 d x ( l 2 ) = 0 d 2 w 2 d x 2 ( l 2 ) = 0
The displacement, slope, and bending moment must be continuous at the touchdown point x = l 1 , while there is a jump in shear force with an amplitude of f t = μ r e s F t at x = l 1 induced by the force F t . Thus, additional conditions at x = l 1 are
{ w 1 ( l 1 ) = w 2 ( l 1 ) d w 1 d x ( l 1 ) = d w 2 d x ( l 1 ) d 2 w 1 d x 2 ( l 1 ) = d 2 w 2 d x 2 ( l 1 ) d 3 w 1 d x 3 ( l 1 ) = d 3 w 2 d x 3 ( l 1 ) + f t E I
With Equations (10) and (11), the nonlinear governing equations are solved numerically by the shooting method [30]. Once the lateral deflections are known, the geometric shortening u 2 is obtained by
u 2 = 1 2 0 l 1 ( d w 1 d x ) 2 d x + 1 2 l 1 l 2 ( d w 2 d x ) 2 d x
The following compatibility condition is employed to link the lateral deflection and the thermal loading induced deflection:
u 1 = u 2
where u 1 is thermal expansion in 0 < x < l s .
We have
u 1 = 0 l s Δ P ¯ ( x ) E A d x
where Δ P ¯ ( x ) = P 0 P ¯ ( x ) .
Thus, this leads to
u 1 = f A ( l s l 1 ) 2 2 E A + ( P 0 P ) l 1 E A
The following formula is obtained by combining Equations (5), (13), and (15).
l s = 1 3 l 1 2 + 2 E A u 2 f A
With Equations (5) and (16), one finally obtains
P 0 = P + f A ( 1 3 l 1 2 + 2 E A u 2 f A 2 3 l 1 )
The bending moment is obtained by
M = E I d 2 w d x 2
where w stands for w 1 or w 2 , and the bending stress σ M is
σ M = M D 2 I
The maximum stress is
σ m = σ P + σ M m
where the stresses σ P and σ M m , induced by axial force P and maximum bending moment M m , respectively, are
{ σ P = P A σ M m = | M m D 2 I |

3. Results

The mathematical model is validated by comparing it with the analytical solution in [26], and the discrepancy between them is discussed. Next, the influence of μ b r k , v o m and μ s is analysed. The results are obtained by employing the analytical formulation developed in Section 2 and taking the parameters in Table 1.
One should note that only the analytical solutions in Figure 4 and Figure 5 come from [26], which is used to validate the numerical results obtained in this study. In [26], the lateral soil resistance is assumed to be constant, while in the present study, nonlinear lateral soil resistance is considered. Moreover, in [26], analytical solutions are obtained due to the assumption of constant lateral soil resistance. In the present study, because the lateral soil resistance is nonlinear, Equation (6) cannot be solved analytically. Thus, the shooting method is used to solve Equation (6) to obtain the numerical results.

3.1. Validation

The solutions obtained in this study were validated by comparing them with the analytical solutions in [26], as shown in Figure 4. An error analysis is shown in Figure 5. In Figure 4, the analytical solutions are obtained by using the formulas derived in Wang et al. [26] with constant lateral soil resistance. To compare with the analytical solutions, the numerical solutions shown in Figure 4 are obtained by assuming μ b r k = 0.5 . For μ b r k = 0.5 , the nonlinear pipe–soil interaction model is reduced to elastic-plastic (see Figure 3). The numerical solutions for μ b r k = 2.0 are also illustrated in Figure 4 to show the influence of the nonlinear pipe–soil interaction.
In Figure 4, there are two branches for each solution, which are denoted as m-b and m-c. The temperature difference at m, i.e., T m , is called the minimum critical temperature difference, since solutions only exist for T 0 > T m .
From Figure 4, the numerical solutions for μ b r k = 0.5 are in good agreement with the analytical solutions, except that there is a slight discrepancy between them around T m . This discrepancy comes from the difference in mobilization distance. For the rigid-plastic model, the resistance is always constant (see Figure 3). For the elastic-plastic model, the lateral soil resistance increases from zero to residual resistance gradually (see Figure 3).
In Figure 4, the discrepancy between the analytical and numerical solutions reduces as the temperature difference increases. The reason for this is that the displacement amplitude increases along with the temperature difference, so that more pipe sections fall into the region of constant lateral soil resistance for numerical solutions.
A more detailed error analysis is illustrated in Figure 5. Figure 6 shows that the mobilization distance is controlled by the parameter a 1 . The mobilization distance is the distance that the lateral resistance reaches μ b r k . In Figure 6, μ b r k = μ r e s = 0.5 , so the mobilization distance in Figure 6 is the distance at which μ reaches 0.5. The elastic-plastic model approaches the rigid-plastic model for larger a 1 , since the mobilization distance becomes smaller. In Figure 5, the discrepancy between the analytical and numerical solutions becomes smaller for larger a 1 . In Figure 4 and Figure 5, the discrepancy in half-buckled length l 2 is larger than in the other parameters. This is because the deflection at the ends of the buckled section is small, and it is affected by the mobilization distance. For the remaining parameters, the discrepancy between the analytical and numerical solutions is small enough.
In Figure 4a, T m becomes larger when the nonlinear pipe–-soil interaction model is considered, which means that lateral buckling can only be triggered at higher temperature differences. At the same temperature difference, both w m and l 2 become smaller when considering the nonlinear pipe–soil interaction model (see Figure 4a,e). Thus, the use of an additional pipe to create lateral deflection, which comes from thermal expansion, also reduces (see Figure 4d), so that l s decreases, as shown in Figure 4f. However, P within the buckled section becomes larger due to the restriction of the breakout resistance. Moreover, at the same temperature difference, the maximum stress σ m becomes larger when nonlinear pipe–soil interaction is considered. This means that the maximum stress is underestimated when assuming the lateral soil resistance to be constant.

3.2. Parametric Study

3.2.1. Influence of μ b r k

The influence of μ b r k on the buckled configuration, post-buckling behaviour, and minimum critical temperature difference T m is shown in Figure 7, Figure 8 and Figure 9, respectively.
In Figure 7, the dashed curves are the pipe sections in contact with the seabed, and the solid curves are the pipe sections suspended due to the existence of the sleeper. In Figure 7a, both the touchdown and suspended pipe segments shrink with larger μ b r k . In Figure 7b, there are two extrema of bending stress in the positive direction. Another three extrema of bending stress occur around the sleeper. One local minimum (absolute value) of bending stress appears at the sleeper, while there are two other local maxima (absolute value) of bending stress close to the sleeper. The occurrence of the local minimum (absolute value) of bending stress at the sleeper is induced by the friction force between the pipeline and the sleeper. For each specific μ b r k , the maximum bending stress is located at the local maxima (absolute value) of bending stress close to the sleeper. For larger μ b r k , all the extrema of bending stress in both positive and negative directions become larger.
In Figure 8a, T m is larger for larger μ b r k . A more detailed analysis on the influence of μ b r k on T m is shown in Figure 9, which shows that T m increases with increasing μ b r k for specific values of v o m and μ s , and the increasing rate of T m reduces with the increase in μ b r k . In Figure 9a, under the same μ b r k , T m is larger for smaller v o m . The increasing rate of T m with increasing μ b r k is also larger for smaller v o m . The reason is that since there are less length of suspended pipeline and larger length of touchdown pipeline with the smaller v o m , the breakout resistance has a larger influence on the initiation of lateral buckling. In Figure 9b, under the same μ b r k , T m becomes larger for larger μ s . The increase in the rate of T m along with the increasing μ b r k remains almost the same for different values of μ s . The reason for this is that the friction force between the pipeline and the sleeper becomes larger for larger μ s ; however, the value of μ s has no influence on the lengths of the suspended or touchdown pipeline segments.
In Figure 8a,e, both the displacement amplitude w m and the half-buckled length l 2 increase with the increasing T 0 , and need more thermal expansion u 1 to form the buckled deflection (see Figure 8d). Thus, larger l s is required for larger T 0 , as shown in Figure 8f. The maximum stress also increases with increasing T 0 (see Figure 8b) since large deflection occurs; however, the axial force P reduces with increasing T 0 (see Figure 8c).
In Figure 8a,e, at a specific temperature difference, both w m and l 2 become smaller for larger μ b r k . The reason for this is that since the breakout resistance is larger for larger μ b r k , the pipeline is subjected to greater lateral soil resistance. The deflection shrinks with larger μ b r k , as shown in Figure 7a. Therefore, both u 1 and l s become smaller with larger μ b r k at the same temperature difference (see Figure 8d,f).
However, P becomes larger with larger μ b r k (see Figure 8c). The reason for this is that the reduction in the axial force reduces, since the greater breakout resistance restricts the deflection of the buckled pipeline. Moreover, σ m along the buckled pipeline becomes larger for larger μ b r k (see Figure 8b). After considering the nonlinear pipe–soil interaction model, both T m and σ m in the pipeline became larger. When the nonlinear pipe–soil interaction model is not included, lateral buckling may fail to be triggered by the sleeper, and the maximum stress along the buckled pipeline may exceed the allowable stress in the design.

3.2.2. Influence of v o m

The influence of v o m on the buckled configuration, post-buckling behaviour, and minimum critical temperature difference T m are shown in Figure 10, Figure 11 and Figure 12, respectively.
From Figure 10a, it is clear that the length of the suspended pipe segment becomes larger with larger sleeper heights, v o m . The deflection of the buckled pipeline enlarges with larger v o m . Because the soil resistance for the suspended pipeline is zero, the buckled pipeline has less restriction from the seabed foundation with larger v o m . The deflection of the buckled pipeline enlarges with larger v o m , but both the local minimum and the local maximum (absolute value) of the bending stress become smaller with larger v o m , as shown in Figure 10b. This is because the deflection of the buckled pipeline is more benign with larger v o m .
In Figure 11a, T m becomes smaller with larger v o m . Figure 12 illustrates the influence of v o m on T m in detail. In Figure 12, T m decreases with increasing v o m for specific values of μ b r k and μ s , and the decreasing rate of T m reduces with increasing v o m . In Figure 12a, under the same v o m , T m becomes larger for larger μ b r k . The influence of μ b r k on T m becomes smaller for larger v o m , since the length of suspended pipeline with zero soil resistance is greater. In Figure 12b, under the same v o m , T m is larger with larger μ s . The decreasing rate of T m with increasing v o m becomes smaller with larger μ s . The influence of μ s on T m is larger with larger v o m . The reason for this is that the concentrated contact force between the pipeline and the sleeper becomes larger with larger v o m , so that the friction force between the pipeline and the sleeper becomes larger with larger v o m . Therefore, an effective way to reduce T m is to increase the sleeper height v o m ; however, the corresponding weakness is that the suspended pipeline will be longer, which may lead to vortex-induced vibration.
In Figure 11a,e, at a specific temperature difference, both w m and l 2 become larger with larger v o m . This is because the length of the suspended pipeline with zero soil resistance increases with increasing v o m , as shown in Figure 10a. There is less restriction from the seabed foundation with larger v o m . Due to the larger deflection with larger v o m , the requirement of additional pipes to feed into the buckled section increases, which creates the need for more thermal expansion (see Figure 11d) and a longer feed-in region (see Figure 11f). The axial force P becomes smaller with larger v o m , since a larger deflection occurs to release more axial force, as shown in Figure 11c. The maximum stress σ m along the buckled pipeline reduces with larger μ b r k (see Figure 11b).

3.2.3. Influence of μ s

The influence of μ s on the buckled configuration, post-buckling behaviour, and minimum critical temperature difference T m are shown in Figure 13, Figure 14 and Figure 15, respectively.
In Figure 13a, the deflection of the buckled pipeline shrinks with larger μ s . This is because the friction force between the pipeline and the sleeper becomes larger with larger μ s , which restricts the deflection of the buckled pipeline. In Figure 13b, the extrema of the bending stress in the positive direction becomes slightly larger with larger μ s . However, both the local minimum and the local maximum (absolute value) of the bending stress close to the sleeper in the negative direction become smaller with larger μ s , as shown in Figure 13b. With larger μ s , the difference between the local minimum and the local maximum of the bending stress close to the sleeper becomes larger. Taking μ s = 0.4 as an example, it is clear that the local minimum (absolute value) of the bending stress is smaller than the local maximum (absolute value) of the bending stress.
In Figure 14a, T m is larger with larger μ s . The effect of μ s on T m is illustrated in Figure 15, with different values of μ b r k and v o m . In Figure 15, T m increases with increasing μ s for specific values of μ b r k and v o m , and the increasing rate of T m slightly reduces with increasing μ s . The friction force between the sleeper and the pipeline becomes larger with larger μ s , which makes it more difficult to trigger the lateral buckling. In Figure 15a, at the same μ s , T m becomes larger with larger μ b r k . The increasing rate of T m with increasing μ s is similar for different values of μ b r k . In Figure 15b, under the same μ s , T m becomes smaller for larger v o m . The increasing rate of T m with increasing μ s is larger for larger v o m . The influence of v o m on T m gradually reduces with increasing μ s . Thus, the friction coefficient between the sleeper and the pipeline μ s should be carefully controlled. When the value of μ s is too large, such as μ s = 0.6 , T m is barely affected by the sleeper height v o m .
In Figure 14a,e, under a specific T 0 , both w m and l 2 reduce with larger μ s . This is because, since the friction force between the sleeper and the pipeline becomes larger with larger μ s , the deflection of the buckled pipeline is restricted by the larger resistance between the sleeper and the pipeline. Thus, the requirements of both u 1 and l s decrease with larger μ s , as shown in Figure 14d,f. Due to the restriction of the larger friction force between the sleeper and the pipeline, the axial force P increases with increasing μ s , as shown in Figure 14c. However, the maximum stress σ m along the buckled pipeline reduces with larger μ s (see Figure 14b), which is induced by the decrease in the maximum bending stress (absolute value) along the buckled pipeline.

4. Conclusions

Through a consideration of the nonlinear pipe–soil interaction model, a mathematical model was proposed to simulate the lateral buckling of subsea pipelines triggered by a sleeper. The model was solved numerically and validated by comparing its predictions with the analytical solutions from [26]. The discrepancy between the numerical and analytical solutions was analysed through the discussion of the mobilization distance. A detailed parametric analysis was presented to show the effect of the breakout resistance, sleeper height, and sleeper friction coefficient on the buckling behaviour of a pipeline laid on a sleeper. The conclusions are:
(i)
The discrepancy between the numerical and analytical solutions comes from the difference between the elastic-plastic and rigid-plastic pipe–soil interaction models, which reduces with decreasing mobilization difference in the elastic-plastic pipe–soil interaction model.
(ii)
When the nonlinear pipe–soil interaction model is taken into account, both the displacement amplitude and the buckled length reduce due to the occurrence of breakout resistance, which decreases further with increasing breakout resistance. However, both the axial force and the maximum stress, along with the buckled pipeline, increase, and increase further with increasing breakout resistance.
(iii)
The deflection of the buckled pipeline enlarges as the sleeper height increases and shrinks as the sleeper friction coefficient increases. The axial force decreases with increasing sleeper height and increases with increasing sleeper friction coefficient. Moreover, the maximum stress along the buckled pipeline decreases with increasing sleeper height and with decreasing sleeper friction coefficient.
(iv)
The minimum critical temperature difference increases with increasing breakout resistance and sleeper friction coefficient, and decreases with increasing sleeper height. The influence of the breakout resistance on the minimum critical temperature difference gradually reduces with increasing sleeper height. Moreover, the sleeper height has little effect on the minimum critical temperature difference when the sleeper friction coefficient is large enough.
In conclusion, it is better to incorporate the nonlinear pipe–soil interaction model into the mathematical model when simulating the lateral buckling of subsea pipelines triggered by a sleeper, since both the minimum critical temperature difference and the maximum stress increase. Moreover, both the sleeper height and the sleeper friction coefficient should be carefully selected and controlled.

Author Contributions

Conceptualization, Z.W. and C.G.S.; methodology, Z.W.; software, Z.W.; validation, Z.W.; investigation, Z.W.; writing—original draft preparation, Z.W.; writing—review and editing, C.G.S.; supervision, C.G.S.; project administration, C.G.S.; funding acquisition, C.G.S. All authors have read and agreed to the published version of the manuscript.

Funding

Zhenkui Wang would like to acknowledge that the work described in this paper was funded by the National Natural Science Foundation of China (grant number: 52001229). This work contributes to the Strategic Research Plan of the Centre for Marine Technology and Ocean Engineering (CENTEC), which is financed by the Portuguese Foundation for Science and Technology (Fundação para a Ciência e Tecnologia—FCT) under contract UIDB/UIDP/00134/2020.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. DNV-RP-F110; Global Buckling of Submarine Pipelines. Det Norske Veritas: Oslo, Norway, 2019.
  2. Hobbs, R.E. In-service buckling of heated pipelines. J. Transp. Eng. 1984, 110, 175–189. [Google Scholar] [CrossRef]
  3. Taylor, N.; Gan, A.B. Submarine pipeline buckling-imperfection studies. Thin-Walled Struct. 1986, 4, 295–323. [Google Scholar] [CrossRef]
  4. Croll, J.G.A. A simplified model of upheaval thermal buckling of subsea pipelines. Thin-Walled Struct. 1997, 29, 59–78. [Google Scholar] [CrossRef]
  5. Karampour, H.; Albermani, F.; Veidt, M. Buckle interaction in deep subsea pipelines. Thin-Walled Struct. 2013, 72, 113–120. [Google Scholar] [CrossRef] [Green Version]
  6. Liu, R.; Xiong, H.; Wu, X.; Yan, S. Numerical studies on global buckling of subsea pipelines. Ocean Eng. 2014, 78, 62–72. [Google Scholar] [CrossRef]
  7. Hong, Z.; Liu, R.; Liu, W.; Yan, S. Study on lateral buckling characteristics of a submarine pipeline with a single arch symmetric initial imperfection. Ocean Eng. 2015, 108, 21–32. [Google Scholar] [CrossRef]
  8. Liu, R.; Wang, X. Lateral global buckling high-order mode analysis of a submarine pipeline with imperfection. Appl. Ocean. Res. 2018, 73, 107–126. [Google Scholar] [CrossRef]
  9. Konuk, I. Coupled lateral and axial soil-pipe interaction and lateral buckling Part II: Solutions. Int. J. Solids Struct. 2018, 132–133, 127–152. [Google Scholar] [CrossRef]
  10. Konuk, I. Coupled lateral and axial soil-pipe interaction and lateral buckling Part I: Formulation. Int. J. Solids Struct. 2018, 132–133, 114–126. [Google Scholar] [CrossRef]
  11. Zhang, X.; Duan, M. Prediction of the upheaval buckling critical force for imperfect submarine pipelines. Ocean Eng. 2015, 109, 330–343. [Google Scholar] [CrossRef]
  12. Zhang, X.; Guedes Soares, C.; An, C.; Duan, M. An unified formula for the critical force of lateral buckling of imperfect submarine pipelines. Ocean Eng. 2018, 166, 324–335. [Google Scholar] [CrossRef]
  13. Liu, R.; Li, C. Determinate dimension of numerical simulation model in submarine pipeline global buckling analysis. Ocean Eng. 2018, 152, 26–35. [Google Scholar] [CrossRef]
  14. Zeng, X.; Duan, M. Mode localization in lateral buckling of partially embedded submarine pipelines. Int. J. Solids Struct. 2014, 51, 1991–1999. [Google Scholar] [CrossRef] [Green Version]
  15. Chee, J.; Walker, A.; White, D. Effects of variability in lateral pipe-soil interaction and pipe initial out-of-straightness on controlled lateral buckling of pre-deformed pipeline. Ocean Eng. 2019, 182, 283–304. [Google Scholar] [CrossRef]
  16. Zhang, X.; Guedes Soares, C. Lateral buckling analysis of subsea pipelines on nonlinear foundation. Ocean Eng. 2019, 186, 106085. [Google Scholar] [CrossRef]
  17. Peek, R.; Yun, H. Flotation to trigger lateral buckles in pipelines on a flat seabed. J. Eng. Mech. 2007, 4, 442–451. [Google Scholar] [CrossRef]
  18. Shi, R.; Wang, L. Single buoyancy load to trigger lateral buckles in pipelines on a soft seabed. J. Eng. Mech. 2015, 141, 1–7. [Google Scholar] [CrossRef]
  19. Wang, Z.; Tang, Y. Antisymmetric thermal buckling triggered by dual distributed buoyancy sections. Mar. Struct. 2020, 74, 102811. [Google Scholar] [CrossRef]
  20. Chee, J.; Walker, A.; White, D. Controlling lateral buckling of subsea pipeline with sinusoidal shape pre-deformation. Ocean Eng. 2018, 151, 170–190. [Google Scholar] [CrossRef] [Green Version]
  21. Silva-Junior, H.C.; Cardoso, C.O.; Carmignotto, M.A.P.; Zanutto, J.C. Reduced Model Device of Solutions to Control Thermal Buckling Effects in HP-HT Subsea Pipelines (OMAE2008-57637). In Proceedings of the International Conference on Ocean, Offshore and Arctic Engineering, Estoril, Portugal, 15 June 2008. [Google Scholar]
  22. de Oliveira Cardoso, C.; Solano, R.F. Performed of triggers to control thermal buckling of subsea pipelines using reduced scale model (ISOPE-I-15-445). In Proceedings of the International Offshore and Polar Engineering Conference, Honolulu, HI, USA, 21–26 June 2015. [Google Scholar]
  23. Bai, Q.; Qi, X.; Brunner, M. Global buckle control with dual sleepers in HP/HT pipelines (OTC-19888-MS). In Proceedings of the Offshore Technology Conference, Houston, TX, USA, 4–7 May 2009. [Google Scholar]
  24. Wang, Z.; Tang, Y. Analytical study on controlled lateral thermal buckling of antisymmetric mode for subsea pipelines triggered by sleepers. Mar. Struct. 2020, 71, 102728. [Google Scholar] [CrossRef]
  25. Hong, Z.; Liu, W. Modelling the vertical lifting deformation for a deep-water pipeline laid on a sleeper. Ocean Eng. 2020, 199, 107042. [Google Scholar] [CrossRef]
  26. Wang, Z.; Tang, Y.; van der Heijden, G.H.M. Analytical study of lateral thermal buckling for subsea pipelines with sleeper. Thin-Walled Struct. 2018, 122, 17–29. [Google Scholar] [CrossRef] [Green Version]
  27. Lagrange, R.; Averbuch, D. Solution methods for the growth of a repeating imperfection in the line of a strut on a nonlinear foundation. Int. J. Mech. Sci. 2012, 63, 48–58. [Google Scholar] [CrossRef]
  28. Wang, Z.; Tang, Y.; van der Heijden, G.H.M. Analytical study of distributed buoyancy sections to control lateral thermal buckling of subsea pipelines. Mar. Struct. 2018, 58, 199–222. [Google Scholar] [CrossRef] [Green Version]
  29. Chatterjee, S.; White, D.J.; Randolph, M.F. Numerical simulations of pipe-soil interaction during large lateral movements on clay. Géotechnique 2012, 62, 693–705. [Google Scholar] [CrossRef]
  30. Wolfram Mathematica, version 11.0; Wolfram Research, Inc.: Champaign, IL, USA, 2016.
Figure 1. Configuration and load distribution. (a) Vertical plane. (b) Horizontal plane.
Figure 1. Configuration and load distribution. (a) Vertical plane. (b) Horizontal plane.
Jmse 10 00757 g001
Figure 2. Axial compressive force distribution.
Figure 2. Axial compressive force distribution.
Jmse 10 00757 g002
Figure 3. Pipe–soil interaction models.
Figure 3. Pipe–soil interaction models.
Jmse 10 00757 g003
Figure 4. Validation. (a) Displacement amplitude w m . (b) Maximum stress σ m . (c) Axial force P . (d) Thermal expansion u 1 . (e) Half-length of buckled region l 2 . (f) Half-length of feed-in region l s . ( v o m = 0.1 m, μ s = 0.1 ) [26].
Figure 4. Validation. (a) Displacement amplitude w m . (b) Maximum stress σ m . (c) Axial force P . (d) Thermal expansion u 1 . (e) Half-length of buckled region l 2 . (f) Half-length of feed-in region l s . ( v o m = 0.1 m, μ s = 0.1 ) [26].
Jmse 10 00757 g004
Figure 5. Error analysis. (a) w m . (b) σ m . (c) P . (d) l 2 ( μ b r k = 0.5 , v o m = 0.1 m, μ s = 0.1 ) [26].
Figure 5. Error analysis. (a) w m . (b) σ m . (c) P . (d) l 2 ( μ b r k = 0.5 , v o m = 0.1 m, μ s = 0.1 ) [26].
Jmse 10 00757 g005
Figure 6. Pipe–soil interaction model with different a 1 .
Figure 6. Pipe–soil interaction model with different a 1 .
Jmse 10 00757 g006
Figure 7. Influence of μ b r k on configurations. (a) Deformed shapes. (b) Bending stresses ( v o m = 0.1 m, μ s = 0.1 , T 0 = 40   ° C ).
Figure 7. Influence of μ b r k on configurations. (a) Deformed shapes. (b) Bending stresses ( v o m = 0.1 m, μ s = 0.1 , T 0 = 40   ° C ).
Jmse 10 00757 g007
Figure 8. Influence of μ b r k on the buckling behaviour. (a) w m . (b) σ m . (c) P . (d) u 1 . (e) l 2 . (f) l s ( v o m = 0.1 m, μ s = 0.1 ).
Figure 8. Influence of μ b r k on the buckling behaviour. (a) w m . (b) σ m . (c) P . (d) u 1 . (e) l 2 . (f) l s ( v o m = 0.1 m, μ s = 0.1 ).
Jmse 10 00757 g008
Figure 9. (a) Influence of μ b r k on T m with different v o m ( μ s = 0.1 ). (b) Influence of μ b r k on T m with different μ s . ( v o m = 0.1 m).
Figure 9. (a) Influence of μ b r k on T m with different v o m ( μ s = 0.1 ). (b) Influence of μ b r k on T m with different μ s . ( v o m = 0.1 m).
Jmse 10 00757 g009
Figure 10. Influence of v o m on configurations. (a) Deformed shapes. (b) Bending stresses ( μ b r k = 2.0 , μ s = 0.1 , T 0 = 40   ° C ).
Figure 10. Influence of v o m on configurations. (a) Deformed shapes. (b) Bending stresses ( μ b r k = 2.0 , μ s = 0.1 , T 0 = 40   ° C ).
Jmse 10 00757 g010
Figure 11. Influence of v o m on the buckling behaviour. (a) w m . (b) σ m . (c) P . (d) u 1 . (e) l 2 . (f) l s ( μ b r k = 2.0 , μ s = 0.1 ).
Figure 11. Influence of v o m on the buckling behaviour. (a) w m . (b) σ m . (c) P . (d) u 1 . (e) l 2 . (f) l s ( μ b r k = 2.0 , μ s = 0.1 ).
Jmse 10 00757 g011
Figure 12. (a) Influence of v o m on T m with different μ b r k ( μ s = 0.1 ). (b) Influence of v o m on T m with different μ s ( μ b r k = 2.0 ).
Figure 12. (a) Influence of v o m on T m with different μ b r k ( μ s = 0.1 ). (b) Influence of v o m on T m with different μ s ( μ b r k = 2.0 ).
Jmse 10 00757 g012
Figure 13. v o m on configurations. (a) Deformed shapes. (b) Bending stresses ( μ b r k = 2.0 , v o m = 0.1 m, T 0 = 40   ° C ).
Figure 13. v o m on configurations. (a) Deformed shapes. (b) Bending stresses ( μ b r k = 2.0 , v o m = 0.1 m, T 0 = 40   ° C ).
Jmse 10 00757 g013
Figure 14. Influence of μ s on the buckling behaviour. (a) w m . (b) σ m . (c) P . (d) u 1 . (e) l 2 . (f) l s ( μ b r k = 2.0 , v o m = 0.1 m).
Figure 14. Influence of μ s on the buckling behaviour. (a) w m . (b) σ m . (c) P . (d) u 1 . (e) l 2 . (f) l s ( μ b r k = 2.0 , v o m = 0.1 m).
Jmse 10 00757 g014
Figure 15. (a) Influence of μ s on T m with different μ b r k ( v o m = 0.1 m). (b) Influence of μ s on T m with different v o m ( μ b r k = 2.0 ).
Figure 15. (a) Influence of μ s on T m with different μ b r k ( v o m = 0.1 m). (b) Influence of μ s on T m with different v o m ( μ b r k = 2.0 ).
Jmse 10 00757 g015
Table 1. Parameters.
Table 1. Parameters.
ParameterValueUnit
External diameter D 323.9mm
Wall thickness t 12.7 mm
Elastic modulus E 206GPa
Steel density ρ 7850 kg / m 3
Coefficient of thermal expansion α 1.1 × 10 5 ° C
Axial friction coefficient μ A 0.5---
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, Z.; Guedes Soares, C. Lateral Buckling of Subsea Pipelines Triggered by Sleeper with a Nonlinear Pipe–Soil Interaction Model. J. Mar. Sci. Eng. 2022, 10, 757. https://doi.org/10.3390/jmse10060757

AMA Style

Wang Z, Guedes Soares C. Lateral Buckling of Subsea Pipelines Triggered by Sleeper with a Nonlinear Pipe–Soil Interaction Model. Journal of Marine Science and Engineering. 2022; 10(6):757. https://doi.org/10.3390/jmse10060757

Chicago/Turabian Style

Wang, Zhenkui, and C. Guedes Soares. 2022. "Lateral Buckling of Subsea Pipelines Triggered by Sleeper with a Nonlinear Pipe–Soil Interaction Model" Journal of Marine Science and Engineering 10, no. 6: 757. https://doi.org/10.3390/jmse10060757

APA Style

Wang, Z., & Guedes Soares, C. (2022). Lateral Buckling of Subsea Pipelines Triggered by Sleeper with a Nonlinear Pipe–Soil Interaction Model. Journal of Marine Science and Engineering, 10(6), 757. https://doi.org/10.3390/jmse10060757

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop