A PRECISE INTEGRATION APPROACH FOR THE DYNAMIC-STIFFNESS MATRIX OF STRIP FOOTINGS ON A LAYERED MEDIUM 1）

提出应用精细积分算法计算多层地基的动力刚度问题.精细积分是计算层状介质中波传播的高效而精确的数值方法.利用傅里叶积分变换将层状地基的波动方程转换为频率-波数域内的两点边值问题的常微分方程组,运用精细积分方法求解格林函数,最后再将得到的频率-波数域内地基表面的动力刚度矩阵转换到频率-空间域内,进而得到刚性条带基础频率域的动力柔度或刚度矩阵.所建议的精细积分算法,可以避免一般传递矩阵计算中的指数溢出问题,对各种情况有广泛的适应性,计算稳定,在高频段可以保障收敛性,并能达到较高的计算精度.     

多层地基条带基础动力刚度矩阵的精细积分算法

提出应用精细积分算法计算多层地基的动力刚度问题. 精细积分是计算层状介质中波传播的高效而精确的数值方法. 利用傅里叶积分变换将层状地基的波动方程转换为频率-波数域内的两点边值问题的常微分方程组, 运用精细积分方法求解格林函数, 最后再将得到的频率-波数域内地基表面的动力刚度矩阵转换到频率-空间域内, 进而得到刚性条带基础频率域的动力柔度或刚度矩阵. 所建议的精细积分算法, 可以避免一般传递矩阵计算中的指数溢出问题, 对各种情况有广泛的适应性, 计算稳定, 在高频段可以保障收敛性, 并能达到较高的计算精度.      

多层地基上轴对称受荷刚性圆板问题

利用刚性圆板表面各点位移相等,并结合刚性圆板与地基表面的位移相容条件与光滑接触条件,经由Hankel变换,推导出了刚性圆板与分层地基表面接触应力的对偶积分方程;求解该对偶积分方程,再由多层地基应力与位移的传递矩阵解,并经Hankel逆变换,得到了多层地基上轴对称受荷刚性圆板问题的解.编制了计算程序,并进行了数值分析与计算.计算结果表明:对均匀地基而言,实际工程的计算分析可只考虑4倍刚性圆板直径以内深度范围内的应力与位移;而地基的分层性对地基的位移和应力有着较大的影响,简单地将均匀地基的结论推广到分层和非均匀地基是不恰当的.     

上覆单相弹性层饱和地基上弹性基础的竖向振动分析

基于Biot动力控制方程,运用Fourier积分变换技术,并按照混合边值条件和连续条件建立了上覆单相弹性层饱和地基上弹性基础竖向振动的对偶积分方程.利用正交多项式将对偶积分方程化简,得到了动力柔度系数随无量纲频率b0的变化关系曲线,从而得到了上覆单相弹性层饱和地基上弹性基础的竖向振动规律.数值分析结果表明,对于弹性基础,当弹性基础的挠曲刚度较大时,发现弹性基础的竖向振动特性与刚性基础的类同,可忽略挠曲刚度对竖向振动的影响,且当无量纲频率较小的时候,动力柔度系数Cv随着无量纲频率b0的变化而发生显著的变化,但当无量纲频率b0较大的时候,动力柔度系数Cv受无量纲频率的影响较小,甚至基本上不受影响.当弹性基础的挠曲刚度较小时,随着挠曲刚度的减小,弹性基础的竖向振动将发生显著的变化,动力柔度系数Cv的实部和虚部的绝对值均变大.     

上覆单相弹性层的饱和地基上刚性圆板的扭转振动分析

采用解析的方法研究了上覆单相弹层的饱和地基上刚性圆板的扭转振动。首先运用积分变换技术分别求解了单相弹性介质和饱和介质情况时的控制方程，然后按混合边值条件建立了上覆单相弹性层的饱和地基上刚性圆板扭转振动的对偶积分方程，并把对偶积分方程化为易于数值求解的第二类Fredholm积分方程，并给出了数值算例。     

下卧刚性基岩的饱和地基上基础的动力分析

采用积分变换与积分方程研究了单层饱和地基上刚性基础的竖向振动问题，首先从饱和土的动力控制方程出发，考虑混合边值条件，获得了一组描述基础振动的对偶积分方程，进而给出了单层饱和地基表面的动力柔度系数。该文的工作可视为对饱和半空间地基上基础振动解答的推广。     

成层非饱和地基上矩形基础稳态响应解析研究

采用积分变换法求解非饱和土的控制方程,得到变换域内土体的位移和应力表达式;建立了单层非饱和土体的刚度矩阵,并组装每一层的刚度矩阵,构成了层状地基的总体刚度矩阵;结合地基上下表面的边界条件,推导了层状非饱和地基表面位移的积分形式解答。另外,选取带有补充项的双重余弦级数解作为矩形基础的振型函数,并令其与地表位移按余弦级数展开的表达式相等,建立了矩形基础与地基的协调条件。最后,联立矩形基础的控制方程、边界条件和基础与地基之间的协调条件,求解得到矩形基础的挠度幅值、弯矩幅值以及基底压力幅值。选取已有文献中的非饱和土参数计算,其结果与文献吻合良好,验证了本文方法的正确性。本文进一步研究了单层地基和双层地基上矩形板的动力响应,分析了土体参数对矩形基础动力响应幅值的影响规律。结果表明:地基层厚与基础尺寸之比为5时,地基就可以看作半空间地基;矩形基础的稳态响应随非饱和土的饱和度的增大而减小;渗透率对矩形基础动力响应的影响不明显。     

层状地基中单桩扭转振动特性分析

采用积分变换和Muki的方法求解了层状地基中单桩的扭转振动问题.在分析过程中,首先对基本控制方程进行Hankel变换,建立了单层地基的初参数解答和刚度矩阵,得到层状地基的递推矩阵;然后利用递推矩阵、边界条件和桩-土变形协调条件建立了层状地基中单桩扭转振动问题的基本积分方程并进行数值求解.文末数值算例给出了退化的层状地基中刚性单桩的扭转变形,其结果与已有经典解答吻合良好.同时,并研究了两层地基中单桩的扭转动力响应,分析了桩-土参数对动力响应的影响,所得结论对工程实践和桩基扭转波检测有一定的指导意义.     

饱和地基上弹性圆板的动力响应

研究弹性圆板在饱和地基上的垂直振动特性,即首先应用Hankel变换方法求解饱和土波动方程,然后按混合边值条件建立饱和地基上圆板垂直振动的对偶积分方程,用一种简便的方法,对偶积分方程可化为易于数值计算的第二类Fredholm积分方程。文末的数值分析得出了板振动的一些规律性,由此表明当板的挠曲刚度D趋于无穷大且不计板的质量时,其结果和无质量刚性圆盘在饱和地基上的振动特性完全一致。     

多层弹性半空间问题解的精确刚度矩阵法

利用二维Fourier积分变换和刚度矩阵法,得到了直角坐标系下多层弹性半空间问题的精确解.在推导的过程中,直接从空间弹性力学的基本方程出发,利用积分变换等数学手段,首先推导出了单层空间弹性问题的刚度矩阵,然后按有限元法组成总体刚度矩阵.通过求解由总体刚度矩阵所构成的代数方程和积分逆变换,就可解出在任意静荷载作用下多层弹性半空间问题的精确解.由于刚度矩阵中不含有正指数项,计算时不会出现溢出现象,从而克服了传递矩阵法的缺点.由于在推导过程中摒弃了应力函数的选择,使得问题的求解更加合理化.最后还给出了计算实例来证明推导结果的准确性.     

Scattering of a plane harmonic SH wave by multiple layered inclusions

Scattering of a plane harmonic SH wave by an arbitrary number of layered inclusions in a half-space is investigated by using a direct boundary integral equation method. The inclusions of arbitrary shape and placement are embedded within an elastic half-space. The effects of multiple scattering, the geometry, and the impedance contrast of the materials for layered inclusions and pipes are considered in detail.     

波动方程的精细逐步积分法

应用现有的波动方程求解方法解决工程实际问题尚存在一定的局限性。本文在结构动力方程精细逐步积分的基础上,提出了波动方程初边值问题的精细逐步积分法,并分别给出了不同边界条件下的精细逐步积分格式。此数值方法虽然是显式积分方法,却是无条件稳定的。分别用精细逐步积分法和其它已有的方法对两个算例进行了计算,一个是有解析解的例子,该例验证了此方法的准确性,另一个例子是求解由波动方程及初始条件和边界条件组成的有杆抽油系统预测模型,此例验证了精细逐步积分法的高效性。     

Contact problems for several transversely isotropic elastic layers on a smooth elastic half-space

The idea of generalized images, first used by the author for the case of crack problems, is applied here to solve a contact problem for n transversely isotropic elastic layers, with smooth interfaces, resting on a smooth elastic half-space, made of a different transversely isotropic material. A rigid punch of arbitrary shape is pressed against the top layer’s free surface. The governing integral equation is derived for the case of two layers; it is mathematically equivalent to that of an electrostatic problem of an infinite row of coaxial charged disks in the shape of the domain of contact. This result is then generalized for an arbitrary number of layers. As a comparison, the method of integral transforms is also used to solve the problem. The main difference of our integral transform approach with the existing ones is in separating of our half-space solution from the integral transform terms. It is shown that both methods lead to the same results, thus giving a new interpretation to the integral transform as a sum of an infinite series of generalized images.     

The axisymmetric mixed boundary-value problem of the vertical vibration of a rigid foundation on saturated layered soil subgrade

Based on the theory of elastic wave propagation in saturated soil subgrade established by the author of this paper, the axisymmetric vertical vibration of a rigid circular foundation resting on partially saturated soil subgrade which is composed of a dry elastic layer and a saturated substratum is studied. The analysis relied on the use of integral transform techniques and a pair of dual integral equations governing the vertical vibration of the rigid foundation is listed under the consideration of mixed boundary-value condition. The results are reduced to the case for saturated half-space. The set of dual integral equations are reduced to a Fredholm integral equation of the second kind and solved by numerical procedures. Numerical examples are given at the end of the paper and plots of the dynamic compliance coefficient C_v versus the dimensionless frequency a₀ are presented. Foundation item: the National Natural Science Foundation of China (59579018) Biography: CHEN Sheng-li (1974-)     

Vertical vibrations of elastic foundation resting on saturated half-space

This paper is mainly concerned with the dynamic response of an elastic foun- dation of finite height bounded to the surface of a saturated half-space.The foundation is subjected to time-harmonic vertical loadings.First,the transform solutions for the governing equations of the saturated media are obtained.Then,based on the assumption that the contact between the foundation and the half-space is fully relaxed and the half- space is completely pervious or impervious,this dynamic mixed boundary-value problem can lead to dual integral equations,which can be further reduced to the Fredholm integral equations of the second kind and solved by numerical procedures.In the numerical exam- ples,the dynamic compliances,displacements and pore pressure are developed for a wide range of frequencies and material/geometrical properties of the saturated soil-foundation system.In most of the cases,the dynamic behavior of an elastic foundation resting on the saturated media significantly differs from that of a rigid disc on the saturated half-space. The solutions obtained can be used to study a variety of wave propagation problems and dynamic soil-structure interactions.     

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

研究横观各向同性饱和土地基上中厚弹性圆板的非轴对称振动问题。基于横观各向同性饱和介质Biot波动方程的一般解,按混合边值问题建立了饱和地基与弹性中厚圆板非轴对称动力相互作用的对偶积分方程,并将对偶积分方程转化为易于计算的第二类Fredholm积分方程;采用数值方法求解该积分方程。数值算例结果表明,当h/a>0.05时,饱和半空间体上中厚度圆板在不同频率下的振动特性与相应频率下的刚性板的振动特性基本相同,当h/a<0.05时,板中心的位移将随h/a的减小而增大。     

Application of generalized images method to contact problems for a transversely isotropic elastic layer on a smooth half-space

The idea, first used by the author for the case of crack problems, is applied here to solve a contact problem for a transversely isotropic elastic layer resting on a smooth elastic half-space, made of a different transversely isotropic material. A rigid punch of arbitrary shape is pressed against the layer’s free surface. The governing integral equation is derived; it is mathematically equivalent to that of an electrostatic problem of an infinite row of coaxial charged disks in the shape of the domain of contact. The case of circular domain of contact is considered in detail. As a comparison, the method of integral transforms is also used to solve the problem. The main difference of our integral transform approach with the existing ones is in separating of our half-space solution from the integral transform terms. It is shown that both methods lead to the same results, thus giving a new interpretation to the integral transform as a sum of an infinite series of generalized images.