三维有限裂纹体分析

基于胡海昌教授提出的一类不含强奇异积分的新型面力边界积分方程,针对含平面裂纹三维有限体问题给出了数值分析和半解析数值分析裂式。文中对含圆裂纹或裂纹系的厚板、圆柱,以及含垂直于表面有的椭圆裂纹的矩形板等有关问题进行了具体求解,求得了相应的应力强度因子。将文中计算结果与已有结果比较表明,该方法具有较高的精度,而计算量较小。     

三维无限体中两平行平片裂纹相互干扰的超奇异积分方程解法

本文利用三维断裂力学的超奇异积分方程求解理论，对三维无限体中两平行平片裂纹在任意载荷作用下的相互干扰问题作了研究。首先导出了以裂纹面移间断（位借）为未知函数的超奇异积分方程组，然后为其建立了有限积分边界元法；在此基础上，讨论用了裂纹面位移间断计算应力强度因子的方法，最后用此计算了两平行平片裂纹的相对位置对裂前沿应力强度因子的影响，其数值结果令人满意。     

饱和土埋置力源的三维动力Lamb问题解答

基于一组弹性土波动方程，应用Fourier级数展开和Hankel积分变换，得到了三维问题饱和土骨架与孔隙水的应力及位移分量在变换域内的积分形式通解．考虑地基表面透水情形，由边界条件导出了半空间饱和土体在埋置力源作用下的三维动力Lamb问题的解答．给出了埋置水平力作用下地基表面竖向位移、径向位移及周向位移的数值解．该研究为运用边界元法求解饱和地基的动力响应课题奠定了理论基础．     

半平面多边缘裂纹反平面问题的奇异积分方程

利用复变函数和奇异积分方程方法,求解弹性范围内半平面多边缘裂纹的反平面问题．提出了满足半平面边界自由的由分布位错密度表示的单边缘裂纹的基本解,此基本解由主要部分和辅助部分组成．将半平面多边缘裂纹问题看作是许多单边缘裂纹问题的叠加,建立了一组Cauchy型奇异积分方程．然后,利用半开型积分法则求解该奇异积分方程,得到了裂纹端处的应力强度因子．最后,给出了几个数值算例．     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

横观各向同性压电材料中裂纹问题的边界元法

采用Somigiliana公式给出了三维横观各向同性压电材料中的非渗漏裂纹问题的一般解和超奇异积分方程,其中未知函数为裂纹面上的位移间断和电势间断．在此基础上,使用有限部积分和边界元结合的方法,建立了超奇异积分方程的数值求解方法,并给出了一些典型数值算例的应力强度因子和电位移强度因子的数值结果,结果令人满意．     

三维有限体平片裂纹的超奇异积分方程与边界元法

利用Ｓｏｍｉｇｌｉａｎａ公式及有限部积分的概念，导出了含任意平片裂纹三维有限体问题的超奇异积分方程组，并联合使用有限部积分与边界元法，建立了数值求解方法．在裂纹前沿附近单元，采用与理论分析一致的平方根位移模型，以提高数值结果的精度．最后计算了若干典型例子的应力强度因子．     

三维非规则非均匀边界元网格的简便的高精度算法

对三维直接边界元中一阶奇异积分、一阶近奇异积分以及非奇异积分进行统一处理，给出了一种提高积分计算精度的简便有效的方法，对非规则非均匀边界元网格可获得比一般方法高得多的计算精度，非常适合边界形状比较复杂的三维实际问题的边界元分析．     

Reissner型板表面裂纹应力强度因子的线弹簧-不连续位移边界积分方程解法

本文由Reissner型板的不连续位移基本解，根据Betti互换定理，导出了Reissuer型板的不连续位移边界积分方程；结合平面问题的不连续位移边界积分方程─—边界元方法和线弹簧模型，给出了Rrissner型板表面裂纹应力强度因子的线弹簧-不连续位移边界积分方程解法。     

三维有限体中两平行裂纹干扰问题的有限部积分与边界元法

利用Ｓｏｍｉｇｌｉａｎａ公式及有限部积分的概念，导出含两平行平片裂纹三维有限体裂纹干扰问题的超奇异积分方程组，联合使用有限部积分与边界元法，建立了数值求解方法，为提高数值计算结果的精度，在裂纹前疝附近单元，采用平方根位移模型，并在此基础雌出直接计算应力强度因子的公式，最后计算若干典型例子裂纹前沿的应力强度因子。     

Penny-shaped cracks

The Somigliana formula is used to reduce an arbitrary elastic crack problem to a system of three integral equations for the components of displacement discontinuity. For the case of a penny shaped crack situated in an infinite isotropic medium with the crack faces subjected to arbitrary tractions, the integral equations are solved explicitly. In particular integral formulae are obtained for the stresses on the plane of the crack beyond the crack-tip, and hence for the stress intensity factors. The special case of uni-directional shear traction on the crack is examined.     

圆盘状裂纹前缘塑性区尺寸及张开位移估计

将Ｄｕｇｄａｌｅ模型推广到三维裂纹问题计算了圆盘状裂纹前缘塑性区尺寸,并结合断裂力学中的Ｂａｒｅｎｂｌａｔｔ－Ｄｕｇｄａｌｅ裂纹模型和三维Ｊ－积分原理计算了圆盘状裂纹前缘张开位移,得到了Ｊ－积分与裂纹张开位移的关系,最后用非线性有限元方法对圆盘状裂纹的前缘塑性区尺寸作了数值分析,确定了公式中的未知常数,并对其正确性作了数值验证,本文的工作推广了Ｄｕｇｄａｌｅ模型的应用范围。     

An axisymmetric problem of an embedded crack in a graded layer bonded to a homogeneous half-space

In an attempt to simulate non-uniform coating delamination, the elasto-static problem of a penny shaped axisymmetric crack embedded in a functionally graded coating bonded to a homogeneous substrate subjected to crack surface tractions is considered. The coating’s material gradient is parallel to the axisymmetric direction and is orthogonal to the crack plane. The graded coating is modeled as a non-homogeneous medium with an isotropic constitutive law. Using Hankel transform, the governing equations are converted into coupled singular integral equations, which are solved numerically to yield the crack tip stress intensity factors. The Finite Element Method was additionally used to model the crack problem. The main objective of this paper is to study the influence of the material non-homogeneity and the crack position on the stress intensity factors for the purpose of gaining better understanding on the behavior of graded coatings.     

A penny shaped crack in a transversely isotropic material under non-axisymmetric impact loads

A solution is given for problems involving non-axisymmetric dynamic impact loading of a penny shaped crack in a transversely isotropic medium. Laplace and Hankel transforms are used to reduce the equations of elasticity to integral equations, and solutions are obtained for the three modes of fracture. The stress intensity factors are determined for a penny shaped crack loaded by concentrated normal impact forces and concentrated radial shear impact forces. The integral equations are solved by numerical methods, and the results are plotted showing how the dynamic stress intensity factors are influenced by the asymmetric loading.     

Dynamic response of an elastic medium containing crack

A fundamental solution for an infinite elastic medium containing a penny-shaped crack subjected to dynamic torsional surface tractions is attempted. A double Laplace–Hankel integral transform with respect to time and space is applied both to motion equation and boundary conditions yielding dual integral equations. The solution of the derived dual integral equations is based on an analytic procedure using theorems of Bessel functions and ordinary differential equations. The dynamic displacements’ field is obtained by inversion of the corresponding Laplace–Hankel transformed variable. Results of a representative example for a crack subjected to pulse surface tractions are obtained and discussed.     

Integration of Singularities in FE/BE Analyses of Soil-foundation Interaction with Non-homogeneous Elastic Soils

This paper addresses the derivation of the boundary integral equations for a non-homogeneous elastic half-space subjected to constant surface tractions on an arbitrarily shaped area on the basis of the respective Green's functions. The type of non-homogeneity considered is a power law variation of Young's modulus with depth below the surface of the half-space. Two different methods, a contour integral and an integration-free approach are presented, applicable for arbitrarily and rectangular shaped boundary elements, respectively. In the first one the divergence theorem is used in order to reduce the integration of a two-dimensional surface element to an integration over the element's confining boundary only. In the second approach no integration at all is needed since the solution is found simply by evaluating functions to be determined at the boundaries of the loaded rectangle.     

反平面剪切裂纹和刚性线问题和带裂纹圆轴扭转问题中的新型积分方程

如果把通常裂纹问题中奇异积分方程中的右端项由应力改为合力,此时积分方程的核也要由奇异核改为对数型奇异核。文中对于反乎面剪切裂纹和刚性线问题和带裂纹圆轴扭转问题,推导出了这种带对数核的积分方程。     

三维间断位移法及强奇异和超奇异积分的处理方法

从积分方程Somigliana等式出发,导出三维状态下单位位错集度的基本解.在此基础上,建立了边界积分方程,并给出了其离散形式.对强奇异和超奇异积分,采用了Hadamard定义的有限部分积分来处理.最后,给出了计算裂纹应力强度因子的算例,并与解析解进行了比较,证实了该方法的有效性.     

Interaction of general plane P wave and cylindrical inclusion partially debonded from its viscoelastic matrix

The interaction of a general plane P wave and an elastic cylindrical inclusion of infinite length partially debonded from its surrounding viscoelastic matrix of infinite extension is investigated. The debonded region is modeled as an arc-shaped interface crack between inclusion and matrix with non-contacting faces. With wave functions expansion and singular integral equation technique, the interaction problem is reduced to a set of simultaneous singular integral equations of crack dislocation density function. By analysis of the fundamental solution of the singular integral equation, it is found that dynamic stress field at the crack tip is oscillatory singular, which is related to the frequency of incident wave. The singular integral equations are solved numerically, and the crack open displacement and dynamic stress intensity factor are evaluated for various incident angles and frequencies. The project supported by the National Natural Science Foundation of China (19872002) and Climbing Foundation of Northern Jiaotong University