www.arpapress.com/volumes/vol17issue1/ijrras_17_1_04.pdf THE STABILITY OF THE PIPELINE LAID ON A PORO-ELASTIC SEABED Ren Yanrong Science School, Beijing Institute of Civil Engineering and Architecture, Beijing 100044, China E-mail: ryr@bucea.edu.cn ABSTRACT The dynamic response of buried submarine pipeline subjected to wave must be take into account in the design. In this paper the two-dimensional finite element model has been established, the poro-elastic model has been adopted, the deformation and force condition of pipeline have been calculated under the load gravity and the combined load of weight and harmonic wave force respectively. The key issue is to simulate the interaction among flow, pipe and soil. The factors such as seabed thickness, water depth, pipe s diameter, elastic modulus, wave period and permeability coeffient have been considered. The numerical results show that pipe s weights, horizontal load, pipe s diameter, water depth, elastic modulus have an important effect, but seabed thickness, wave period have little effect. The research results can provide scientific gist for the pipeline s design. Key words:wave loading, seabed pipeline, pipe and soil interaction, numerical analysis. 1. INTRODUCTION Because of the ever increasing engineering activities in offshore and coastal regions, phenomena such as wave attenuation caused by percolation in seabed, seafloor instability, and wave induced soil response around submarine structures have attracted the attention from coastal and marine geotechnical engineers. Particularly, since inadequate design of marine structures and foundations can lead to costly failures such as jacking and floatation of pipeline, it is necessary to improve our knowledge on the interaction among wave, seabed and submarine structures. In general, the fluctuations of wave pressure at the surface of the seabed exert excess pore pressures and effective stresses, which have been recognized as dominant factors in analysing the instability of a seabed. Once the seabed under the pipeline becomes unstable, the instability of the pipeline can occur. It has also been reported that the soil supporting the pipeline may fail due to the liquefaction of the seabed, resulting in the self-burial of the entire pipeline(bijker and leeuwenstein [1],1984; Sumer and Fredsoe [2],1997);this is the reason of causing the instability of the pipeline, but another reason is that the lateral resistance of pipes on seabed can not offset against the horizontal component of the hydrodynamic force, pipe breakout occurs. Although the importance of the wave-soil-structure interaction phenomenon has been addressed in the literature, most previous investigations have only concerned with the soil response owing to two-dimensional progressive wave or three-dimensional short-crested waves(jeng [43],1997).Only few investigations have been carried out for the wave-soil-pipe interaction problem in recent years(mei and Foda [44],1981).However, their research only considered the seabed response under the pipeline; Jeng [54] et al(1998)carried out the seabed response around a pipe laid on the seabed, in their research, soil permeabilty and shear modulus varied with the seabed depth, however, they couldn t consider the contact surface between the pipe and seabed and didn t describe the behavior of the pipeline on the seabed. Gao fuping [6] has proposed an improved analysis method for the on-bottom stability of a submarine pipeline which is based on the relationship between Um/gD 0.5 and Ws/D2. The proposed analysis method may provide a helpful tool for the engineering practice of pipeline on-bottom stability design. Gao fuping [7] has employed a hydrodynamic loading method in a flow flume for simulating ocean currents induced submarine pipeline stability on a sandy seabed. The pipeline stability in currents is compared with that in waves, which indicates that the pipeline laid directly upon the sandy seabed is more laterally stable in currents than in waves. Liu Jing [8] has simulated the interaction among flow, pipe and soil. The factors such as contact effect, frictional coefficient between pipe and soil, buried depth, pipe radius. Based on the numerical results the vertical displacement and hoop stress should are underestimated dramatically without considering the contact effect in the case of smaller buried depth. Meanwhile, porous water pressure in coarse sand attenuates slower than that in fine sand, so pipe embedded in fine sand is more stable and safer than that in coarse sand. The object of the paper is to describe the movement behavior of the pipeline on the seabed by means of the contact surface prescribed in the general program ABAQUS.A two-dimensional finite element model has been established. Based on the model, the mechanism of the wave-induced seabed response in the vicinity of a pipeline and the behavior of the pipeline has been discussed in detail. 27
2.THE ESTABLISHMENT OF THE FINITE MODEL In this study,we consider a soil column of a finite thickness(h)beneath a pipeline(with a radius R),as depicted in fig.1. Rigid impermeable bottom Fig.1. Definition of wave-seabed-pipe interaction problem In this study, the wave crests are assumed to propagate in the positive x-direction, the seabed is simulated by the poro-elastic model, and the pipe is considered to be fully elastic. The rectangular elements with 8 nodes and 232 elments will be used to model the pipe,the coupled pore pressure elements with 4 nodes and 2680 elements will be used to model the seabed. The finite element model is as follows: Fig.2 Finite element model 3.BOUNDARY VALUE PROBLEM 3.1. Governing equations In this study, the seabed is modelled by a poro-elastic theory. The consolidation equation(biot,1941),extending from Terzaghi s theory, is generally accepted as the governing equation for flow of a compressible pore fluid in a compressible porous medium. The consolidation equation is: 2 n w p w p K t K t In which: is water density, n is porosity of soil, p is pore pressure, t is time, volumetric strain of the soil and w porosity is as follows: 1 1 S u v, K w P wo x y The soil skeleton is assumpted to be ideal and isotropic, The equation is as follows under the plane strain: 28
u u x 2G( ) x 1 2u v u y 2G( ) y 1 2u u v xy G( ) yx y x 2 u 1 v v x ( x z ) x E 1 v 2 w 1 v v z ( z x ) x E 1 v The body force is ignored, so the equilibrium equation has changed to be : x xy p x y x xy y p x y y 2 G p G u, 1 2u x x 2 G p G v. 1 2u y y 3.2. Boundary conditons 3.2.1. Defining the Displacement boundary conditions The displacement of the bottom is zero, and there is no normal flow. P u w 0, y h Z 3.2.2. Defining the force boundary conditions The dynamic pore pressure P has been produced because of the wave loadings, the expression is: p( x,0; t) p0 cos( kxt) wh p0 cosh kd 2 The dynamic pressure P is Harmonic oscillation, the p 0 is the amplitude. The direction of force is positive at the wave crest, negative at the trough of wave. The stress and pore water pressure are produced because of the dynamic wave loadings. The strong dynamic water pressure can produce and assumulate superpore pressure from the medium dense soil to loose soil. The effective stress and tangent stress disappear to zero at the seabed. o, p( x,0; t) p0 cos( kxt) when y 0, y The wave pressure of the surface of the seabed is: wh p0 cosh kd 2 In which:d -water depth. H-wave height k-wave number, -angular speed, xy 2 l 2 T Because of the pipeline s exsitence, the wave will produce the turbulence. The wave force has been computated according to the Morsion Equation, the horizontal force will be imposed on the left of the pipeline, the lift force will be imposed on the bottom of the pipeline. 29
z(cm) 4. NUMERICAL RESULTS 4.1. Parameters needed Some parameters needed are as follows: Table1 Parameters Wave parameters Wave period 1(s),15.0(s) Water depth 20m Wave length 121.1m Wave height 5.0m Water bulk module 210 9 N/m 2 Sea water density 1030 kg/m 3 Soil parameters Soil thickness Degree of saturation 0.975 Soil s porosity 0.4 Possion s ratio 0.4 Shear stiffness 5 10 6 N/m 2 Soil s coefficient of permeability 10-4 m/s(coarse sand)10-5 m/s (fine sand) 10-2 m/s (gravel) The bulk modulus of soil particles 410 10 N/m 2 The bulk modulus of soil skeleton 10010 6 N/m 2 Pipe s diameter 0.4m 4.2. Numerical results 4.2.1. The effect of permeability coefficient 1.6 1.4 1.2 1.0 0.8 0.6 0.4 0.2 1E-5 1E-4 1E-3 permeability coefficient Fig 3 The relationship between pipe s penetration and permeability coefficient From the figure, it can be found that the pipe s penetration will increase with the permeability coefficient s decreasing, so the pipe is stable. It is because that the pipe have the penetration under the effect of pipe s weights and horizontal load, another reason is that the water loss in the fine sand is larger than that of coarse sand under the same horizontal load, so the soil particle under the pipe is easy to be loose and it is easy to cause the penetration for the pipe, which is advantageous for the pipe s stability. 4.2.2. The effect of the pipe s weight The computation has been carried out by changing pipe s weight when pipe s diameter is constant. 30
lateral coefficient horizontal displacement(mm) pipe's penetration(mm) 15.0 14.5 14.0 13.5 13.0 12.5 12.0 11.5 0.6 0.8 1.0 1.2 1.4 1.6 pipe's weight Ws(kN/m) 13 12 D=0.4m 11 10 9 8 0.6 0.8 1.0 1.2 1.4 1.6 pipe's subweights Ws(kN/m) Figure 4 The relationship between pipe s penetration and pipe s weight Figure 5 The relationship between pipe s horizontal displacement and pipe s weight 1.0 0.8 0.6 0.4 0.2 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 pipe's penetration Ws(kN/m) Figure 6 The relationship between lateral coefficient and pipe s weight 31
horizontal displacement(mm) pipe's penetration(mm) From the above figure, it can be found that the pipe s penetration will increase with pipe s weight increasing, but horizontal displacement will be decrease,which can explain that the pipe s stability will be increase. The average value of lateral coefficient is 0.887. 4.2.3. The effect of environmental load In order to consider the effect of environmental load, when pipe s diameter is 0.4, pipe s weight is 0.656kN, the results will be figured are as follows by changing pipe s environmental load. 12 10 8 6 4 2 0 10 15 20 25 30 35 environmental load(kn) 14 12 10 8 6 4 2 0 10 15 20 25 30 35 horizontal load(kn) Figure 7The relationship between pipe s penetration and environmental load Figure 8The relationship between horizontal displacement and environmental load From figure 7 and 8,it can be concluded that the pipe s horizontal displacement and penetration will increase with the horizontal load s increasement. 32
4.2.4. The effect of elastic modulus P/P 0 E=5e5 E=2.5e5 E=7.5e5 0.2 0.4 0.6 0.8 1.0 1.2 1.4 Figure 9 The relationship between pore pressure and elastic modulus xx '/P 0 E=5e5 E=2.5e5 E=7.5e5 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 Figure 10 The relationship between x-stress and elastic modulus yy '/P 0 E=5e5 E=2.5e5 E=7.5e5 0.2 0.4 0.6 0.8 1.0 1.2 Figure11 The relationship between normal stress and elastic modulus 33
penetration(mm) xy '/P 0 E=5e5 E=2.5e5 E=7.5e5 0.1 0.2 0.3 0.4 0.5 0.6 Figure 12 The relationship between tangent stress and elastic modulus 20 18 16 14 12 10 8 6 4 200000 300000 400000 500000 600000 700000 800000 elastic modulus(n/m 2 ) Figure13The relationship between pipe s penetration and elastic modulus From the above figures, it can be concluded that the value of elastic modulus has decisive effect on the distribution the pore pressure. The pore pressure is decreasing with the elastic modulus around the seabed, but increasing when is bigger than 0.5. The effective stress is increasing with the decreasing of the elastic modulus. The pipe s penetration will decrease with the increase of elastic modulus. 4.2.5. The effect of period P/P0 T=12.5s T=15s T=17.5s 0.2 0.4 0.6 0.8 1.0 Figure14 The relationship between pore pressure and wave period 34
xx '/P 0 T=12.5 T=15 T=17.5 0 2 4 6 8 0.10 0.12 0.14 Figure 15The relationship between x stress and wave period yy ' /P 0 T=12.5s T=15s T=17.5s 0 5 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 Figure 16 The relationship between normal stress and wave period xy '/P 0 T=12.5s T=15s T=17.5s 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 Figure 17 The relationship between tangent stress and wave period 35
pipe's penetration(mm) 2.6 2.4 2.2 2.0 1.8 1.6 1.4 12.5 13.0 13.5 14.0 14.5 15.0 15.5 16.0 16.5 17.0 17.5 18.0 18.5 19.0 19.5 2 wave period(s) Figure 18 The relationship between pipe s penetration and wave period From the above figures, it can be found that wave period had not obvious effect on the soil s stress and pore water pressure. The normal stress and pipe s penetration had little increase with the increase of the wave period. The reason is that the pore pressure will increase because of the increase with wave load. 4.2.6. The effect of water depth P/P 0 d=40m d=60m d=100m 0 1 2 3 4 5 Figure19 The relationship between pore pressure and water depth xx '/P 0 d=40m d=60m d=100m 0 1 2 3 4 5 Figure 20 The relationship between x-stress and water depth 36
yy '/P 0 d=40m d=60m d=100m 0 1 2 3 4 5 6 7 8 Figure 21 The relationship between normal stress and water depth xy '/P 0 d=40m d=60m d=100m 0 5 10 15 20 25 Figure 22 The relationship between tangent stress and water depth From the above figures, it can be found that the pore pressure, x-stress will increase with the increase of water depth but tangent stress on the contrary.the normal stress will increase with the decrease of water depth, but when is larger than 0.4,the normal stress will increase with the water depth. 4.2.7. The effect of seabed thickness The results are as follows: p /P 0 h=10m h=20m h=40m 0.2 0.4 0.6 0.8 1.0 Figure 23 The relationship between pore pressure and seabed thickness 37
' xx /p 0 h=10m h=20m h=40m 0.5 1.0 1.5 2.0 2.5 3.0 Figure 24 The relationship between x-stress and seabed thickness ' yy /p 0 h=10m h=20m h=40m 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 3.2 Figure 25 The relationship between normal stress and seabed thickness ' xy /p 0 h=10m h=20m h=40m 0 1 2 3 4 5 6 7 8 9 10 11 12 Figure 26The relationship between tangent stress and seabed thickness 38
horizontal displacement(mm) pipe's penetration(mm) 3.0 2.5 2.0 1.5 1.0 4 6 8 10 12 seabed thickness(m) Figure 27 The relationship between pipe s penetration and seabed thickness 19 18 17 16 15 14 13 12 11 4 6 8 10 12 seabed thickness(m) Figure 28 The relationship between horizontal displacement and seabed thickness 4.2.8. The effect of pipe s diameter The pipe s diameter is an important parameter in pipeline s design. The computation results are as follows: p /p 0 D=0.4m D=0.8m D=1.0m 0.2 0.4 0.6 0.8 1.0 Figure 29 The relationship between pore pressure and pipe s diameter 39
' xx /p 0 D=0.4m D=0.8m D=1.0m 0 2 4 6 8 10 Figure 30 The relationship between x-stress and pipe s diameter ' yy /P 0 D=0.4m D=0.8m D=1.0m 0 2 4 6 8 10 Figure31 The relationship between normal stress and pipe s diameter ' xy /p 0 D=0.4m D=0.8m D=1.0m 0 1 2 3 4 5 Figure 32 The relationship between tangent stress and pipe s diameter 40
pipe's penetration(mm) 5 4 3 2 1 0.4 0.5 0.6 0.7 0.8 0.9 1.0 pipe's diameter(m) Figure33 The relationship between pipe s penetration and pipe s diameter From the above figures, it can be found that when is less than 0.4,pore pressure will increase with the pipe s decrease, but on the contrary when is larger than 0.4. Tangent stress will decrease and the pipe s penetration will increase with the increase of pipe s diameter. 5. CONCLUSIONS (1)The pipe's penetration,horizontal displacement,the distribution of pore pressure can be obtained by changing the parameters of soil and horizontal load. (2)In these parameters, pipe s weights, horizontal load, pipe s diameter, water depth, elastic modulus have an important effect, but seabed thickness, wave period have little effect. ACKNOWLEDGMENT This paper is supported by The Importation and Development of High-Caliber Talents Project of Beijing Municipal Institutions. The Fund number is 21301413105. 6. REFERENCES [1]. Bijker E W and Leeuwestein W. Interaction between pipelines and the seabed under the influence of waves and currents. Seabed mechanics, Proc. Int. Union of Theoretical and Appl. Mech. Conf., B. Denness,ed., Graham and Trotman, Gettysburg, Md.,1984;235-242. [2]. SummerB.M.,Fredsoe,J.,1997. Hydrodynamics Around Cylindrical Structures, World Scientific, Singapore. [3]. Jeng D S. Wave-induced Seabed Instability in front of a Breakwater [J]. Ocean Engineering, 1997, 24(10):887-917. [4]. Mei CC.,Foda M F. Wave-induced Stresses around a Pipe laid on a Poro -elastic Seabed.Geotechnique,1981b,31: 509-517 [5]. Jeng, Dong-Sheng,Wave-induced seabed response in a cross-anisotropic seabed in front of a breakwater: an analytical solution. Ocean Engineering,vol. 25 issue 1 January, 1998. p. 49-67 [6]. Gao fuping,jengds,wu YX, Improved Aanlysi Method for Wave-Induced Pipeline Stability on Sandy Seabed, Journal of Transportation Engineering-ASCE,132(7):590-596,2006. [7]. Gao fuping,yan SM,Yang B,WU YX,Ocean currents-induced pipeline lateral stability on sandy seabed, Journal of Engineering Mechanics-Asce,133(10):1086-1092,2007. [8]. Liu Jing,L I Minggao,, Study on Interaction between Submarine Buried Pipeline and Porous Elastic Seabed, Water Resources and Power,27(5),159-162,2009. 41