横观各向同性含液饱和多孔介质黏弹性瞬态反应

在考虑横观各向同性含液饱和多孔介质固体骨架和流体可压缩性以及固体骨架的黏弹性特征下,基于横观各向同性含液饱和多孔介质u-w形式的三维动力控制方程,以固相位移u、液相相对位移w为基本未知量,综合运用Laplace变换、双重Fourier变换等方法,在直角坐标系下通过引入中间变量,将六元2阶动力控制方程组化为两组各含4个未知变量的常微分方程组,给出了直角坐标系下横观各向同性含液饱和多孔介质三维黏弹性动力反应的积分形式一般解;作为理论推导的验证,通过引入初始条件和边界条件,对横观各向同性含液饱和多孔介质半空间黏弹性瞬态反应问题进行了求解。解答的退化验证表明,所推导的理论解是正确的。     

三维水质预测问题的Laplace变换边界元法及其推广

本文用边界元法解决非稳定系统中水质预测的计算问题,建立了三维水质预测问题的边界元法计算模型,为简化计算手续、提高计算精度,讨论了Laplace变换边界元法,推广了Laplace变换边界元法中的Schapery方法。     

横观各向同性饱和土体三维粘弹性动力分析

采用针对横观各向同性饱和土体u-w形式三维粘弹性动力方程,考虑土骨架的粘弹性性质且基于粘弹性理论,通过运用Fourier 展开、Laplace 和Hankel 积分变换方法和引入中间变量,将含有粘弹性参数的六元二阶偏微分运动控制方程组,化为2组各含4个未知变量的常微分方程组,从而给出了柱坐标系下粘弹性横观各向同性饱和土体在非轴对称动力荷载作用下的瞬态反应的土骨架位移分量、孔隙流体相对于土骨架的位移分量瞬态反应一般解。在此基础上,引入初始条件和边界条件,对垂直向和水平向动力荷载作用下半空间边值问题进行了求解。根据动力时域解答的一般解,利用Laplace和Hankel 数值逆变换技术,编制了相应的数值计算程序。并进行了实例验证和弹性、粘弹性解的对比分析。结果表明,在进行横观各向同性饱和土体动力分析时,考虑土骨架的粘弹性是必要的。     

Laplace方程Cauchy问题的一种数值解法

研究了Laplace方程Cauchy问题的数值求解，该问题是一个典型的病态问题。利用格林（Green）公式将Laplace方程的Cauchy问题转化为Hausdorff矩问题。本文利用一种新方法，即矩问题的积分方程方法，求解矩问题，设计了二维Laplace方程Cauchy问题稳定的算法，给出了近似解的误差估计，并对二维Cauchy问题进行了数值模拟。     

瓦斯抽采过程中孔壁的动态响应分析

将含瓦斯煤岩体简化为气饱和的各向同性均质体,基于混合物理论建立了含瓦斯煤岩体的物理方程、几何方程、连续方程和动力控制方程。瓦斯抽采孔的长度远大于直径,将瓦斯抽采过程中抽采负压引起的孔壁的动态响应问题简化为二维平面应变问题,借助Laplace变换和Laplace逆变换分别建立了相关问题的频域和时域解答。以河南省某煤矿的煤岩体瓦斯参数为例,对孔边的动态响应进行了数值计算,分析了无量纲径向位移、径向应力和环向应力等的变化规律,结果表明：径向位移和径向应力等动态响应主要集中在孔边的一定区域内,且该区域随着抽采时间的增大而扩大;孔壁产生波动较大的环向应力,且该环向应力随着抽采负压的增大而增大。     

多层地基非轴对称Biot固结的理论解

提出了一种新方法来推导多层地基非轴对称Biot固结问题的传递矩阵,相应的理论推导工作量较少。基于柱坐标系下非轴对称Biot固结的基本方程,通过引入中间变量,并对坐标 进行Fourier级数展开,对时间 和坐标 进行Laplace-Hankel变换,得到了6×6阶和2×2阶的两组常微分方程;然后,两组常微分方程进行关于 的Laplace变换和逆变换,得到了单层地基非轴对称Biot固结问题的传递矩阵;结合边界条件和连续条件,运用传递矩阵法得到多层地基非轴对称Biot固结问题在Laplace-Hankel变换域内的解,通过Laplace-Hankel逆变换得到了该问题物理域的真实解答。编制了计算程序,并进行了数值计算与分析。该方法具有计算速度较快的特点,方便工程应用。     

渗透破坏若干问题的解析解及机理分析

讨论了柱坐标下不同介质有限区域内用分离变量法求解Laplace 方程的问题。通过一系列计算试验,给出了柱坐标下Laplace 方程在不同介质有限区域内用分离变量方法求的解,并用该解对渗透破坏现象进行了讨论和实例验证。认为渗透系数分布不均匀是发生渗透破坏的根本原因,渗透破坏朝着渗透系数差异较大的方向发展,并对管涌的复杂机理进行了引申阐述。该工作对于流土、管涌机理的更深入研究具有重要的意义。     

粘弹性准饱和土中球空腔的动力响应

从工程实际出发,采用粘弹性两相介质模型,考虑土骨架的粘性以及流体与固体之间的耦合作用,利用Laplace变换求解了粘弹性准饱和土中球空腔的动力响应问题,得到了变换域内的解析解。借助数值Laplace反变换,数值分析了粘弹性准饱和土中球空腔动力响应的位移、应力及孔压的变化规律。为分析地下结构动力响应提供了一种有效的方法,模型符合工程实际,有一定的工程应用价值。     

半空间埋置动点源荷载问题的传递矩阵法

引用了流体饱和两相多孔介质的动力控制方程分析半空间埋置动点源荷载问题的位移和变形。经过Laplace Hankel变换 ,控制方程化成常微分方程组。利用数学软件mathmatic对上述方程组求解 ,可以得到单层砂土的传递矩阵。分析过程中 ,假设在两层面上 ,位移与应力相互连续 ,可以借鉴有限元的思想进行耦合计算。这样就获得了在饱和砂土中施加竖向动荷载问题的Laplace Hankel变换解 ,其最终的解还需要通过Laplace Hankel逆变换得到     

饱和土中单桩水平瞬态响应研究

基于Biot动力固结方程,应用Novak薄层方法研究了饱和土中单桩在水平冲击荷载作用下的动力响应问题。首先引入势函数对方程解耦,再通过算子分解和分离变量法,结合初始边界条件,并联立桩基振动微分方程,推导了Laplace变换域内桩身位移函数及内力表达式。采用Laplace逆变换的优化模型求得了时域内瞬态响应的封闭解。将该解退化到单相介质条件下的桩顶位移-时间曲线与已有的边界元方法的结果基本吻合,验证了解答的正确性。参数分析结果表明：桩-土模量比、长径比、渗透系数对桩身位移均影响显著,而同条件下长径比超过一特定值后则影响较小,且桩土模量比是影响桩身弯矩大小及其分布的重要参数。     

Boundary Element Analysis (Laplace Transform Solution) of Groundwater Unsteady Flow to a Multiple Well System in a Confined Aquifer1

Abstract The calculations of unsteady flow to a multiple well system with the application of boundary element method (BEM) are discussed. The mathematical model of unsteady well flow is a boundary value problem of parabolic differential equation. It is changed into an elliptic one by Laplace transform to eliminate time variable. The image function of water head H can be solved by BEM. We derived the boundary integral equation of the transformed variable H and the discretization form of it, so that there is no need to discretize the boundaries of well walls and it becomes easier to solve the groundwater head H by numerical inversion.     

Thermo-mechanical performance of layered transversely isotropic media around a cylindrical/tubular heat source

We develop a new numerical model based on a precise integration method to investigate the coupled thermo-mechanical performance of layered transversely isotropic media around a cylindrical/tubular heat source. To obtain the relational matrices of the extended precise integration method, we first convert the governing equations of the problem into ordinary differential matrix equations through the Laplace–Hankel transform. Then, the cylindrical heat source is divided into a series of plane heat sources, and the plane temperature load term is added to the state vector between layer elements. By combining the layer elements, we build a layered transversely isotropic numerical model containing a cylindrical heat source in the transformed domain. Finally, we solve the model in the transformed domain and obtain the solution of the problem in the real domain through the Laplace–Hankel transform inversion. The accuracy of this method is verified by comparing the solutions with the results of the analytical method and the finite element method. Then, we study the influence of the anisotropy of thermal parameters, the embedded depth, the length/radius ratio, the type of heat source and the stratification of the medium on the thermo-mechanical coupled performance.

     

Boundary element analysis of contaminant transport in fractured porous media

A two-dimensional boundary integral method to analyse the flow of contaminant in fractured media having a two- or three-dimensional orthogonal fracture network is presented. The method assumes that the fractures provide the paths of least resistance for transport of contaminants while the matrix, because of its low permeability, acts as ‘storage blocks’ into which the contaminant diffuses. Laplace transform is used to eliminate the time variable in the governing equation in order to facilitate the formulation of a boundary integral equation in the Laplace transform space. Conventional boundary element techniques are applied to solve for the contaminant concentrations at specified locations in the spatial domain. The concentration in the time domain is then obtained by using an efficient inversion technique developed by Talbot. The method is able to analyse the behaviour of waste repositories which have diminishing concentration due to the mass transport of the contaminant into the surrounding fractured media.     

A method of computing the consolidation behaviour of layered soils using direct numerical inversion of Laplace transforms

A method is presented for obtaining the consolidation behaviour of a layered soil subjected to strip, circular, or rectangular surface loadings, or subjected to fluid withdrawal due to pumping. The solution method involves applying a Fourier or Hankel transform to the field quantities along with a Laplace transformation. The effect of the Fourier or Hankel transform is to reduce a two- or three-dimensional problem or one involving axial symmetry, to one involving only a single spatial dimension. In cases where the soil is horizontally layered, this has great advantages over conventional methods, such as finite element or finite difference methods, since very little computer storage and data preparation time is required. Solution of the time dependent problem is achieved by applying a Laplace transformation to the field variables, obtaining solutions in Laplace transform space, and then numerically inverting the transformed solutions to obtain the real time behaviour. This eliminates the need for ‘marching type’ schemes where a solution is found from one at a previous time. By direct inversion of the Laplace transform, a solution may be obtained directly at any given time.     

Extended precise integration method for axisymmetric thermo‐elastic problem in transversely isotropic material

This paper presents a numerical solution for the analysis of the axisymmetric thermo‐elastic problem in transversely isotropic material due to a buried heat source by means of extended precise integral method. By virtue of the Laplace–Hankel transform applied into the basic governing equations, an ordinary differential matrix equation is achieved, which describes the relationship between the generalized stresses and displacements in transformed domain. An extended precise integration method is introduced to solve the aforementioned matrix equation, and the actual solution in the physical domain is acquired by inverting the Laplace–Hankel transform. Numerical examples are carried out to demonstrate the accuracy of the proposed method and elucidate the influence of the character of transverse isotropy, the anisotropy of linear expansion coefficient, the anisotropy of thermal diffusivity, and medium's stratification on the thermo‐elastic response. Copyright © 2015 John Wiley & Sons, Ltd.     

Semi-analytical solutions of two-dimensional plane strain consolidation in unsaturated soils subjected to the lateral semipermeable drainage boundary

The study presents semi-analytical solutions of two-dimensional plane strain consolidation problem in unsaturated soils incorporating the lateral semipermeable drainage boundary by adopting Fourier sine series and Laplace transform. The two-dimensional plane strain consolidation equations in the form of two-order partial differential equations with three variables are firstly converted to two-order partial differential equations with two variables, which are similar to those of one-dimensional consolidation problem. The four-order ordinary differential equations about excess pore-air and excess pore-water pressures are got by applying Laplace transform and the substitution method. Then, the solutions of excess pore pressures and settlement are achieved in the Laplace transform domain. Afterwards, on the basis of Crump's method, the inverse Laplace transform is conducted to obtain the analytical solutions in time domain. The comparison is conducted to verify the exactness of the obtained solutions, and the two-dimensional plane strain consolidation property with the lateral semipermeable drainage boundary is illustrated and discussed. Parametric studies are demonstrated for the excess pore pressures and normalized settlement with the change of the boundary parameters, air-water and lateral-vertical permeability coefficients, and the distance and depth. It can be found that the lateral semipermeable drainage boundary impedes the consolidation rate obviously, and when different investigated parameters are adopted, the consolidation property is similar to each other under the later permeable and semipermeable drainage boundary conditions.     

频率域中有限元素法的探讨

本文指出了用拉普拉斯变换处理有限元方程和用富里叶变换所得到的结果是等价的。从物理方面讨论了频率域有限元的局限性。     

Numerical analysis of rafts on visco-elastic media using eigenvector expansions

Problems of a raft on a visco-elastic continuum may be converted to equivalent elastic problems by means of a Laplace transformation. In this paper a numerical solution of the equivalent problem is obtained by finite element or other techniques. This solution is converted into the form of an eigenvector expansion and then transformed into a numerical solution to the original problem by inverting the Laplace transform. This form of solution has the advantage of applying to any visco-elastic model, and of requiring little additional computation as a result of changing the visco-elastic model, the relative raft-soil stiffness or the load pattern. The application of the method is illustrated by results for circular rafts, strip footings of finite length and rectangular rafts, and particular attention is paid to a realistic soil creep function which is asymptotic to a linear function of log time.     

结构体系的一种新型半解析动力算法——地下结构动力分析的拉普拉斯积分变换解法

得到了地下结构在水平地震作用下的动力运动方程,利用拉普拉斯变换和结构的对接条件,得到了求解待定系数的超越方程。利用拉普拉斯逆变换得到了结构的位移表达式和确定拉普拉斯参数极点的方程。进一步研究发现,拉普拉斯变换参数s与结构频率ω之间存在某种特殊的关系:s2=?ω2。为了避免确定待定系数及拉普拉斯参数极点时求解复杂超越方程的困难,建议通过数值的办法先把结构体系的频率求出来,再通过拉普拉斯逆变换留数定理求解结构的位移和内力。     

Semi‐analytical solution to one‐dimensional consolidation for unsaturated soils with semi‐permeable drainage boundary under time‐dependent loading

This paper presents semi‐analytical solutions to Fredlund and Hasan's one‐dimensional consolidation of unsaturated soils with semi‐permeable drainage boundary under time‐dependent loadings. Two variables are introduced to transform two coupled governing equations of pore‐water and pore‐air pressures into an equivalent set of partial differential equations, which are easily solved by the Laplace transform. The pore‐water pressure, pore‐air pressure and settlement are obtained in the Laplace domain. Crump's method is adopted to perform the inverse Laplace transform in order to obtain semi‐analytical solutions in time domain. It is shown that the present solutions are more general and have a good agreement with the existing solutions from literatures. Furthermore, the current solutions can also be degenerated into conventional solutions to one‐dimensional consolidation of unsaturated soils with homogeneous boundaries. Finally, several numerical examples are provided to illustrate consolidation behavior of unsaturated soils under four types of time‐dependent loadings, including instantaneous loading, ramp loading, exponential loading and sinusoidal loading. Parametric studies are illustrated by variations of pore‐air pressure, pore‐water pressure and settlement at different values of the ratio of air–water permeability coefficient, depth and loading parameters. Copyright © 2017 John Wiley & Sons, Ltd.