Comparison of 3D Solid and Beam–Spring FE Modeling Approaches in the Evaluation of Buried Pipeline Behavior at a Strike-Slip Fault Crossing

: Validated 3D solid ﬁnite element (FE) models offer an accurate performance of buried pipelines at earthquake faults. However, it is common to use a beam–spring model for the design of buried pipelines, and all the design guidelines are ﬁtted to this modeling approach. Therefore, this study has focused on (1) the improvement of modeling techniques in the beam–spring FE modeling approach for the reproduction of the realistic performance of buried pipelines, and (2) the determination of an appropriate damage criterion for buried pipelines in beam–spring FE models. For this paper, after the veriﬁcation of FE models by pull-out and lateral sliding tests, buried pipeline performance was evaluated at a strike-slip fault crossing using nonlinear beam–spring FE models and nonlinear 3D solid FE models. Material nonlinearity, contact nonlinearity, and geometrical nonlinearity effects were considered in all analyses. Based on the results, pressure and shear forces caused by fault movement and pipeline deformation around the high curvature zone cause local conﬁnement of the soil, and soil stiffness around the high curvature zone locally increases. This increases the soil–pipe interaction forces on pipelines in high curvature zones. The beam–spring models and design guidelines, because of the uniform assumption of the soil spring stiffness along the pipe, do not consider this phenomenon. Therefore, to prevent the underestimation of forces on the pipeline, it is recommended to consider local increases in soil spring stiffness around the high curvature zone in beam–spring models. Moreover, drastic increases in the strain responses of the pipeline in the beam–spring model is a good criterion for a damage evaluation of the pipeline.


Introduction
The pipeline network has been spread all over the world to provide the essential needs of human societies (e.g., for transmission of gas, water, oil, wastewater, and chemical products). Hence, there are a great deal of pipelines crossing seismic hazardous areas, such as high curvature zones [1]. Earthquakes are the greatest threat to structures [2]. In the case of buried pipelines, most damage arises due to permanent ground deformation (PDG), for example, fault dislocations, liquefaction, and landslides. However, small regions within the pipeline network and earthquake fault zones are prone to PGD; because they cause large deformation of the pipelines, the damage potential is very high. On the other hand, there are very few pipeline cases that are damaged by wave propagations [3,4]. Accordingly, this paper has focused on the stability analysis of buried pipelines subjected to the strike-slip fault movement during PGD. Earthquake fault movements cause plastic deformation and even rupture to buried pipelines. The damage of buried pipelines due to previous earthquakes has been reported and has caused severe health and environmental issues [5][6][7][8][9][10][11][12][13][14][15][16]. Damage to the pipeline system not only affects urban infrastructures' serviceability after earthquakes, but also, owing to the leakage of ecologically dangerous materials (e.g., chemicals, natural gas, fuel, or liquid waste), they can result in health and environmental issues [17]. Therefore, it is essential to develop analysis methods for buried pipelines to ensure the reliable resistance and economic design of these pipelines. Remarkably, because PGD can cause severe pipeline damage, the construction of buried pipelines at earthquake fault intersections is a key problem in engineering design tasks [1,18]. The effect of active fault crossings on buried pipelines has been investigated through several experimental studies [19,20]. Palmer et al. [21] described a large-scale testing facility at Cornell University and its working principles. O'Rourke and Bonneau [22] then performed large-scale tests to evaluate the effects of ground rupture on high-density polyethylene (HDPE) pipelines and the performance of steel gas pipelines distributed with 90 • elbows. Lin et al. [23] performed small-scale tests to analyze the performance of buried pipelines under strike-slip faults. For the first time, the centrifuge-based approach was proposed by O'Rourke et al. [20] to model the ground faulting effects on buried pipelines, and several centrifuge tests have been performed to investigate the response of buried HDPE pipeline subjected to faulting displacement [24][25][26][27].
In 1975, the first analytical research on buried pipelines at a fault crossing was conducted by Newmark and Hall [28]. They evaluated a simplified analytical pipeline model by assuming cable-like behavior for a pipeline in small displacement, which was later extended by Kennedy et al. [29,30]. In those studies, the bending stiffness of the buried pipeline at the high curvature zone was neglected, which caused an overestimation of the bending stress-strain and increased the axial forces. Wang and Yeh [31] developed Kennedy et al. [30] model to strike-slip fault crossings by the extension of a simplified pipeline bending stiffness in their calculation. Wang and Yeh included the lateral soil yielding and introduced four segments to the buried pipeline at a fault crossing by partitioning. Two segments were adjacent to the fault crossing (the "high curvature zone"), and two were farther from the high curvature zone. However, the soil yielding starting point and partitioning assumptions were not realistic. Karamitros et al. [32] developed a method for strike-slip faults, wherein the pipeline is partitioned into four segments that are analyzed based on the beam-on-elastic-foundation and elastic beam theories. After the analytical solution, the longitudinal soil-pipe interaction was determined in addition to the steel pipe material's non-linearity using a bilinear stress-strain relationship. The Karamitros et al. [32] model was developed by Trifonov and Cherniy [33] for normal fault crossings. They removed the symmetry conditions and estimated the axial elongation. Karamitros et al. [34], with the same assumption, extended their study [32] to normal-slip faults; however, this model has the same shortcomings as their previous study [32]. In 2012, an analytical model based on the stability model of Trifonov and Cherniy [33], including the operational loads (internal pressure and temperature gradient), was developed for the stress-strain analysis of buried pipelines at a fault crossing by Trifonov and Cherniy [35]; however, their study had the same shortcomings as their previous governing differential equation [33]. In 2020 and 2021, Talebi and Kiyono [36,37] introduced a novel nonlinear governing equation that includes the longitudinal sliding behavior of a pipe within soil during large PGDs, lateral elastoplastic soil-pipe interaction springs, and longitudinal forces made by geometrical nonlinearity effects. They removed the unrealistic assumptions and remarkably increased the accuracy and application area of the analytical methods for the problem of buried pipelines at active strike-slip fault crossings.
The efficacy of FEM-based analysis to assess the behavior of a buried pipeline crossing an active fault has been proved in the literature. FEM has been used to evaluate the buried pipeline performance with the assessment of criteria such as local buckling, ovalization, and tensile damages [18,[38][39][40][41]. Vazouras et al. [39] modeled a hybrid (shell and solid elements beside the equivalent springs) pipeline buried in solid soil. Liu et al. [40] modeled a buried pipeline at a reverse fault crossing using the FE commercial code ABAQUS, in which the pipe was modeled as shell elements and the soil-pipe interaction was modeled as non-linear soil springs. In their study, the pipe and the soil-pipe interaction were modeled as shell elements and as non-linear soil springs, respectively. Besides, they had an investigation on the buckling of buried pipelines influenced by the yield strength and strain hardening parameters. Demirci et al. [41] studied the behavior of a continuous buried pipeline subjected to a reverse fault motion through a new experimental centrifuge modeling of a pipeline crossing a reverse fault. Additionally, 3D FEM analyses were used for more detailed results. Literature reviews show that various modeling approaches, including beam, shell, hybrid I (beam+shell), and hybrid II (spring+shell) for a pipe and springs, and soil continuum solid elements for soil modeling have been employed to evaluate pipeline performance against fault movement.
A 3D FE (3D solid)-based analysis, the most detailed approach for modeling a pipeline at a high curvature zone, including the shell elements and solid elements to simulate the buried pipeline and surrounding soil, respectively, can provide the most realistic evaluation of buried pipeline performance, including the local buckling, ovalization, and tensile damages. Typically, the 3D solid model is used for research purposes, owing to the modeling complexity. The beam-spring model is a simplified FE modeling approach, which uses beam and spring elements to model the pipe and soil-pipe interactions, respectively. Although it is pretty common to use a beam-spring model for the design purpose of buried pipelines at a fault crossing and all the design guidelines are extended based on this modeling approach for the analysis of buried pipelines, there is a gap in the literature for a detailed comparison of the beam-spring model and the 3D solid model results (a 3D solid model has almost realistic results). This comparison provides helpful conclusions for the development of the beam-spring modeling approach and the improvement of buried pipeline design guidelines.
This paper intends to present a detailed comparison between a nonlinear beam-spring modeling approach and a nonlinear 3D solid modeling approach to evaluate the performance of buried pipelines at strike-slip fault crossings to improve the beam-spring modeling approach and pipe damage criteria in this approach. In this regard, firstly, to verify the validity of the longitudinal and lateral soil-pipe interaction of the 3D solid model, a pipe and soil box model of pull-out and lateral sliding tests, based on Vazouras et al. [39], was simulated. The 3D solid model calibrated based on Vazouras et al. [39] and highly verified the results reported in the literature. Force-displacement diagrams of the equivalent soil-pipe interaction springs in the longitudinal, lateral, and vertical directions were extracted by displacement control FE analyses. Secondly, FE-base beam-spring models for identical pull-out and lateral sliding tests in previous 3D-solid models were created, and soil-pipe interaction springs were created based on the 3D solid model results and the results of beam-spring models were verified against the 3D-solid models for pull-out and lateral sliding tests. Finally, a 3D solid and a beam-spring FE model with identical properties for the problem of a buried pipeline at a strike-slip fault crossing were created. The responses of the buried pipeline modeled by both modeling approaches were compared to evaluate the pipe performance in both modeling approaches, and to extract any requirements in the development of the beam-spring modeling approach, and the development of the damage parameters and criteria.

Analytical Evaluation of Buried Pipeline Behavior
Details of the strike-slip fault angle (ψ) and fault displacement (δ) are schematically shown in Figure 1, in which δ x and δ y are the longitudinal and lateral components of the fault displacement, respectively.
As mentioned in Section 1, Talebi and Kiyono [36,37] have extended the linear and nonlinear analytical methods for the stability analysis of buried pipelines at strike-slip faults with high accuracy results. Since, in this paper, the analytical method was just used for the parametric interpretation of the FE analyses results instead of the complex analytical approaches in the papers by [36,37], a simple linear analytical method mentioned by Talebi and Kiyono [42] was employed. A simplified differential equilibrium equation for the pipeline crossing the strikeslip is expressed in Equation (1). By a closed-form solution, Equation (1) yields to Equation (4) [42].
where w y is the transverse displacement of the fault and β is In analytical methods, L c is the soil yielding zone at the high curvature zone of the pipeline. In fact, L c is the distance between the pipe-fault intersection point and the first point with zero deflection on the pipeline (shown in Figure 2). L c is the first point after the fault line where deflection reaches zero. Based on Equation (4), where x = L c and w = 0, Regarding Equation (6), the high curvature zone length (L c ) has a direct relationship with EI and an inverse relationship with soil stiffness (k) [36].

Modeling Procedure
The comparison of buried pipeline performance at a strike-slip fault crossing using 3D solid and beam-spring FE modeling approaches was the main focus of the current study. Both the mentioned FEM modeling approaches were applied to investigate the problem of a buried pipeline at a strike-slip fault. The soil-pipe interaction modeling is the most complicated and important part of creating an FE model for the mentioned problem. Soilpipe interaction modeling parameters can severely affect the pipe performance responses during faulting. Therefore, to have accurate models and results using both modeling approaches, in the first step, the soil-pipe interaction properties and modeling of the buried pipeline were validated using pull-out tests and lateral sliding test simulations. In the 3D solid model, the soil material model and contact properties are determinants for the calibration of the soil-pipe interaction. In the beam-spring model's soil-pipe interaction, the spring properties are determinants. To validate the longitudinal soil-pipe interaction of the 3D solid FE model, three pull-out tests were simulated based on the model by Vazouras et al. [39] and the FE models' results were verified against their results. Based on the 3D solid model's pull-out tests, an identical beam-spring model was created, and its soil-pipe interaction results were verified against the 3D solid model. After verification of the 3D solid model soil-pipe interaction forces, three lateral sliding tests were created using the same 3D solid models, and lateral soil-pipe interaction forces were extracted. In the same manner, an identical beam-spring model for a lateral sliding test was created, and its lateral soil-pipe interaction forces were verified against the 3D solid model's results. After calibration of the soil-pipe interaction parameters for the 3D solid and beam-spring models, we attempted to create FE models for the problem of a buried pipeline at a strikeslip fault crossing. All the FE models were created in the multi-purpose finite element program ABAQUS [43]. Both mentioned approaches include both geometrical and material nonlinearity effects.

Longitudinal and Lateral Test Modeling
Two 3D FE models were created to study the longitudinal and lateral soil-pipe interaction and calibration of the 3D solid FE model soil-pipe interaction force-displacement relationship. The models were longitudinal pull-out and lateral sliding tests of the buried pipeline to show the force-displacement relationship of the longitudinal and lateral soilpipe interaction, respectively. All the FE models were created in the multi-purpose finite element program ABAQUS [43].
Both the 3D solid and beam-spring models' results were obtained for an X65 steel 36" pipeline with an outer diameter of 0.914 m, a thickness of 0.0095 m, Young's modulus of 21 Tpa, a Poisson ratio of 0.3, and a density of 7850 kg/m 3 . The Young's modulus for pipe material was assumed to be 100 times stiffer than X65 steel for decreasing the pipeline deformation effect in the soil-pipe interaction evaluation. The pipeline was assumed to be buried in undrained clay. The soil had a density of 2000 kg/m 3 , a Young's Modulus of 25 MPa, a Poisson ratio of 0.5 cohesion of 50 kPa, and a friction angle of 0 • . As with real cases, it was assumed that the buried pipeline was surrounded by a thin layer of sand. Thus, a frictional soil-pipe interaction was employed. In addition, the soil box was modeled with the dimensions of 20 m × 10 m × 5 m. Table 1 reports the pipe and soil properties.
In the 3D solid model, the soil material was defined as an elastic-perfectly plastic Mohr-Coulomb constitutive model. The pipe elements were 4-node shell S4R element type, and the soil elements were 8-node linear brick, reduced integration C3D8R element type. The geometrical nonlinearity effect was taken into account for all analyses by the Nlgeom method, conducted by the FE program of ABAQUS.
In the beam-spring model, B31, RB3D2, and CONN3D2 elements were employed for the pipe elements, rigid bodies, and soil, respectively. The soil-pipe interaction was modeled through equivalent nonlinear soil springs in the longitudinal, lateral (horizontal), and vertical directions extracted from the 3D FE simulation results based on the pull-out and sliding tests. The Nlgeom method was the same as the 3D solid model's and was implemented into the beam-spring model to have geometrical nonlinearity effects.

Longitudinal Pull-Out Test Analyses
As shown in Figure 3, three 3D solid FE cases for longitudinal pull-out tests were created to evaluate the longitudinal soil-pipe interaction. The FE analysis results were verified with the pull-out test results of Vazouras et al. [39] in Figure 4. As illustrated, good compatibility can be seen between the results of the current study and the results reported by Vazouras et al., 2015. Figure 5 shows the 3D solid models' longitudinal forcedisplacement diagram of the unit length of the pipeline for the pull-out test results in the three cases with friction coefficients of 0.2, 0.3, and 0.4, in which the longitudinal soil-pipe interaction curves are elastic-perfectly plastic curves for all the cases. Therefore, the pipeline slides when the longitudinal force reaches the corresponding displacement of the yielding force in each case. After verification of the 3D solid model for the pull-out test for the 0.2, 0.3, and 0.4 friction coefficient cases, the case with the 0.3 friction coefficient (τ max = 10 kPa) was selected as the main case for this paper. To compare the beam-spring model versus the 3D solid model, it was required to create a FE model with identical properties and geometry and verify that their soil-pipe interaction behavior was almost identical. To achieve this goal, a beam-spring model with equivalent properties to a 3D solid pull-out test model (µ = 0.3 and τ max = 10 kPa) was created. The soil-pipe interaction spring in the beam-spring model was created based the on 3D solid force-displacement curve in Figure 5, and the modeling details of the beam-spring model are clarified in Section 4. As presented in Figure 6, the beam-spring model's longitudinal soil-pipe interaction force-displacement curve is almost the same as the 3D solid model's, and it was verified against the 3D solid model's results.

Lateral Sliding Test Analyses
After verification of the 3D solid model's longitudinal soil-pipe interaction results, a 3D solid FE model with identical properties with the longitudinal pull-out 3D solid FE model was extended to simulate a lateral sliding test for the reproduction of the lateral soil-pipe interaction force-displacement relationship. The displacement contours for the lateral sliding test model are illustrated in Figure 7. Lateral soil-pipe interaction forcedisplacement curves for cases with friction coefficients of 0.2, 0.3, and 0.4 are shown in Figure 8. Based on Figure 8, all the cases with different friction coefficients have identical behavior in the elastic range and their discrepancy in a nonlinear range is almost negligible. In addition, it shows the soil material properties' effect on the lateral soil-pipe interaction is predominant in comparison with the contact properties. On the contrary, based on Figure 5, the contact properties have a higher impact on the longitudinal soil-pipe interaction force-displacement relationship in comparison with the soil material.
The lateral soil-pipe interaction spring properties of the verified beam-spring model in the pull-out test were extended based on the 3D solid lateral sliding test model results for the µ = 0.3 and τ max = 10 kPa case. As shown in Figure 9, the lateral soil-pipe interaction force-displacement curve of the beam-spring model in the lateral sliding test is almost the same as the 3D solid model's for the same case. Therefore, the beam-spring lateral soil-pipe interaction is verified to have almost identical behavior with the 3D solid model in the case with µ = 0.3 and τ max = 10 kPa.

FE Modeling
Both the verified 3D solid and verified beam-spring models were extended for the problem of a buried steel pipeline subjected to a 60 • strike-slip fault movement, based on the verified pull-out test and lateral sliding test models. In addition, both the 3D solid and the beam-spring modeling properties for the soil, pipe, springs, and contacts were identical with the FE models for the pull-out and sliding tests (except the pipe material properties). The soil box dimensions were 60 m × 10 m × 5 m; the 3D solid model is shown in Figure 10, and the beam-spring model is shown in Figure 11. A friction coefficient of 0.3 was assumed for contact modeling, which is equivalent to the demonstrated soil with a maximum shear stress equal to 10 kPa. The fault movements and boundary conditions were also applied to the soil box's faces. The properties of the pipeline modeled by the 3D solid and beam-spring methods were the same (see Table 1), and the soil-pipe interaction curves are shown in Figure 12. The fault displacement components were applied to the ends of the rigid elements at the end of the soil spring elements and the pipeline was free to move in the axial direction on both ends.  The pipeline's steel material plasticity was modeled based on the Ramberg-Osgood equation, shown in Equation (7) (Figure 13). The Ramberg-Osgood parameters for X65 steel material are referenced in Table 2.
in which a and r are equal to 38.31 and 31.51, respectively.

Analyses Results
The 3D solid and beam-spring FE models were analyzed with six cases with 60 • strike-slip fault displacements of (δ) 0.17D, 0.5D, 1D, 2D, 3D, and 4D, where D is the pipe diameter. Figure 14 shows the Mises stress outputs of the buried pipeline deformation modeled by the 3D solid method and its buckling for different fault movements. As seen, the pipeline starts buckling when the fault displacement is almost over 1D. Lateral displacements of the pipeline for both modeling approaches at the strike-slip fault crossing on a crown/invert of a pipe section (neutral axis) are compared in Figure 15. In general, good compatibility is seen between the lateral displacement response of the buried pipeline for the 3D solid and beam-spring models.
However, due to a shorter high curvature zone of the 3D solid model than that of the beam-spring model, a gap can be seen around the faulting zone between the two model lateral displacement responses at the high curvature zone. Based on Equation (6), by increasing the soil stiffness and decreasing the pipe cross-section bending stiffness (EI), the high curvature zone length (L c ) decreases. Since the material properties are identical in the 3D solid model and the beam-spring model, E is the same for both models. Therefore, a shorter high curvature zone length in the 3D solid models took place because of the local increasing soil stiffness at the faulting zone (k) and the decreasing of the inertia moment (I) of the pipe-section at the buckled locations. More precisely, the local stiffening of the soil in the 3D solid model is due to the local confinement of the soil at the fault line during the strike-slip fault movement. The beam-spring model cannot include the effects of soil confinement at the high curvature zone and local buckling of the buried pipe in the analysis, which are the main weakness of this modeling approach. Both modeling approaches present steady lateral displacement outside the high curvature zones. Furthermore, the lateral displacement results of the 3D solid model are lower than those of the beam-spring models at a steady-state range (further from the high curvature zone) for large fault movements (more than 2D). This happened because of the buckling of the buried pipeline at maximum bending moment locations (4-5 m from the fault plane) in the 3D solid model. Therefore, the strain energy release rate at the buckling location increased, and the strain rate in other locations along the pipeline decreased.
Stress and strain outputs are shown for the springline (point B in Figure 10) for the pipeline at both sides of the fault. The springline at the right side of the fault line has tensile stress, and the left side is in compression owing to the bending of the pipeline during faulting. As shown in Figures 16 and 17, the distance between the maximum tensile/compression stress and the strain of the buried pipeline and fault line for the case of the 3D solid model is shorter than that of the beam-spring model. This is again because of the shortening of the high curvature zone due to the local stiffening of the soil and the local weakening of the I at the buckled zones in the 3D solid models.    Figure 10) of the buried pipeline. Because of the symmetry of the problem, the results are shown only for one side of the pipe cross-section. The left springline of the pipe section has tensile stress and the right springline is in compression owing to the bending of the pipeline. As already discussed, local buckling in the pipeline occurs when the fault movement is almost over 1D. This local buckling affects the stress and strain responses. As shown in Figures 16 and 17, the strain and stress responses of the buried pipeline before the appearance of local buckling in the 3D solid model (δ ≤ 1D) were similar to the beam-spring model. However, the strain responses in the 3D solid model after buckling (δ > 1D) were less than the beam-spring model. As shown in Figures 16 and 17, the distance between maximum tensile/compression stress and strain of the buried pipeline and the faultline in the 3D solid models are shorter than the beam-spring model. This is again because of the shortening of the high curvature zone due to the local stiffening of the soil at the faulting zone and the local weakening of I at the buckled zones in the 3D solid models.
Because of the local stiffening of the soil around the high curvature zone and the local weakening of I at the buckled zones in the 3D solid models, the pipeline experienced higher stresses in the elastic range in comparison with the beam-spring model, and the maximum stress of the pipeline reached the yielding stress earlier than the beam-spring models. As shown in Figure 16, in the 3D solid model cases with (δ > 1D), in locations where the longitudinal stress of the springlines on the pipe cross-sections have excided the yielding stress of the pipe material (490 MPa) at the local buckling hinges' locations because of the dropping of the pipe's bending stiffness, the stress response of the pipeline locally drops, and the stress oscillates because of wrinkles on the pipe. By increasing δ, the length of this local stress drop-zone increases. In a similar manner (Figure 17), in the 3D solid model cases, with the strain in cases with (δ > 1D), the longitudinal strains of the pipeline in the local bending locations dropped and oscillated because of wrinkles along the pipe in those locations. Indeed, the pipeline strain on the springline in the beam-spring model drastically increased after the occurrence of local buckling in the 3D solid model, which represents that the pipeline was damaged. Therefore, large plastic strain in the beam-spring model is a criterion for damage detection in this modeling approach. Moreover, the pipeline experiences high stress at the elastic zone, and the pipeline yields in the 3D solid model earlier than in the beam-spring model owing to the local stiffening of the soil material around the high curvature zone.
Most of the buried pipeline designs are performed by beam-spring model analyses, and pipeline design guidelines are fitted for this modeling approach. However, none of the existing design guidelines have included the effect of local stiffening of the soil at the high curvature zone and the local weakening of the inertia moment of the pipe at buckled zones on buried pipelines' force-displacement and stress-strain behavior. Local stiffening of the soil in high curvature zones of the pipe and local weakening of the inertia moment of the pipe at buckled zones can increase the soil-pipe interaction forces and, correspondingly, increase stresses and damage vulnerability to the buried pipelines. In other words, it causes the underestimation of forces on the buried pipeline at earthquake faults crossing. Therefore, it is recommended to assign stiffer nonlinear soil-pipe interaction spring properties to the high curvature zone of the pipeline than the rest of the pipeline in the beam-spring modeling approach.

Conclusions
In this study, six cases of buried pipelines subjected to a 60 • strike-slip fault have been evaluated through beam-spring and 3D solid modeling approaches. Finally, the following have been found:

1.
In the 3D solid model, because of the pressure and shear forces caused by the fault movement on the soil and pipeline around the high curvature zone, local confinement takes place, and soil stiffness surrounding the pipeline increases locally around the high curvature zone.

2.
Because of the locally stiffening of the soil and the local weakening of the inertia moment of the pipe at buckled zones in the 3D solid models, the high curvature zone length of the pipeline is shorter than in the beam-spring models.

3.
Due to the locally stiffening of the soil and the local weakening of the pipe crosssection moment of inertia in the 3D solid models, the soil-pipe interaction forces increase. Therefore, the buried pipeline in the 3D solid models before the occurrence of local buckling experiences higher stress and strain in comparison with the beamspring models.

4.
In the 3D solid models, damage to the pipeline is visually detectable. Whereas, in the beam-spring models, observation of a local large strain on the pipeline is the best criterion for damage (e.g., buckling and ovalization) evaluation of the buried pipelines.

5.
In the 3D solid model, because of the local softening of the bending stiffness of the pipe cross-section, the pipe stress decreases at the buckled zone, which decreases the lateral deflection of the pipeline at further distances from the high curvature zone. 6.
Creating the 3D solid model is much more complex than the beam-spring model and it is easy to make mistakes in modeling for an amateur analyst. Moreover, modeling and analyzing take much time and cost. However, it can reproduce detailed and accurate results and cover all phenomena. 7.
Beam-spring models cannot include the effects of soil confinement at the high curvature zone and the local buckling of the buried pipe in the analysis, which are the main weaknesses of this modeling approach. 8.
Existing buried pipeline design guidelines have not included the effect of the local stiffening of the soil and the local softening of the pipe's moment of inertia in the analysis recommendation. Therefore, to prevent the underestimation of the forces on the pipeline in the beam-spring modeling approach, it is recommended to assign stiffer nonlinear soil-pipe interaction spring properties to the high curvature length of the pipeline than the rest of the pipeline.