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.     

内部压力作用下矩形板中源于椭圆孔的分支裂纹应力强度因子的一种数值分析

应用一种边界元方法来研究内部压力作用下矩形板中源于椭圆孔的分支裂纹。该边界元方法由Crouch与Starfied建立的常位移不连续单元和笔者最近提出的裂尖位移不连续单元构成。在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界。本数值结果进一步证实这种数值方法对计算有限大板中复杂裂纹的应力强度因子的有效性,同时该数值结果可以揭示裂纹体几何对应力强度因子的影响。     

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.     

A numerical analysis of perpendicular cracks under general in-plane loading with a hybrid displacement discontinuity method

In this paper, a numerical analysis of perpendicular cracks under general in-plane loading is performed by using a hybrid displacement discontinuity method which consists of the non-singular displacement discontinuity element presented by Crouch and Starfied and the crack tip displacement discontinuity elements by the author. In the boundary element implementation the left or the right crack tip displacement discontinuity element is placed locally at corresponding left or right crack tip on top of the ordinary non-singular displacement discontinuity elements that cover the entire crack surface and the other boundary. The present numerical results show that the numerical approach is simple, yet very accurate for calculating numerically stress intensity factors for perpendicular cracks under general in-plane loading.     

The hybrid crack-tip element approach to thermo-elastic cracks

The stress intensity factors for a thermal crack problem are computed directly from the crack tip singular field as it was embedded in a specialized crack-tip element. The method gives accurate results without modeling a very refined mesh near the crack tip and without evaluating additional line and/or area integrals for the “modified” path independent integrals as suggested by most current methods for thermoelastic cracks.     

双轴载荷作用下源于椭圆孔的分支裂纹的一种边界元分析

利用一种边界元方法来研究双轴载荷作用下无限大板中源于椭圆孔的分支裂纹．该边界元方法由Crouch与Starfied建立的常位移不连续单元和笔者提出的裂尖位移不连续单元构成．在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界,文中算例说明本数值方法对计算平面弹性裂纹的应力强度因子是非常有效的。该文对双轴载荷作用下无限大板中源于椭圆孔的分支裂纹的数值结果进一步证实本数值方法对计算复杂裂纹的应力强度因子的有效性,同时该数值结果可以揭示双轴载荷及裂纹体几何对应力强度因子的影响。     

A Self-Similar Crack Expansion method for three-dimensional cracks in an infinite or semi-infinite medium

The Self-Similar Crack Expansion (SSCE) method is used to calculate stress intensity factors for three-dimensional cracks in an infinite medium or semi-infinite medium by the boundary integral element technique, whereby, the stress intensity factors at crack tips are determined by calculating the crack-opening displacements over the crack surface. For elements on the crack surface, regular integrals and singular integrals are precisely evaluated based on closed form expressions, which improves the accuracy. Examples show that this method yields very accurate results for stress intensity factors of penny-shaped cracks and elliptical cracks in the full space, with errors of less than 1% as compared with analytical solutions. The stress intensity factors of subsurface cracks are in good agreement with other analytical solutions.     

层状磁电复合材料界面共线裂纹平面问题分析

本文研究了面内电磁势载荷作用下双层压电压磁复合材料中共线界面裂纹问题．考虑了压电材料的导磁性质和压磁材料的介电性质,引入了界面电位移和磁感强度的连续性条件．利用Fourier 变换得到一组第二类Cauchy 型奇异积分方程．进一步导出了相应问题的应力强度因子、电位移强度因子和磁感强度强度因子的表达式,给出了应力强度因子的数值结果．结果表明电磁载荷会导致界面裂纹尖端I、II 混合型应力奇异性,同时还伴随着电位移和磁感强度的奇异性．比较了双裂纹左右端的应力强度因子,发现在面内极化方向上施加面内磁势载荷时共线裂纹内侧尖端区域的两个法向应力场发生互相干涉增强．     

准脆性材料断裂过程区中的圆饼微裂纹损伤计算

准脆性工程材料及结构在外力作用下,不仅引起内部缺陷变化和微裂纹的出现及发展,且使得其结构承载能力降低或性能劣化．在其材料失效过程中常存在裂缝与断裂损伤过程区．为研究材料细观缺陷或微裂纹与宏观破坏的规律,通过细观力学方法,对于代表性体积单元RVE中的圆饼型微裂纹的尺寸与密度变化,探讨其宏观断裂过程区力学参量与损伤之间的量化关系．借助宏观断裂过程区的黏聚裂纹模型,将损伤单元RVE嵌入到宏观裂缝端部的断裂过程区中,对其进行联接细观损伤到宏观破坏的力学多尺度研究．文中也通过实验数据,对其理论计算结果进行了算例的讨论与分析．     

Brittle fracture dynamics with arbitrary paths I. Kinking of a dynamic crack in general antiplane loading

We present a new method for determining the elasto-dynamic stress fields associated with the propagation of anti-plane kinked or branched cracks. Our approach allows the exact calculation of the corresponding dynamic stress intensity factors. The latter are very important quantities in dynamic brittle fracture mechanics, since they determine the crack path and eventual branching instabilities. As a first illustration, we consider a semi-infinite anti-plane straight crack, initially propagating at a given time-dependent velocity, that changes instantaneously both its direction and its speed of propagation. We will give the explicit dependence of the stress intensity factor just after kinking as a function of the stress intensity factor just before kinking, the kinking angle and the instantaneous velocity of the crack tip.     

用光弹性法测双层复合板的应力强度因子

本文用光弹性法分析了双层复合材料结构的裂纹问题，研究了两种不同材料的存在对应力强度因子的影响，以及讨论了由于裂纹不断加深，应力强度因子的变化规律，并与计算的结果作了比较，在裂纹不接近界面情况下，两者结果吻合的较好。     

压电薄板板边导电裂纹尖端的力电场分析

采用有限元方法，分析了压电薄板板边不同长度导电裂纹尖端的力电场分布规律，发现导电裂纹尖端的应力场和电场强度存在明显的集中和奇异现象，集中和奇异的程度与裂纹长度有关。而且，在裂纹延长线上分别存在两点，这里的应力和电场对裂纹长度不太敏感，总等于无裂纹时薄板的均匀应力和均匀电场强度；同时，还研究了导电裂纹尖端的应力强度因子和电场强度因子对裂纹长度的依赖关系，发现在线性本构的前提下，导电裂纹尖端的应力强度因子与电场强度因子之间具有近似的线性关系。     

基于权函数法的核主泵飞轮应力强度因子研究

采用权函数法确定了含裂纹飞轮在离心力和接触压力作用下的应力强度因子,计算了在同步转速工况下裂纹尖端应力强度因子的值,并与解析法和有限元法计算结果进行了对比。结果表明:权函数法与解析法的误差在3%以内,与有限元法的误差在1%以内,验证了权函数法计算应力强度因子的准确性高;在结构不变的情况下,权函数法可以快速求解不同载荷条件、不同长度裂纹的应力强度因子。通过控制变量法研究了不同参数对应力强度因子的影响,结果表明,飞轮裂纹尖端总应力强度因子随着裂纹尺寸、旋转转速、飞轮尺寸外径与内径比值的增大而增大。     

Variation of the stress intensity factor along the crack front of interacting semi-elliptical surface cracks

Summary  In this study, the interaction between two semi-elliptical co-planar surface cracks is considered when Poisson's ratio ν = 0.3. The problem is formulated as a system of singular integral equations, based on the idea of the body force method. In the numerical calculation, the unknown density of body force density is approximated by the product of a fundamental density function and a polynomial. The results show that the present method yields smooth variations of stress intensity factors along the crack front very accurately, for various geometrical conditions. When the size of crack 1 is larger than the size of crack 2, the maximum stress intensity factor appears at a certain point, β₁=177^∘, of crack 1. Along the outside of crack 1, that is at β₁=0∼90^∘, the interaction can be negligible even if the two cracks are very close. The interaction can be negligible when the two cracks are spaced in such a manner that their two closest points are separated by a distance exceeding the small crack's major diameter. The variations of stress intensity factor of a semi-elliptical crack are tabulated and charted. Received 30 August 1999; accepted for publication 22 February 2000     

SELF-SIMILAR CRACK EXPANSION METHOD FOR THE STRESS INTENSITY FACTOR ANALYSIS

The Self-Similar Crack Expansion (SSCE) method is proposed to evaluate stress intensi-ty factors at crack tips, whereby stress intensity factors of a crack can be determined by the crackopening displacement over the crack, not just by the local displacement around the crack tip. The crackexpansion rate is estimated by taking advantage of the crack self-similarity. Therefore, the accuracy ofthe calculation is improved. The singular integrals on crack tip elements are also analyzed and are pre-cisely evaluated in terms of a special integral analysis. Combination of these two techniques greatly in-creases the accuracy in estimating the stress distribution around the crack tip. A variety of two-dimen-sional cracks, such as subsurface cracks, edge cracks, and their interactions are calculated in terms ofthe self-similar expansion rate. Solutions are satisfied with errors less than 0.5% as compared with theanalytical solutions. Based on the calculations of the crack interactions, a theory for crack interactionsis proposed such that for a group of aligned cracks the summation of the square of SIFs at the right tipsof cracks is always equal to that at the left tips of cracks. This theory was proved by the mehtod ofSelf-Similar Crack Expansion in this paper.     

On the higher order asymptotic analysis of a non-uniformly propagating dynamic crack along an arbitrary path

Transient mixed-mode elastodynamic crack growth along arbitrary smoothly varying paths is considered. Asymptotically, the crack tip stress field is square root singular with the angular variation of the singular term depending weakly on the instantaneous values of the crack tip speed and on the mode-I and mode-II stress intensity factors. However, for a material particle at a small distance away from the moving crack tip, the local stress field will depend not only on the instantaneous values of the crack tip speed and stress intensity factors, but also on the past history of these time dependent quantities. In addition, for cracks propagating along curved paths the stress field is also expected to depend on the nature of the curved crack path. Here, a representation of the crack tip fields in the form of an expansion about the crack tip is obtained in powers of radial distance from the tip. The higher order coefficients of this expansion are found to depend on the time derivative of crack tip speed, the time derivatives of the two stress intensity factors as well as on the instantaneous value of the local curvature of the crack path. It is also demonstrated that even if cracks follow a curved path dictated by the criterion K ₁₁ ^d =0, the stress field may still retain higher order asymmetric components related to non-zero local curvature of the crack path.     

A crack perpendicular to the bimaterial interface in finite solid

The dislocation simulation method is used in this paper to derive the basic equations for a crack perpendicular to the bimaterial interface in a finite solid. The complete solutions to the problem, including the T stress and the stress intensity factors are obtained. The stress field characteristics are investigated in detail. It is found that when the crack is within a weaker material, the stress intensity factor is smaller than that in a homogeneous material and it decreases when the distance between the crack tip and interface decreases. When the crack is within a stiffer material, the stress intensity factor is larger than that in a homogeneous material and it increases when the distance between the crack tip and interface decreases. In both cases, the stress intensity factor will increase when the ratio of the size of a sample to the crack length decreases. A comparison of stress intensity factors between a finite problem and an infinite problem has been given also. The stress distribution ahead of the crack tip, which is near the interface, is shown in details and the T stress effect is considered.     

The interaction of doubly periodic cracks

In this paper, an extremely accurate and efficient method for computing the interaction of a set of or multiple sets of general doubly periodic cracks has been presented on the basis of superposition principle, pseudo-traction method, and isolating analysis technique. A great number of typical examples are given in this paper. The stress intensity factors (SIF), the minimum strain energy density factors (SED) of crack tip and the critical stress (CRS) of crack growth are calculated with the accuracy of six significant digits for the rectangularly distributed periodic cracks and five significant digits for the general doubly periodic cracks. The relation of the interaction effect of the double periodic cracks with the periods and the ratio of crack length to crack spacing is analyzed. Also in this paper, the key technique problems for this method are discussed.     

A numerical analysis of cracks emanating from an elliptical hole in a 2-D elasticity plate

This paper is concerned with the stress intensity factors (SIFs) of cracks emanating from an elliptical hole in an infinite or a finite plate under biaxial loads by using a boundary element method, which consists of the non-singular displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements due to the author. In the boundary element implementation the left or the right crack-tip element is placed locally at the corresponding left or right crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. A few numerical examples are included to show that the present approach is very efficient and accurate for the calculating the SIFs of crack problems in an infinite or a finite plate. The present numerical results of cracks emanating from an elliptical hole under biaxial loads can reveal the effect of the elliptical aspect ratio and the transverse load on the SIFs.     

Special approach for the determination of fracture mechanical parameters at an interface crack tip

An interface crack of finite length is considered between two semi-infinite planes with an artificial contact zone at one of the two crack tips. A transcendental equation and certain simple asymptotic formulas are established for the real contact zone (in the Comninou-Dundurs sense) in terms of the stress intensity factors (SIFs) of the considered model. In these terms analytical expressions are also provided for the energy release rate and for the SIF of the classical interface crack model with an oscillating singularity at the crack tip. The appropriate length of the artifical contact zone is shown to be attainable on the basis of the analysis of the stresses at the crack tip. The use of the proposed model is suggested for integrity assessment of inhomogeneous structural elements of composites containing interface cracks. Received 26 March 1997; accepted for publication 12 September 1997