Received 2 June 2016; accepted 21 June 2016; published 24 June 2016
Pipeline transportation is the main means of energy transportation, and it is the lifeline engineering of modern society. China is an earthquake prone country   , buried pipeline as an important part of the long pipeline project that is inevitable to be built in the earthquake prone area. Now the experts and scholars of domestic and foreign on the seismic response laws of buried pipeline have been done a large number of studies  -  . The commonly used pipe soil mutual interaction models are: elastic foundation beam model  , soil spring model  and nonlinear contact model  . The spring soil analysis model of buried pipeline in tension and compression on viscoelastic boundary was established by being used ADINA finite element analysis software on the basis of previous research in this paper  -  . The buried pipeline seismic response was analyzed, and it was compared with the simple boundary model or entity model of soil to be found the optimal simplified model, which had certain guiding significance for the actual pipeline seismic design.
2. Finite Element Implementation of Buried Pipeline Model with Tension and Compression Spring
2.1. The Selection of the Unit
The pipe body was adopted the four node shell unit because the buried pipeline was a three-dimensional modeling. The soil around the buried pipeline was replaced by soil spring, the spring unit was chosen when being modelled in ADINA. The nonlinear analysis was carried out, and the tension spring was used to simulate the soil spring. The nonlinear properties were needed to open the panel, and the appropriate force displacement relation and the parameters of the nonlinear elastic material were set up.
2.2. Determination of Material Parameters
The Q235 was selected as the pipeline steel in this paper, the curves of experiment of simplified bilinear stress-strain relationship was used in ADINA finite element software. The interaction between in soil and pipe was simulated through three directional nonlinear spring, and the soil spring was included three direction (axis, horizontal, vertical). The pipe body was adopted four node shell elements, and was used steel pipes, pipe diameter was equal to 0.764 meter, pipeline wall thickness was equal to 0.01 meter, the burial depth of the pipe was equal to 2 meter, Young’s modulus was equal to E = 2.1e11pa, density was equal to ρ = 7800 kg/m3, poisson’s ratio was equal to μ = 0.3, it was above underground water level, and it was soil type in IV site.
2.3. Seismic Wave
The seismic waves of El Centro were used, the duration time of the earthquake was 10 second. The time interval was 0.02 second. As shown in Figure 1 and Figure 2, the seismic waves were horizontally input in order to be
Figure 1. El Centro acceleration record chart.
Figure 2. Shell unit.
considered the effect of ground motion of different intensity. This wave was had a high magnitude, and it was more suitable to be used as a design basis than other seismic waves. The input was adjusted according to the requirement of the 8 degree earthquake, seismic acceleration records were from north to south.
3. The Determination of Effective Area of Viscoelastic Boundary Model
3.1. The Simulation Method of Viscoelastic Boundary in ADINA
The viscoelastic boundary could be implemented in ADINA, which was a linear single degree of freedom spring damping unit in ADINA. In this spring, the stiffness matrix and damping matrix of the ground spring were all one-dimensional.
3.2. The Determination of Effective Area of Viscoelastic Boundary Model
Simple boundary: Fixed boundary was used in both ends of pipeline;
Viscoelastic boundary: Spring damper unit was adopted in the two ends of the pipe.
The pipeline was used four node Shell unit, as shown in Figure 3. The pipeline was divided into unit grids along the pipeline circumferential and axial direction, as shown in Figure 4. The circumferential direction of the pipe was divided into 16 units.
Viscoelastic boundary: The seismic response displacement and response stress of the pipeline length of 30 m, 40 m, 60 m and 80 m were calculated and analyzed.
As shown in Table 1, the maximum horizontal displacement and the maximum compressive stress of the buried pipeline were increased with the increase of the pipe length. With the increase of the calculation range, it was tended to be stable. When the range of calculating was increased to 60 m, it was basically stable, and was approximately 60d. So we were in the calculation area of the viscoelastic boundary, the pipe length was 60 m.
3.3. The Comparative Analysis of Calculation Results of Simple Boundary and Viscoelastic Boundary Pipeline Model
Centro El acceleration wave in horizontal direction was introduced, and the difference between the simple
Figure 3. Viscoelastic boundary of soil spring model of underground.
Figure 4. The solid soil model.
Table 1. Comparative analysis of the maximum horizontal displacement force effect of different length of buried pipelines under earthquake and the maximum compressive stress.
boundary and the viscoelastic boundary pipeline and the soil mass was observed.
3.4. The Solid Soil Model
The solid modeling was used in the soil, as shown in Figures 4-6. Although it was relatively time consuming, the relative efficiency was relatively low, but the calculation result was relatively reliable. In this paper, the ADINA software was used to establish the solid model of the buried pipeline, and the pipeline length was 60 meters. The purpose of this paper was to compare the numerical results of the finite element model of the buried pipeline with the finite element model of the buried pipeline.
3.5. The Comparative Analysis of Calculation Results of Solid Soil Model, Simple Boundary and Viscoelastic Boundary Pipeline Model
As shown in Table 2, the calculation results of solid soil model, simple boundary and viscoelastic boundary pipeline model are compared and analyzed.
Figure 5. Finite element model of buried pipeline.
Figure 6. Finite element model of soil.
Table 2. Horizontal pipe 27,694 node displacement and stress.
4. Results and Discussions
The finite element calculation model of the buried pipeline based on the viscoelastic boundary was established. By being based on the finite element analysis software ADINA, the calculation area of the buried pipeline with tension and compression spring based on the viscoelastic boundary under the action of Centro El acceleration wave was analyzed. It was concluded that the simplified calculation area was 60 times the diameter.
By being used ADINA finite element analysis software, the ground motion response of the pipeline under the action of the IV field and the seismic action was analyzed and compared with the simple boundary, the viscoelastic boundary and the solid pipeline model. It was found that it was more practical to simulate the effect of the tension spring and the soil of the solid earth. But the solid earth model was time consuming. The maximum compressive stress was much greater than the simple boundary entity model of viscoelastic soil boundary horizontal maximum compressive stress, which was larger error, and the viscoelastic boundary was achieved, and theoretical basis for seismic design of pipeline was provided.
 Wang, X.Z., Dong, F. and Liu, Y. (2012) Research Progress of Long Distance Pipeline Safety under the Action of Earthquake Disaster. Safety and Environmental Engineering, 18, 28-33.
 Jing, L., Yi, X.D., et al. (2006) 3D Viscous-Spring Artificial Boundary in Time Domain. Earthquake Engineering and Engineering Vibration, 5, 93-102.