Controlling parameter of the stress intensity factors for a planar interfacial crack in three-dimensional bimaterials

Numerical solutions of singular integral equations are discussed in the analysis of a planar rectangular interfacial crack in three-dimensional bimaterials subjected to tension. The problem is formulated as 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 density functions are chosen to express singular behavior along the crack front of the interface crack exactly. The calculation shows that the present method gives smooth variations of stress intensity factors along the crack front for various aspect ratios. The present method gives rapidly converging numerical results and highly satisfied boundary conditions throughout the crack boundary. The stress intensity factors are given with varying the material combination and aspect ratio of the crack. It is found that the stress intensity factors K_I and K_II are determined by the bimaterial constant ε alone, independent of elastic modulus ratio and Poisson’s ratio.     

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.     

Variation of the stress intensity factor along the front of a 3-D rectangular crack subjected to mixed-mode load

Summary  The singular integral equation method is applied to the calculation of the stress intensity factor at the front of a rectangular crack subjected to mixed-mode load. The stress field induced by a body force doublet is used as a fundamental solution. The problem is formulated as a system of integral equations with r ⁻³-singularities. In solving the integral equations, unknown functions of body-force densities are approximated by the product of polynomial and fundamental densities. The fundamental densities are chosen to express two-dimensional cracks in an infinite body for the limiting cases of the aspect ratio of the rectangle. The present method yields rapidly converging numerical results and satisfies boundary conditions all over the crack boundary. A smooth distribution of the stress intensity factor along the crack front is presented for various crack shapes and different Poisson's ratio. Received 5 March 2002; accepted for publication 2 July 2002     

Stress intensity factors of a rectangular crack meeting a bimaterial interface

Using the hypersingular integral equation method based on body force method, a planar crack meeting the interface in a three-dimensional dissimilar materials is analyzed. The singularity of the singular stress field around the crack front terminating at the interface is analyzed by the main-part analytical method of hypersingular integral equations. Then, the numerical method of the hypersingular integral equation for a rectangular crack subjected to normal load is proposed by the body force method, which the crack opening dislocation is approximated by the product of basic density functions and polynomials. Numerical solutions of the stress intensity factors of some examples are given.     

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.     

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

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

Stress intensity factors for interface crack in finite size specimen using a generalized variational approach

A generalized variational approach together with eigenfunction expansion is applied to determine the stress intensity factors for interface crack in finite size specimen. Application is also made of the complex potentials such that a complex stress intensity factor with components corresponding to the Mode I and II stress intensity factors can be identified with one of the leading coefficients in the eigenfunction expansion. Obtained are the numerical values of the stress intensity factors for an interface edge crack in a bimaterial rectangular specimen. The outside boundary is subjected to uniform stress normal and parallel to the crack. Solutions are also obtained for the same crack aand specimen geoinetry is subjected to a pair of equal and opposite concentrated forces along the open end away from the edge crack. The third example pertains to the case of three-point bending where the centre concentrated load is directed along the interface dividing the two materials. Numerical results are obtained for four different combinations of the bimaterial specimen with an interface edge crack.     

双材料界面裂纹应力强度因子的边界元分析

采用双材料基本解建立边界元法基本方程,计算双材料界面裂纹尖端附近的应用力和位移场。不离散界面,并设置面力奇异四分之一点裂尖单元以提高计算精度。数值结果表明,本文的方法具有较高的精度和效率。     

TENSION OF A CLAMPED RECTANGULAR PLATE CONTAINING A CENTRAL CRACK

Using the fundamental solution of a single crack and the Fourier transform solution of an infinite strip, the tension problem of a clamped rectangular plate containing a central crack is reduced to solve a system of singular integral equations. Then, the normal stress on clamped side and the stress intensity factars of the central cruck are carried out by means of Gauss-Jacobi integration formulas. The comparison of numerical results is shown in the "table of stress intensity factars".     

An effective boundary element method for analysis of crack problems in a plane elastic plate

A simple and effective boundary element method for stress intensity factor calculation for crack problems in a plane elastic plate is presented. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements proposed by YAN Xiangqiao. In the boundary element implementation the left or the right crack-tip displacement discontinuity element was placed locally at the corresponding left or right each crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Test examples (i. e. , a center crack in an infinite plate under tension, a circular hole and a crack in an infinite plate under tension) are included to illustrate that the numerical approach is very simple and accurate for stress intensity factor calculation of plane elasticity crack problems. In addition, specifically, the stress intensity factors of branching cracks emanating from a square hole in a rectangular plate under biaxial loads were analysed. These numerical results indicate the present numerical approach is very effective for calculating stress intensity factors of complex cracks in a 2-D finite body, and are used to reveal the effect of the biaxial loads and the cracked body geometry on stress intensity factors.     

Some basic problems in anisotropic bimaterials with a viscoelastic interface

This research is devoted to the study of anisotropic bimaterials with Kelvin-type viscoelastic interface under antiplane deformations. First we derive the Green’s function for a bimaterial with a Kelvin-type viscoelastic interface subjected to an antiplane force and a screw dislocation by means of the complex variable method. Explicit expressions are derived for the time-dependent stress field induced by the antiplane force and screw dislocation. Also presented is the time-dependent image force acting on the screw dislocation due to its interaction with the Kelvin-type viscoelastic interface. Second we investigate a rectangular inclusion with uniform antiplane eigenstrains embedded in one of the two bonded anisotropic half-planes by virtue of the derived Green’s function for a line force. The explicit expressions for the time-dependent stress field induced by the rectangular inclusion are obtained in terms of the simple logarithmic and exponential integral functions. It is observed that in general the stresses exhibit the logarithmic singularity at the four corners of the rectangular inclusion. Our results also show that when one side of the rectangular inclusion lies on the viscoelastic interface, the interfacial tractions are still regular at the two corners of the inclusion which are located on the interface. Last we address a finite Griffith crack normal to the viscoelastic interface by means of the obtained Green’s function for a screw dislocation. The crack problem is formulated in terms of a resulting singular integral equation which is solved numerically. The time-dependent stress intensity factors at the two crack tips are obtained and some interesting features are discussed.     

Generalized stress intensity factors in the interaction within a rectangular array of rectangular inclusions

Summary To evaluate the mechanical strength of fiber-reinforced composites, it is necessary to consider singular stresses at the end of fibers because they cause crack initiation, propagation, and final failure. A square array of rectangular inclusions under longitudinal tension is considered in this paper. The body-force method is applied to a unit cell region. Then, the problem is formulated as a system of singular integral equations, where the unknown functions are the densities of body forces distributed in infinite plates having the same elastic constants as those of the matrix and inclusions. The unknown functions are expressed as piecewise-smooth functions using power series and two types of fundamental densities which express singular stresses. Generalized stress intensity factors at the corners of inclusions are systematically calculated with varying the shape and spacing of a square array of square and rectangular inclusions.     

Tension of a clamped rectangular plate containing a central crack

Using the fundamental solution of a single crack and the Fourier transform solution of an infinite strip, the tension problem of a clamped rectangular plate containing a central crack is reduced to solve a system of singular integral equations. Then, the normal stress on clamped side and the stress intensity factors of the central crack are carried out by means of Gauss-Jacobi integration formulas. The comparison of numerical results is shown in the “table of stress intensity factors”. This work was supported by the Science Fund of the Academy of Sciences of China.     

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

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

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

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

Full field analysis of an anti-plane interface crack subjected to dynamic body forces

In this study, the transient full field response of an interface crack between two different media subjected to dynamic body force at one material is investigated. For time t < 0, the bimaterial medium is stress free and at rest. At t = 0, a concentrated anti-plane dynamic point loading is applied at the medium as shown in Fig. 1. The total wave field is due to the effect of this point loading and the scattering of the incident waves by the interface crack. An alternative methodology that is different from the conventional superposition method is used to construct the reflected, refracted and diffracted wave fields. A useful fundamental solution is proposed in this study and the full field solution is determined by superposition of the fundamental solution in the Laplace transform domain. The proposed fundamental problem is the problem of applying an exponentially distributed traction (in the Laplace transform domain) on the interfacial crack faces. The Cagniard–de Hoop method of Laplace inversion is used to obtain the transient solution in time domain. Exact transient closed form solutions for stresses and stress intensity factors are obtained. Numerical results for the time history of stresses and stress intensity factors during the transient process are discussed in detail.     

Intensity of singular stress fields causing interfacial debonding at the end of a fiber under pullout force and transverse tension

In this study, singular stress fields at the ends of fibers are discussed by the use of models of rectangular and cylindrical inclusions in a semi-infinite body under pullout force. Those singular stresses have not been discussed yet in the previous studies for pullout problems although they are important for causing interfacial initial debonding. The body force method is used to formulate those problems as a system of singular integral equations where unknowns are densities of the body forces distributed in a semi-infinite body having the same elastic constants as those of the matrix and inclusions. In order to compare the results with the previous solutions, tension problems of a fiber in a semi-infinite body are also considered. Then, generalized stress intensity factors at the corner of rectangular and cylindrical inclusions are systematically calculated for various geometrical conditions with varying the elastic ratio, length, and spacing of the location from edge to inner of the body. The effects of elastic modulus ratio and aspect ratio of inclusion upon the stress intensity factors are discussed for pullout problems.     

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

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

Concentrated impact loading on one face of a semi-infinite crack in an anisotropic elastic material

The response of an unbounded anisotropic elastic body containing a semi-infinite crack subjected to a concentrated impact force on one of the crack faces is studied. An exact solution of the dynamic stress intensity factors is obtained from a linear superposition of the solution of Lamb’s problem and a solution of a dislocation emitting from the crack tip. The stress intensity factors exhibit square-root singularity upon the arrival of the Rayleigh wave at the crack tip. As the Rayleigh wave passes through the crack tip, the stress intensity factors either instantaneously assume the static values or gradually approach to zero. Several numerical examples are given for isotropic, cubic and orthotropic materials.     

The influence of elastic waves on dynamic stress intensity factors (three-dimensional problems)

Summary The fundamental solutions of the displacement discontinuity for three-dimensional problems in Laplace space are deduced in thsi paper. The displacement discontinuity method and the equivalent stress method were combined and used to determine dynamic stress intensity factors for three-dimensional time-dependent crack problems. The stress intensity factors were calcualted for dynamically loaded cracks with rectangular, circular, and elliptical crack fronts. The influence of elasticity waves (in particular surface waves) on the magnitude of the stress intensity factor and on the displacement of the crack surfaces was analysed.On leave from the Central-South University of Technology, Changsha, Hunan Province, P. R. China.