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.     

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.     

3D dynamic stress intensity factors at three rectangular cracks in an infinite elastic medium subjected to a time-harmonic stress wave

Summary Dynamic stresses around three coplanar cracks in an infinite elastic medium are determined in the paper. Two of the cracks are equal, rectangular and symmetrically situated on either side of the centrally located rectangular crack. Time-harmonic normal traction acts on each surface of the three cracks. To solve the problem, two kind of solutions are superposed: one is a solution for a rectangular crack in an infinite elastic medium, and the other one is that for two rectangular cracks in an infinite elastic medium. The unknown coefficients in the combined solution are determined by applying the boundary conditions at the surfaces of the cracks. Finally, stress intensity factors are calculated numerically for several crack configurations. Received 14 July 1998; accepted for publication 2 December 1998     

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

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

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

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

Curve cracks lying along a parabolic curve in anisotropic body

ＣＵＲＶＥＣＲＡＣＫＳＬＹＩＮＧＡＬＯＮＧＡＰＡＲＡＢＯＬＩＣＣＵＲＶＥＩＮＡＮＩＳＯＴＲＯＰＩＣＢＯＤＹＨｕＹｕａｎ－ｔａｉ（胡元太）ＺｈａｏＸｉｎｇ－ｈｕａ（赵兴华）（ＳｈａｎｇｈａｉＵｎｉｖｅｒｓｉｔｙ;ＳｈａｎｇｈａｉＩｎｓｔｉｔｕｔｅｏｆＡ...     

Closed-form weight function for edge crack problems

An accurate crack surface opening displacement formula for edge cracks in finite bodies has been derived, and was subsequently used to develop wide-range closed-form weight functions for edge cracks. The analytical weight function is general and can be easily used to determine accurate stress intensity factors and crack surface opening displacements for arbitrarily loaded edge cracks with high efficiency. Examples have been given for problems of edge crack(s) with a number of typical load conditions for illustration.     

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

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

A computational study of local stress intensity factor solutions for kinked cracks near spot welds in lap-shear specimens

In this paper, the local stress intensity factor solutions for kinked cracks near spot welds in lap-shear specimens are investigated by finite element analyses. Based on the experimental observations of kinked crack growth mechanisms in lap-shear specimens under cyclic loading conditions, three-dimensional and two-dimensional plane-strain finite element models are established to investigate the local stress intensity factor solutions for kinked cracks emanating from the main crack. Semi-elliptical cracks with various kink depths are assumed in the three-dimensional finite element analysis. The local stress intensity factor solutions at the critical locations or at the maximum depths of the kinked cracks are obtained. The computational local stress intensity factor solutions at the critical locations of the kinked cracks of finite depths are expressed in terms of those for vanishing kink depth based on the global stress intensity factor solutions and the analytical kinked crack solutions for vanishing kink depth. The three-dimensional finite element computational results show that the critical local mode I stress intensity factor solution increases and then decreases as the kink depth increases. When the kink depth approaches to 0, the critical local mode I stress intensity factor solution appears to approach to that for vanishing kink depth based on the global stress intensity factor solutions and the analytical kinked crack solutions for vanishing kink depth. The two-dimensional plane-strain computational results indicate that the critical local mode I stress intensity factor solution increases monotonically and increases substantially more than that based on the three-dimensional computational results as the kink depth increases. The local stress intensity factor solutions of the kinked cracks of finite depths are also presented in terms of those for vanishing kink depth based on the global stress intensity factor solutions and the analytical kinked crack solutions for vanishing kink depth. Finally, the implications of the local stress intensity factor solutions for kinked cracks on fatigue life prediction are discussed.     

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.     

Interaction integrals for fracture analysis of functionally graded piezoelectric materials

This paper presents domain form of the interaction integrals based on three independent formulations for computation of the stress intensity factors and electric displacement intensity factor for cracks in functionally graded piezoelectric materials. Conservation integrals of J-type are derived based on the governing equations for piezoelectric media and the crack tip asymptotic fields of homogeneous piezoelectric medium as auxiliary fields. Each of the formulation differs in the way auxiliary fields are imposed in the evaluation of interaction integral and each of them results in a consistent form of the interaction integral in the sense that extra terms naturally appears in their derivation to compensate for the difference in the chosen crack tip asymptotic fields of homogeneous and functionally graded piezoelectric medium. The additional terms play an important role of ensuring domain independence of the presented interaction integrals. Comparison of the numerically evaluated intensity factors through the three consistent formulations with those obtained using displacement extrapolation method is presented by means of two examples.     

Limiting velocity of crack propagation in elastic materials

A crack is represented as a continuous set of linear dislocations. Simple analytical expressions are obtained for the potential and kinetic energies of the environment of moving cracks and the attached mass of cracks for an arbitrary form of the stress applied to the crack P(x). It is shown that the indicated analytical expressions are bilinear integrals of the functions P(x) and ∂P(x)/∂x. These integrals are calculated in explicit form for a constant stress over the entire crack length and the stress due to the action of molecular adhesion forces in a narrow region near the crack openings. It is shown that the calculation results does not depend on the form of molecular adhesion forces. The proposed approach to describing cracks and calculations based on it has made it possible for the first time to obtain a complete analytical expression for the limiting crack propagation velocity in elastic materials as a function of the main mechanical characteristics of such materials. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 4, pp. 158–166, July–August, 2009.     

3D dynamic stress intensity factors around two parallel square cracks in an infinite elastic medium under impact load

Summary  Transient stresses around two parallel cracks in an infinite elastic medium are investigated in the present paper. The shape of the cracks is assumed to be square. Incoming shock stress waves impinge upon the two cracks normal to tzheir surfaces. The mixed boundary value equations with respect to stresses and displacements are reduced to two sets of dual integral equations in the Laplace transform domain using the Fourier transform technique. These equations are solved by expanding the differences in the crack surface displacements in a double series of a function that is equal to zero outside the cracks. Unknown coefficients in the series are calculated using the Schmidt method. Stress intensity factors defined in the Laplace transform domain are inverted numerically to the physical space. Numerical calculations are carried out for transient dynamic stress intensity factors under the assumption that the shape of the upper crack is identical to that of the lower crack. Received 2 February 2000; accepted for publication 10 May 2000     

Dynamic behavior of two unequal parallel permeable interface cracks in a piezoelectric layer bonded to two half piezoelectric materials planes

IntroductionDuetotheintrinsicelectro_mechanicalcouplingbehavior,piezoelectricmaterialsareveryusefulinelectronicdevices.However,mostpiezoelectricmaterialsarebrittlesuchasceramicsandcrystals.Therefore ,piezoelectricmaterialshaveatendencytodevelopcriticalcracksduringthemanufacturingandthepolingprocesses.So ,itisimportanttostudytheelectro_elasticinteractionandfracturebehaviorsofpiezoelectricmaterials.Theincreasingattentiontothestudyofcrackproblemsinpiezoelectricmaterialshasledtoalotofsignificantw…     

THE BEHAVIOR OF TWO PARALLEL SYMMETRIC PERMEABLE CRACKS IN PIEZOELECTRIC MATERIALS

ＩｎｔｒｏｄｕｃｔｉｏｎＩｔｉｓｗｅｌｌ＿ｋｎｏｗｎｔｈａｔｐｉｅｚｏｅｌｅｃｔｒｉｃｍａｔｅｒｉａｌｓｐｒｏｄｕｃｅａｎｅｌｅｃｔｒｉｃｆｉｅｌｄｗｈｅｎｄｅｆｏｒｍｅｄａｎｄｕｎｄｅｒｇｏｄｅｆｏｒｍａｔｉｏｎｗｈｅｎｓｕｂｊｅｃｔｅｄｔｏａｎｅｌｅｃｔｒｉｃｆｉｅｌｄ .Ｔｈｅｃｏｕｐｌｉｎｇｎａｔｕｒｅｏｆｐｉｅｚｏｅｌｅｃｔｒｉｃｍａｔｅｒｉａｌｓｈａｓａｔｔｒａｃｔｅｄｗｉｄｅａｐｐｌｉｃａｔｉｏｎｓｉｎｅｌｅｃｔｒｉｃ＿ｍｅｃｈａｎｉｃａｌａｎｄｅｌｅｃｔｒｉｃｄｅｖｉｃｅｓ,ｓｕｃ…     

Effects of the curvature radius of a cylindrical shell on the stress intensity factors of surface cracks

The effects of the curvature radius of a cylindrical shell on stress intensity factors are investigated in circumferential (inner and outer) semielliptical surface cracks in a cylindrical shell. What is new in this paper is to have given: (1) The stress intensity factors for surface cracks in a cylindircal shell are determined by photoelastic technique. (2) By a special method photoelastic slices are handled for obtaining a clear caustic curve, and the stress intensity factors for surface cracks in a cylindrical shell are determined by the caustic method. (3) An approximate equation of curvature correction factor F_c is proposed. (4) Effects of the curvature radius R of a cylindrical shell on the stress intensity factors of surface cracks are obtained. The results of this paper are in fair agreement with already existing analytical results. The approximate equation of curvature correction factor F_c can be widely used for engineering purposes.     

拉压循环载荷下正交异性钢桥面板肋-面板焊缝裂纹扩展分析

采用四步法计算了考虑循环载荷中压应力影响的正交异性钢桥面板的肋-面板焊缝表面裂纹扩展。第一步是基于正交异性钢桥面板的疲劳分析模型,计算肋-面板焊缝处的应力,第二步是通过肋-面板焊缝的三维局部模型,用Schwartz-Neumann交替法计算焊缝表面裂纹的应力强度因子分布,第三步是用二维断裂力学模型和增量塑性损伤模型,计算循环载荷中的压应力对裂纹扩展的影响,第四步是用第二步中的三维裂纹分析结果和第三步中的二维断裂力学模型得到的裂纹扩展公式,计算钢桥面板的肋-面板焊缝表面裂纹扩展。计算结果表明,对应于正交异性钢桥面板肋-面板焊缝处的循环应力,本文所用模型的裂纹尖端反向塑性区导致裂纹扩展率增加50%以上。研究结果为正交异性钢桥面板肋-面板焊缝裂纹的疲劳寿命分析提供了研究基础。     

Fatigue crack growth behavior of Al7050-T7451 attachment lugs under flight spectrum variation

This paper discusses an analytical and experimental investigations of the fatigue crack growth behavior in attachment lugs subjected to a randomized flight-by-flight spectrum. In the analysis, the stress intensity factors for through-the-thickness cracks initiating from lug holes were compared by weight function method, boundary element method (BEM), the interpolation of Brussat’s solution. The stress intensity factors of a corner crack at a transition region were obtained using two parameter weight function method and correction factors. Fatigue life under a load spectrum was predicted using stress intensity factors and Willenborg retardation model considering the effects of a tensile overload. Experiments were performed under a load spectrum and compared with the fatigue life prediction using the stress intensity factors by different methods. Changes of fatigue life and aspect ratio according to the clipping level of the spectrum were discussed through experiment and prediction. Effect of the spectrum clipping level on the fatigue life was experimentally evaluated by using beach marks of fractured surface.     

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

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

两种各向异性材料界面共线裂纹的反平面问题

本文研究两种各向异性材料界面共线裂纹的反平面剪切问题。利用复变函数方法，提出了一般问题公式和某些实际重要问题的封闭形式解。考察了裂纹尖端附近的应力分布并给出了应力强度因子公式。从本文解签的特殊情形，可以直接导出两种各向同性材料界面裂纹，均匀各向异性材料共线裂纹以及均匀各向同性材料共线裂纹的相应问题公式，其中包括已有的经典结果。