横观各向同性层状地基动应力响应的求解与分析

动力响应问题的求解对于地基在外荷载作用下引起的弹性波动问题研究有重要的意义。本文提出了一种求解横观各向同性层状地基在施加时间简谐荷载作用下任意点的应力响应的算法。此算法利用傅里叶变换将广义平面应变问题频率-空间域的动力方程转化到频率-波数域内,结合对偶变量的引入,利用高精度的精细积分算法对状态方程进行求解,在得到频率-波数域内的位移响应的基础上,利用傅里叶逆变换得到任意点的动应力响应。简谐荷载不仅可以施加在地基表面,而且可以施加在地基内部。对比算例验证了本文算法的准确性,同时对地基各向异性特性、激励频率和阻尼比对动应力响应的影响进行了参数分析,为工程实际提供可靠的数值依据。     

墙下条形基础与层状横观各向同性地基共同作用

运用对偶积分方程来求解层状横观各向同性地基与墙下条形基础的共同作用问题。从直角坐标平面应变问题控制方程出发,通过傅里叶（Fourier）变换和层间连续性条件,可以得到层状横观各向同性地基的传递矩阵解。基于该传递矩阵解,并利用条形基础与地基接触的混合边值条件,推导出一组关于基础挠度和地基反力的对偶积分方程。考虑墙下条形基础受到竖向集中荷载的情况,利用弹性薄板理论先求解出条形基础挠度;随后应用雅可比（Jacobi）正交多项式和级数展开的方法,将对偶积分方程转化为线性代数方程组进行求解。编制了相应的计算程序,其计算结果与有限元软件ABAQUS的结果基本吻合,从而验证了所提理论的正确性。算例分析表明,板土相对刚度与地基成层性对地基反力、地表沉降和沿z轴竖向正应力有很大的影响。     

横观各向同性饱和地基三维瞬态Lamb问题

研究了横观各向同性饱和土地基在地表动力荷载作用下的三维瞬态响应。基于饱和多孔介质的三维Biot波动理论,利用Laplace变换,建立圆柱坐标系下横观各向同性饱和土的波动方程;解耦波动方程后,根据算子理论,并借助Fourier展开和Hankel变换技术,得到瞬态荷载作用下,饱和土介质的土骨架位移和应力、孔隙水相对位移和孔隙水压力的一般解;利用一般解,给出横观各向同性饱和地基在地表集中荷载激励下的瞬态Lamb问题的解答。数值算例结果表明,采用各向同性饱和介质的动力学模型,不能准确描述具有明显各向异性特性的饱和土地基的瞬态动力特性。     

横观各向同性地基中管桩扭转振动响应解析解

考虑地基沉积过程中产生的竖向和水平向力学性质的差异,对横观各向同性地基中管桩扭转振动频域响应进行了理论研究。基于横观各向同性材料的本构关系以及桩-土耦合扭转振动,建立了桩土系统定解问题,通过Laplace变换和分离变量法求得了桩周土和桩芯土扭转振动位移形式解。通过桩-土接触面的连续条件,求得了管桩扭转频域响应解析解,并得到了桩顶复动刚度和速度导纳的表达式。将所得解退化到横观各向同性地基中实心桩解以及均匀地基中管桩解,并与已有文献进行了对比,验证了解的合理性。通过数值算例,分析了桩周土和桩芯土的横观各向同性力学参数对桩顶扭转复刚度及速度导纳的影响。     

横观各向同性地基模型的理论分析及其应用

根据横观各向同性体各向同性面水平时本构方程推导出了横观各向同性面倾斜时的本构方程并得出了横观各向同性地基在上述两种情况下的自重应力计算公式。运用本论文的理论 ,从自重应力角度能很好的解释水平应力大于竖直应力和水平应力的各向异性这些重要的实测现象     

三维横观各向同性成层地基的传递矩阵解

通过解耦变换推导出三维直角坐标系下横观各向同性地基的非耦合状态方程;利用双重Fourier变换以及Cayley-Hamilton定理得到了单层地基的传递矩阵;结合边界条件和层间连续条件进而得其传递矩阵解。编制了相应程序并进行了数值计算与分析,结果表明：数值结果与已有文献结果十分吻合,地基的横观各向同性性质与成层性质对受荷地基中竖向位移和应力的影响较为显著。     

横观各向同性地基上刚性矩形基础分析

对横观各向同性地基上刚性矩形板进行了求解。首先,利用表面受矩形均布荷载作用下的层状横观各向同性地基的位移解答,获得地基的柔度矩阵;然后,通过刚性矩形基础与层状横观各向同性地基的协调条件,建立刚性矩形基础与横观各向同性地基共同作用的方程,进而求得基础的地基反力。通过编制相应的程序,确定了合理的网格划分值;最后,进行算例分析,分析了地基横观各向同性性质、矩形刚性基础的长宽比以及地基分层性对地基反力的影响。分析表明：以上3种因素对地基反力有重要影响。     

多层横观各向同性地基轴对称固结的传递矩阵解

从横观各向同性地基轴对称Biot固结的基本方程出发,通过关于t的Laplace变换和关于r的Hankel变换,得到关于z的一阶常微分方程组。然后,对变换域内的基本未知量进行线性化处理,建立了变换域内的基本状态变量在z = 0处和任意深度处z的显式关系。利用传递矩阵法,结合层间连续性条件和边界条件,得到了多层横观各向同性地基的Biot固结轴对称问题的解答。该解答能避免随着层数增加而需要求解大型方程组的困难,明显地提高了计算效率。     

列车荷载下轨道系统-层状横观各向同性饱和地基动力响应

基于Biot波动理论,构建列车荷载-轨道系统-双层状横观各向同性饱和地基模型,将模型分为上覆路轨系统和地层系统。对上覆路轨系统和地层系统处理,并利用双重Fourier变换技术,在变换域中将横观各向同性饱和地基动力响应的求解简化为求解一个6阶控制方程的特征值问题,进而得到了列车荷载作用下双层横观各向同性饱和地基力响应的解析结果。利用离散Fourier逆变换得到数值计算结果,重点分析了上下土层的刚度和泊松比对位移和孔隙水压力和剪切应力响应的影响,结果表明,上、下土层刚度差异对地基动力响应有较大影响,土层各向异性参数中模量的影响较泊松比大。计算结果可为软土路基加固深度的确定提供理论依据。     

水平条形均布荷载下横观各向同性弹性地基的位移分析

分析了水平条形均布荷载下横观各向同性地基位移规律,归纳了位移随弹性模量比n之不同而变化的规律,提出了水平条形均布荷载下横观各向同性地基位移的近似简化解,以供岩土工程设计人员参考。     

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.     

Response of multilayered transversely isotropic medium due to axisymmetric loads

A novel procedure associated with the precise integration method (PIM) and the technique of dual vector is proposed to effectively calculate the magnitude and distribution of deformations in a homogeneous multilayered transversely isotropic medium. The planes of transverse isotropy are assumed to be parallel to the horizontal surface of the soil system. The linearly elastic medium is subjected to four types of vertically acting axisymmetric loads prescribed either at the external surface or in the interior of the soil medium. There are no limits for the thicknesses and number of soil layers to be considered. By virtue of the governing equations of motion and the constitutive equations of the transversely isotropic elastic body, and based on the Hankel integral transform and a dual vector formulation in a cylindrical coordinate system, the partial differential motion equations can be converted into first‐order ordinary differential matrix equations. Applying the approach of PIM, it is convenient to obtain the solutions of ordinary differential matrix equations for the continuously homogeneous multilayered transversely isotropic elastic soil in the transformed domain. The PIM is a highly accurate algorithm to solve the sets of first‐order ordinary differential equations, which can ensure to achieve any desired accuracy of the solutions. What is more, all calculations are based on the standard method with the corresponding algebraic operations. Computational efforts can be reduced to a great extent. Finally, numerical examples are provided to illustrate the accuracy and effectiveness of the proposed approach. Some more cases are analyzed to evaluate the influences of the elastic parameters of the transversely isotropic media on the load‐displacement responses. Copyright © 2015 John Wiley & Sons, Ltd.     

3D dynamic response of a transversely isotropic multilayered medium subjected to a moving load

The dynamic problem of a transversely isotropic multilayered medium is reducible to quasi‐static problem by introducing a moving system that travels synchronously with the load. Based on the governing equations in the moving system, the ordinary differential equations in the Fourier transformed domain are deduced. An extended precise integration method is adopted to solve the ordinary differential equations, and the solution in the physical domain is recovered by the inverse Fourier transform. Numerical examples are performed to verify the accuracy of the presented method and to analyze the influence of material properties and the load characteristic.     

Static analysis of vertically loaded pile group in multilayered transversely isotropic soils

In this paper, a coupling approach is presented to study the static responses of vertically loaded pile group embedded in multilayered transversely isotropic soils. The individual pile in pile group is modeled by the finite element method, while the analytical layer-element method is applied to represent the soil's behavior. Then, the interaction equation of piles and soils is obtained by considering the force equilibrium and displacement compatibility conditions and solved by a FORTRAN program. The results computed by the proposed approach compare favorably with those from some existing solutions and field test. Some typical parametric analysis cases are investigated to study the effect of soil anisotropy, pile stiffness ratio, and pile spacing on the behavior of vertically loaded pile group.     

Coupled consolidation and heat flow analysis of layered soils surrounding cylindrical heat sources using a precise integration technique

The engineering applications of energy piles, geological radioactive waste disposals, and mining wells of geothermal and petroleum are usually associated with strong coupled behaviour of consolidation and heat flow. This paper aims to present an efficient precise integration technique (PIT) for the analysis of such behaviour within layered saturated soils surrounding cylindrical heat sources (ie, with a cross section as a point, ring, or disc). Each soil layer, together with its embedded part of heat source, is divided into 2^N layer elements with equal thickness. Then any pair of adjacent two layer elements are merged into a heat source on the interface. With the aid of Taylor series expansion and recurrence formula for adjacent layer elements combination, such problems can be solved by means of an improved PIT. Typical examples are performed to examine the effects of heat source type and soils layered properties on the coupled consolidation and heat flow responses. The elevation of the clay surface increases with time because of thermal expansion and reaches a peak value before showing a tendency of getting stabilised because of the dissipation of pore pressure becoming dominant. Such a peak cannot be achieved in sand case because of no accumulation of pore pressure. The influencing area of the heat source was found to be limited to near the source. These quantitative results serve as good verification of the presented technique, which proves to be remarkably efficient and several orders more accurate than traditional numerical techniques in that it ideally reaches the accuracy limit of the hardware of the computers used.     

Analytical layer element solutions for deformations of transversely isotropic multilayered elastic media under nonaxisymmetric loading

This paper presents the analytical layer element solutions for deformations of transversely isotropic elastic media subjected to nonaxisymmetric loading at an arbitrary depth. The state vectors for the nonaxisymmetric problem are deduced through the substitution of the Hu Hai‐chang solutions into the basic equations for the transversely isotropic elastic media. From the state vectors, the analytical layer element of a single layer is obtained in the Hankel transformed domain. The analytical layer element is an exact and symmetric stiffness matrix whose elements are without positive exponential functions, which can not only simplify the calculation but also improve the stability of computation. On the basis of the continuity conditions between adjacent layers, the global stiffness matrix is obtained by assembling the interrelated layer elements. The solutions for the multilayered elastic media in the transformed domain are obtained by solving the algebraic equation of the global stiffness matrix, which satisfies the boundary conditions. The actual solutions in the physical domain are further obtained by inverting the Hankel transform. Finally, some cases are analyzed to verify the solutions and evaluate the influences of the transversely isotropic character and stratified character of the media on the load–displacement responses. The numerical results show that the variations of the elastic properties between layers have a great effect on the displacements of the multilayered media. Copyright © 2014 John Wiley & Sons, Ltd.     

Fundamental solutions for a fluid‐saturated,transversely isotropic,poroelastic solid

The fundamental solutions were obtained for step‐like point forces acting in three orthogonal directions and an instantaneous fluid point source in a fluid‐saturated, porous, infinite solid of transversely isotropic elasticity and permeability. After expressing the governing equations in the form of matrix in the Laplace space, we employed Kupradze's method together with the triple Fourier transforms. This method reduces the simultaneous partial differential equations with respect to three displacement components and a pore fluid pressure to a differential equation in terms of only one potential scalar function, which can be operationally solved in the transformed space. After the Laplace inversion of the potential, the residue theorem was applied to its Fourier inverse transform with respect to one of the transformation variables. The Fourier transforms with respect to two other variables were rewritten into the Hankel transforms. Copyright © 2002 John Wiley & Sons, Ltd.