波浪与新型双排开孔圆筒防波堤相互作用三维数值模拟

双排开孔圆筒防波堤是基于圆筒、板式结构的一种复合式新型结构型式;基于不可压缩两相流模型建立三维数值波浪水槽,通过RNG k-ε湍流模型进行湍流封闭,并采用TruVOF方法捕捉自由液面,开展波浪与双排开孔圆筒防波堤相互作用数值模拟,探究相对排间距、开孔率对新型双排开孔圆筒防波堤消浪性能的影响,分析了后排开孔圆筒防波堤附近的复杂水动力现象和流动特性.结果表明,在本文研究工况范围内,沿程平均波高随相对排间距的增大先增大后减小,随开孔率的增大而增大,周期对沿程平均波高的影响没有明显规律;当B/D=9, e=23.11%时,新型双排开孔圆筒防波堤消浪效果最优,反射系数在0.4～0.46之间,透射系数在0.3～0.35之间,耗散系数在0.8～0.85之间;自由液面破碎、水气掺混、环状涡运动演化是新型双排开孔圆筒防波堤紊动耗能消波的主要原因;相对排间距会引起后排防波堤附近涡量分布以及剪切层形态的变化,从而导致不同的紊动特性,影响双排开孔圆筒防波堤消浪特性.研究结果可以为新型双排开孔圆筒防波堤工程设计与消浪机理研究提供理论支撑.     

多孔材料中声波的传播与演化

采用两相多孔介质的拉格朗日模型来描述一种理论流体充填的多孔弹性固体材料,其中孔隙度的变化满足一个附加的平衡方程。     

Diffraction of waves by semi-infinite breakwater using finite and infinite elements

A finite and infinite element model is derived to predict wave patterns around a semi-infinite breakwater in water of constant depth. Both circular and square meshes of elements are used. The wave theory used is that of Berkhoff. The appropriate boundary conditions for finite and infinite boundaries are described. The singularity in the velocity at the breakwater tip is modelled effectively using the technique of Henshell and Shaw originally developed in elasticity. The results agree well with the analytical solution. In addition the problem of waves incident upon a semi-infinite breakwater and parabolic shoal, where both diffraction and refraction are present, is solved. There is no analytical solution for this case. The combination of finite and infinite elements is found to be an effective and accurate technique for such problems.     

基于卷积神经网络的涵洞式直立堤波浪透射预测

涵洞式直立堤是一种具有特殊用途的海岸工程结构物,对其透浪特性的研究具有重要工程意义.然而,目前众多学者对于涵洞式直立堤波浪透射问题的研究主要以理论分析、实验模拟及数值计算为主.随着机器学习技术的发展,传统水动力学问题迎来了新的求解理念.机器学习算法可根据训练数据集自主学习相应的规律,以数据映射的方式建立水动力学特征预测...     

A numerical model for the isolation of moving-load induced vibrations by pile rows embedded in layered porous media

A numerical model is developed to analyse the isolation of moving-load induced vibrations using pile rows embedded in a layered poroelastic half-space. Based on Biot’s theory and the transmission and reflection matrices (TRM) method, the free wave field solution for a moving load applied on the surface of a layered poroelastic half-space and the fundamental solution for an harmonic circular patch load are determined. Using Muki and Sternberg’s method, the second kind of frequency domain Fredholm integral equations for the dynamic interaction between pile rows and the layered poroelastic half-space are derived. The time domain solution is recovered via inverse Fourier transform in order to obtain the amplitude reduction ratio and thus assess the vibration isolation efficiency of pile rows. A special case of the present model shows good agreement with an existing solution. Numerical results of this study show that the speed of moving loads has an important influence on the isolation of vibrations by pile rows: the same pile rows can achieve better isolation efficiencies for higher speed loads than for lower speed loads. Pile rows embedded in a two-layered poroelastic half-space with a softer overlying layer usually generate better vibration isolation effects than those with a stiffer overlying layer. Finally, better isolation vibration may be realized by increasing the pile length and decreasing the net spacing between neighboring piles in a pile row.     

Poromechanical response of a finite elastic cylinder under axisymmetric loading

Drained or undrained cylindrical specimens under axisymmetric loading are commonly used in laboratory testing of soils and rocks. Poroelastic cylindrical elements are also encountered in applications related to bioengineering and advanced materials. This paper presents an analytical solution for an axisymmetrically-loaded solid poroelastic cylinder of finite length with permeable (drained) or impermeable (undrained) hydraulic boundary conditions. The general solutions are derived by first applying Laplace transforms with respect to the time and then solving the resulting governing equations in terms of Fourier–Bessel series, which involve trigonometric and hyperbolic functions with respect to the z-coordinate and Bessel functions with respect to the r-coordinate. Several time-dependent boundary-value problems are solved to demonstrate the application of the general solution to practical situations. Accuracy of the numerical solution is confirmed by comparing with the existing solutions for the limiting cases of a finite elastic cylinder and a poroelastic cylinder under plane strain conditions. Selected numerical results are presented for different cylinder aspect ratios, loading and hydraulic boundary conditions to demonstrate the key features of the coupled poroelastic response.     

成层饱和介质平面波斜入射问题的一维化时域方法

地震波斜入射下自由场的输入是大型结构抗震分析中亟待解决的问题之一,尤其是成层饱和多孔介质自由场问题,由于问题的复杂性,目前研究甚少. 本文基于Biot提出的饱和多孔介质动力方程,建立了一种新的求解平面波斜入射下基岩上覆饱和多孔介质成层场地自由场分析的一维化时域计算方法. 该方法首先根据Snell定律将饱和多孔介质二维空间问题转化为一维时域问题,通过对深度方向的有限元离散,得到饱和多孔介质波动问题的一维化有限元方程,然后采用单相弹性介质精确人工边界条件模拟基岩半空间的波动辐射和输入特征,通过考虑基岩与饱和多孔介质间透水或不透水边界条件以及不同饱和多孔介质交界面边界条件,形成基岩上覆成层饱和介质系统的整体有限元方程,最后采用中心差分法与Newmark平均加速度近似格式相结合的方法对时间进行离散,得到节点的动力时程的显式表达. 典型场地的地震反应分析表明,本文方法的计算结果与传递矩阵法结合傅里叶变换的计算结果完全吻合,证明了其有效性.      

Fluid pressure response in poroelastic materials subjected to cyclic loading

When cyclic loading is applied to poroelastic materials, a transient stage of interstitial fluid pressure occurs, preceding a steady state. In each stage, the fluid pressure exhibits a characteristic mechanical behavior. In this study, an analytical solution for fluid pressure in two-dimensional poroelastic materials, which is assumed to be isotropic, under cyclic axial and bending loading is presented, based on poroelasticity. The obtained analytical solution contains transient and steady-state responses. Both of these depend on three dimensionless parameters: the dimensionless stress coefficient; the dimensionless frequency; and, the axial-bending loading ratio. We focus particularly on the transient behavior of interstitial fluid pressure with changes in the dimensionless frequency and the axial-bending loading ratio. The transient properties, such as half-value period and contribution factor, depend largely on the dimensionless frequency and have peak values when its value is about 10. This suggests that, under these conditions, the transient response can significantly affect the mechanical behavior of poroelastic materials.     

Wave interaction with a new type perforated breakwater

In this study examined is the wave interaction with a new modified perforated breakwater, consisting of a perforated front wall, a solid back wall and a wave absorbing chamber between them with a two-layer rock-filled core. The fluid domain is divided into three sub-domains according to the components of the breakwater. Then by means of the matched eigenfunction expansion method, an analytical solution is obtained to assess the hydrodynamic performance of the new structure. An approach based on a step approach method is introduced to solve the complex dispersion equations for water wave motions within two-layer porous media. Numerical results of the present model are compared with previous limiting cases. The effects of rock fill on the reflection coefficient and the horizontal wave force are discussed. The project supported by the Program for Changjiang Scholars and Innovative Research Teams in Universities (IRT0420). The English text was polished by Keren Wang.     

地震波倾斜入射下两相饱和土层响应数值模拟

基于Biot两相饱和多孔介质理论建立了频域内饱和土层的刚度矩阵,对从基岩倾斜入射地震波下饱和土层的波动响应进行了数值研究,入射波采取修正的Kanai-Tajimi谱。通过土层地表面处的水平以及竖直放大系数研究了在不同工况条件下土层的波动响应,研究结果表明,土层的波动响应在入射角接近临界入射角时达到峰值;土层的波动响应随着基岩埋深基以及土层阻尼比的增大显著减小;而土层与基岩的波速比对土层的波动响应的影响存在临界波速比,使土层响应达到极值。对于含有夹层的地质环境,研究了波动响应沿深度的分布,通过与均质土层下响应对比,显示非均质夹层能显著减小土层地表处的波动响应,尤其软夹层能较好地阻止入射波向上的传播。     

Validation of the Slow Compressional Wave in Porous Media: Comparison of Experiments and Numerical Simulations

We perform numerical simulation of ultrasonic experiments on poroelastic samples, in which Biot's slow compressional wave had been observed. The simulation is performed using OASES modeling code, which allows to compute elastic wave fields in layered poroelastic media. Modeled were the experiments of Plona (1980), Rasolofosaon (1988), and our own measurements. In all the three situations, a good agreement between experiment and simulations has been observed. This further confirms the fact that Biot's theory of poroelasticity, on which the simulations were based, adequately describes the behavior of the porous materials under investigations at ultrasonic frequencies.     

Modeling the interaction of solitary waves and semi-circular breakwaters by using unsteady reynolds equations

A vertical 2 -D numerical wave model was developed based on unsteady Reynolds equations. In this model, the k-epsilon models were used to close the Reynolds equations, and volume of fluid( VOF) method was used to reconstruct the free surface. The model was verified by experimental data. Then the model was used to simulate solitary wave interaction with submerged, alternative submerged and emerged semi-circular breakwaters. The process of velocity field, pressure field and the wave surface near the breakwaters was obtained. It is found that when the semi-circular breakwater is submerged, a large vortex will be generated at the bottom of the lee side wall of the breakwater ; when the still water depth is equal to the radius of the semi-circular breakwater, a pair of large vortices will be generated near the shoreward wall of the semi-circular breakwater due to wave impacting, but the velocity near the bottom of the lee side wall of the breakwater is always relatively small. When the semi-circular breakwater is emerged, and solitary wave cannot overtop it, the solitary wave surface will run up and down secondarily during reflecting from the breakwater. It can be further used to estate the diffusing and transportation of the contamination and transportation of suspended sediment.     

Transient response of fluid pressure in a poroelastic material under uniaxial cyclic loading

Poroelasticity is a theory that quantifies the time-dependent mechanical behavior of a fluid-saturated porous medium induced by the interaction between matrix deformation and interstitial fluid flow. Based on this theory, we present an analytical solution of interstitial fluid pressure in poroelastic materials under uniaxial cyclic loading. The solution contains transient and steady-state responses. Both responses depend on two dimensionless parameters: the dimensionless frequency Ω that stands for the ratio of the characteristic time of the fluid pressure relaxation to that of applied forces, and the dimensionless stress coefficient H governing the solid-fluid coupling behavior in poroelastic materials. When the phase shift between the applied cyclic loading and the corresponding fluid pressure evolution in steady-state is pronounced, the transient response is comparable in magnitude to the steady-state one and an increase in the rate of change of fluid pressure is observed immediately after loading. The transient response of fluid pressure may have a significant effect on the mechanical behavior of poroelastic materials in various fields.     

潜堤对波浪在直立式防浪墙上爬高的影响

对直立式防浪墙前潜堤的空间布置对波浪在直墙上爬高的影响进行了研究.建立简化模型,即直立式防浪墙前海床设计为平底,潜堤设计为直立薄板. 改变潜堤高度及其与防浪墙之间距,研究不同的潜堤布置对波浪在直墙上爬高的影响. 针对线性波浪场,利用数值波浪水槽模拟了潜堤作用下直墙上波浪的爬高现象. 同时建立了理论模型,系统地分析了潜堤布置形式对直墙上波浪爬高的影响. 研究结果表明潜堤相对于水深的高度越大,对波浪在直墙上爬高的影响越大,而在一定的相对高度条件下,潜堤与直立墙之间距对波浪爬高的影响呈现出一定的规律性.      

高速荷载下多孔饱和地基的动力响应

研究高速荷载作用下梁与多孔饱和半空间的动力响应。由Fourier变换求解多孔饱和固体的动力基本方程，根据梁与半空间的接触条件得出多孔饱和半空间上梁的垂直位移的表达式。文中的数值算例考虑了荷载移动速度对梁的动力位移的影响，并与相应的弹性半空间问题作了对比。从算例中可以发现荷载移动速度对动力位移有很大的影响，当移动速度与半空间的表面波速相近时，地面会当产生很大的振动，同时还发现当速度大于介质的剪切波速时，多孔饱和半空间上梁的动力响应与弹性半空间上梁的动力响应有很大的差别。     

Dynamic response of an axially loaded rigid sphere embedded in a saturated poroelastic medium

In this paper, an analytical solution for the response of a rigid sphere embedded in a full space poroelastic medium subjected to a dynamic lateral load is derived. The solution is obtained using Biots theory for acoustic waves. In this solution, the displacements of the solid skeleton and the pore pressure are expressed in terms of three scalar potentials. These potentials correspond to the wave velocities of the slow and fast compressional waves and to the shear wave. The governing equation for the dynamic motion is expressed in the frequency domain using Fourier transformation. Different boundary and load conditions were investigated. Curves showing variation in the fluid pressure and solid displacements with the loads frequency were plotted in non-dimensional forms.     

Three-dimensional steady-state closed form solution for multilayered fluid-saturated anisotropic finite media due to surface/internal point source

A three-dimensional(3 D)steady-state solution of fluid saturated anisotropic finite media is presented.The eigenequation method and the pseudo-Stroh formalism are used to obtain the exact solution for homogeneous saturated finite media.The propagator matrix method is introduced to deal with the corresponding multilayered poroelastic media.The poroelastic solutions due to surface or internal point fluid source are obtained.The comparison of the results of the saturated isotropic media in a half space and those obtained by the finite element method is given to illustrate the accuracy of the solution in a finite domain.Numerical solutions of a sandwich poroelastic medium are presented to analyze its hydromechanical behaviors.Two ratios of the horizontal permeability to vertical permeability and different source positions are investigated.The results show that the fluid parameters and source positions have great influence on the hydromechanical behaviors of the layered media.     

Stress and pressure fields around two wellbores in a poroelastic medium

The problem of two circular wellbores of different size in a poroelastic medium is considered in the present work. The constitutive behaviour of the poroelastic medium is assumed to comply with the classical Biot model for isotropic porous materials infiltrated by compressible fluid. The wellbores are assumed infinitely long and the fluid flow is taken stationary, thus making it possible to perform a plane strain analysis. Owing to the geometrical layout of the system, bipolar cylindrical coordinates have been adopted. Three different sets of BCs on the pressure field and on the fluid flux have been considered, founding the corresponding forms of the pressure field. Based on Helmholtz representation, a displacement potential has been introduced, and the corresponding stress field in the poroelastic medium has been assessed. However, such a solution does not satisfy the BCs at the edges of the wells. Then, an auxiliary stress function, which allows accomplishing the BCs, is introduced, leading to the complete solution of the problem. The cases of two coaxial wellbores (eccentric annulus), a single hole bored in a poroelastic half plane and two intersecting holes have been considered also. The proposed approach allows evaluating the pore pressure and the stress and strain fields in the system varying the amplitude of the wells and the physical parameters of the porous material. In particular, the evaluation of the peak values of the stress components around the circular boreholes plays a key role in a variety of engineering contexts, with particular reference to the stability analysis of wellbores and tunnels and failure of vascular vessels in biological tissues.     

Complex variable function method for the scattering of plane waves by an arbitrary hole in a porous medium

Based on the complex variable function method, a new approach for solving the scattering of plane elastic waves by a hole with an arbitrary configuration embedded in an infinite poroelastic medium is developed in the paper. The poroelastic medium is described by Biot's theory. By introducing three potentials, the governing equations for Biot's theory are reduced to three Helmholtz equations for the three potentials. The series solutions of the Helmholtz equations are obtained by the wave function expansion method. Through the conformal mapping method, the arbitrary hole in the physical plane is mapped into a unit circle in the image plane. Integration of the boundary conditions along the unit circle in the image plane yields the algebraic equations for the coefficients of the series solutions. Numerical solution of the resulting algebraic equations yields the displacements, the stresses and the pore pressure for the porous medium. In order to demonstrate the proposed approach, some numerical results are given in the paper.     

Rayleigh Waves on an Exponentially Graded Poroelastic Half Space

In this paper we consider the propagation of seismic waves in isotropic poroelastic half spaces with continuously varying elastic properties, namely with an exponentially decaying depth profile. The present paper shows that the problem leads naturally to a bicubic equation. We obtain explicit inhomogeneous plane wave solutions in an exponential evanescent form with respect to the depth of half space. Further, these solutions are used to solve the boundary value problem of a Rayleigh surface wave and the secular equation is established. The results obtained theoretically are exemplified for numerical data and represented graphically for a representative poroelastic material.