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

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

横观各向同性饱和地基与中厚圆板的非轴对称动力相互作用

研究横观各向同性饱和土地基上中厚弹性圆板的非轴对称振动问题,即首先利用Fourier展开和Hankel变换技术,求解了简谐激励下横观各向同性饱和土地基的非轴对称Biot波动方程,然后按混合边值问题建立地基与弹性中厚圆板非轴对称动力相互作用的对偶积分方程,并将对偶积分方程化为易于数值计算的第二类Fredholm积分方程．文末给出了算例．数值结果表明,在一定频率范围内,地基表面的位移幅值随激振频率增加而增大,随距离的增大以振荡形式衰减变化．     

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

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

上覆单相弹性层的饱和地基上刚性基础竖向振动的轴对称混合边值问题

运用作者提出的饱和土弹性波动方程,从理论上研究了上覆单相弹性土层的饱和地基上刚性基础的竖向振动轴对称问题,即采用Hankel积分变换技术并按混合边值条件建立了部分饱和地基上刚性基础竖向振动的对偶积分方程,并将其蜕化为完全饱和地基的情形;该对偶积分方程可化为易于数值计算的第二类Fredholm积分方程。文末的算例给出了地基表面动力柔度系数C_v随无量纲频率a₀的变化曲线。     

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

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

无限横观各向同性弹性圆柱中刚性圆盘的扭转振动

分析了无限横观各向同性弹性圆柱中，刚性圆盘的扭转振动问题．绘制的图形显示了材料的各向异性性质对应力强度因子的影响．     

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

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

无限非均质横观各向同性黏弹性介质中具有圆柱孔洞时的振动

在无限介质中,研究了横截面为圆的柱形孔洞表面上瞬时径向力或扭转引起的扰动,讨论了高阶黏弹性和横观各向同性弹性参数的非均匀性对扰动产生的影响．根据高阶黏弹性Voigt模型,将非零应力分量简化为径向位移分量项表示,这对横观各向同性和高阶黏弹性固体介质是合宜的．导出了含有弹性和黏弹性参数以幂律变化时的应力方程．在瞬时力和扭转边界条件下,求解该方程,求得径向位移分量以及和它相关的应力分量,用修正的Bessel函数项来表示．对瞬时径向力作用问题进行了数值分析,并给出了不同阶的黏弹性和非均质性时的位移和应力变化图形．扭转作用时扰动的数值解可以用类似的方法研究,这里不再深入讨论．     

二维横观各向同性厚板在空间变化热源及体力作用下的动力学响应

研究热源和体力作用下的横观各向同性厚板的二维问题,板的上表面无应力作用,但有规定的表面温度作用;板的下表面置于刚性基础之上,并处于绝热状态.采用Green和Naghdi提出的广义热弹性理论,通过Laplace和Fourier双重变换,在Laplace-Fourier变换域中,得到位移和温度场的控制方程.数值求解双重变换的逆变换,采用一个基于Fourier级数展开的方法,数值地求解Laplace变换的逆变换.对材料镁(Mg)进行数值计算,并用图形表示其结果.推演出各向同性材料铜(Cu)的数值结果,并用图形与横观各向同性材料镁进行比较.同时研究了体力的影响.     

横观各向同性电磁弹性介质中裂纹对SH波的散射

研究横观各向同性电磁弹性介质中裂纹和反平面剪切波之间的相互作用．根据电磁弹性介质的平衡运动微分方程、电位移和磁感应强度微分方程,得到SH波传播的控制场方程．引入线性变换,将控制场方程简化为Helmholtz方程和两个Laplace方程A·D2通过Fourier变换,并采用非电磁渗透型裂面边界条件,得到了柯西奇异积分方程组．利用Chebyshev多项式求解积分方程,得到应力场、电场和磁场以及动应力强度因子的表达,并给出了数值算例．     

3D dynamic responses of a multi-layered transversely isotropic saturated half-space under concentrated forces and pore pressure

Dynamic Green's function plays an important role in the study of various wave radiation, scattering and soil-structure interaction problems. However, little research has been done on the response of transversely isotropic saturated layered media. In this paper, the 3D dynamic responses of a multi-layered transversely isotropic saturated half-space subjected to concentrated forces and pore pressure are investigated. First, utilizing Fourier expansion in circumferential direction accompanied by Hankel integral transform in radial direction, the wave equations for transversely isotropic saturated medium in cylindrical coordinate system are solved. Next, with the aid of the exact dynamic stiffness matrix for in-plane and out-of-plane motions, the solutions for multi-layered transversely isotropic saturated half-space under concentrated forces and pore pressure are obtained by direct stiffness method. A FORTRAN computer code is developed to achieve numerical evaluation of the proposed method, and its accuracy is validated through comparison with existing solutions that are special cases of the more general problems addressed. In addition, selected numerical results for a homogeneous and a layered material model are performed to illustrate the effects of material anisotropy, load frequency, drainage condition and layering on the dynamic responses. The presented solutions form a complete set of Green's functions for concentrated forces (including horizontal load in x(y)-direction, vertical load in z-direction) as well as pore pressure, which lays the foundation for further exploring wave propagation of complex local site in a layered transversely isotropic saturated half-space by using the BEMs.     

Boussinesq indentation of a transversely isotropic half-space reinforced by a buried inextensible membrane

The normal indentation of a rigid circular disk into the surface of a transversely isotropic half-space reinforced by a buried inextensible thin film is addressed. By virtue of a displacement potential function and the Hankel transform, the governing equations of this axisymmetric mixed boundary value problem are represented as a dual integral equation, which is subsequently reduced to a Fredholm integral equation of the second kind. Two important results of the contact stress distribution beneath the disk region as well as the equivalent stiffness of the system are expressed in terms of the solution of the Fredholm integral equation. When the membrane is located on the surface or at the remote boundary, exact closed-form solutions are presented. For the limiting case of an isotropic half-space the results are verified with those available in the literature. As a special case, the elastic fields of a reinforced transversely isotropic half-space under the action of surface axisymmetric patch loads are also given. The effects of anisotropy, embedment depth of the membrane, and material incompressibility on both the contact stress and the normal stiffness factor are depicted in some plots.     

A transversely isotropic medium containing a penny-shaped crack subjected to a non-uniform axisymmetric loading via an anchored smooth rigid disk

A smooth rigid circular anchor disk encapsulated by a penny-shaped crack is embedded in and unbounded transversely isotropic medium. The lamellar rigid disk exerts a nonuniform axisymmetric loading to the upper face of the crack. With the aid of an appropriate stress function and Hankel transform, the governing equations are converted to a set of triple integral equations which in turn are reduced to a Fredholm integral equation of the second kind. For some transversely isotropic materials the normalized stiffness of the system falls well outside of the envelope pertinent to isotropic media. It is shown that mode I stress intensity factor is independent of the material properties and solely depends on the ratio of the radius of the rigid disk to that of the crack; moreover, for the cases where this ratio is less than about 0.9 a simple explicit approximate expression for the mode I stress intensity factor is derived. In contrast, the normalized mode II stress intensity factor is independent of the mentioned geometrical parameters but depends on the elastic properties of the material; depending on the material properties, the normalized mode II stress intensity factor can vary between 0 to ∞ for transversely isotropic materials and between 0 to π/4 for isotropic materials.     

Interaction of a harmonic torsional wave with ring-shaped defects in an elastic body

We have solved the problem of determination of the stressed state in an isotropic elastic body near ring-shaped defects (a crack or a thin rigid inclusion) as a result of the action of a harmonic torsional wave. The method of solution is based on the use of discontinuous solutions of the equation of torsional vibrations and lies in the reduction of the initial boundary-value problems to integral equations for the unknown jumps of angular displacement or tangential stress.