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

本文首先从弹性力学的基本方程出发，利用Hankel积分变换等数学手段，推导出了单层弹性半空问轴对称问题的刚度矩阵，然后按传统的有限元方法组成总体刚度矩阵。通过求解由总体刚度矩阵所构成的代数方程和Hankel积分逆变换就可解出静荷载作用下多层弹性半空间轴对称问题的精确解。由于刚度矩阵的元素中不含有正指数项，计算时不会出现溢出的现象，从而克服了传递矩阵法的缺点。由于在推导过程中摒弃了应力函数的选择，使得问题的求解更加理论化和合理化，同时也为进一步研究这类问题如温度场，动力学等方向奠定了理论基础。最后，文中还给出了计算实例来证明推导结果的准确性。     

求解动荷载作用下多层粘弹性半空间轴对称问题的精确刚度矩阵法

利用Laplace—Hankel联合积分变换，推导出了单层粘弹性半空间轴对称问题在动荷载作用下，层间完全接触情况的刚度矩阵，然后按传统的有限元方法组成总体刚度矩阵。通过求解由总体刚度矩阵所构成的代数方程就可解出动荷载作用下多层粘弹性半空间轴对称问题的矩阵。由于刚度矩阵的元素中只含有负指数项，计算时不会出现溢出的现象。本文还成功地应用了Durbin的Laplace逆变换的数值方法，求解出了多层粘弹性体的时域解。最后，文中还给出了路面动弯沉的计算结果与实测结果的对比来证明推导结果的准确性。     

弹性地基上四边自由Reissner矩形中厚板的有限积分变换法

将弹性地基视为Winkler模型,利用二维有限积分变换的方法推导出了弹性地基上四边自由矩形中厚板位移和内力的精确解.由于在求解过程中不需要预先人为选取位移函数,而是从弹性地基上中厚板的基本方程出发,直接利用有限积分变换的数学方法求出可以完全满足四边自由边界条件,弹性地基上矩形中厚板问题的精确解,使得问题的求解更加合理.最后通过计算实例验证了所采用方法及所推导出的公式的正确性.     

轴对称横观各向同性层状弹性半空间问题受力分析

本文从柱坐标系弹性力学基本方程出发,将位移场和应力场在径向进行Ｈａｎｋｅｌ变换,利用常微分主程求解原理,直接得出在轴对称荷载作用下横观各向同性半无限弹性空间的位移场,利用此结果推导出单层元的刚度矩阵。     

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

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

层状地基任意形状刚性基础动力响应求解

提出了基于积分变换、对偶方程与精细积分算法求解多层地基任意形状刚性基础的动力刚度问题. 首先在频率波数域内圆柱坐标体系中利用圆形微元的对称与反对称特性建立多层地基中格林影响函数的波动方程,然后将应力和位移关系表示成对偶形式进行精细积分求解以提高计算精度和稳定性. 再将任意形状刚性基础与地基的交界面离散化为一系列圆形微元,利用格林影响函数建立其平动与转动动力刚度的矩阵方程. 该求解方法高效、准确并且计算稳定,适于任意复杂多层地基任意形状基础动力刚度的计算.      

三维粘弹性层状半空间埋置集中荷载动力格林函数求解—修正刚度矩阵法

针对三维粘弹性层状半空间埋置集中荷载作用下动力响应问题,在柱面坐标下,结合径向Hankel积分变换,提出了一种新的求解方法—修正刚度矩阵法。方法基于位势函数理论,将三维问题分解为平面内反应(P-SV波型)和平面外反应（SH波型）两个二维问题的叠加;借鉴结构力学中超静定结构的位移法原理,首先固定荷载所在层的上下界面,通过对波动方程的特解和齐次解叠加得到“固端”反力。进而放松两“固端约束”,利用直接刚度法求得各层面位移,荷载作用层内反应另需叠加上该“固定层”内解,并将特解部分积分（直达波）由全空间解析解代替,解决了当接收点和源点作用水平面接近时的积分收敛问题。算例分析表明,对于低频（可退化为静力状态）和高频问题,本文方法均具有很高的计算效率和精度。     

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

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

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

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

波数-频率域内地基土表面位移Green函数的理论分析

建立了柱面坐标系下分层弹性半空间地基土模型。利用钟阳刚度矩阵法和Haskell-Thomson传递矩阵法推导出所有分层土体之间的振动传递关系;根据Helmholtz定理将土体的位移向量分解成势函数的形式,推导出弹性半空间表面应力与位移之间的关系;再将分层土体和半空间地基土通过位移与应力之间的关系进行耦合,得到分层弹性半空间地基土模型表面位移与应力之间的关系。结合单位脉冲荷载作用下地基土表面的边界条件,推导出波数-频率域内地基土表面位移Green函数的解析解,用Matlab程序语言对理论进行实现并通过算例对地基土表面位移Green函数的特征进行了分析和总结。