波浪作用下海床动力反应的数值分析

近海和离岸建筑物的基通常处于连续不断的小风浪作用之下，可将其变形视为弹性。基于二维广义Biot理论，提出了线性或非线性波浪作用下饱和弹性海床动力应应的时域有限元数值解法。静力平衡条件和Biot方程组成的边值方程可视为其特例，在比较算例中，数值计算得到的孔压和有效应力幅值沿海床深度的分布与解析解十分吻合。土骨架和孔隙流体的加速度对海床动力反应的影响很小。具体算例表明，线性波沿缓坡海床传播时，土层中超静孔压和有效应力幅值随之增大，有可能发生滑动坡坏。所提出的数值解法能够灵活地处理非线性波浪荷载，海底复杂地质条件和波浪沿缓坡传播等复杂情况。     

波浪作用下孔隙海床-管线动力相互作用分析

波浪作用下海床中的孔隙水压力与有效应力是影响海底管线稳定性的主要因素。然而,在目前的海床响应分析中一般将管线假定为刚性,并不能合理地考虑海床与管线之间的相互作用效应,同时也没有考虑土体和管线加速度对海床动力响应的惯性影响,从而无法确定由此所引起的管线内应力。为此考虑管线的柔性,分别采用饱和孔隙介质的Biot动力固结理论和弹性动力学理论列出了海床与管线的控制方程,进而采用摩擦接触理论考虑海床与管线之间的相互作用效应,基于有限元方法建立了海床-管线相互作用的计算模型及其数值算法。通过变动参数对比计算讨论了管线几何尺寸、海床土性参数对波浪所引起的管线周围海床孔隙水压力和管线内应力的影响。     

随机波作用下埋管海床动态响应及液化研究

基于广义Biot动力理论和Longuet-Higgins线性叠加模型,构建波浪-海床-管线动态响应的有限元计算模型,求解随机波作用下,多层砂质海床中管线周围土体孔隙水压力和竖向有效应力的分布。采用基于超静孔隙水压力的液化判断准则,得出液化区的最大深度及横向范围,从而判断海床土体液化情况。考虑海洋波浪的随机性,将海床视为多孔介质,海床动态响应计算模型采用u-p模式,孔隙水压力和位移视为场变量。并考虑孔隙水的可压缩性、海床弹性变形、土体速度、土体加速度以及流体速度的影响,忽略孔隙流体惯性作用。参数研究表明:土体渗透系数、饱和度以及有效波高等参数对海床土体孔隙水压力、竖向有效应力和液化区域分布有显著影响。     

波流共同作用下海底人工边坡动态响应分析

为了研究波流共同荷载作用下开挖基槽附近海床动态响应和液化破坏情况,提出一个二维耦合计算模型,采用雷诺时均纳维-斯托克斯（RANS）方程描述波浪运动情况,通过设定侧边界条件实现稳定流场。海床部分通过求解Biot固结方程,得到波流荷载下海床中的应力和位移情况。将模型计算结果与水槽试验数据和解析解进行比较,验证了波流模型和海床模型的有效性。在此模型基础上,分析得到了开挖之后海床新的应力和固结状态。同时,通过参数分析得到了波流耦合情况下波浪形态的变化,以及海流对海床液化情况和孔压情况的影响。最后,通过线性回归计算得到最大液化深度与流速的拟合关系曲线。计算结果可用于判断基槽开挖后不规则海床的液化情况,对相关研究和实际工程具有一定参考意义。     

波流作用下黄河三角洲多层粉质土孔压累积及液化分析

在波浪和海流作用下,多层粉质土海床可能会发生液化破坏,进而影响工程构筑物的安全.以黄河三角洲为背景,选取垦东地区三种典型粉质土海床,基于Terzaghi固结方程和Biot渗流理论,验证了波流作用下多层土海床累积孔压响应模型,探讨了三种不同典型海床累积孔压与液化特征的差异,以及波流参数对海床动力响应的影响.结果表明,上覆...     

浅水区波浪非线性效应对砂质海床动力响应的影响

以广义Biot动力固结理论为基础,运用一阶椭圆余弦波和二阶Stokes波等非线性波浪理论考虑浅水区波浪荷载的非线性效应,在时域上采用有限元方法对非线性波浪力作用下饱和砂质海床的动力响应进行了数值求解,并与线性波浪作用下海床动力响应特性进行了对比分析。结果表明,随着波长与水深之比L／d及无量纲参数T(g/d)^1/2的增大,非线性波浪对海床动力响应的影响增大。与线性波浪理论相比,孔隙水压力与有效应力幅值的增大效应非常显著。因此在近海海洋建筑物设计与工程场地评价中,波浪力的非线性特性必须引起注意。     

浅水非线性波作用下沙质底床孔隙水压力响应数值分析

近岸水深较浅,波浪具有较强的非线性,海床破坏与波浪作用下孔隙水压力的分布有着密切的关系。波浪场控制方程采用雷诺时均方程和k-ε紊流模型,入射波采用椭圆余弦波,采用PLIC-VOF法追踪自由表面;海床域以Biot动力固结理论为基础,建立了非线性波浪与海床相互作用的弱耦合数学模型,获得椭圆余弦波作用下沙质海床中孔隙水压力响应规律。计算结果表明,与线性波浪相比,浅水非线性波作用下沙质海床中孔隙水压力幅值增大非常显著。     

非线性海床土对钢悬链式立管触地点动力响应和疲劳损伤影响分析

基于钢悬链式立管(SCR)动力分析程序CABLE3D,采用大挠度柔性梁理论建立SCR的运动方程,将线性海床模型扩展为考虑海床土吸力的非线性海床模型,采用非线性有限元方法对控制方程进行离散,时域内积分采用Newmark-β法,开发出新的计算程序。通过算例分析上部浮体垂荡运动幅值、海床土剪切强度、海床土剪切强度梯度对SCR触地点区域动力响应和疲劳损伤的影响。分析结果表明:SCR触地点区域动力响应和疲劳损伤对上部浮体垂荡运动幅值和海床土剪切强度的变化较为敏感,疲劳损伤在触地点区域最大,远大于悬垂段和流线段,在设计过程中应采取一定的加强措施。     

随机波作用下海床动态响应分析

真实的海洋波浪是随机的,而前人对海床的动态响应分析大都是选用线性波或者Stokes波理论,对海床的模拟大都采用Biot拟静力模型,忽略了流体速度及土体位移加速度的影响。联合使用Longuet-Higgins随机波模型(采用Jonswap谱)以及动力u-p形式的海床响应计算模型,使用COMSOL Multiphysics多场耦合软件的PDE模块输入方程进行有限元计算,得到随机波作用下整体海床动态响应结果。将随机波结果与一阶Stokes波和椭圆余弦波结果进行对比,并对渗透系数和饱和度进行参数分析,研究表明渗透系数和饱和度对于随机波作用下海床动态响应影响显著。     

波浪作用下单桩基础周围海床液化机制研究

建立波浪作用下单桩周围三维海床动力响应模型,考虑自重影响下的海床长时间固结过程。采用已有物理模型试验数据对模型进行验证,证实其具有较好的适用性。模拟波浪作用下单桩周围三维海床液化区域,通过定量分析超孔隙水压力和土体初始有效应力的变化,讨论单桩插入深度对海床液化的影响机制。研究表明,单桩插入深度发生变化时,土体初始有效应力对海床液化的影响要大于超孔隙水压力,且影响程度随着插入深度的增加而逐渐增大。     

Spatial and temporal variation process of seabed dynamic response induced by the internal solitary wave

Internal solitary wave(ISW) is often accompanied by huge energy transport, which will change the pore water pressure in the seabed. Based on the two-dimensional Biot consolidation theory, the excess pore water pressure in seabed was simulated, and the spatiotemporal distribution characteristics of excess pore water pressure was studied. As the parameters of both ISW and seabed can affect the excess pore water pressure, the distribution of pore water pressure showed both dissipation and phase lag...     

Numerical simulation of the nonlinear wave-induced dynamic response of anisotropic poro-elastoplastic seabed

Abstract

In this paper, a 2D poro-elastoplastic model for wave-induced dynamic response in an anisotropic seabed is derived analytically. The seabed is treated as a porous medium and characterized by Biot’s consolidation equations. The soil plasticity and wave non-linearity are included in the model and both the pore fluid and the soil skeleton are assumed to be compressible. The nonlinear ocean waves are respectively considered as progressive and standing waves. The previous experimental data is used to validate the proposed model. Numerical results demonstrate that the influence of nonlinear wave components should not be ignored without committing substantial error. A significant difference between progressive and standing waves is also observed for the development of residual pore pressure, as well as the distribution of liquefied zone. A detailed parametric investigation reveals that the nonlinear wave-induced seabed response is also affected significantly by cross-anisotropic soil parameters.     

Numerical study for waves propagating over a porous seabed around a submerged permeable breakwater: PORO-WSSI II model

The phenomenon of the wave, seabed and structure interactions has attracted great attentions from coastal geotechnical engineers in recent years. Most previous investigations have based on individual approaches, which focused on either flow region or seabed domain. In this study, an integrated model (PORO-WSSI II), based on the Volume-Averaged/Reynolds-Averaged Navier-Stokes (VARANS) equations and Biot's poro-elastic theory, is developed to investigate the mechanism of the wave-permeable structure-porous seabed interactions. The new model is verified with the previous experimental data. Based on the present model, parametric studies have been carried out to investigate the influences of wave, soil and structure parameters on the wave-induced pore pressure. Numerical results indicated: (i) longer wave period and larger wave height will obviously induce a higher magnitude of pore pressure at the leading edge of a breakwater; (ii) after a full wave-structure interaction, the magnitude of pore pressure below the lee side of a breakwater decreases with an increasing structure porosity while it varies dramatically with a change of structure height; and (iii) the seabed thickness, soil permeability and the degree of saturation can also significantly affect the dynamic soil behaviour.     

Parametric study of the wave-induced residual liquefaction around an embedded pipeline

Wave-induced liquefaction in a porous seabed around submarine pipeline may cause catastrophic consequences such as large horizontal displacements of pipelines on the seabed, sinking or floatation of buried pipelines. Most previous studies in relation to the wave and seabed interactions with embedded pipeline dealt with the wave-induced instaneous seabed response and possible resulting momentary liquefaction (where the soil is liquefied instantaneously during the passage of a wave trough), using theory of poro-elasticity. Studies for the interactions between a buried pipeline and a soil undergoing build-up of pore pressure and residual liquefaction have been comparatively rare. In this paper, this complicated process was investigated by using a new developed integrated numerical model with RANS (Reynolds averaged Navier–Stokes) equations used for governing the incompressible flow in the wave field and Biot consolidation equations used for linking the solid–pore fluid interactions in a porous seabed with embedded pipeline. Regarding the wave-induced residual soil response, a two-dimensional poro-elastoplastic solution with the new definition of the source term was developed, where the pre-consolidation analysis of seabed foundation under gravitational forces including the body forces of a pipeline was incorporated. The proposed numerical model was verified with laboratory experiment to demonstrate its accuracy and effectiveness. The numerical results indicate that residual liquefaction is more likely to occur in the vicinity of the pipeline compared to that in the far-field. The inclusion of body forces of a pipeline in the pre-consolidation analysis of seabed foundation significantly affects the potential for residual liquefaction in the vicinity of the pipeline, especially for a shallow-embedded case. Parametric studies reveal that the gradients of maximum liquefaction depth with various wave and soil characteristics become steeper as pipeline burial depth decreases.     

Stability and liquefaction analysis of porous seabed subjected to cnoidal wave

Cnoidal wave theory is appropriate to periodic wave progressing in water whose depth is less than 1/10 wavelength. However, the cnoidal wave theory has not been widely applied in practical engineering because the formula for wave profile involves Jacobian elliptic function. In this paper, a cnoidal wave-seabed system is modeled and discussed in detail. The seabed is treated as porous medium and characterized by Biot's partly dynamic equations (u–p model). A simple and useful calculating technique for Jacobian elliptic function is presented. Upon specification of water depth, wave height and wave period, Taylor's expression and precise integration method are used to estimate Jacobian elliptic function and cnoidal wave pressure. Based on the numerical results, the effects of cnoidal wave and seabed characteristics, such as water depth, wave height, wave period, permeability, elastic modulus, and degree of saturation, on the cnoidal wave-induced excess pore pressure and liquefaction phenomenon are studied.