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

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

波浪作用下海床孔压累积过程离散元数值模拟

海床土层在波浪的循环荷载作用下会逐渐累积孔压,降低土层的稳定性,并威胁海上工程。为了研究孔隙水压力的累积机制,本文提出离散元孔隙密度流方法,并改进研发离散元分析软件MatDEM,实现了海床沉积物孔压的累积过程模拟。基于现场试验装置及土体力学参数建立离散元模型,通过对比试验和数值模拟结果发现:对海床沉积物施加波浪荷载后,表层土体中产生较高孔压,并逐渐向深层传递;在循环波浪荷载作用下,土颗粒间孔压累积范围逐渐增加;当孔压累积时间足够长时,土层中孔压收敛于所施加最大荷载与最小荷载的平均值,此时若孔压达到初始有效应力,土体将发生液化,内部土颗粒成为再悬浮沉积物;在周期性波浪荷载作用下,土颗粒液化悬浮后发生移动,浅层颗粒位移量大,土体整体表现为圆弧形移动。     

波浪作用下沙床中的孔隙水压力响应模型实验研究

本文介绍了在波浪水槽中进行的波浪作用下沙床中的孔隙水压力量测实验，并给出和分析了部分实验数据，实验中的波浪采用和规则波和不规则波。沙床采用了均质沙床和分层沙床，同时对沙床上有直立墙的模型也进行了实验。     

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

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

波浪作用下斜坡上单桩周围海床孔隙水压力响应特性试验研究

单桩基础周围斜坡海床中的波致孔隙水压力响应与纯斜坡海床存在较大差异。为了解不同波高、波周期条件下,单桩基础周围波浪传播变形及其对斜坡海床孔压振荡响应的影响,在波浪水槽末端铺设了长6 m、坡度1∶16的斜坡砂床进行试验。通过改变桩身位置和波浪参数,测量斜坡段各处波面形态,采集单桩周围孔隙水压力,分析了桩身位置及波浪参数对斜坡海床孔压响应的影响。结果表明：相同入射波条件下,随距坡脚水平距离增加,波高、近底流速和桩周孔隙水压力幅值都随之增大;桩周孔隙水压力幅值分布规律为：桩前孔压幅值明显大于桩侧与桩后孔压幅值。当Keulegan-Carpenter数大于6时,随着波高和波周期增大,马蹄涡产生的负压区使得桩侧海床孔隙水压力与纯斜坡海床孔隙水压力差值迅速增加。     

波浪作用下粉土中的孔压响应及其在粉土海床稳定性评价中的应用

粉土在波浪等动荷载作用下极易发生液化破坏,而孔隙水压力在粉土动力学行为中扮演了一个很重要的角色,其发展变化会直接影响到土体的稳定性。因此,通过室内波浪作用下的粉土孔压响应模型试验探讨了孔压与波浪之间的响应情况,发现波浪能量的影响沿土层深度递减,水深条件相同时,响应的孔压随波高的增大而增加,当波浪作用足够长时间后粉土发生液化破坏,此时粉土内累积的孔压小于上覆土体的自重应力。根据结果提出了1种评价粉土海床稳定性的方法。     

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

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

波致液化粉质土有效应力变化的试验数据分析

海床不稳定性的现象很多都是海底液化引起的。根据之前的研究,粉土在液化之后有效应力会依然存在,关于波浪作用下海床液化之后有效应力变化的研究很少。采用波浪水槽实验,在未液化和液化2种情况下,分别施加不同波高的波浪,对底床各层位土体孔隙水压力进行采集,并对比研究。实验结果表明,液化后有效应力的相对值相对于液化前有很明显的减小,并且在深度上表现出显著的差异性,这种差异性随着波高的增大而减小。当相同深度处同一种波高作用一定时间时,有效应力会出现极大值,然后有效应力会减小。     

波浪诱发松散海床渐进式液化的数值分析

近海区域广泛分布着第四纪新沉积的松散海洋土,波浪荷载作用下松散海床会发生液化进而对近海结构物的稳定性存在巨大威胁。本文采用中国科学院流体-结构-海床相互作用数值计算模型FSSI-CAS 2D,选用Pastor-Zienkiewicz-Mark Ⅲ（PZⅢ）弹塑性本构研究了波浪诱发的松散海床液化问题。分析了波浪荷载引起的松散海床内超孔隙水压力、有效应力以及应力角的时程变化特性,并预测了松散海床的渐进液化过程。计算结果表明,波浪荷载作用下松散海床内残余孔压会累积增长,海床表面最先发生液化,然后逐渐向下发展至液化最大深度。同时指出海床内超孔隙水压力的竖向分布特征和应力角的变化时程均可以作为判断海床液化的间接参数。最后,通过应力状态分析,讨论了海床渐进式液化的发展过程和趋势。     

非线性波浪作用下埋置管道上波浪力简化计算

由于浅水区波浪的非线性影响显著,浅埋管道受非线性波浪荷载的影响大,为了保证管道长期运行的稳定性,在管道设计过程中需要充分考虑由非线性波浪引起的波浪力。考虑孔隙水和海床土的压缩性,基于Biot固结理论和一阶近似椭圆余弦波理论,利用分离变量法推导了非线性波浪作用下浅水区埋置管道周围海床的渗流压力,进而给出了埋置管道上的波浪力压力解析解,并与已有的文献结果进行比较。计算结果表明,在椭圆余弦波的作用下,海底管道周围海床内的渗流压力呈正弦分布,且管道所受的波浪力随着管径的增大而增大。     

The numerical study of wave-induced pore water pressure response in highly permeable seabed

The coupling numerical model of wave interaction with porous medium is used to study waveinduced pore water pressure in high permeability seabed.In the model,the wave field solver is based on the two dimensional Reynolds-averaged Navier-Stokes(RANS) equations with a k-ε closure,and Forchheimer equations are adopted for flow within the porous media.By introducing a Velocity-Pressure Correction equation for the wave flow and porous flow,a highly efficient coupling between the two flows is implemented.The numerical tests are conducted to study the effects of seabed thickness,porosity,particle size and intrinsic permeability coefficient on regular wave and solitary wave-induced pore water pressure response.The results indicate that,as compared with regular wave-induced,solitary wave-induced pore water pressure has larger values and stronger action on seabed with different parameters.The results also clearly show the flow characteristics of pore water flow within seabed and water wave flow on seabed.The maximum pore water flow velocities within seabed under solitary wave action are higher than those under regular wave action.     

An integrated three-dimensional model of wave-induced pore pressure and effective stresses in a porous seabed: II. Breaking waves

The evaluation of the wave-induced liquefaction potential is particularly important for coastal engineers involved in the design of marine structures. Most previous investigations of the wave-induced liquefaction have been limited to two-dimensional non-breaking waves. In this paper, the integrated three-dimensional poro-elastic model for the wave-seabed interaction proposed by [Zhang, H., Jeng, D.-S., 2005. An integrated three-dimensional model of wave-induced pore pressure and effective stresses in a porous seabed: I. A sloping seabed. Ocean Engineering 32(5/6), 701–729.] is further extended to simulate the seabed liquefaction potential with breaking wave loading. Based on the parametric study, we conclude: (1) the liquefaction depth due to breaking waves is smaller than that of due to non-breaking waves; (2) the degree of saturation significantly affects the wave-induced liquefaction depth, and no liquefaction occurs in full saturated seabed, and (3) soil permeability does not only significantly affect the pore pressure, but also the shear stresses distribution.     

Effects of dynamic soil behavior and wave non-linearity on the wave-induced pore pressure and effective stresses in porous seabed

Most previous investigations for the wave-induced soil response have only considered the quasi-static soil behavior under linear wave loading. However, it is expected that the dynamic soil behavior and wave non-linearity will play an important role in the evaluation of wave-induced seabed response. In this paper, we include dynamic soil behavior and wave non-linearity into new analytical models. Based on the analytical solution derived, the effects of wave non-linearity on the wave-induced seabed response with dynamic soil behavior are examined. Numerical results demonstrate the significant effects of wave non-linearity and dynamic soil behavior on the wave-induced effective stresses. The applicable range of dynamic and quasi-static approximations is also clarified for engineering practice.     

Analysis of wave-induced liquefaction in seabed deposits of silt

The dynamic stress introduced in half elastic space by wave loading is characterized by the equation between the magnitude of half cyclic axial stress and cyclic torsion shear stress and the principal stress,whose direction rotates continuously and compression stress on seabed can be calculated by the use of small amplitude wave theory.With relationship curves of saturated silt of liquefaction cycles and cyclic stress ratios obtained by cyclic triaxial-torsional coupling shear tests and curve fitting method to different data points of relative density,it is suggested that the cyclic stress ratio corresponding to constant liquefaction impedance be taken as the critical cyclic stress ratio which implies liquefaction.There exists a linear relationship between critical cyclic stress ratio and relative density under different relative densities.Empirical formula for critical cyclic stress ratios of seabed liquefaction induced by wave loading under different relative densities is established.The possibility of seabed silt liquefaction and its influence factors are analyzed based on the small-amplitude wave theory and the data acquired in laboratory tests.     

不同破碎波对沙质海床作用的实验研究

破碎波对近海海岸地形以及海岸建筑物影响强烈,通过物理模型实验对孤立波、规则波作用下破碎带的床面形态以及孔隙水压力进行分析。破碎波冲击海床,破碎处床面上形成沙坝和沙坑,与规则波相比,孤立波破碎时对床面的冲刷更加剧烈,床面形成的沙坝和沙坑尺度更大,且土体内孔隙水压力幅值也较大。同时研究了波面变化对孔隙水压力的影响,发现波面变化历时曲线与孔隙水压力历时曲线相似,与孔隙水压力梯度历时曲线更为相似,说明波面变化更能反映海床内部孔隙水压力梯度的变化。通过探讨波浪与海床之间相互耦合作用,发现破碎带地形变化使得波浪出现不同破碎类型,分析得出卷破波比崩破波作用下孔隙水压力幅值大。     

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

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

An experimental study on non-linear progressive wave-induced dynamic stresses in seabed

In this study, the dynamic stresses within the seabed induced by non-linear progressive waves were explored through a series of hydraulic model tests on a movable bed within a wave flume. By comparing Stokes’ 2nd-order wave theory with the theory of wave-induced dynamic stresses within the seabed as proposed by Yamamoto et al. [1978. On the response of a poro-elastic bed to water waves. Journal of Fluid Mechanics 87 (1), 193–206.] and Hsu and Jeng [1994. Wane-induced soil response in an unsaturated anisotropic seabed of finite thickness. International Journal for Numerical and Analytical Methods in Geomechanics 18, 785–807], the experimental results show that the pressure on the seabed surface, the pore water pressure within the seabed as well as the vertical and the horizontal stresses are all smaller than their theoretical values. If we were to obtain the characteristics of seabed soil, the analytical solution of Hsu and Jeng [1994. Wane-induced soil response in an unsaturated anisotropic seabed of finite thickness. International Journal for Numerical and Analytical Methods in Geomechanics 18, 785–807] might agree to the simulation of the wave-induced effective stresses and shear stress in the sandy seabed. A different phase shift exists among all the three soil stresses. Their influences on the three dynamic stresses within seabed soil are important for seabed stability, and can be used in the verification of numerical models. In the whole, the non-linear progressive waves and the naturally deposited seabed are found to have a strong interaction, and the behavior of the induced dynamic stresses within the seabed is very complicated, and should be investigated integrally.     

波致海床滑动稳定性计算方法及滑动失稳特征研究

波浪引起的海床不稳定性是海洋工程中需要考虑的重要问题。在对现有波致海床滑动稳定性计算方法进行分析的基础上,提出了一种波致海床滑动稳定性计算的全应力状态法,将其与现有计算方法进行了对比分析,并进一步研究了波致砂土海床和软土海床的滑动失稳特征。结果分析表明,全应力状态法在波致海床滑动稳定性分析中具有较好的适用性。对于砂土海床,其滑动稳定性受饱和度的影响较大,且当海床计算厚度约为0.2倍波长时对应的滑动深度最大。波浪作用下坡度不超过2°的均质软土海床,其最危险滑动面的位置仅与波长有关,其滑动深度约为0.21倍波长,滑动面半弦长约为0.33倍波长;海床表面的波压力数值只影响其安全系数的大小,而不影响其滑动深度。