双材料中平片裂纹问题的超奇异积分方程解法

利用三维断裂力学的超奇异积分方程方法,对双材料空间中重直于界面的平片裂纹Ⅰ型问题进行了研究。首先根据双材料空间的弹性力学基本解,使用边界积分方程方法,在有限部积分的意义下导出了以裂纹面位罗间断为未知函数的超奇异积分方程,并为其建立了数值法。在此基础上,讨论了用裂纹面位移问题计算应力强度因子的方法。最后用此计算了几个典型的Ⅰ型下片裂纹问题的应力强度因子,其数值结果令人满意。     

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

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

垂直磁压电材料界面三维裂纹的超奇异积分法

应用有限部积分概念和广义位移基本解,垂直于磁压电双材料界面三维复合型裂纹问题被转 化为求解一组以裂纹表面广义位移间断为未知函数的超奇异积分方程问题. 进而,通过主部 分析法精确地求得裂纹尖端光滑点附近的奇性应力场解析表达式. 然后,通过将裂纹表面 位移间断未知函数表达为位移间断基本密度函数与多项式之积,使用有限部积分法对超奇异 积分方程组建立了数值方法. 最后,通过典型算例计算,讨论了广义应力强度因子的变化规 律.     

各向异性平面含斜裂纹的奇异积分方程方法

本文应用平面弹性复变方法，将无限各向异性平面中的任意斜裂纹问题归结为求解一组解析函数边值问题，通过构造适当的积分变换将边值问题转化为奇异积分方程，进而应用Lobotto-Chebyshev数值求积公式，求出该奇异积分方程的数值解，并得到了应力强度因子的近似表达式，最后，给出了一些实例的数值结果，对特例的数值结果与精确结果进行比较，吻合的很好。     

含圆形嵌体弹性平面中径向裂纹问题的超奇异积分方程方法

根据含圆形嵌体平面问题在极坐标下的弹性力学基本解,使用Betti互换定理,在有限部积分意义下将问题归结为两个以裂纹岸位移间断为基本未知量、对于Ⅰ型和Ⅱ型问题相互独立的超奇异积分方程,对含圆形嵌体弹性平面中的径向裂纹问题进行了研究.根据有限部积分原理,建立了问题的数值算法.计算结果表明,嵌体半径、裂纹位置及材料剪切弹性模量等都对裂纹应力强度因子具有较为明显的影响.     

椭圆平片裂纹前沿应力强度因子解析解的超奇异积分方程方法

对无限大三维均质弹性体中任意平片裂纹的超奇异积分方程,巧妙地引入椭球坐标系和利用裂纹表面位移间断人有平方根的特性,获得了受任意方向均布压力作用下椭圆平片裂纹问题的超奇异积分方程的解析解。运用这些解析解和应力强度因子的定义,得到了裂纹前沿Ⅰ型、Ⅱ型和Ⅲ型和应力强度因子的精确表达式,所得结果与现有精确解完全一致。     

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

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

磁电热弹耦合材料三维多裂纹超奇异积分法

用超奇异积分方程法将多场耦合载荷作用下磁电热弹耦合材料内含任意形状和位置三维多裂纹问题转化为求解一以广义位移间断为未知函数的超奇异积分方程组问题,退化得到内含任意形状平行三维多裂纹问题的超奇异积分方程组;推导出平行三维多裂纹问题的裂纹前沿广义奇异应力场解析表达式、定义了广义(应力、应变能)强度因子和广义能量释放率;应用有限部积分概念及体积力法,为超奇异积分方程组建立了数值求解方法,编制了FORTRAN程序,以平行双裂纹为例,通过典型算例,研究了广义(应力、应变能)强度因子随裂纹位置、裂纹形状及材料参数变化规律,得到裂纹断裂评定准则. 最后,分析了裂纹间干扰、屏蔽作用及其在工程实际中的应用.      

与界面垂直相触的三维裂纹的超奇异积分方程和奇性指数

本文使用有限部积分原理和两相材料空间弹性力学问题的点力基本解导出了与界面垂直相触的三维平片解纹的超奇异积分方程组;     

Analysis of multiple interfacial cracks in three-dimensional bimaterials using hypersingular integro-differential equation method

By using the concept of finite-part integral, a set of hypersingular integro-differential equations for multiple interracial cracks in a three-dimensional infinite bimaterial subjected to arbitrary loads is derived. In the numerical analysis, unknown displacement discontinuities are approximated with the products of the fundamental density functions and power series. The fundamental functions are chosen to express a two-dimensional interface crack rigorously. As illustrative examples, the stress intensity factors for two rectangular interface cracks are calculated for various spacing, crack shape and elastic constants. It is shown that the stress intensity factors decrease with the crack spacing.     

三维动态裂纹问题的超奇异积分方程法

基于弹性材料的动态基本方程,结合广义Betti-Rayleigh互易等式与时域下的边界积分方程,推导得到时域下的超奇异积分方程组。引入Laplace域下的动态基本解,将经过主部分析的积分核函数分解为静态和动态部分,其中动态积分核不具有奇异性。在裂纹前沿附近单元,采用与理论分析一致的平方根位移模型。结合Lubich时间卷积实现拉氏变换,采用配置点法计算超奇异积分,获得问题的数值解。并针对椭圆裂纹算例编写Fortran程序,得到冲击荷载作用下张开型裂纹的动态应力强度因子变化规律,数值结果稳定且收敛速度快。     

INTERFACIAL CRACK ANALYSIS IN THREE-DIMENSIONAL TRANSVERSELY ISOTROPIC BI-MATERIALS BY BOUNDARY INTEGRAL EQUATION METHOD

The integral-differential equations for three-dimensional planar interfacial cracks of arbitrary shape in transversely isotropic bimaterials were derived by virtue of the Somigliana identity and the fundamental solutions, in which the displacement discontinuities across the crack faces are the unknowns to be determined. The interface is parallel to both the planes of isotropy. The singular behaviors of displacement and stress near the crack border were analyzed and the stress singularity indexes were obtained by integral equation method. The stress intensity factors were expressed in terms of the displacement discontinuities. In the non-oscillatory case, the hyper-singular boundary integral-differential equations were reduced to hyper-singular boundary integral equations similar to those of homogeneously isotropic materials.     

A crack perpendicular to and terminating at a bimaterial interface

Using dislocation simulation approach, the basic equation for a finite crack perpendicular to and terminating at a bimaterial interface is formulated. A novel expansion method is proposed for solving the problem. The complete solution to the problem, including the explicit formulae for theT stresses ahead of the crack tip and the stress intensity factors are presented. The stress field characteristics are analysed in detail. It is found that normal stresses {ie27-1} and {ie27-2} ahead of the crack tip, are characterised byQ fields if the crack is within a stiff material and the parameters |p _T | and |q _T | are very small, whereQ is a generalised stress intensity factor for a crack normal to and terminating at the interface. If the crack is within a weak material, the normal stresses {ie27-3} and {ie27-4} are dominated by theQ field plusT stress. This work was supported by the Swedish Research Council for Engineering Sciences.     

Finite-part integral and boundary element method to solve flat crack problems

Using the Somigliana formula and the concepts of finite-part integral, a set of hypersingular integral equations to solve the arbitrary fiat crack in three-dimensional elasticity is derived and its numerical method is then proposed by combining the finitepart integral method with boundary element method. In order to verify the method, several numerical examples are carried out. The results of the displacement discontinuities of the crack surface and the stress intensity factors at the crack front are in good agrernent with the' theoretical solutions.     

Fracture analysis of mode-II crack perpendicular to imperfect bimaterial interface

The problem of a mode-II crack close to and perpendicular to an imperfect interface of two bonded dissimilar materials is investigated.The imperfect interface is modelled by a linear spring with the vanishing thickness.The Fourier transform is used to solve the boundary-value problem and to derive a singular integral equation with the Cauchy kernel.The stress intensity factors near the left and right crack tips are evaluated by numerically solving the resulting equation.Several special cases of the mode-II crack problem with an imperfect interface are studied in detail.The effects of the interfacial imperfection on the stress intensity factors for a bimaterial system of aluminum and steel are shown graphically.The obtained observation reveals that the stress intensity factors are dependent on the interface parameters and vary between those with a fully debonded interface and those with a perfect interface.     

Hypersingular integral equation method for a three-dimensional crack in anisotropic electro-magneto-elastic bimaterials

Using Green’s functions, the extended general displacement solutions of a three-dimensional crack problem in anisotropic electro-magneto-elastic (EME) bimaterials under extended loads are analyzed by the boundary element method. Then, the crack problem is reduced to solving a set of hypersingular integral equations (HIE) coupled with boundary integral equations. The singularity of the extended displacement discontinuities around the crack front terminating at the interface is analyzed by the main-part analysis method of HIE, and the exact analytical solutions of the extended singular stresses and extended stress intensity factors (SIFs) near the crack front in anisotropic EME bimaterials are given. Also, the numerical method of the HIE for a rectangular crack subjected to extended loads is put forward with the extended crack opening dislocation approximated by the product of basic density functions and polynomials. At last, numerical solutions of the extended SIFs of some examples are obtained.     

Method of strongly singular integral equation for the torision of cylinder with edge crack

From the dislocation type solution of the torsion of single crack, by using the concept of finite part integrals, we reduce the torsion problem of cylinder with a single crack into an integral equation with strong singularity. The numerical method is also obtained and several numerical examples are calculated successfully at the end of this paper. Project supported by P. H. D. Foundation of the State Education Commission of China     

Some particular solutions for penny-shaped crack problem by using hypersingular integral equation or differential-integral equation

Summary A hypersingular integral equation or a differential-integral equation is used to solve the penny-shaped crack problem. It is found that, if a displacement jump (crack opening displacement COD) takes the form of (a ²−x ²−y ²)^1/2 x ^m y ⁿ , where a denotes the radius of the circular region, the relevant traction applied on the crack face can be evaluated in a closed form, and the stress intensity factor can be derived immediately. Finally, some particular solutions of the penny-shaped crack problem are presented in this paper. Received 1 July 1997; accepted for publication 13 October 1997