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倍左右;近断层脉冲地震动作用下,管道上浮破坏和横向破坏两种破坏模式同时存在,且由于速度脉冲效应,管道横向破坏风险大于上浮破坏风险。     

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

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

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

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

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

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.     

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.     

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.     

Two-dimensional effective stress analysis of liquefaction of irregular ground including soil–structure interaction

An effective stress method is presented for analysis of seismic response and liquefaction of irregular ground including soil–structure interaction, based on an implicit–explicit finite element method. A pore water pressure is computed with iteration from the total stress considering an undrained condition. The simulated pore water pressure is in reasonably good agreement with the experimental data. The proposed method of analysis is compared with other well-known methods for a one-dimensional model, which is in good agreement. The present effective stress method is also applied to liquefaction problems involving a two-dimensional soil–structure model. The structure is modelled by not only a rigid model but also as a multi-degree-of-freedom system with bi-linear springs. The numerical results are considered to be significant from the viewpoint of earthquake engineering.     

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

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

跨越断层埋地管线地震反应数值分析

跨越断层埋地管线在地震中的破坏是非常严重的,地震本身和管土相互作用体系中都存在很多不确定性因素,所以管线在断层运动过程中反应比较复杂。本文利用有限元理论和数值模拟手段,建立了管土作用模型,采用非线性接触问题研究方法详细地分析了管线由断层运动而产生的反应,对影响管线的各种因素进行了分析,包括位错量、跨越角度、断层运动形式、埋设深度、初始轴向力、断层裂缝宽度、填覆土质和管径。通过研究,得到一些初步结论。