Determination of stress intensity factor with direct stress approach using finite element analysis

In this article,a direct stress approach based on finite element analysis to determine the stress intensity fac-tor is improved.Firstly,by comparing the rigorous solution against the asymptotic solution for a problem of an infinite plate embedded a central crack,we found that the stresses in a restrictive interval near the crack tip given by the rigorous solution can be used to determine the stress intensity fac-tor,which is nearly equal to the stress intensity factor given by the asymptotic solution.Secondly,the crack problem is solved numerically by the finite element method.Depending on the modeling capability of the software,we designed an adaptive mesh model to simulate the stress singularity.Thus, the stress result in an appropriate interval near the crack tip is fairly approximated to the rigorous solution of the corre-sponding crack problem.Therefore,the stress intensity factor may be calculated from the stress distribution in the appro-priate interval,with a high accuracy.     

Equivalence of the notch stress intensity factor,tip opening displacement and energy release rate for a sharp V-notch

For an infinite elastic plane with a sharp V-notch under the action of symmetrically loading at infinity, the length of crack initiation ahead of the V-notch’s tip is estimated according to a modified Griffith approach. Applying a new conservation integral to the perfectly plastic strip (Dugdale model) ahead of the V-notch’s tip, the relationship between notch stress intensity factor (NSIF) and notch tip opening displacement (NTOD) is presented. Also, the relationship between NSIF and perfectly plastic strip size (PPSS) is found. Since there are three fracture parameters (NSIF, NTOD, and PPSS) with changeable notch opening angle in two basic relationships, it is necessary to select one critical parameter with changeable notch opening angle or two independent critical parameters, respectively. With the help of a characteristic length, it is found by this new conservation integral that the NSIF, NTOD and energy release rate are equivalent in the case of small-scale yielding. Especially, the characteristic length possesses clear physical meaning and, for example, depends on both the critical NSIF and SIF or both the NTOD and CTOD, respectively, in which SIF and CTOD are from the tip of a crack degenerated from the sharp V-notch. The dependence of NSIF on NTOD and PPSS is presented according to the equivalence, and the critical NSIF depending on the critical NTOD with a notch opening angle is also predicted.     

A singular element for calculating the stress intensity factor

This paper proposes a new type of special element (sectorial singular element) for calculating the linear elastic stress intensity factor. The shape of element not only accords with the demands of finite element analysis, but also coincides with the theory of linear elastic fracture mechanics. The accuracy and economy of the result in this paper are satisfactory.     

Stress field at a tip of a prefabricated spiral V-notch

Based on the tranditional V-notched blasting, a technique of spirally V-notched blasting to loosen earth and rock was presented. Fracture mechanics and Westergaard stress function were adopted to build a complex stress function to derive the plane stress and strain fields at one tip of the crack under a quasi-static pressure. An expression was formulated to define the stress intensity factor of spiral V-notch loosen blasting. Factors that have effects on the stress intensity factor were studied. It is demonstrated that spiral V-notch loosen blasting is an extension of vertical V-notch blasting, straight cracking, and alike theories.     

Plastic-strain distribution at the root of a sharp notch

A new experimental technique employing a pattern of microhardness indentations has been developed to permit an accurate measurement of the plastic-strain distribution at the root of a mechanically introduced notch. The method allows measurements to be made to within 0.001 in. from the crack tip. The longitudinal strain distribution was determined for notched 2024-T3 aluminum-alloy specimens, and both the longitudinal and the transverse distributions were obtained for notched specimens of high-purity aluminum. The results from the experiments indicated that a relationship exists between the grain size and both the strain magnitude and distribution. The specimens having the largest grain size developed the highest strain gradients, and exhibited the greatest degree of symmetry about the center line. The symmetry of the strain distribution around the notch improved as the applied load approached that required to produce crack propagation whereas, at low loads, there was little evidence of any symmetry.     

Photoelastic analysis of the singular stress field in a bimaterial wedge

This paper presents an experimental investigation of the singular stress field near the vertex of a bimaterial wedge using a digital photoelastic technique. Special attention is given to the casting of bimaterial wedge specimens and analysis technique for extracting stress intensity factors from photoelastic samples. Different bimaterial wedge specimens are made of two different photoelastic materials bonded through a special casting procedure and loaded in simple tension. A new multiple-parameter method is developed to obtain the stress intensity factor reliably from the isochromatic fringe patterns and the series representation of the stress field at the vertex of the wedge. Experimental results are compared with finite element predictions, and good agreement is observed.     

A study on stress intensity factors and singular stress fields of three-dimensional interface crack

By using the finite-part integral concepts and limit technique, the hypersingular integrodifferential equations of three-dimensional (3D) planar interface crack were obtained; then the dominant-part analysis of 2D hypersingular integral was further used to investigate the stress fields near the crack front theoretically, and the accurate formulae were obtained for the singular stress fields and the complex stress intensity factors. After that, a numerical method is proposed to solve the hypersingular integrodifferential equations of 3D planar interface crack, and the problem of elliptical planar crack is then considered to show the application of the method. The numerical results obtained are satisfactory. Project supported by the Foundation of Solid Mechanics Open Research Laboratory of State Education Commission at Tongji University and the National Natural Science Foundation.     

Crack initiation at the rectangular step notch surface under combined loading

Summary  Under external forces acting on the face of a notch, cracks originate at corners, and the system is liable to fail. An analysis is presented of the stress field in the neighborhood of the notch tips, based on the integral representation of the biharmonic solution and on numerical methods. Computations were performed for constant loading or constant displacement distributed along one face of the notch. The coefficients in the principal terms of the asymptotic formulae for the circumferential and shear stresses depend on the angle and height of the notch face and on the boundary conditions. The maximal values of these coefficients determine the stress intensity factors for the opening and shear modes. The angles corresponding to the maximal values of the intensity factors indicate the directions of initiation of opening and sliding cracks. Received 30 May 2000; accepted for publication 3 April 2001     

复合材料中矩形夹杂角端部力学行为分析

提出了一种分析矩形夹杂角端部奇异应力场的新型杂交有限元方法,该方法在分析矩形夹杂角端部奇异应力场时,需要在夹杂端部构造一个超级单元。超级单元的刚度矩阵可以通过夹杂端部特征问题数值解建立。我们用这种方法计算了单向载荷作用下无限大均质板中单个矩形夹杂角端部奇异应力场,并与现有的数值解进行了比较。比较结果表明:本文提出的方法是可行的、有效的,而且数值结果精度高。为说明本文方法适用范围更广,文章最后讨论了各向异性弹性材料和横观各向同性压电材料中矩形夹杂角端部电弹性场行为。     

Strain energy approach to compute stress intensity factors for isotropic homogeneous and bi-material V-notches

A strain energy approach (SEA) is developed to compute the general stress intensity factors (SIFs) for isotropic homogeneous and bi-material plates containing cracks and notches subject to mode I, II and III loading conditions. The approach is based on the strain energy of a control volume around the notch tip, which may be computed by using commercial finite element packages. The formulae are simple and easy to implement. Various numerical examples are presented and compared to corresponding published results or results that are computed using different numerical methods to demonstrate the accuracy of the SEA. Many of those results are new, especially for the cases of bi-material notches where the problem is quite complicated.     

A unified approach to problems of stress concentration near V-shaped notches with sharp and rounded tip

The paper proposes a unified approach to problems of stress concentration near notches with sharp and rounded tip based on the method of singular integral equations. A solution for an elastic region having a V-shaped notch with rounded tip of large curvature is first found. Then, the stress intensity factor at the tip of a sharp-tipped notch is calculated by passing to the limit. Numerical results are obtained for a slit and a square hole in an elastic plane and an edge notch in a half-plane __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 2, pp. 70–87, February 2007. For the centenary of the birth of G. N. Savin.     

On the asymptotic crack-tip stress fields in nonlocal orthotropic elasticity

The crack-tip stress fields in orthotropic bodies are derived within the framework of Eringen’s nonlocal elasticity via the Green’s function method. The modified Bessel function of second kind and order zero is considered as the nonlocal kernel. We demonstrate that if the localisation residuals are neglected, as originally proposed by Eringen, the asymptotic stress tensor and its normal derivative are continuous across the crack. We prove that the stresses attained at the crack tip are finite in nonlocal orthotropic continua for all the three fracture modes (I, II and III). The relative magnitudes of the stress components depend on the material orthotropy. Moreover, non-zero self-balanced tractions exist on the crack edges for both isotropic and orthotropic continua. The special case of a mode I Griffith crack in a nonlocal and orthotropic material is studied, with the inclusion of the T-stress term.     

On the stress state of a piezoceramic body with a flat crack under symmetric loads

A static-equilibrium problem is solved for an electroelastic transversely isotropic medium with a flat crack of arbitrary shape located in the plane of isotropy. The medium is subjected to symmetric mechanical and electric loads. A relationship is established between the stress intensity factor (SIF) and electric-displacement intensity factor (EDIF) for an infinite piezoceramic body and the SIF for a purely elastic material with a crack of the same shape. This allows us to find the SIF and EDIF for an electroelastic material directly from the corresponding elastic problem, not solving electroelastic problems. As an example, the SIF and EDIF are determined for an elliptical crack in a piezoceramic body assuming linear behavior of the stresses and the normal electric displacement on the crack surface __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 11, pp. 67–77, November 2005.     

Elastic stress distributions for hyperbolic and parabolic notches in round shafts under torsion and uniform antiplane shear loadings

Closed-form solutions are developed for the stress fields induced by circumferential hyperbolic and parabolic notches in axisymmetric shafts under torsion and uniform antiplane shear loading. The boundary value problem is formulated by using complex potential functions and two different coordinate systems, providing two classes of solutions. It is also demonstrated that some solutions of linear elastic fracture and notch mechanics reported in the literature can be derived as special cases of the general solutions proposed herein.Finally the analytical frame is used to link the Mode III notch stress intensity factor to the maximum shear stress at the notch tip, as well as to give closed-form expressions for the strain energy averaged over a finite size volume surrounding the notch root.     

压应力场中爆生裂纹分布与扩展特征实验分析

采用动静组合加载实验装置和数字激光焦散线实验系统,进行了0、3、6、9 MPa等4种压应力场中PMMA试件的爆破致裂实验,分析了沿静态主应力方向扩展的裂纹运动学和力学行为。实验结果表明:首先,静态竖向载荷在预制炮孔周围产生应力集中,在炮孔壁上下端部处出现最大拉应力;随后,在动态爆炸载荷的叠加作用下,裂纹优先在炮孔壁上最大拉应力位置处起裂,并沿最大主应力方向扩展;裂纹扩展过程中,静态竖向载荷越大,裂纹扩展速度越大,且裂纹尖端应力强度因子值越大。     

Mixed-mode thermal stress intensity factors from the finite element discretized symplectic method

A finite element discretized symplectic method is introduced to find the thermal stress intensity factors (TSIFs) under steady-state thermal loading by symplectic expansion. The cracked body is modeled by the conventional finite elements and divided into two regions: near and far fields. In the near field, Hamiltonian systems are established for the heat conduction and thermoelasticity problems respectively. Closed form temperature and displacement functions are expressed by symplectic eigen-solutions in polar coordinates. Combined with the analytic symplectic series and the classical finite elements for arbitrary boundary conditions, the main unknowns are no longer the nodal temperature and displacements but are the coefficients of the symplectic series after matrix transformation. The TSIFs, temperatures, displacements and stresses at the singular region are obtained simultaneously without any post-processing. A number of numerical examples as well as convergence studies are given and are found to be in good agreement with the existing solutions.     

Stress concentration around a spheroidal crack caused by a harmonic wave incident at an arbitrary angle

A method is proposed to study the stress concentration around a shallow spheroidal crack in an infinite elastic body. The stress concentration is due to the diffraction of a low-frequency plane longitudinal wave by the crack. The direction of wave propagation is established in which the combined concentration of mode I and mode II stresses is maximum __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 1, pp. 70–77, January 2006.     

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     

The expression of stress and strain at the tip of three-dimensional notch

ＴＨＥＥＸＰＲＥＳＳＩＯＮＯＦＳＴＲＥＳＳＡＮＤＳＴＲＡＩＮＡＴＴＨＥＴＩＰＯＦＴＨＲＥＥ－ＤＩＭＥＮＳＩＯＮＡＬＮＯＴＣＨＱｉａｎＪｕｎ（钱俊）（ＥａｓｔＣｈｉｎａＩｎｓｔｉｔｕｔｅｏｆＴｅｃｈｍｏｌｏｇｙ）ＮａｎｊｉｎｇＬｏｎｇＹｕ－ｑｉｕ（龙驭...