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.     

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.     

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.     

辽东湾近海海底土小应变动力特性试验研究

开展场地地震动反应分析是众多海洋工程活动顺利实施和长期安全运行的重要保障,其中查明海底土小应变动力特性是十分关键的。文章借助海上原位钻孔取样和室内共振柱仪对辽东湾近海海底土小应变动力特性开展试验研究,探讨前人总结提出的滨海海底土初始动剪切模量G_max预测经验公式在辽东湾海域的适用性,对比文章试验结果与周边海域已有海底土剪切波速测试结果,分析小应变范围内海底土动剪切模量的衰减特征和阻尼比发展演变趋势,并同时进行定量化模拟预测,评价海底土动剪切模量衰减曲线的归一化特征。研究结果表明：(1)共振柱试验得到的原状海底土剪切波速与周边海域已有海底土剪切波速测试结果吻合良好;(2)相比起陆地土,海洋土在小应变范围内归一化动剪切模量G/G_max衰减得更慢,阻尼比增长得也更慢;(3)通过引入临界剪应变,不同土类和埋深对应的海底土小应变动剪切模量衰减曲线可以进行归一化。研究成果对合理开展海洋工程场地地震动反应分析评价具有一定的指导意义。     

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的动力固结理论和弹性动力学理论建立波浪荷载作用下孔隙海床与海底隧道的动力分析模型,并同时考虑海底隧道与海床之间的接触效应、边界效应等对海底隧道内力的影响。最后通过变换海床的变形模量、渗透系数和海底隧道的半径、埋深等,观察其对海床土孔隙水压力和海底隧道内力的影响,为海底隧道的设计提供依据。     

水平向不均匀土对箱基动力阻抗的影响

采用薄层元素法和有限单元法,建立了地基-箱形基础动力相互作用的三维分析模型,该模型可考虑箱形基础周围土的不均匀性影响.利用该模型,分析了水平方向不均匀土对箱形基础的水平和摇摆阻抗的影响.分析了水平方向不均匀土的厚度、剪切弹性模量以及材料阻尼比的影响.分析结果表明:水平方向不均匀土使得箱形基础的各种阻抗均降低,特别是阻抗的虚部.水平方向不均匀土对摇摆阻抗的影响较大.水平方向不均匀土的剪切弹性模量的影响较大而其材料阻尼的影响较小.随着水平方向不均匀土的厚度的增加和基础埋深的增加,其影响增大.     

Transient dynamic response of fluid-saturated soil under a moving cyclic loading

Based on the Theory of Porous Media (TPM), a mathematical model of a two-dimensional incompressible fluid-saturated elastic soil is established, and the periodic boundary conditions are presented to analyze the transient dynamic response of this soil under a moving cyclic loading. The differential quadrature method (DQM) and the second-order backward difference scheme are applied to discretize the governing equations on the spatial and temporal domains, respectively. As application, a typical two-dimensional wave-induced transient problem with a seabed of finite thickness is analyzed, and the numerical results are compared with the analytical results presented in the present work. In addition, a transient dynamic response of fluid-saturated soil under limit moving vehicle loadings is studied. The effects of the velocity of vehicle and the volume fraction on the settlement and the pore water pressure are studied.     

土层对地震的随机反应分析

利用随机振动理论,本文研究了非均匀土层对地震的动力反应问题,假定土层的剪切模量随深度呈指数函数形式分布,关于基岩输入地震动加速度的功率谱密度函数,考虑了两种形式：白噪声谱和有色谱,为了比较起见,文中考虑了均匀土层的情形,数值计算表明;（１）在基岩输入地震动加速度的功率谱密度为白谱的情况下,土层的最大期望反应（相对位移,绝对加速度,剪应变）大于基岩输入地震动的功率谱为有色谱的情形;（２）在土层的上部     

浪作用下海洋底床动态响应的研究

有关波浪作用下的底床动态响应越来越引起人们的重视。本文从海洋土的特点出发，针对各向同性底床和各向异性底床，详细论述了在线性波加载下，波浪衰减和底床动态响应这两方面的研究现状，在分析和比较已有研究成果的基础上，对今后的研究方向提出了自己的看法。     

水域隧道地震响应分析

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

埋深对地铁车站地震反应的影响分析

关于埋深对地下结构地震反应的影响的研究对象多见于地下隧道,对地铁车站地震反应受埋深影响变化规律缺乏深入研究。本文基于ANSYS有限元软件,采用改进的简化方法建立三种不同埋深的地铁车站结构有限元模型,以两种基岩波的水平向和竖向地震动作为激励,求解各模型中地铁车站结构重要部位的地震反应。分析不同埋深时地铁车站结构惯性作用、侧面土体和上部土体三个因素对地铁车站地震反应的影响情况。分析结果表明:在双向地震作用下,地铁车站侧壁弯矩、剪力、轴力和中柱轴力随埋深的增加而增加,中柱剪力和弯矩随埋深增加而减少。埋深越深,侧面土体对地铁车站地震反应影响越大;上部土体使中柱轴力不断增加;结构自身的惯性作用对其地震反应的贡献逐渐减小。     

Transient mechanical wave propagation in semi-infinite porous media using a finite element approach

In this study, we propose a numerical investigation in the time domain of the mechanical wave propagation of an impulsional load on semi-infinite soil. The ground is modelled as a porous saturated viscoelastic medium involving complete Biot theory. All the couplings and a hysteretic Rayleigh damping are taken into consideration. An accurate and efficient finite element method using a matrix-free technique and an expert multigrid system are applied. Our results present the displacements of the fluid and solid particles over the surface and in depth. The arrival times of body and surface waves are studied. Particularly, the compressional wave of the second kind is highlighted. The influence of the different couplings and more specifically, the influence of the permeability on the response of the soil are analyzed.     

Thin layer element method for dynamic soil-structure interaction analysis of axi-symmetric structure in submerged soil

The three-dimensional thin layer element method is formulated for the dynamic response analysis of an axi-symmetric structure in submerged soil. Biot's wave equation for fluid-filled porous medium is used in the formulation. The three-dimensional thin layer element method computes the wave numbers and their associated mode shapes, for both Rayleigh waves and Love waves in submerged soil, which define the characteristics of the waves. The submerged condition affects the characteristics of the Rayleigh waves in soil. As a result, it alters substantially the soil-structure interaction stresses if the permeability of the soil is relatively large and, to less extent, the response of the structure. The thin layer element method is far more efficient than the finite element method for analyzing the fluid-filled porous medium, yet capable of taking into account a multi-layered inhomogeneous soil.     

Modeling shear rigidity of stratified bedrock in site response analysis

Where a distinct soil-rock interface exists, the bedrock medium is commonly treated as elastic half-space and the bedrock surface as the lower boundary of the soil-column model for site response analyses (or the lower boundary of the finite element model for soil-structure interaction analyses). While shear wave velocity in bedrock varies with depth, there has been no consensus amongst scientists and practitioners over the value of “effective depth” into bedrock at which the “half-space” shear wave velocity value should be taken for modeling purposes. This paper reports an interesting and important observation that the effective depth into bedrock is sensitive to the shear wave velocity profile of the overlying soil sediments. A simple and heuristic method, namely Resonant Period Equivalence (RPE) Method, is proposed herein for representing a stratified elastic bedrock of inhomogeneous properties by an equivalent homogeneous elastic half-space medium, which is characterized by a single equivalent shear wave velocity (V_R) value. The proposed calculation method has been verified by extensive comparative analyses involving the use of programs SHAKE and NERA and employing the complete shear wave velocity models of both the soil sediments and the underlying stratified bedrock.     

基于流固耦合统一计算框架的三维大规模海域场地地震反应分析?以东京湾为例

海域场地地震响应分析是确定海洋工程结构抗震设计地震动输入的重要环节。然而,针对海水、饱和土、基岩之间的流固耦合分析,目前一般通过对3种介质方程进行离散,然后整体求解或分区耦合求解的方式进行,过程复杂而低效。因此,大规模海域场地地震反应分析仍是一个挑战性问题。本文基于流固耦合统一计算框架求解海域近场波动问题,采用透射边界模拟无限域,通过将海水和基岩视为孔隙率分别等于1和0的广义饱和多孔介质,使得海水、饱和土、基岩之间的相互耦合可在统一计算框架中实现,避免不同介质求解器之间的数据交换。采用集中质量显式有限元并行计算,不同进程之间采用MPI进行数据交换,提高计算效率;采用逐元技术,按单元类别存储单元刚度,大大节省了内存,便于大规模计算。通过自编程,输入界面高程数据和材料参数,实现建模-自由场-三维地震动模拟全流程自动化。以东京湾为例,使用该方法和程序在超级计算机上模拟SV波垂直入射时的地震响应,证实了该方法用于三维大规模海域地震波场模拟的高效性和可行性。     

Earthquake analysis with generalized plasticity model for saturated soils

A plasticity based constitutive model for anisotropic behaviour of soils is implemented in a finite element procedure based on the generalized Biot theory for the dynamic non-linear response of porous materials. The model represents a version in the hierarchical approach of constitutive modelling and allows for inelastic response during loading, unloading and reloading. The procedure has been verified previously with respect to closed-form solutions for wave propagation in porous media. In this paper, it is used to predict the behaviour of a realistic structure-saturated porous soil system subjected to earthquake loading. Both linear and non-linear analyses have been performed. It has been found that the predicted responses from the two analyses are significantly different; for example, in comparison with the linear analysis the non-linear response shows increased magnitudes and zones of concentration of pore water pressures, increased magnitudes of horizontal displacements, decreased magnitudes of vertical displacements and increased magnitudes of shear stresses.     

Response of a randomly inhomogeneous layered media to harmonic excitations

We present a unified formulation of an analytical method for evaluating the response of a random soil medium to surface or earthquake excitations. Specifically, we are interested in this study to the case of a horizontally stratified layered soil profile. Soil properties, mass density and shear modulus of each layer are modeled as spatial random fields. The soil profile is laying on a homogeneous half-space. Integro-differential equations are formulated and solved using the Laplace transform method. Numerical results are, firstly, obtained in terms of effects of soil properties on distributions of mean dimensionless displacement and stress with depth. Then, the amplification function of the surface layer of a soil profile is formulated and a parametric study is conducted to examine the effects of stochastic variations of mass density and shear modulus on the amplification function. We, lastly, computed the seismic response of a site located in the basin of Metidja in Algeria to base rock accelerations which were recorded at Keddara's station during the May 21, 2003 Boumerdes earthquake, in terms of free surface Fourier spectrum amplitudes and accelerations. Comparison between computed accelerations and recorded ones at Dar Elbeida site has proved the effectiveness of the approach.     

土体参数对地表加速度峰值和反应谱的影响

运用一维土层地震反应的等效线性化波动方法,研究了单层均质土密度、动剪切模量比和阻尼比、土层厚度和土层剪切波速等参数变化对土层地表加速度峰值和反应谱的影响,所得到的结论对于指导地震安全性评价中的勘察试验工作、提高土层地震反应的计算精度具有一定的参考价值。     

Dynamic responses of coastal structures during earthquakes including sediment–sea–structure interaction

The effects of the interaction among sea water, sediment, backfill-soil and coastal structures (embankments) were included in the present study. The formulation is derived from fundamental theories in various fields, including marine hydrodynamics, flow in porous medium, and structural dynamics. The hybrid finite-difference and finite element methods were used in the analysis. The finite-difference method was used to calculate the nonlinear hydrodynamic pressures of sea water as well as the pore water in the sediment acting on the coastal embankment faces by seismic-wave actions. The fluid-filled solid mixture was used to model sediment and back-fill soil and the corresponding dynamic responses were also evaluated by finite difference method. The dynamic response of the coastal structures was calculated by finite element method. The numerical results are presented for various water depths and ground motion intensities. The significant dynamic forces on coastal structures were calculated during earthquakes and the possible sliding of the coastal embankment will occur and the special foundation treatment should be made.