Three-dimensional poro-elasto-plastic model for wave-induced seabed response around submarine pipeline

The evaluation of the wave-induced excess pore pressure around a buried pipeline is particularly important for pipeline engineers involved in the design of offshore pipelines. Existing models for the wave-induced seabed response around submarine pipeline have been limited to poro-elastic soil behavior and de-coupled oscillatory and residual mechanisms for the rise in excess pore water pressure. To overcome the shortcoming of the existing models, in this study a three-dimensional poro-elasto-plastic soil model with submarine pipeline is established, in which both oscillatory and residual mechanisms can be simulated simultaneously. With the proposed model, a parametric study is conducted to investigate the relative differences of the predictions of the wave-induced pore pressure with poro-elasto-plastic model. Based on numerical examples, it can be concluded that the poro-elasto-plastic behaviors of soil have more significant influence on wave-induced pore pressure of seabed around submarine pipeline. As the seabed depth increases, the normalized pore pressures decrease rapidly at the upper part of seabed, and then change slightly at the lower part of the seabed. Soil permeability and wave period have obvious influence on the wave-induced normalized pore pressure.     

海底悬跨管线的地震响应和振动控制

渤海湾地区属于地震高发地区,在此处铺设海底管线必须考虑地震作用的影响,本文基于Ｍｏｒｉｓｏｎ方程建立了悬跨海底管线在地震作用下的动力方程,该方程考虑了海水对管道反应的影响,并对影响悬跨海底管线地震反应的因素进行了分析,在此基础上建议了控制管道反应的方法。     

Effects of cross-anisotropic soil behaviour on the wave-induced residual liquefaction in the vicinity of pipeline buried in elasto-plastic seabed foundations

In this paper, a two-dimensional integrated numerical model is developed to examine the influences of cross-anisotropic soil behaviour on the wave-induced residual liquefaction in the vicinity of a pipeline buried in a porous seabed. In the wave model, the RANS (Reynolds Averaged Navier–Stokes) equation is used to govern the wave motion. In the seabed model, the residual soil response in the vicinity of an embedded pipeline is considered with the 2-D elasto-plastic solution, where the phase-resolved shear stress is used as a source for the build-up of residual pore pressure. Classical Biot׳s consolidation equation is used for linking the solid-pore fluid interaction. The validation of the proposed integrated numerical model is conducted by the comparisons with the previous experimental data. Numerical examples show that the pore pressures can accumulate to a large value, thus resulting in a larger area of liquefaction potential in the given anisotropic soil compared to that with isotropic solution. The influences of anisotropic parameters on the wave-induced residual soil response in the vicinity of pipeline are significant. A high rate of pore pressure accumulation and dissipation is observed and the liquefaction potential develops faster as the anisotropic parameters increase. Finally, a simplified approximation based on a detailed parametric investigations is proposed for the evaluation of maximum liquefaction depth (z_L) in engineering application.     

波浪作用下成层海床孔隙水压力响应的简化分析

采用不排水条件下孔隙水压力发展模式,作为Terzagh i一维固结方程中考虑波浪循环作用所引起的孔隙水压力源项,对于成层海床建立了推广的一维动力固结方程,运用数理方程中的分离变量法与G reen函数求解了成层海床在波浪作用下残余孔隙水压力的发展规律,进而对成层海床的液化势进行了评判。对比计算与分析表明,海床表层土的渗透性及其厚度对于海床的整体抗液化性能具有显著的影响,低渗透性的表层导致海床孔隙水压力的显著积累,此时表层置换法是防治液化的有效途径。     

地震作用下可液化场地管道上浮动力响应及影响因素分析

地震作用下土体发生液化之后,由于超静孔隙水压力的产生和土体抗剪强度的降低,管道易发生上浮破坏。为研究管道上浮动力反应的影响因素,基于OpenSees有限元软件,通过目标反应谱和谱匹配等方法选取地震波,考虑不同管土特性和地震动特性,对地震作用下管道上浮动力反应进行了二维数值模拟。结果表明:土体相对密度、管径和管道埋深对管道上浮反应的影响较大,分别给出了土体相对密度、管径、管道埋深对管道上浮位移的影响规律及对应拟合公式;长持时地震动作用下,超静孔隙水压力消散较慢,管道上浮位移可达短持时地震动作用下管道上浮位移的2倍左右;近断层脉冲地震动作用下,管道上浮破坏和横向破坏两种破坏模式同时存在,且由于速度脉冲效应,管道横向破坏风险大于上浮破坏风险。     

Dynamic response of a porous seabed around pipeline under three-dimensional wave loading

The evaluation of the wave-induced pore pressure around a buried pipeline is particularly important for pipeline engineers involved in the design of offshore pipelines. Most previous investigations of the wave-induced dynamic response around an offshore pipeline have limited to two-dimensional cases. In this paper, a three-dimensional model including buried pipeline is established, based on the existing DYNE3WAC models. Based on the proposed numerical model and poro-elastic soil material assumption, the effects of wave and soil characteristics, such as wave period, water depth, shear modulus and permeability, and configuration of pipelines, such as pipeline radius and pipeline buried depth, on the wave-induced excess pore pressure will be examined. Numerical results indicated that the normalized excess pore pressures versus z/h near the pipeline increase as the obliquity angle, wave period and water depth increase, and they decrease as the burial depth and radius of pipeline increase above the pipeline. Soil permeability has obvious influence on the wave-induced normalized excess pore pressure, and different soil material will result in distinct computation results.     

地基液化导致沉箱码头破坏及地基加固方法的非线性数值分析

地震作用引发的地基液化,往往导致沉箱基础的破坏。本文基于Biot两相饱和多孔介质动力耦合理论,采用FE-FD耦合数值分析方法,对液化海床沉箱基础的地震反应进行非线性有效应力分析。在数值分析过程中,建立了以土骨架位移和超静孔隙水压力表达的us-pw动力固结方程和循环弹塑性本构模型,该方法能够很好地模拟地震作用下沉箱码头的动力特性及液化破坏的影响。通过数值模拟计算,分析了采用碎石桩进行置换砂区域的防液化加固方法,并就碎石桩处理区域的选择提出了建议。     

地震荷载下饱和砂土孔压增长时程计算方法

本文采用逐波累计方法,依据适于不同固结条件的孔压增量模型,计算每一应力循环的孔压变化,可给出非均等固结复杂情况下孔压实时增长过程。采用不同固结比两种砂进行若干典型地震荷载动三轴液化试验,试验结果与所提方法的计算结果有良好对应关系,表明本文方法可有效地描述地震荷载作用下非均等固结饱和砂土孔压增长过程。     

砂土液化引起大位移对地下管道影响的非线性分析

地下管线是生命线工程的主要部分,已经成为现代工农业生产和城镇生活的大动脉。已有震害调查表明,饱和砂土液化引起的地基大变形(侧向变形和沉降)是导致强震区生命线工程震害的主要原因。采用三维非线性有限差分分析方法来研究砂土液化引起的大位移对地下管道的破坏特征,分析砂土液化的斜坡变形特征、孔隙水的演化过程。结果表明,砂土液化引起的大位移对地下管道有破坏作用,导致管道变形规律与其斜坡的位移规律相同,地下管线的变形随着振动频率和幅值的增加其非线性增大。     

A practical numerical method for large strain liquefaction analysis of saturated soils

Accurate prediction of the liquefaction of saturated soils is based on strong coupling between the pore fluid phase and soil skeleton. A practical numerical method for large strain dynamic analysis of saturated soils is presented. The u–p formulation is used for the governing equations that describe the coupled problem in terms of soil skeleton displacement and excess pore pressure. A mixed finite element and finite difference scheme related to large strain analysis of saturated soils based on the updated Lagrangian method is given. The equilibrium equation of fluid-saturated soils is spatially discretized by the finite element method, whereas terms associated with excess pore pressure in the continuity equation are spatially discretized by the finite difference method. An effective cyclic elasto-plastic constitutive model is adopted to simulate the non-linear behavior of saturated soils under dynamic loading. Several numerical examples that include a saturated soil column and caisson-type quay wall are presented to verify the accuracy of the method and its usefulness and applicability to solutions of large strain liquefaction analysis of saturated soils in practical problems.     

Response analysis of a submarine cable under fault movement

Based on the performance of submarine cables in past earthquakes, an analytical method to determine cable performance under seabed fault movement is proposed in this paper. First, common types of earthquake damage to submarine cables are summarized, which include seabed displacement induced by fault movement, submarine landslides and seabed soil liquefaction, etc. The damage is similar to damage observed to buried pipelines following land earthquakes. The Hengchun earthquake of Dec. 26, 2006 is used as a case study. The M7.2 earthquake occurred in the South China Sea at 20:26 Beijing Time, and caused 14 international submarine cables to sever and break. The results show that the proposed method predicts damage similar to that observed in the Hengchun earthquake. Based on parametric studies of the influence of the water depth and the magnitude of the submarine earthquake, countermeasures to prevent damage to submarine cables are proposed.     

遮帘式板桩码头地基地震液化破坏机理

遮帘式板桩码头作为一种新型的板桩结构型式,其抗震性能研究是设计建造过程中的重要环节。在FEM-FDM水土耦合计算的平台上引入循环弹塑性本构模型,借助FORTRAN编程软件形成饱和砂土动力液化分析的数值方法,可有效模拟饱和砂土在地震动力作用下的非线性及大变形特性,同时也可模拟砂土液化流动对遮帘桩和前墙的动土压力。研究表明:地震作用下可液化土层超孔隙水压力比增长并发生较大的水平流动变形,对前墙的水平破坏大于竖向破坏;前墙剪力最大值位于海床与前墙交界处;遮帘桩剪力最大值位移与前墙底平行的位置;后拉杆拉力逐渐变大,前拉杆拉力逐渐变小。通过对板桩码头地震液化灾害的分析,可为抗震和抗液化设计提供参考依据。     

横向非一致地震激励下埋地管道的动力响应分析

为研究埋地管道在地震激励时管-土相互作用的动力响应问题,研发双向层状剪切连续体模型土箱,建立管G土相互作用有限元分析模型,对横向非一致地震激励下埋地管道地震响应进行数值模拟分析,并与试验结果进行对比.结果表明:数值模拟和振动台试验结果中的管道应变峰值均呈现出沿管道中间大两端小的现象,管道中间应变峰值最小达到两端的1．6倍左右;管道加速度、 土体加速度峰值均随着加载等级的提高而增大,涨幅愈加明显,多峰频率由0~10Hz逐渐向10~ 20Hz频域扩散,管道运动更为自由;土体位移随着加载等级的提高呈现逐级增大的现象,在加载等级增加到0．4g 时位移曲线斜率减小,土体非线性表现明显.数值模拟和振动台试验对比分析的结论表明数值模拟分析的合理性和试验结果的可靠性,为研究横向非一致激励对埋地管道地震响应的影响提供了依据.     

Evaluation of sheet pile-ring countermeasure against liquefaction for oil tank site

The evaluation of a countermeasure against liquefaction which uses a sheet pilering for oil tank sites is presented. The simulation of earthquake responses observed at tank sites with and without sheet pile-ring is first performed to validate the three-dimensional finite element numerical model. Using the numerical model, liquefaction analysis is performed and the excess pore water pressure generated in the soil and the settlement of tank are investigated. The comparison of two- and three-dimensional models is also conducted to assess the applicability of two-dimensional analysis. The results show that the numerical model could simulate the observed earthquake responses of tank-ring-soil system, and that the excess pore water pressure and the settlement of the tank could be significantly reduced using a sheet pile-ring. The two-dimensional analysis proves to be capable of representing the main features of the dynamic response of the three-dimensional tank-ring-soil system.     

Investigation of the dynamic interaction between a high-pressure pipeline and the moving liquid inside under seismic loading

The response of an open high-pressure pipeline containing moving liquid to horizontal seismic loading is investigated. The problem of the interaction between the pipe and the liquid inside it is solved after thorough kinematic analysis, involving two-stage application of the Coriolis theorem. The basic axial loading on the pipeline is also taken into account. As a result, the equation of motion of an arbitrary one-span part of the pipeline under seismic loading is obtained. The harmonic balance is used to find the solution. As an illustration to the theory a numerical example is given with the data of a real water power station pressure pipeline under a past earthquake record.     

Uplift mechanism for a shallow-buried structure in liquefiable sand subjected to seismic load: centrifuge model test and DEM modeling

Based on a centrifuge model test and distinct element method(DEM), this study provides new insights into the uplift response of a shallow-buried structure and the liquefaction mechanism for saturated sand around the structure under seismic action. In the centrifuge test, a high-speed microscopic camera was installed in the structure model, by which the movements of particles around the structure were monitored. Then, a two-dimensional digital image processing technology was used to analyze the microstructure of saturated sand during the shaking event. Herein, a numerical simulation of the centrifuge experiment was conducted using a two-phase(solid and fl uid) fully coupled distinct element code. This code incorporates a particle-fl uid coupling model by means of a "fi xed coarse-grid" fl uid scheme in PFC3D(Particle Flow Code in Three Dimensions), with the modeling parameters partially calibrated based on earlier studies. The physical and numerical models both indicate the uplifts of the shallow-buried structure and the sharp rise in excess pore pressure. The corresponding micro-scale responses and explanations are provided. Overall, the uplift response of an underground structure and the occurrence of liquefaction in saturated sand are predicted successfully by DEM modeling. However, the dynamic responses during the shaking cannot be modeled accurately due to the restricted computer power.     

某核电站导流堤地基液化及地震永久变形分析

在地震反应分析中常采用总应力分析法,但总应力法没有考虑孔隙水压力的变化规律和液化随时间的发展过程。基于二维有效应力动力分析方法,结合Biot动力固结理论,采用自行开发的有效应力动力分析程序对某核电站导流堤地基进行液化分析,给出在地震作用下砂土层的液化范围,并计算出永久变形。所得结论可以给类似工程提供理论指导。     

Numerical modeling for wave–seabed–pipe interaction in a non-homogeneous porous seabed

An evaluation of the wave-induced pore pressures and effective stresses has been recognized by marine geotechnical engineers as an important factor in the design of marine pipelines. Most previous investigations for such a problem have considered the pipeline as a rigid material. Thus, the internal stresses within the pipeline have not been examined in the wave–seabed–pipe interaction problem. In this paper, we consider the pipeline itself to be an elastic material, and link the analysis of the pipeline with the wave–seabed interaction problem. Based on the numerical model presented, the effects of pipe geometry and variable soil characteristics on the wave-induced pore pressure and internal stresses will be discussed in detail. It is found that the internal normal stresses in the angular direction (σ_pθ) and shear stress (τ_p) within the pipe are much larger than the amplitude of wave pressure at the surface of the seabed.     

水域隧道地震响应分析

本文基于Biot动力固结理论和弹性动力学理论,考虑海床（土壤）的两相性、黏弹性人工边界及流（水）-固耦合作用,建立了隧道-土-流体相互作用的力学模型,讨论了P波作用下有无水的情况以及水深、水域隧道埋深、海床土性质和地震波入射角等因素对隧道及其周围海床应力的影响。结果表明:隧道周围海床土的孔隙水压力和隧道内应力随着水深的增加而增加;地震波特性和海床土特性对隧道的内应力和海床土的孔隙水压力均有较大的影响;海床土的渗透性和隧道埋深对隧道的内应力影响较小,而对隧道周围海床土的孔隙水压力影响较大;地震动的入射角对隧道的内应力和隧道附近土层的孔隙水压力均有较大影响。