Abstract
Due to the increasing number of vehicles in urban areas and the degradation of pavement performance, the failure rate of polyethylene pipeline, across roads, has increased rapidly. Based on relevant data, the traffic loads are regarded as a random variable; also its distribution is determined. First, the nonlinear contact interaction model was established, among traffic loads, soil and buried polyethylene pipeline, using finite element analysis software. Second, considering the decline of the pavement performance, the mechanical characteristics of buried gas pipeline suffering from different traffic loads were analyzed. Based on statistics of vehicle load, the distribution law of its maximum equivalent stress was determined by numerical simulation. Particularly, the strength distribution law of polyethylene pipeline was obtained through experiments. The reliability index of polyethylene pipeline was calculated using central point method and Monte Carlo method. Finally, based on the degradation curve of pavement performance, the remaining life of buried polyethylene pipeline was assessed.
Introduction
As an important transportation method, the pipeline is widely used in lifeline projects, such as water supply, drainage, and gas delivery. Due to the convenient installation and maintenance, polyethylene (PE) pipeline is widely used in the transportation of urban gas. As we know, pipeline transportation uses buried form generally. There are two types of pipeline loading forms. One is internal loads (internal pressure) and the other is external load divided into natural loads and artificial loads. Buried pipelines were attacked by natural loads frequently, such as ground settlement, landslide, and debris flow.1–3 The common artificial loads include the repeated load (traffic load), 4 static load (surface load), 5 impact load, 6 and so on. Apparently, the buried pipeline in cities or towns is mainly affected by artificial loads. 7 Especially, in recent years, with the increasing number of vehicles in urban areas and the degradation of pavement performance, the rate of failure accidents has increased rapidly. Therefore, the influence of traffic load on buried pipelines gradually becomes the research focus. The relative research on the loads, applied to pipes, has been developing since the Marston theory was put forward. At present, most researches on traffic load are focused on the response of buried pipelines or the failure of pipelines. These studies involve different types of pipes, such as concrete pipes,8–10 steel pipes, 11 or flexible pipes.12,13 In addition, these studies can also be roughly divided into experimental studies,14–17 numerical simulation study,18,19 and their combination. 20
However, the traditional design of the pipeline employs the method based on strength design criterion, such as the allowable stress. Due to the improvement of the public’s awareness of risk, more frequent natural disasters as well as the improvement of economic competitiveness, in order to meet social expectations of pipeline safety, the design and evaluation method based on the reliability must be developed. 21 At present, the design and evaluation method based on the reliability has received a great deal of attention even if the acquisition of target reliability and relevant variables exist many problems. Many scholars have studied the reliability of pipelines, such as the study on the reliability of corroded metal pipelines,22–24 PE pipelines with pressure 25 as well as buried pipelines under the external load of landslide or longitudinal ground movement.26,27
Since the complexity in studying special material characteristics and the mechanical response of buried pipe suffering from external load, the reliability analysis of buried PE pipeline under traffic loads has not been reported. Moreover, the probabilistic assessment of an engineering system performance might contain a substantial number of uncertainties in system behavior due to its functioning mode, 25 which makes the reliability analysis more difficult. For buried PE pipeline under traffic loads, this uncertainty is affected by many factors, such as vehicle axle speed, coupling between vehicle, and underground pavement structure. Partial difficulties are also reflected as follows: (1) relationship assessment among the random variables of the traffic load–soil–pipeline system, (2) finding methods (or relationships) to reduce number of involved random variables, and (3) difficulty in collecting relevant data. 28
The aim of this work is to conduct the reliability analysis for a gas PE pipeline subject to traffic loads. Actually, once the pipeline is buried, the layered pavement structure and the form of pipe–soil contact are determined. Thus, the reliability of the pipeline is only related to external loads.29,30 Based on the stress–strength interference model and the idea of reducing random variables, the reliability calculation of pipeline is simplified to the problem for finding the maximum stress distribution of pipeline and the strength distribution of material. Using tensile test, the characteristics of PE materials were analyzed and the probability distribution of yield strength was calculated. Numerical simulations using finite element method are performed to determine pipe–soil system response under traffic loads and obtain the relationship between traffic load and pipeline stress. Based on the statistical data of traffic loads, the probability distribution rule of pipeline stress was obtained. Based on the von Mises yield criterion, the ultimate state equation was established. Then the reliability index of pipeline was calculated using the center point method and Monte Carlo method. Finally, a graphical method of remaining life prediction based on the degradation of pavement performance is proposed to estimate the remaining life of buried PE pipeline subject to traffic loads.
Uniaxial tensile test
The real stress–strain curves of most PE pipeline are shown in Figure 1. The curve is divided into (a), (b), and (c) stages: 1 (a) the stage shows an elastic stage similar to steel; (b) this stage is the plastic softening stage or necking stage, at which the stress decreases gradually with the increase of strain; (c) this stage is the rising stage. At the end of stable necking, the stress will increase for a short time, and then, the specimen will snap. This article focuses on the change of the pipeline in the elastic stage only. Once the pipe yields, irreversible plastic deformation will occur and the pipeline will be damaged.

Stress–strain relationship of polyethylene material.
For the sake of uniaxial tensile test of PE pipeline materials, the sample was produced according to ISO 6259-3:2015. 31 The tensile rate has been set as 1, 3, 10, 30, 100, and 300 mm/min, respectively. The test apparatus is shown in Figure 2. The test data are shown in Figure 3. This article only cares the data near yield. The yield strength of PE pipeline at different tensile rates (v) was calculated, as shown in Figure 4. It can be seen from the figure that the yield strength of PE pipeline has obvious rate correlation. The greater the tensile rate or strain rate, the greater the yield strength. When the strain reaches a certain degree, the greater the tensile rate, the faster the stress drops in the contraction stage. In conclusion, the relationship between yield strength and tensile rate was obtained.

Uniaxial tensile test apparatus.

Engineering stress–strain relationship of uniaxial tensile tests at different velocities.

Yield stress at different velocities of tensile.
Polyethylene material is a special kind of material. Its material characteristics vary with temperature and time. Suleiman and Coree 32 established a hyperbola constitutive model considering the strain rate by a more systematic method, and the curves obtained are in good agreement with the curves in Figure 3. Since this article only considers its stage before yield, this model is chosen as the constitutive model of PE pipeline materials. The general formula of hyperbola constitutive model is as follows
where

The relationship between constants a, b and the velocity of tensile.
The initial elastic modulus of PE pipeline is inversely proportional to its constant a. The relationship between them is shown as follows
Calculation of maximum pipeline stress under traffic load
Conditional hypothesis
Pipe fittings and joints of pipeline were not considered in this work.
Drucker–Prager yield criterion was used to deal with the plastic deformation of soils.
The interaction behavior between soil and pipe was a limited sliding contact.
The material properties of PE were assumed to not change over temperature and time. In other words, the creep and relaxation properties of PE would not be considered.
Material parameters
The PE80 grade pipeline was used in this work, the standard size is standard size ratio (SDR), D/t = 17.6), and the pipe diameter was 400 mm. Drucker–Prager model was adopted for soil and rigid basement soil, which can simulate the elastic–plastic state of soil under impact loads well and avoid instability in the calculation of large deformation. 21 Because the purpose of this article is not the mechanical response of pavement structure, the cement stabilized macadam layer and asphalt concrete surface layer are regarded as homogeneous linear elastic materials. The material parameters of each layer of soil are shown in Tables 1–3.33,34 The density of PE pipeline was 950 kg/m3 and its Poisson’s ratio was 0.4.
Material parameters.
Soil modeling parameters (D-P model).
Hardening parameters of Drucker–Prager soil.
The failure criterion is the basis for judging the failure, which depends on the failure mode. In this article, the criterion of the von Mises yield criterion was adopted, that is, when the equivalent stress of pipeline reaches the yield stress of pipe material, the pipeline failure is considered. The equivalent stress is von Mises equivalent stress, 35 and its expression is
where
Model establishment
The whole road model was simplified into four layers, including rigid basement, filling soil layer, cement stabilized macadam layer, and asphalt concrete layer. The pipeline was away from the pipeline side boundary more than six times the pipe diameter. The length and depth of the foundation of the model were set large enough to minimize the impact of boundary conditions. The calculation model is shown in Figure 6.

Model sketch.
The entire model was a cuboid with a scale of 8 m × 6.5 m × 3.5 m, which simulates a common single lane with a width of 3.5 m. 33 Due to the symmetry of loads, semi-model was adopted for modeling. The bottom of the soil remained fixed in all direction (Ux = Uy = Uz = 0). The upper surface of the soil was free surface without any constraints. Symmetrical boundary was imposed on the other surfaces of the soil and the both ends of the pipeline, as shown in Figure 7. At the moment, the pipeline–soil contact model is the most accurate method for evaluating the interaction of the pipeline–soil at the large deformation segment and can explain the friction and separation of the pipeline–soil interface. The contact between the outer surface of the PE pipeline and the surrounding soil was set as a limited-slip contact. A contact algorithm was employed to simulate the interface between the outer surface of the pipeline and the surrounding soil. Interface friction had been taken into account in the contact algorithm in the form of defining an appropriate friction coefficient. The friction coefficient equals to 0.18. 26 Eight-node and hexahedron linear shrinkage integral unit (C3D8R) were employed for modeling both pipeline and surrounding soil. The mesh of the soil around the pipeline was a local mesh refinement, as shown in Figure 7.

Finite element model of (a) the refined mesh around pipeline, (b) the boundary conditions, and (c) the meshing of pipeline.
Numerical calculation
Relevant regulations stipulate that the buried depth of PE pipeline under the road should not be less than 0.9 m.36,37 When the buried depth is set as 1 m and the change range of vehicle loads (T) is 0.4–1.8 MPa, the mechanical characteristics of PE pipeline under traffic loads are analyzed; also the distribution law of its maximum equivalent stress is determined.
It can be seen from Figure 8(a) and (b) that the road is obviously deformed under the action of vehicle. By comparing to all the results, corresponding stress concentrations appear in the middle of the model, and meanwhile, the deeper the dent is, the greater the equivalent stress value it will be as the loads increase. As shown in Figure 8(c), a large deformation and stress concentration emerge in the middle of the PE pipeline, and the maximum equivalent stress appears at the lower side of the inner surface of the middle pipeline. By comparing all the calculated results, it can be found that as the action increases, the deformation of pipeline grows and the stress values strengthen and concentrated.

Stress nephogram of model and pipeline (traffic load = 1.8 MPa, E1 = 120 MPa, E2 = 100 MPa): (a) whole model, (b) soil and pavement, and (c) pipe.
It is notable that the complex change of pavement performance is simplified into the reduction of the elastic modulus of road material in this article. As for obtaining the relationship between the stress or strain of PE pipeline and traffic loads, the calculated results are plotted into Figure 9, of which E1 and E2 represent the elastic modulus of asphalt concrete layer and cement stabilized macadam layer, respectively. Figure 9 shows that the maximum equivalent stress of the pipeline increases along with the raise of the vehicle loads and the aging of the pavement materials. However, due to the pipeline internal pressure, the maximum equivalent stress will not be lower than a certain value.

Stress versus vehicle load.
The curve of buried pipe stress and load under traffic load was fitted to obtain the mathematical expression of the relationship between traffic load and stress. The curve was fitted with the following function, and the relevant fitting parameters were listed in Table 4
where ST is pipe stress (MPa); T are traffic loads (MPa); a, b, and c are coefficients.
Fitting of the relationship between pipe stress and traffic loads.
Calculation of pipeline reliability
Distribution law of vehicle loads
In terms of probability distribution characteristics of vehicles, a lot of studies on vehicle loads have been made domestically, including normal, exponential, and lognormal distribution. The multi-peak model can be simplified into a combination of normal distribution for different types of vehicle loads. 38 The characteristic parameters of traffic loads distribution, which obtained by Wang based on actual measured data, are used in this article. 39 Specifically, China’s freight cars at present mainly use tires of 9.00-20 (about 15%), 10.00-20 (about 25%), 11.00-20 (about 40%), and 12.00-20 (about 15%). These four types of tires are selected as vehicles; their specific statistical data are shown in Table 5, and the statistical data of wheel pressure of large- and medium-sized buses are shown in Table 6.
Statistical table of freight car load.
Statistical table of bus load.
The data in Tables 5 and 6 are calculated by the numerical calculation tool MATLAB. The t-test method is used to determine whether the tire pressure is normally distributed. If it is normally distributed, its characteristic parameters are given, as shown in Table 7.
Characteristic parameters of tire pressure distribution.
The judgment parameter of normal distribution is h. When h = 0, it means that the original hypothesis cannot be rejected at the set significant level—the sample is subject to normal distribution. When h = 1, it means that the original hypothesis can be rejected at the set significance level—the sample is not subject to normal distribution. s is the probability when the hypothesis is true in the test. When it is small, the original hypothesis could be problematic. Therefore, the front wheel of 11.00-20 is subject to N∼ (1.10, 0.112), the rear wheel of 11.00-20 is subject to N∼ (1.29, 0.032), a large bus is subject to N∼ (0.94, 0.122), and a medium bus is subject to the normal distribution of N∼ (0.82, 0.12).
Yield strength distribution
Considering the rate correlation of PE material, when at a temperature of 23°, stretch rate is 50 mm/min and sample size is 214; pipeline tensile tests are conducted according to ISO 6259-3:2015; also the yield strength of PE80 pipeline material is obtained by statistics, as shown in Figure 10. According to the final experimental data, its yield strength follows the normal distribution of N∼ (21.13, 2.172).

Histogram of yield stress frequency distribution of PE pipeline materials.
Reliability calculation
The reliability of PE pipeline can be defined as the probability that the pipeline meets demands for gas supply under specified environmental conditions. The stress–strength interference model insists that if the stress interferes with the allowable strength of the system, the system might fail. 40 As for fracture failure, the system reliability R is the probability that the strength is greater than the stress, which can be expressed as
In the formula, P stands for the failure probability; S and R represent stress and strength distribution, respectively, and all of them are random variables with certain statistical distribution laws. The central point method in the first-order second-moment method (FOSM) was used to calculate the reliability index to establish the function of buried PE pipeline structure under traffic loads
where
Using the linearization rule of random variables, the function was expanded to Taylor series at the mean point
The linearized limit equation of state was as follows
Final reliability index was as follows
where
The reliability index
PE: polyethylene.
Monte Carlo sampling number = 10,000.

The reliability index
It can be seen from Table 8 that under the different actions of vehicles, the reliability of PE pipeline also performs differently. The larger the vehicle load, the smaller is the value of reliability index. The reliability index of PE pipeline is about 3.15 under the action of 11.00-20 front wheel, 2.68 under the action of 11.00-20 rear wheel, 3.65 under the action of large bus, and 4.15 under the action of medium bus. With the decline of pavement performance, the reliability index of the pipeline shows a sharply downward trend in Figure 11. In other words, the reliability decreases accordingly.
It is assumed that the degradation curve of road performance is shown in the green curve in Figure 12. The life of PE pipeline can be obtained by combining the degradation curve of road performance (E–n curve) and the changing curve of PE pipeline reliability with road performance (r–E curve). For example, when the target reliability index of PE pipeline is rlimit, the pavement performance Elimt can be obtained from the r–E curve. Then, the limit service life nlimit can be obtained from the E–n curve. In addition, when service time n of PE pipeline is known, the corresponding residual life of PE pipeline is nlimit– n years.

Life prediction diagram of PE pipeline.
Conclusion
In this article, pipeline mechanical model is established and a simple analysis carried out. Stress and stress distribution of pipeline are calculated by single random variables of traffic load. According to von Mises yield criterion, the limit state equation is set up. The reliability index of PE pipeline which under the action of traffic was calculated using central point method and Monte Carlo method.
After the actual tensile test, the yield strength of PE pipeline material at the tensile rate of 50 mm/min is generally in accordance with the normal distribution of N∼ (21.13, 2.172).
The maximum equivalent stress of the PE pipeline increases along with the raise of traffic loads and the aging of pavement materials. There is an exponential relationship between stress and traffic loads.
The limit state equation of buried PE pipeline under traffic load is established, and the reliability index of pipeline is calculated using the central point method and Monte Carlo method. The larger the vehicle load, the smaller is the value of reliability index. With the decline of pavement performance, the reliability index of the pipeline shows a sharply downward trend.
Based on the studies about the aging performance of pavement materials, the remaining service life of the PE pipeline, under a certain target reliability index, can be calculated, which provides a new prospect for the life assessment of PE pipeline crossing the road.
The stress environment of the buried PE pipeline under traffic load is simplified in the current work, considering the more independent random variables. Future work should focus on forming an entire pipeline life assessment method through a comprehensive analysis of these factors: the aging of actual pavement materials, the number of vehicles, and the reliability of pipelines.
Footnotes
Appendix 1
Acknowledgements
This work was supported by the China Scholarship Council and Science & Technology Program of Sichuan Province, China (No.2019YFS0075).
Handling Editor: Vesna Spasojevic Brkic
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received the following financial support for the research, authorship, and/or publication of this article: This project was supported by the China Scholarship Council and Science & Technology Program of Sichuan Province, China (No.2019YFS0075).
