轴对称环形片状界面裂纹问题分析

讨论受拉伸载荷作用的轴对称环形片状界而裂纹问题．该问题归结为求解一组超奇异积分-微分方程．方程中的未知位移间断近似表示为基本密度函数与多项式之积,其中基本密度函数考虑到问题的对称性用二维界面裂纹精确解表示．在圆形片状裂纹的情况下,数值结果与现有理论解作比较的结果表明,数值结果与相应界面圆形片状裂纹和均质体圆形片状裂纹的精确解均吻合得很好．文中以图表形式给出应力强度因子与材料组合和几何条件之间的关系．     

平行于异材界面的三维椭圆平片裂纹分析

讨论了拉伸载荷作用下平行于两相材料界面的椭圆平片裂纹问题．首先,使用有限部积分概念和两相材料界面完全接合时的点力基本解导出了一组以裂纹表面位移差为未知函数的超奇异积分方程组．该组方程表明,此时三种裂纹模型同时存在;其次,在数值求解该组方程的过程中,未知函数裂纹表面位移差被近似为位移差的基本密度函数与多项式之积．基本密度函数反映了裂纹前沿应力奇性性态;最后,以拉伸载荷为例,讨论了椭圆平片裂纹与界面的距离、裂纹形状比和不同材料组合对应力强度因子的影响,并以图表形式给出。     

Multi-interface cracks in a piezoelectric layer bonded to two half-piezoelectric spaces under anti-plane shear loading

The behavior of four parallel symmetry permeable interface cracks in a piezoelectric layer bonded to two half-piezoelectric spaces under anti-plane shear loading is investigated. By using the Fourier transform, the problem can be solved with the help of two pairs of triple integral equations. These equations are solved by the Schmidt method. This process is quite different from that papers adopted previously. The normalized stress and electrical displacement intensity factors are determined for different geometric and property parameters for permeable crack surface conditions. Numerical examples are provided to show the effect of the geometry of the interacting cracks, the thickness and the materials constants of the piezoelectric layer upon the stress and electric displacement intensity factors of the cracks. It is found that the electric displacement intensity factors for the permeable crack surface conditions are much smaller than the results for the impermeable crack surface conditions.     

Numerical solutions of singular integral equations for planar rectangular interfacial crack in three dimensional bimaterials

Stress intensity factors for a three dimensional rectangular interfacial crack were considered using the body force method. In the numerical calculations, unknown body force densities were approximated by the products of the fundamental densities and power series; here the fundamental densities are chosen to express singular stress fields due to an interface crack exactly. The calculation shows that the numerical results are satisfied. The stress intensity factors for a rectangular interface crack were indicated accurately with the varying aspect ratio, and bimaterial parameter.     

Interaction of three interfacial Griffith cracks between bonded dissimilar orthotropic half planes

The paper deals with the interaction of a pair of outer cracks on a central crack situated at the interface of two dissimilar orthotropic half-planes. The mixed boundary value problem is reduced to solving a pair of simultaneous singular integral equations which have finally been solved numerically by using Jacobi polynomials. The analytical expressions for stress intensity factors at the central crack tip and the expression of the strain energy release rate have been derived for general loading. Numerical values of the interaction effects of the outer cracks on the central crack have been calculated through stress magnification factors. It is seen that the interaction effects are either shielding or amplification depending on the size of the outer cracks and their spacing from the central crack.     

Interaction of elastic waves with a periodic array of collinear inplane cracks

A boundary integral equation method is applied to the study of the interaction of plane elastic waves with a periodic array of collinear inplane cracks. Numerical results are presented for the dynamic stress intensity factors. The effects of the wave type, wave frequency, wave incidence angle, and crack spacing on the dynamic stress intensity factors are analyzed in detail. The project supported by the Committee of Science and Technology of Shanghai and Tongji University     

Isotropic/orthotropic two-layered strips containing cracks terminating at and going through an interface

Crack problems for isotropic/orthotropic two-layered strips have been investigated. A system of two singular integral equations can be derived by using Fourier integral transformation and boundary conditions of crack problems. After stress singularities at crack tips or other special points are determined for internal and edge cracks, and for cracks terminating at and going through the interface, the system of singular integral equations is solved numerically by Gauss-Jacobi or Gauss-Chebyshev integration formulas for stress intensity factors at the tips and other singular points of cracks. Finally, possible crack growth behavior for cracks approaching and going through the interface is discussed.     

The solution of integral equations with strongly singular kernels applied to the torsion of cracked circular cylinder

In this paper,the functions of warping displacement interruption defined on the crack lines are taken for the fundamental unknown functions.The torsion problem of cracked circular cylinder is reduced to solving a system of integral equations with strongly singular kernels.Using the numerical method of these equations,the torsional rigidities and the stress intensity factors are calculated to solve the torsion problem of circular cylinder with star-type and other different types of cracks.The numerical results are satisfactory.     

界面裂纹问题中的权函数方法

本文将Ｐａｒｉｓ等确定均匀材料中裂纹尖端应力强度因子的权函数方法推广应用到界面裂纹问题，给出了界面裂纹尖端附近或无限大体半无限界面裂纹问题的权函数的显式表达式。利用此权函数表达式可以很简便地求解界面裂纹尖端附近一些外来作用引起的应力强度因子，比如任意分布力、相变应变、位错和热等。作为一个算例，本文计算了界面一侧一个刃型位错引起的应力强度因子。     

三维体内含垂直于双材料界面混合型裂纹的超奇异积分解法

利用双材料位移基本解和Somigliana公式,将三维体内含垂直于双材料界面混合型裂纹问题归结为求解一组超奇异积分方程。使用主部分析法,通过对裂纹前沿应力奇性的分析,得到用裂纹面位移间断表示的应力强度因子的计算公式,进而利用超奇异积分方程未知解的理论分析结果和有限部积分理论,给出了超奇异积分方程的数值求解方法。最后,对典型算例的应力强度因子做了计算,并讨论了应力强度因子数值结果的收敛性及其随各参数变化的规律。     

Analysis of stress intensity factors of a ring-shaped interface crack

In this paper, numerical solutions of singular integral equations are discussed in the analysis of axi-symmetric interface cracks under torsion and tension. The problems of a ring-shaped interface crack are formulated in terms of a system of singular integral equations on the basis of the body force method. In the numerical analysis, unknown body force densities are approximated by the products of the fundamental density functions and power series, where the fundamental densities are chosen to express a two-dimensional interface crack exactly. The accuracy of the present analysis is verified by comparing the present results with the results obtained by other researchers for the limiting cases of the geometries. The calculation shows that the present method gives rapidly converging numerical results for those problems as well as for ordinary crack problems in homogeneous material. The stress intensity factors of a ring-shaped interface crack are shown in tables and charts with varying the material combinations and also geometrical conditions.     

Dynamic stress intensity factors for two parallel interface cracks between a nonhomogeneous bonding layer and two dissimilar elastic half-planes subject to an impact load

Some composite materials are constructed of two dissimilar half-planes bonded by a nonhomogeneous elastic layer. In the present study, a crack is situated at the interface between the upper half-plane and the bonding layer of such a material, and another crack is located at the interface between the lower half-plane and the bonding layer. The material properties of the bonding layer vary continuously from those of the lower half-plane to those of the upper half-plane. Incoming shock stress waves impinge upon the two interface cracks normal to their surfaces. Fourier transformations were used to reduce the boundary conditions for the cracks to two pairs of dual integral equations in the Laplace domain. To solve these equations, the differences in the crack surface displacements were expanded in a series of functions that are zero-valued outside the cracks. The unknown coefficients in the series were solved using the Schmidt method so as to satisfy the conditions inside the cracks. The stress intensity factors were defined in the Laplace domain and were inverted numerically to physical space. Dynamic stress intensity factors were calculated numerically for selected crack configurations.     

Investigation of the behavior of a griffith crack at the interface between two dissimilar orthotropic elastic half-planes for the opening crack mode

The behaviors of an interface crack between dissimilar orthotropic elastic halfplanes subjected to uniform tension was reworked by use of the Schmidt method. By use of the Fourier transform, the problem can be solved with the help of two pairs of dual integral equations, of which the unknown variables are the jumps of the displacements across the crack surfaces. Numerical examples are provided for the stress intensity factors of the cracks. Contrary to the previous solution of the interface crack, it is found that the stress singularity of the present interface crack solution is of the same nature as that for the ordinary crack in homogeneous materials. When the materials from the two half planes are the same, an exact solution can be otained.     

Scattering of SHPlease define abbreviation SH waves by arc-shaped interface cracks between a cylindrical magneto-electro-elastic inclusion and matrix: near fields

Summary In this paper, the scattering of SH waves by a magneto-electro-elastic cylindrical inclusion partially debonded from its surrounding magneto-electro-elastic material is investigated by using the wavefunction expansion method and a singular integral equation technique. The debonding regions are modeled as multiple arc-shaped interface cracks with non-contacting faces. The magneto-electric impermeable boundary conditions are adopted. By expressing the scattered fields as wavefunction expansions with unknown coefficients, the mixed boundary-value problem is firstly reduced to a set of simultaneous dual-series equations. Then, dislocation density functions are introduced as unknowns to transform these dual-series equations to Cauchy singular integral equations of the first type,which can be numerically solved easily. The solution is valid for arbitrary number and size of the arc-shaped interface cracks. Finally, numerical results of the dynamic stress intensity factors are presented for the cases of one debond. The effects of incident direction, crack configuration and various material parameters on the dynamic stress intensity factors are discussed. The solution of this problem is expected to have applications in the investigation of dynamic fracture properties of magneto-electro-elastic materials with cracks.The work was supported by the National Natural Science Fund of China (Project No. 19772029) and the Research Fund for Doctors of Hebei Province, China (Project No. B2001213).     

Uniformly moving screw dislocation interacting with interface cracks in anisotropic bimaterials

This paper attempts to investigate the problem for the interaction between a uniformly subsonic moving screw dislocation and interface cracks in two dissimilar anisotropic materials. Using Riemann–Schwarz’s symmetry principle integrated with the analysis singularity of complex functions, we present the general elastic solutions of this problem and the closed form solutions for interface containing one and two cracks. The expressions of stress intensity factors at the crack tips and image force acting on moving dislocation are derived explicitly. The results show that the stress intensity factors at the crack tips decrease with increasing velocity of dislocation, and larger dislocation velocity leads to the equilibrium position of dislocation leaving from crack tips. The presented solutions contain previously known results as the special cases.     

Dynamic stress intensity factors of two collinear mode-III cracks perpendicular to and on the two sides of a bi-FGM weak-discontinuous interface

The mechanical model was established for the anti-plane dynamic fracture problem for two collinear cracks on the two sides of and perpendicular to a weak-discontinuous interface between two materials with smoothly graded elastic properties, as opposed to a sharp interface with discontinuously changing elastic properties. The problem was reduced as a system of Cauchy singular integral equations of the first kind by Laplace and Fourier integral transforms. The integral equations were solved by Erdogan's collocation method and the dynamic stress intensity factors in the time domain were obtained through Laplace numerical inversion proposed by Miller and Guy. The influences of geometrical and physical parameters on the dynamic stress intensity factors were illustrated and discussed, based on which some conclusions were drawn: (a) to increase the thickness of the FGM strip on either side of the interface will be beneficial to reducing the DSIF of a crack perpendicular to a bi-FGM interface and embedded at the center of one of the FGM strips; (b) To increase the rigidity of the FGM strip where the crack is located will increase the DSIF. However, when the material in one side of the interface is more rigid, the DSIF of the interface-perpendicular embedded crack in the other side will be reduced; (c) To decrease the weak-discontinuity of a bi-FGM interface will not necessarily reduce the stress intensity factor of a crack perpendicular to it, which is different from the case of interfacial crack; (d) For two collinear cracks with equal half-length, when the distance between the two inner tips is less than about three times of the half-length, the interaction of them is intensified, however, when the distance is greater than this the interaction becomes weak.     

界面下主微裂纹间干涉的研究

涉及两相正交各向异性体界面干涉问题的研究,多裂纹问题被分解为只含单裂纹的子问题,利用位错理论和裂面应力自由条件,列出一组可数值求解位错密度函数的奇异积分方程,从耐 注得应力强度因子。     

Two mode III internal cracks located within two bonded functionally graded piezoelectric half planes respectively

This paper studies the mode III crack problem of two bonded functionally graded piezoelectric half planes which contain a crack respectively. These two cracks are located normal to the interface. All the material properties are assumed to vary along the direction of the crack line. A system of singular integral equations for electrically impermeable and permeable cracks is derived and solved numerically by using the Gauss–Chebyshev integration formula. The influence of the nonhomogeneous parameters and the dependence of the crack interactions on the stress and electric displacement intensity factors are investigated.     

Analysis of three-dimensional edge cracks under tensile loading

Three-dimensional edge cracks are analyzed using the Self-Similar Crack Expansion (SSCE) method with a boundary integral equation technique. The boundary integral equations for surface cracks in a half space are presented based on a half space Green's function (Mindlin, 1936). By using the SSCE method, the stress intensity factors are determined by the crack-opening displacement over the crack surface. In discrete boundary integral equations, the regular and singular integrals on the crack surface elements are evaluated by an analytical method, and the closed form expressions of the integrals are given for subsurface cracks and edge crakcs. This globally numerical and locally analytical method improves the solution accuracy and computational effort. Numerical results for edge cracks under tensile loading with various geometries, such as rectangular cracks, elliptical cracks, and semi-circular cracks, are presented using the SSCE method. Results for stress intensity factors of those surface breaking cracks are in good agreement with other numerical and analytical solutions.     

横观各向同性材料的三维断裂力学问题

从三维横观各向同性材料弹性力学理论出发, 使用Hadamard有限部积分概念, 导出了三维状态下单位位移间断(位错)集度的基 本解. 在此基础上, 进一步运用极限理论, 将任意载荷作用下, 三维无限大横观各向 同性材料弹性体中, 含有一个位于弹性对称面内的任意形状的片状裂纹问题, 归结为求 解一组超奇异积分方程的问题. 通过二维超奇异积分的主部分析方法, 精确地求得了裂纹前沿光滑点附近的应力奇异指数和奇异应力场, 从而找到了以裂纹表面位移间断表示的应力强度因子表达式及裂纹局部扩展所提供 的能量释放率. 作为以上理论的实际应用,最后给出了一个圆形片状裂纹问题 的精确解例和一个正方形片状裂纹问题的数值解例. 对受轴对称法向均布载荷作用下圆形片状裂纹问题, 讨论了超奇异积分方程的精确求解方法, 并获得了位移间断和应力强度因子的封闭解, 此结果与现有理论解完全一致.