饱和地基上弹性基础的竖向振动特性研究

采用解析的方法研究了饱和地基上受一简谐竖向荷载作用下弹性基础的动力响应．在分析中,首先利用积分变换技术获得了饱和介质基本控制方程的变换解,然后基于基础-半空间完全放松接触、半空间表面完全透水或不透水的假设,建立了该动力混合边值问题的对偶积分方程,并把该对偶积分方程进一步化为易于数值求解的第二类Fredholm积分方程A·D2文末数值算例给出了动力柔度系数、位移和孔隙水压力随振动频域和土-基础体系物理力学参数特性的变化曲线．结果表明:饱和地基上弹性基础的动力响应完全不同于饱和地基上刚性圆板的动力响应．所用方法可用于研究波的传播、土-结构动力相互作用等许多问题．     

伴有排水的两相饱和介质动力问题的LAMB积分公式

饱和土在动力作用下伴有排水的反应是工程中常见的现象.由于饱和土以两相饱和介质模拟,在已经求得的集中力作用下,两相饱和介质基本解Green函数和三项流相Green函数基础上,由Betti定理出发,推出排水问题的Lamb积分公式和应力公式,并代入相关参数,通过计算得到了中心扩张源问题的位移场、排水状态、孔隙压力等解答.最后给出了当两相介质蜕化为单相介质后,与经典解答比较的结果.     

横观各向同性饱和地基的三维动力响应

首先引入位移函数,将直角坐标系下横观各向同性饱和土Biot波动方程转化为2个解耦的六阶和二阶控制方程;然后基于双重Fourier变换,求解了Biot波动方程,得到以土骨架位移和孔隙水压力为基本未知量的积分形式的一般解,并用一般解给出了饱和土总应力分量的表达式．在此基础上系统研究了横观各向同性饱和半空间体的稳态动力响应问题,考虑表面排水和不排水两种情况,得到了半空间体在任意分布的表面谐振荷载作用下,表面位移的稳态动力响应,文末给出了算例．     

流体饱和多孔隙介质波动方程小波有限差分法

研究流体饱和多孔隙介质中波动方程的数值模拟．针对求解二维弹性波方程问题,提出小波有限差分法．该方法综合了小波多分辨分析计算灵活、计算效率高特性和有限差分易于实现的优点．数值模拟的结果显示,此方法对于求解流体饱和多孔隙介质方程的数值模拟是有效稳定的．     

饱和多孔弹性杆流固耦合动力响应的保结构算法

研究了不可压饱和多孔弹性杆的流固耦合动力响应问题.基于多孔介质理论,根据多孔介质流固混合物动量方程、孔隙流体动量方程及体积分数方程,建立流固耦合不可压饱和多孔弹性杆的轴向振动方程;引入正则变量,构造饱和多孔弹性杆轴向振动方程的广义多辛保结构形式、广义多辛守恒律及广义多辛局部动量误差;采用中点Box离散方法得到轴向振动方程的广义多辛离散格式、广义多辛守恒律数值误差及局部动量数值误差;数值模拟不可压饱和多孔弹性杆的轴向振动过程及流相渗流速度分布,考察了流固两相耦合系数对轴向振动过程及广义多辛守恒律误差和局部动量误差的影响.结果表明,已构造的广义多辛保结构算法具有很高的精确性和长时间的数值稳定性.     

不可压饱和多孔弹性杆动力响应的多辛方法

研究了不可压饱和多孔弹性杆的一维动力响应问题.基于多孔介质理论,在流相和固相微观不可压、固相骨架小变形的假定下,建立了不可压流体饱和多孔弹性杆一维轴向动力响应的数学模型.利用Hamilton空间体系的多辛理论,构造了不可压饱和多孔弹性杆轴向振动方程的多辛形式及其多种局部守恒律.采用中点Box离散方法得到轴向振动方程的多辛离散格式和局部能量守恒律以及局部动量守恒律的离散格式;数值模拟了不可压饱和多孔弹性杆的轴向振动过程,记录了每一时间步的局部能量数值误差和局部动量数值误差.结果表明,已构造的多辛离散格式具有很高的精确性和较长时间的数值稳定性,这为解决饱和多孔介质的动力响应问题提供了新的途径.     

结构非线性风振的路径积分方法

通过引入相关脉动风速滤波,将结构非线性风振方程转变为Ito随机微分方程的形式;该方程的解过程具有Markov性质．在时域内将状态方程展开,利用其瞬时线性化随机方程的解析解,基于路径积分给出了结构非线性风振响应概率密度的形式解,得到了一种分析结构非线性风振响应的新方法．对桅杆算例的数值分析表明,该方法较线性频域分析方法和非线性时域积分方法具有更好的准确性和有效性．     

多孔介质瞬态分析中非分裂PML及时域有限元实现

在波场的数值模拟中,完全匹配层(perfectly matched layer,PML)已经被证明是一种十分有效的吸收技术,并得到了广泛的应用.为了解决具有无限域的多孔介质中2阶弹性波动方程数值模拟中的吸收边界问题,提出了一种非分裂格式的PML(non-splitting perfectly matched layer,NPML).首先,基于Biot多孔介质波动理论,建立了以固相和流相位移表示的2阶动力控制方程,其中考虑了固体颗粒和孔隙流体的可压缩性、惯性以及孔隙流体的粘性.其次,根据复伸展坐标变换的定义,通过Laplace变换获得了非分裂格式PML的频域表达式.然后,借助辅助函数将该方程变换到时间域内,得到了一种有效的非分裂PML.最后,基于Galerkin近似方法,给出了其时域有限元计算格式.通过数值算例分析了该非分裂格式的PML在饱和介质动力响应分析中的有效性.     

饱和粘土非线性渗流规律与径向固结

以多孔介质等效渗流概念得出粘土平均孔隙微尺度数量级范围在0.01 μm至0.1μm,与测试结果一致．实验结果表明饱和粘土微尺度孔隙渗流为非线性流．理论推导表明固液界面作用与渗透率平方根或孔隙半径成反比,固液界面相互作用是导致饱和粘土非线性渗流的重要原因．提出了精确描述饱和粘土微尺度孔隙非线性渗流基本规律的数学模型,其参数量纲明确,物理意义清楚．建立了饱和粘土非线性渗流新定律,统一描述了从较低到较高水力梯度全过程渗流曲线特征,达西定律是其特例．基于新定律,建立了饱和粘土非线性渗流定流量径向固结数学模型．以粘性边界层思想与稳态依序替换法及积分方程法,导出了粘土非线性渗流平均质量守恒方程及活动边界运动方程,给出了饱和粘土非线性渗流超孔隙水压力分布公式与平均固结度计算公式,获得了粘土层压力分布规律和平均固结度随时间变化规律．结果表明:饱和粘土非线性渗流使活动边界运动速度减小．研究结果为粘土地质工程与岩土工程应用提供了新的科学依据．达西渗流径向固结计算是新的非线性渗流固结计算的特例．     

层状横观各向同性饱和地基上圆板的非轴对称振动

研究了层状横观各向同性饱和地基上弹性圆板的非轴对称振动问题．首先,通过方位角的Fourier变换,将圆柱坐标系下横观各向同性饱和土的三维动力方程转化为一阶常微分方程组,基于径向Hankel变换,建立问题的状态方程,求解状态方程后得到传递矩阵;其次,利用传递矩阵,结合层状饱和地基的边界条件、排水条件及层间接触和连续条件,给出了任意简谐激振力作用下层状横观各向同性饱和地基动力响应的通解;然后,按混合边值问题建立层状饱和地基上弹性圆板非轴对称振动的对偶积分方程,并将对偶积分方程化为易于数值计算的第二类Fredholm积分方程,并给出了算例．     

Extended finite element modeling of deformable porous media with arbitrary interfaces

In this paper, an enriched finite element method is presented for numerical simulation of saturated porous media. The arbitrary discontinuities, such as material interfaces, are encountered via the extended finite element method (X-FEM) by enhancing the standard FEM displacements. The X-FEM technique is applied to the governing equations of porous media for the spatial discretization, followed by a generalized Newmark scheme used for the time domain discretization. In X-FEM, the material interfaces are represented independently of element boundaries and the process is accomplished by partitioning the domain with some triangular sub-elements whose Gauss points are used for integration of the domain of elements. Finally, several numerical examples are analyzed, including the dynamic analysis of the failure of lower San Fernando dam, to demonstrate the efficiency of the X-FEM technique in saturated porous soils.     

横观各向同性饱和弹性多孔介质非轴对称动力响应

应用Fourier展开和Hankel变换求解了简谐激励下横观各向同性饱和弹性多孔介质的非轴对称Biot波动方程,得到了一般解。用一般解给出了多孔介质总应力分量的表达式。最后对求解横观各向同性饱和弹性多孔介质非轴对称动力响应边值问题的方法作了系统说明,并且给出了数值分析特例。     

Poroelastodynamic potential method for transversely isotropic fluid-saturated poroelastic media

By introduction of two scalar potentials, an analytical method is developed for the solution of poroelastodynamic boundary value problems in transversely isotropic fluid-saturated poroelastic media. The governing equations of motion are considered in the framework of Biot's complete model without any assumption or simplification. As a case of application, solutions in three dimensions for a transversely isotropic fluid saturated porous half space loaded by an arbitrary distribution of time harmonic tractions at the free surface is derived. The free surface of the half space may be considered either permeable or impermeable. As a particular solution, Green's functions for uniform vertical and horizontal circular patch loads are presented as semi-infinite integrals which may be evaluated by means of an appropriate numerical method proposed. The accuracy of the solutions is verified both analytically and numerically against the preceding solutions. Some numerical results are also presented to clarify the influence of different degrees of anisotropy and frequency of excitation on the response of the medium.     

Solution of Poisson’s equation by analytical boundary element integration

The solution of Poisson’s equation is essential for many branches of science and engineering such as fluid-mechanics, geosciences, and electrostatics. Solution of two-dimensional Poisson’s equations is carried out by BEM based on Galerkin Vector Method in which the integrals that appear in the boundary element method are expressed by analytical integration. In this paper, the Galerkin vector method is developed for more general case than presented in the previous works. The integrals are computed for constant and linear elements in BEM. By employing analytical integration in BEM computation, the numerical schemes and coordinate transformations can be avoided. The presented method can also be used for the multiple domain case. The results of the analytical integration are employed in BEM code and the obtained analytical expression will be applied to several examples where the exact solution exists. The produced results are in good agreement with the exact solution.     

粘性流体两相多孔介质非线性动力问题的罚有限元法

采用基于混合物理论的多孔介质模型,给出粘性流体饱和两相多孔介质非线性动力问题的控制场方程以及相应边值和初值问题的提法,用Galerkin加权残值法导出罚有限元公式,并给出该非线性方程组的迭代求解方法。考虑了体积分数和渗透率与变形相关的情况。用编制的有限元程序计算分析了一维多孔柱体在脉冲载荷作用下的瞬态响应,数值结果表明文中方法正确有效。     

A Boundary Integral Formulation for Poroelastic Materials

Wave propagation in porous media is an important topic, e.g. in geomechanics or the oil-industry. We formulate a linear system of coupled partial differential equations based on Biot's theory with the solid displacements and the pore pressure as the primary unknowns. To solve this system of coupled partial differential equations in a semi-infinite homogeneous domain the BEM (Boundary element method) is especially suitable. Starting from a representation formula a system of two boundary integral equations is derived. These boundary integral equations are used to solve related boundary value problems via a direct approach. Coercivity of the resulting bilinear form is shown, from which unique solvability of the variational formulation follows from injectivity. Using these results we derive the unique solvability of the related boundary integral equations. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Convolution quadrature time-domain boundary element method for 2-D fluid-saturated porous media

For large-scale wave analyses of fluid-saturated porous media, a conventional time-domain boundary element method (BEM) cannot be applied because of the following reasons: (1) no time-domain fundamental solutions are known for some problems, (2) the method sometimes suffers from instability, and (3) the analyses require large amounts of computational time and memory. In this study, an innovative time-domain BEM is developed for a fluid-saturated porous medium. The formulation presented herein overcomes the above disadvantages using a convolution quadrature method (CQM), first proposed by Lubich, and hybrid-parallelization with both MPI and OpenMP. Problems involving the scattering of an incident plane wave by cavities in a 2-D poroelastic medium are solved as a means of validating the proposed method.     

A dynamic boundary value problem arising in the ecology of mangroves

We consider an evolution model describing the vertical movement of water and salt in a domain splitted in two parts: a water reservoir and a saturated porous medium below it, in which a continuous extraction of fresh water takes place (by the roots of mangroves). The problem is formulated in terms of a coupled system of partial differential equations for the salt concentration and the water flow in the porous medium, with a dynamic boundary condition which connects both subdomains.We study the existence and uniqueness of solutions, the stability of the trivial steady state solution, and the conditions for the root zone to reach, in finite time, the threshold value of salt concentration under which mangroves may live.