Stress Analysis of Buried Pipeline Installed by Horizontal Directional Drilling Using ANSYS Finite Element Software

In designing a horizontal directional drilling (HDD) pipeline project, designers face the challenge of determining the regions of maximum and minimum stresses on pipelines, ensuring the stability of the bore-hole from collapse and minimizing the stresses induced on the pipeline due to the bore-profile. This study analyses the stress induced on an HDD pipeline system using the ANSYS Version 18, mechanical APDL finite element (FE) software. The pipeline used as the case study was a gas transmission pipeline installed in south-west Nigeria. A macro-file for ANSYS Version 18, mechanical APDL used to model the pipeline was developed. The results showed that the maximum and minimum stresses induced on the HDD pipeline were at the top and bottom of the pipe, respectively; while the stresses on the sides were uniform (≈888 kg/cm) all through the pipeline, irrespective of element number. The maximum stress occurred at the curvature point with the highest entry angle (10 ̊), resulting in a maximum deflection at this point. The model stress validation performed by comparing results with theoretical solutions, both with respect to radius of curvature and internal pressure, showed percentage difference (errors) less than 10%. The cross sectional area validation showed a percentage difference of 0.059%.


Stresses on Buried Pipeline
The loads acting on the pipe include its weight, the internal operating pressure and the external pressure due to the overhead soil load [1]. Underground loads on a buried pipeline induced by the soil pressure can be divided into two parts.
Firstly, there is the soil pressure which increases with the depth of cover. Then, there is the additional force caused by ground loads which decreases with increase in depth of cover [2]. The depth of embedment of pipe, type of backfill, the pipe thickness and surcharge loads are the prime factors that affect the behavior of buried pipes [3].
HDD pipelines are subjected to a combination of tension, bending, and external pressure. Drilled path design and pipe specification may be governed by these installation loads, as they could be more severe than operational loads [4].
The borehole profile and pipe properties must be considered in selecting the suitable pipe materials for an HDD installation. These two factors should be considered together to ensure that the pipeline can be installed and operated without risk of damage [5].

Finite Element Method (FEM) for Pipeline Stress Analysis
In the analysis and design of underground structures, the finite element method (FEM) has been used extensively as an effective numerical technique. The results more reasonably reflect the actual situation, providing an improved design procedure to avoid pipe wear and breakage [1].
Its application ranges from the study of aircraft or automobile structural framework, complicated thermal system (nuclear power plant), to the analysis of a fluid flowing through a duct, over a weir, or through the earth [6] [7] [8] developed a 3D FE program using MATLAB to examine the performance of buried pipeline subjected to large ground motion; crossing active fault.
Spectrum technique with 1D beam units was used by Wu et al. [9] to undertake dynamic earthquake analysis on oil pipelines using CAESAR II stress analysis software. Hence, critical sections and maximum stress region were identified visibly. It was concluded that reducing the bend angle could minimize the stress and displacement of underground pipelines during an earthquake. However, this study was not extended to stress analysis of deeply buried HDD pipelines. Also, Huang, et al. [10] utilized CAESAR II software to investigate the stress on underground oil pipelines in steep and high slope regions. Comparing the different simulation results showed that the movement and stress of underground oil pipelines in steep and high slope regions were drastically reduced with enhanced safety by proposing an improved laying process.
Gabar and Bilgin [11] used ABAQUS FE software package to model buried pipeline behavior under surface-induced actions, and described the pipe-soil interactions and pipe deformation. They concluded that vertical pipe deformation increased with increasing loads for various soil layers, which may cause buckling.
Xiao, et al. [1] performed analysis of glass reinforced plastics (GRP) pipeline using ANSYS and the results were generally consistent with that of the traditional CECS190 standard, based on the maximum long term vertical deformation. It Nirmala and Rajkumar [3] also used ANSYS as FEA tool on buried UPVC pipe and suggested that the depth of backfill increases the deflection of buried pipes, and that the behaviour of buried pipes varied due to the type of backfill material. Hence, it can be observed that ANSYS has not been used as an FEA tool for HDD pipeline stress analysis.

Materials and Methods
The model was created using ANSYS, Version 18, Mechanical APDL following the flow diagram in Figure 1. The co-ordinates for the key points, exit angle, entry angle, element type and the geometric properties (such as, pipe outer diameter, thickness, length of crossing etc.); and the material properties of the steel pipe (such as, the specified minimum yield stress, poisson ratio and young modulus) were all inputted into the model at the pre-processing stage.
At the solution stage, the loads and boundary conditions such as the degree of freedom and model solution setting were clearly specified. The final stage of the model, post-processing stage, was where the deflection plots and the stress contour diagrams of the pipeline system were plotted.

Design Input Parameters
The pipeline used to create the model on ANSYS, Version 18, Mechanical APDL, was a gas pipeline installed by Horizontal Directional Drilling (HDD) in south-west Nigeria. The gas pipeline design parameters are shown in Table 1.
These parameters were input into the ANSYS platform to give the section outline in Figure 2. The ANSYS, Version 18, Mechanical APDL solver was used to model the pipeline. A macro file for non-linear analysis was prepared for modeling the pipeline system.

Modeling and Meshing
Four (4) co-ordinates points were identified on the pipeline design profile to  Figure 2. Section outline.
F. T. Olumoko, C. V. Ossia World Journal of Engineering and Technology serve as the key-points. The pipeline was assumed to move along the x-y direction and therefore no co-ordinate was given for the z-direction. The radius of curvature for the sag-bends was modeled as fillet between lines 1, 2 and lines 2, 3 respectively. For this model, the radius of curvature for both sag-bends was expected to be same, with reference to the design profile. However, the exit and entry angles were different. The number of pipe section mesh was defined for the model. The "SECNUM" command was used to identify the particular section for meshing. This command reads the section division as the amount of mesh sizes. Figure 3 shows the pipe section divisions used for the pipe meshing.
The "LMESH" command was used to generate nodes and line elements along lines. For this model, meshing was done on all the lines by selecting "ALL".

Loads and Boundary Conditions
The "ANTYPE" command was used to specify the type of analysis and the restart status. For this model, the static analysis was required. The "DK" command was utilized to state the degree of freedom at the key points. For this model, at all the key points, only one degree of freedom was permitted which was in the x-direction. The y-and z-directions were constrained. The "FK" command was used to define the applied load in the x-direction at the key points. This applied load is the force of the rig which pulls the pipe in the x-direction. The "ALLSEL" command was used to select all entities.

Theoretical Equivalent Resultant Stress
After running the ANSYS model using the relevant input parameters, the maximum equivalent stress obtained was S ANSYS . The percentage error when comparing the resultant equivalent stress values from both theoretical calculation and ANSYS model was evaluated using Equation (5).
For the ANSYS model validation and investigate the effect of the HDD radius of curvature on the pipeline stresses, the internal pressure was kept constant at P = 127.5 kg/cm 2 and different R-values (90,000 cm ≤ R ≤ 120,000 cm) were applied to obtain different S theory , S ANSYS and percentage (%) difference which results were shown in Figure 5. Similarly, to investigate the effect of internal pressure on the HDD pipeline stresses, the pipe radius of curvature was kept constant at R = 100,000 cm while the internal pressure applied was varied in the range 90 kg/cm 2 ≤ P ≤ 127.5 kg/cm 2 to obtain different S theory , S ANSYS and percentage (%) difference and the results shown in Figure 6. Figure 5 showed the radius of curvature and the corresponding maximum equivalent stress obtained from theoretical method and ANSYS, respectively. From the Figure 5 it can be observed that the percentage difference between the theoretically calculated maximum resultant stress and the maximum stress obtained from the ANSYS model was less than 10% with a maximum difference of 3.49% at a radius of curvature of 90,000 cm (900 m) and a minimum % difference of 0.22% at a radius of curvature of 110,000 cm (1100 m). Furthermore, the percentage error followed a parabolic distribution with respect to the radius of curvature, with a minima at R = 110,000 cm as shown in Figure 5. Figure 6 showed the variation of Von-Mises (maximum equivalent stress) with respect to internal fluid pressure of the pipeline. The results in Figure 6 showed a comparison between the theoretically calculated and ANSYS equivalent resultant stress. From Figure 6, it can be observed that the percentage difference between the theory and ANSYS results was less than 10% with a maximum difference of 9.81% at an internal pressure of 90 kg/cm 2 and a minimum difference of 1.44% at an internal pressure of 127.5 kg/cm 2 .

Cross Sectional Area Validation
A comparison was also done between the theoretically calculated cross-section area of steel pipe and the ANSYS model cross-sectional area. The theoretically calculated cross-sectional area of the pipe is given as Equation (6).
where e A and i A are the external and internal cross-sectional area of the pipe, respectively; and A s was calculated as A s = 579.12 cm 2 . World Journal of Engineering and Technology From Figure 3, the cross-sectional area of the ANSYS model cross sectional area was obtained as 578.439 cm 2 . Therefore, the percentage difference in the cross-sectional area was obtained as: 579.12 578.439 Percentage Difference *100% 0.059% 579.12 578.439

Element Numbering
The pipeline was discretized into 24 elements with 3 node points each by the ANSYS program for FEA to be done on each element. However, the element numbering was not arranged sequentially by APDL, but as shown in Figure 7.

Stress Contour Diagram and Maximum Stress Distribution
The stress contour diagram for the equivalent stress or Von Mises stress was plotted to investigate the stress distribution on the pipe.    it can be seen that this was not the case. Figure 10 showed the minimum stress distribution around the pipe for each element which also indicates that the stress distribution along the sides of the pipe all through the pipeline was relatively uniform regardless of the location.

Minimum Stress Distribution
However the top and bottom minimum stress varies considerably at the curved sections. The minimum stress at the curved section indicates a lower stress level on the elements with maximum stresses.

Deflection Plot
A nodal solution of the pipe deflection was plotted to show the deflection at each section of the pipe. From Figure 11 it can be observed that the maximum and minimum deflection experienced by the pipe was 413.942 cm and 0.000486 cm, respectively which were displayed on the top left-hand side of the Figure 11. The exact point of maximum deflection can also be observed in Figure 11. The location point of the maximum and minimum deflection point for the pipeline system was shown in Figure 12.
The plot of the pipeline deflection was shown on Figure 11. This plot in Figure 11    indicates that the maximum deflection on the pipeline was experienced at the point of maximum stress level. While, Figure 12 showing the points of the maximum and minimum deflection showed that the minimum deflection was obtained at the entry point of the HDD installation.

Conclusion
Hence, ANSYS can be used as veritable tool for the analyses and design of a pipeline system to be installed by horizontal directional drilling (HDD). The results obtained from this study using the ANSYS model were corroborated those of theoretical methods. The point of curvature with the higher angle (that is, if both points of curvature are the same and the entry and exit angles are different) experiences the maximum induced stress. Maximum deflection of the pipeline should be expected at the point with maximum stress as this was due to the high bending moment experienced at this point. On a particular pipeline system, the maximum and minimum stresses were induced at the pipe top-side and bottom-side, respectively; while the stresses at the sides of the pipeline remained relatively uniform with slight variation.