1. Introduction
The underground pipeline system is crucial for meeting essential material needs in daily life and industrial production. Ensuring the structural safety of buried pipelines, including detection, maintenance, and management for long-term safe operation, is a significant challenge that primarily depends on two factors: structural integrity and effective soil-pipe interaction [
1]. Furthermore, the two factors are interdependent. Uneven foundation conditions can lead to reduced support for buried pipelines, resulting in structural damage to the pipes [
2,
3]. Additionally, pipe damage, such as leakage and cracking, can erode the underlying soil, further compromising soil stability [
4]. It has been reported that failures in underground pipelines are a leading cause of urban road collapse accidents [
5]. However, substantial attention has traditionally been focused on the structural strength and integrity of pipelines, often based on the assumption of uniform bedding [
6].
During the service life, pipelines often experience non-uniform longitudinal soil support due to various factors. In fact, the non-uniform bedding support can induce cracks or local subsidence of buried pipe, thereby diminishing structural-bearing capacity [
7,
8]. Consequently, there is growing attention to analyzing the damage and risks to pipeline structures posed by non-uniform soil-pipe interaction including local voids [
9,
10]. Voids around underground pipelines are generally caused by internal factors or non-structural external factors [
11]. External causes include groundwater or engineering activities that wash away finer particles from loosely filled soil, leading to erosion voids as water infiltrates and transports these particles downstream. Internal causes primarily stem from improper installation or structural damage due to aging, overload, or adverse conditions, manifesting as drips, leaks, and seepage. These issues cause continuous infiltration and erosion of soil particles, forming voids around the pipeline. For rigid pipelines, the loss of soil support due to such voids reduces resistance and affects pipe-soil interaction. Talesnick and Baker [
12] reported a void at the pipe invert during excavation, observing a significant gap. Additionally, economical and effective repair technologies for the defects around the pipe have been proposed to prevent greater failures and consequences [
13,
14,
15].
Prestressed Concrete Cylinder Pipe (PCCP), shown in
Figure 1, is an indispensable component in water supply networks and is extensively utilized in trans-regional water diversion projects. Due to the characteristics of large diameter and pressure water delivery, the consequences will be disastrous and catastrophic once failure occurs. Moreover, it is a composite segmented pipeline consisting of concrete, mortar, steel cylinder, and prestressing wires (see
Figure 2), with the steel wires serving as the primary load-bearing component. Once cracks appear in the mortar, it can provide conditions for external corrosive substances to contact the prestressed steel wire, posing a serious threat to operational safety and service life [
16]. As long-span buried pipelines, PCCPs may experience uneven pipe-soil interaction due to local voids caused by internal or external factors during their service [
17]. Case studies also have demonstrated that a reduction in longitudinal pipe-soil interaction in segmented buried pipelines can lead to more severe failure consequences [
18,
19,
20,
21]. However, current PCCP design and research often assume a uniform foundation, primarily focusing on the structural response of pipes buried in uniform bedding [
22]. The response under uneven support conditions has received relatively little attention. Additionally, the spigot-and-bell joint is particularly vulnerable to damage due to uneven support, which can generate shear forces and angular stresses at the joint, leading to cracks or leaks. This can result in the pipeline reaching its service life prematurely, even before the structural strength is compromised. Therefore, it is essential to analyze the response and damage to segmented PCCP under uneven foundation conditions to inform prospective and preventive technology development.
In this study, the cause and potential harm of a local void around the buried pipeline are analyzed first. It then conducts a comparative qualitative analysis to assess the influence of the void position on shear and rotation angle using the beam-on-elastic-spring approximation model. Based on the results, the finite element model of the PCCP with a local soil void is developed to quantitatively evaluate the pipe response and material damage. This study aims to evaluate the impact of uneven pipe-soil interaction (localized void) on buried PCCP, and the findings could provide a foundation for assessing structural safety throughout their lifecycle.
2. Mechanical Analysis
The presence of soil voids results in localized loss of support and alters the original stress distribution in the pipeline-soil interaction. Specifically, voids beneath pipe joints can induce shear forces and rotational effects. In this section, the expressions of shear force and rotation angle with unilateral and bilateral voids are derived, respectively, using the mechanical beam-on-elastic-spring approximation model [
23]. The model considers two rigid pipes with rigidity significantly greater than the stiffness of the surrounding soil, connected by the moment-release joint. Shear forces develop to ensure that the vertical deformations of the two pipes at the joint remain consistent. A qualitative analysis investigates the impact of parameters such as void location, length, soil stiffness, and pipe diameter.
2.1. Analysis Model
- (1)
Unilateral void
The support and response diagram of a unilateral local soil void on the pipe joint is shown in
Figure 3. The physical balance and moment balance equation of the pipe with a unilateral soil void are established as follows:
where
k is soil stiffness;
l and
d is the pipe length and outside diameter, respectively;
,
is vertical displacement under vertical force load (
,
);
(
=
−
) is void length.
Then, the vertical displacement at the center of the pipe can be obtained as follows:
According to the moment balance equation, the respective joint angle (
,
), and the respective vertical displacement (
,
) of the joint without taking into account the interactions can be obtained as follows:
If the response of two beams is independent, the relative vertical movement of the joint is as follows:
- (2)
Void spanning the joint
When the void spans the joint (the void longitudinal length is both
), the schematic of soil local support and the mechanical response is shown in
Figure 4. The physical balance and moment balance equations are established as Equations (7) and (8), and the derivation of displacement and rotation angle is similar to the unilateral void mentioned above.
2.2. Shear Force and Rotation
When the relative displacement occurs between the spigot and bell end, the shear forces and the relative rotation angle of the joint occur. Regardless of the gasket deformation, the shear force is to eliminate most of the net shear displacement of the joint, that is, . The incremental deformations of the beam resulting from shear can be calculated using the following Equations (3)–(6), where .
Therefore, the shear force with a unilateral soil void can be evaluated from
where
.
Similarly, the relative rotation angle
is the difference between load angle and shear angle.
Considering the uniform load is caused by backfill over the pipes,
, the values for joint shear and rotation for uniform length under the case of the unilateral soil void become
where
,
.
Repeating the process above, the expressions with void spanning the joint can be derived by
According to Equations (13) and (14), the shear of the buried rigid pipeline with unilateral void is related to three factors, namely, the void length, the pipe length, and the overlying load (see
Figure 5a). Additionally, the relative rotation angle of the joint is also influenced by pipe diameter and foundation stiffness, besides the above factors. Both shear force and rotation angle are positively correlated with the overlying load. When the void spans the joint symmetrically, the shear force is zero, but the rotation angle is greater compared to a unilateral void and negatively correlated with the pipe diameter (as shown in
Figure 5b).
Since the relative rotation angle of the joint serves as a critical safety indicator for segmented pipelines, longitudinal non-uniform support caused by a local void spanning the joint is relatively unfavorable. PCCP pipelines are segmented with rigid pipes connected by flexible joints; therefore, the qualitative assessment results above serve as a reasonable reference. To quantitatively analyze the impact of a localized soil void on stress response and materials damage to PCCP, an FE simulation model including pipes and the surrounding soil was developed in the following section.
3. Numerical Analysis
In this section, the size effect of local void is investigated by establishing a buried PCCP nonlinear finite element model considering the pipe-soil interaction. Based on the pipe damage and soil failure results, the response characteristics of the buried PCCP under local soil voids are further summarized.
3.1. Finite Element Modeling
As known above, the adjacent buried pipes with a void spanning the joint can lead to a larger relative rotation angle, representing a more unfavorable condition. Moreover, the joint is prone to seepage and even leakage in adverse conditions, forming a void. Because of the difficulty of a void geometry measurement, a simple shape that is easy to simulate is used to study the influence based on the assumptions of previous scholars in numerical modeling [
3,
24,
25]. In this study, void angle (VA), void depth (VD), and void length (VL) are used as the dimension control parameters, as shown in
Figure 6. To study parameter sensitivity, the void angles are set as 20°, 40°, 60°, and 90°, with void depths of 50 and 100 and void lengths of 62, 125, and 250. For the convenience of expression, the format of the void size defined in this section is VA (°)-VD (mm)-VL (cm). The birth-death element technology in finite element is a common method in structural engineering simulation for the simulation of drilling, excavation, backfilling, and other processes, which deactivate or reactivate specific regions/sets through model changes. Therefore, the formation and evolution process of the void is simulated by the birth-death element method.
Based on the nonlinear analysis software ABAQUS2019, the finite element model of two adjacent PCCPs, including the bell-and-spigot joint, and their surrounding soil were assembled, as depicted in
Figure 7. The key parameters for the embedded PCCP were as follows: burial depth H = 2.5 m, inside diameter D = 3.8 m, work pressure P = 0.6 MPa, and single pipe length L = 5 m. The section was modeled with five layers, and the interactions between adjacent layers make them work together. The concrete core, mortar, and sounding soil were modeled using an eight-node reduced-integration “brick” element (C3D8R), while a four-node reduced-integration shell element (S4R) was used for the steel cylinder. To account for the stiffness contribution of the steel wire, a two-node truss element (T3D2) was employed for the prestressed steel wire, and the equivalent cooling method was used to simulate the prestress [
26]. The stress-strain relationships of PCCP are derived from the current AWWA C304 standards [
27], and the strain damage threshold is 1155με for mortar. The material parameters of the PCCP and the surrounding soil are shown in
Table 1 and
Table 2. The pipe-pipe interaction was modeled using the contact mechanics method, which employs a contact algorithm to describe interface friction. A fully bonded interface was assumed for the pipe-soil interaction, achieved by defining a tie constraint. [
7]. In addition, the other finite element implementation adheres to the process outlined in the published paper [
28]. The whole model comprises about 240,000 elements and 270,000 nodes.
Twelve steps were performed to simulate the state of the pipeline during operation (gravity of the backfills, pipe, and fluid, operating internal pressure, and the effect of groundwater), and then the soil void was introduced. Three types of boundary conditions were implemented as follows: the bottom of the model was fixed in three directions, the top surface was left unconstrained to act as a free surface, and the remaining four surfaces were restrained against displacement in their normal directions.
3.2. Pipe Stress Analysis
The local void in the bedding results in uneven support, altering the interaction between the pipe and soil and consequently affecting pipe stress.
Figure 8 shows the influence of the void on the maximum principal stress of the mortar layer at various characteristic positions along the lateral longitudinal path. Under operational loads, the overall PCCP deformation is “elliptical.” The maximum principal stress at the pipe waist is markedly higher than that at the crown and invert, and the value is about 3.5 MPa, 1.0 MPa, and 0.2 MPa, respectively, with minimal variation along the longitudinal path. The introduction of a local void (20-100-125) causes the redistribution of pipe stress and has an obvious effect on the pipe invert. Its influence gradually decreases along the longitudinal path, while the maximum principal stress increases by about 1.1 MPa. In other words, the increase can lead to the tensile state of the mortar and concrete, which may result in damage to the pipe material.
As seen in
Figure 8, the most vulnerable location in a local void is the pipe invert. In order to evaluate the impact of void evolution on pipe stress, the maximum principal stress is extracted for varying void sizes, as illustrated in
Figure 9. The presence of a local void in the bedding results in an increase in the maximum principal stress at the pipe end (above the void), with the increment gradually decreasing along the longitudinal direction away from the void. The size of the void has a significant effect on the maximum principal stress, with void length exerting a greater impact compared to the depth and angle. A longer void length results in a larger longitudinal response range. When the VL reaches 2.5 m, the maximum principal stress increases by approximately 8 times compared to the case with no void.
Figure 10 shows directional stresses in the mortar layer of the outer middle section, comparing conditions without and with a soil void (void size: 20-100-125). The axial stress is uniformly distributed along the circumferential path and less affected by the void. The circumferential stress is basically equal to the maximum principal stress; that is, the direction of maximum principal stress is circumferential. Under the void condition, the maximum principal stress at the pipe invert increases by approximately 0.78 MPa, resulting in tensile stress. Therefore, it is essential to further investigate the material damage to the mortar layer. Additionally, the presence of a void significantly increases the circumferential stress, which is consistent with the field test results on concrete pipes [
15].
Differing from the mortar, the maximum principal stress of the concrete core is basically equal to the axial stress, indicating that the direction of the maximum principal stress is axial and evenly distributed along the path.
Figure 11 and
Figure 12 illustrate the impact of a void on concrete stress at the spigot and bell ends, respectively. The concrete core experiences circumferential compression due to prestressing. The compressive stress resulting from prestressing is greater than the circumferential stress under the operation loads, maintaining a compressive state in circumferential stress.
The void primarily affects the circumferential stress of the concrete core, reducing the circumferential compressive stress by approximately 0.9 MPa at the outer side of the pipe invert (see
Figure 11a) and by 0.6 MPa at the spigot and bell ends (see
Figure 12a). The circumferential distribution on the inner is opposite to that on the outer path. The circumferential compressive stress at the invert of the inner concrete increases because of the void, with rises of about 0.7 MPa at the spigot (see
Figure 11b) and 0.8 MPa at the bell end (see
Figure 12b). Notably, the maximum principal stress of the inner spigot end (
Figure 11b) and the outer bell end (
Figure 11a) is tensile. The value of the inner spigot end is about 2.4 MPa, but it does not reach the tensile strength of the concrete. In conclusion, the soil void affects the circumferential stress of the concrete core, which remains compressive. Namely, the void does not cause crack damage to the concrete but rather alters its stress distribution. The potential damage to the material is mainly the crack of the mortar protective layer.
The soil void spanning the joint causes a change in the stress state of the pipe, including a redistribution of the bending moment in the concrete core. The formation of local voids leads to a decrease in the bending moment at critical positions, which is consistent with the findings in the published literature (see
Figure 13 below).
Figure 14 further compares the influence of void size (length, depth, and angle) on the bending moment along the circumferential path. The greatest impact is observed at the pipe invert, where the maximum reduction is approximately 70%. Furthermore, the void depth and angle have a lesser influence on the bending moment than the void length. As the length of the void increases, the bending moment of the concrete gradually increases at the pipe invert.
3.3. Joint Angle and Material Damage
The vertical displacement of the pipeline is about 1.3 cm under operational loads. During service, adverse operational conditions may create localized ground voids, causing unsupported soil and subsequent pipe displacement.
Figure 15 shows the vertical displacement of pipe with different soil void sizes. It is evident that the void causes vertical displacement, and the closer to the void, the greater the displacement. When the void is 20-100-125, the maximum vertical displacement of joint reaches about 1.33 cm. With the continuous expansion, the vertical displacement increases, and the displacement increases to about 1.38 cm when the void evolves to 90-100-125. Consequently, the local soil void in the bedding can cause vertical displacement of the pipe, while the displacement due to void is quite small and can be ignored for large-diameter PCCP.
The joint is the relatively weak position of the segmented pipeline. A local soil void spanning the joint can cause the vertical displacement of the pipe on longitudinal path, resulting in the relative angle of the pipe joint.
Figure 16 illustrates the effect of void size on the relative angle. As void evolution (increasing in angle, length, and depth), the relative angle of the pipe joint tends to increase continuously. However, the maximum angle observed is approximately 0.02°, which maintains a significant safety margin relative to the permissible angle limit of 0.5° [
29]. Thus, a small local void does not significantly alter the joint angle, indicating a low risk of leakage due to the localized void.
PCCP is a composite material pipe, and the cracks on the mortar layer and concrete are the initial signs of pipe damage. The maximum principal strain of the mortar with (the void size is 20-50-250) and without (operation condition) void is shown in
Figure 17. The PCCP exhibits an “elliptic” deformation under operational loads, with the maximum principal strain of the mortar layer measuring approximately 160 με at the pipe springline and about 10 με at the pipe invert. The local void changes the stress state of the pipe, and the maximum principal strain slightly increases to about 165 με and 60 με at the waist and invert, respectively. The values remain significantly below the specified limit for mortar crack damage [
27].
To sum up, due to the design safety factor and the large rigidity of the large-diameter PCCP, the local void (with a maximum size of 90-100-250) has little influence on the vertical displacement and relative angle for large PCCP. Consequently, the risk of leakage from excessive angle changes or significant material cracking induced by local voids is low. This further indicates that large PCCP has good structural safety and can withstand local ground inhomogeneity in a certain range.
3.4. Influence of Local Void on Soil Stability
To analyze the stability state of soil around the void, the Molar-Coulomb nonlinear constitutive relationship model is employed, and its shear failure criterion can be expressed as follows:
where
is shearing stress,
is the critical shear strength of soil is determined by cohesion
and internal friction angle
.
According to the relationship between normal stress and shear stress, the strength theory expressed by the maximum principal stress can be derived as follows:
Then, a constant
for shear failure is defined in post-processing,
, that is, the soil is shear failure when
.
Figure 16 shows the nephogram of bedding soil under three cases: operational load (no void), void formation (20-50-62), and expansion (20-100-125). Under operational loads (
Figure 18a), the bedding soil is evenly stressed with no shear failure occurring. When the soil void appears, the bedding can no longer provide uniform support, resulting in a change in stress state and approaching the shear critical value at the pipe waist (
Figure 18b). The soil at the pipe waist and on both sides of the void has shear failure when the void continues to enlarge in size, as shown in
Figure 18c. Therefore, the soil state around the void is influenced by its size, with larger voids posing a greater risk of collapse. Additionally, loose soil is prone to loss, which can exacerbate void expansion, weaken the interaction between the pipe and soil, and jeopardize pipeline safety.
In summary, the large stiffness and high safety factor of a PCCP ensure that a small local void in the bedding does not typically lead to severe material damage to the pipe itself. However, it does significantly alter the stress distribution within both the pipe and the surrounding soil. As voids evolve, there is an increased risk of the mortar protective layer cracking. Moreover, the support loss due to void and the secondary damage from soil collapse are the main potential threats to the safe operation of the pipeline. Therefore, timely measures, such as grouting, are essential to prevent more serious consequences.
4. Conclusions and Recommendations
Ensuring stable pipe-soil interaction is crucial for the safe operation of buried pipelines. Soil voids can lead to uneven support, which is an unavailable condition for buried segmented PCCPs. This paper derives expressions for the shear force and rotation angle of a joint under a local void using a beam-on-elastic-spring approximation model. It then qualitatively assesses the influence of parameters such as soil stiffness, pipe length, pipe diameter, and void length. Based on the results, a nonlinear three-dimensional numerical model was developed to investigate the effects of bedding voids on the stress state and material damage to large buried PCCPs in service. The evaluation facilitates the prevention of functional failures and contributes to the longevity of the pipeline. The main conclusions are as follows:
(1) For rigid pipelines, a unilateral local void induces shear force and rotation angle, and the shear force is related to void length, pipe length, and overburden load. Symmetric bilateral void results in zero shear force but significantly increases the rotation angle, which is positively correlated with overburden load and negatively correlated with the pipe diameter.
(2) The local void alters the pipe-soil interaction, primarily impacting circumferential stress while having a lesser effect on axial stresses in the mortar layer of PCCP. Compared to void depth and angle, void elongation in the longitudinal direction has a more pronounced effect on pipe stress. When the void length reaches 2.5 m, the maximum principal stress increases about eight times compared with the no-void condition. Additionally, the local void leads to a reduction in the bending moment on the concrete, with a maximum decrease of approximately 70% at the pipe invert. As the void evolves, shear failure of the surrounding soil occurs, leading to local displacement nearby and potential collapse towards the void, leaving the pipe waist unsupported. Furthermore, the loose backfilled soil at the pipe invert may contribute to the formation of a larger-scale void.
(3) For large-diameter PCCP, a localized soil void in the bedding is unlikely to cause severe damage to the pipe material or increase the risk of joint leakage due to the high stiffness and safety factors incorporated in the design. However, it significantly alters the stress conditions between the pipe and the soil. As the void continues to evolve, collapse of the surrounding soil may occur, leading to a loss of foundation-bearing capacity and subsequent damages that pose a primary potential threat to pipeline safety. Therefore, timely grouting repairs at the locations of voids are essential to mitigate these potential hazards.