Strain energy density prediction of crack initiation and direction in cracked T-beams and pipes

In this paper, the S-theory is applied to determine crack initiation and direction for cracked T-beams and circumferentially cracked pipes. It makes use of a parameter called strain energy density factor, S, which is a function of the stress intensity factors. The strain energy density theory provides a more general treatment of fracture mechanics problems by virtue of its ability in describing the multiscale feature of material damage and in dealing with mixed mode crack propagation problem. A simple method for obtaining approximate stress intensity factors is also applied. It takes into account the elastic crack tip stress singularity while using the elementary beam theory. Some basic loading conditions in beams and pipes are studied.     

Modelling of crack growth by fracture damage zone

The strain energy density (SED) criterion is applied for analyzing the full range of mixed mode fracture from tensile to shear loading. A fracture damage zone (FDZ) local to the crack tip is defined and discussed in connection with the influence of crack geometry, loading and local material property. The size of FDZ tends to change continuously from statically to cyclically applied load conditions. It can be estimated from the uniaxial mechanical properties of the material. Both experimental and analytical results are examined for subcritical crack growth under static loading that depends on the type steel structures the fracture behavior of which could be represented by a single curve for the given specimen geometry.     

Cracking characteristics of mixed mode dislocations near a lip-like mode crack

This work is concerned with the cracking characteristics of mixed mode dislocations near a lip-like mode crack, stress intensity and strain energy density factor are obtained by using conformal mapping, singularity analysis and Cauchy integrals. Shielding effect generated by screw dislocation near a lip-like mode crack decreases with the increment of the distance between screw dislocation and crack tip. Larger distance between two faces of the crack leads to the shielding effect waning. The strain energy density factor of mode III decreases with the increment of the distance between dislocation and crack tip. Larger distance between two faces of lip-like mode crack also leads to the strain energy density factor waning and encourages crack initiation; the shielding effects generated by edge dislocation near the crack decrease with the increment of the distance between edge dislocation and crack tip.     

Evaluation of crack initiation angle under mixed mode loading at diverse strain rates

Crack initiation angle, under mixed mode loading at several strain rates, is analysed using an experimental–numerical approach. The physical phenomenon for the problem at hand is influenced by the local and global conditions. One of such factors is the strain rate at the crack tip. For this purpose, PMMA plates with centred angled cracks under mixed mode loading were tested. The strain rate at the neighbourhood of the crack tip before crack propagation was evaluated. Considering that this material is strain rate sensitive, the numerical models were calibrated with the modulus of elasticity measured in tension tests at the observed strain rates. Numerical evaluations were performed with the finite element method in conjunction with the volume energy density criterion. An improvement in the evaluation of the crack propagation angle was observed. In order to complete the analysis, the crack initiation angle was also evaluated with the strain energy density factor S, considering the mechanical properties of PMMA, as evaluated at the observed strain rates, and the stress intensity factors k₁ and k₂. Results are in agreement with those observed experimentally.     

Fracture behavior of a nickel-base single crystal superalloy as predicted by the strain energy density criterion

Results from compact tension specimens and scanning electronic microscopy are used to study the fracture characteristics of nickel-base single crystal superalloy DD3. Finite element analyses are also made to better understand the outcome. The experiments showed that the initiation and growth of cracks were dependent on the crystallographic orientation of specimens, and the fracture surfaces coincided with the slip planes. Three-dimensional stresses near the crack tip are obtained. The strain energy density criterion is applied to predict the fracture of the superalloy. Agreement with experiments is good.     

Asymptotic character of mixed mode in plane deformation of crack in rubber-like material

The asymptotic stress and strain distribution near a crack tip in rubber-like materials is determined by finite element for in-plane mixed mode loading. For large strain, the crack tip field is always in a state of uniaxial tension. The shear load affects only the orientation of the deformed near tip field in the space. A good agreement is obtained between the theoretical and numerical results.     

Mixed mode fatigue crack growth with a sudden change in loading direction

With a sudden change in the maximum load level, there will be a corresponding change in the crack driving force regardless of whether the load is applied monotonically or cyclically. The effective strain energy density factor range ΔS_p,eff has been used to correlate mixed mode fatigue crack propagation where the crack growth direction is not known as an a priori. Examined in this work is a sudden change of load direction on fatigue crack growth while the load level remains unchanged. Yielding is assumed to be localized near the crack tip such that the crack growth behavior can be described adequately by the elastic stress field. Under the conditions investigated, minimal change on crack growth rates is observed. No firm conclusion could be drawn on deviation of crack path for the case considered.     

Stress multiaxiality factor for crack growth prediction using the strain energy density theory

The multiaxiality factor defined as the ratio of the von-Misses equivalent stress to the volumetric stress has been reported to be related to the initiation and progression of failure in structures. It is demonstrated in the present paper that the location around the crack tip where the multiaxiality factor obtains minimum value is an indicator of the direction of minimum material fracture resistance for crack propagation. It is also proposed that the location along the direction of crack propagation path where multiaxiality factor obtains minimum value is considered as the critical distance away from the crack tip, where the strain energy density should be evaluated and compared to its critical value. Theoretical predictions correlate well with the test results for the investigated cases.     

Review of brittle fracture criteria in case of static and cyclic mixed mode loading

The paper gives in the first part in pressed form a survey of brittle fracture criteria using a reference intensity factor in case of static mixed mode loading. Criteria (expressed in terms of different quantities such as stress, deformation and strain energy) usually refer to a parameter that is characteristic of the material response at fracture. Criteria include information on two basic hypotheses (crack propagation direction and unstable crack growth). In the second part a generalized method is suggested for application of cyclic reference intensity factor in case of cyclic mixed mode loading. Three basic hypotheses describe crack growth direction, stable crack growth steps and unstable crack growth.     

压缩载荷作用下岩石的细观损伤和塑性研究

建立岩石微裂纹扩展的细观力学模型,研究了岩石的细观损伤和塑性性质.压缩载荷下微裂纹尖端翼裂纹稳定扩展表征岩石的细观损伤,采用应变能密度准则求解复合型断裂的翼裂纹扩展长度,微裂隙统计的二参数Weibull函数模型反映绝对体积应变对微裂纹分布数目影响,进而用翼裂纹扩展所表征的应力释放体积和微裂纹数目来表示含有微裂隙的岩石损伤演化变量;宏观塑性屈服函数采用Voyiadjis等的等效塑性应变的硬化函数,反映了塑性内变量对硬化函数的影响;建立岩石模型的本构关系及其数值算法,并用回映隐式积分算法编制了模型的本构程序.分析弹塑性损伤模型的围压对岩石损伤的影响,并从围压和短微裂隙长度等因素分析模型的岩石的损伤和宏观塑性特性.     

Inclined semi-elliptical crack for predicting crack growth direction based on apparent stress intensity factors

Linear elastic criterion of the inclined semi-elliptical crack growth direction is elaborated on the basis of the strain energy density theory. Stress and displacement fields are presented for higher order terms asymptotic expansion. Solutions for elastic stress intensity factors are accounting for the function describing of the crack tip fields near the free surface of plate. The mixed mode behavior of crack growth direction angle along the semi-elliptical crack front for different combination of biaxial loading, inclination crack angle and surface flaw geometry is determined.     

复合型韧性断裂实验和控制参数

通过对不同韧性材料在各种平面复合载荷形式下裂纹启裂阶段裂端变形、启裂位置和扩展方向的系统的宏微观实验验证及计算分析,考察了韧性断裂参数空穴扩张比的分布特征和裂纹启裂及扩展方向的关系．得到：对于不同韧性的材料,在裂端的钝化变形区域,空穴扩张比的极大值区对应干裂纹的启裂位置,裂纹启裂时钝化裂端前缘空穴扩张比的临界值不敏感于复合比的变化．而对于裂纹启裂后的扩展方向,则需根据具体材料在相应的特定区域中比较空穴扩张比参数极大值的分布特征,需经进一步的分析,从而确定裂纹的扩展方向．实验及计算结果表明,尽管复合型断裂时裂纹启裂及扩展的机理极其复杂,用于韧性材料复合型断裂的空穴扩张比参数仍能很好地预测裂端的启裂及扩展方向,可作为复合型韧性断裂过程的控制参数．     

Mixed mode fatigue crack growth and the strain energy density factor

The strain energy density factor S was first proposed by Sih for the prediction of the critical of the load and failure direction under monotonic, mixed mode loading condition. It seems a natural extension to apply the same concept to fatigue crack propagation. However, a close examination of the existing theory indicates that the Strain Energy Density Factor cannot logically account for the phenomena of the R-ratio effect and crack arrest. Thus, modification is necessary before the concept can be applied successfully for the prediction of mixed mode fatigue crack propagation.Based on the concept of hysteresis energy dissipation, an effective strain energy density factor range, ΔS_p,eff, is proposed for the correlation of fatigue crack growth data. ΔS_p,eff is consistent with the concept of crack closure. Experimental investigation indicates that it could predict the crack growth rates and trajectories.     

Stress intensity factor range and propagation mode of surface cracks under rolling–sliding contact

Finite element analyses were conducted in order to evaluate the mode I and mode II stress intensity factors for inclined edge cracks under cyclic contact load under rolling and rolling–sliding condition. The SIF range depends on crack orientation, crack length to Hertzian contact zone half-width ratio, friction between the crack faces and friction on the contact surface. The results were combined in two compact functions that determine the ΔK_I and ΔK_II values. The crack propagation mode and direction were investigated using both the maximum stress criterion and the minimum strain energy density criterion. The results are displayed in graph form, which allows a fast evaluation of the crack growth condition.     

Crack initiation behavior of functionally graded piezoelectric material: prediction by the strain energy density criterion

Transient response of a functionally graded piezoelectric medium is considered for a through crack under the mixed-mode in-plane mechanical and electric load. Integral transforms and dislocation density functions are employed to reduce the problem to singular integral equations. The energy density factor criterion is applied to obtain the maximum of the minimum energy density factor. This determines the direction of crack initiation. Numerical results display the effects of material constants, loading combination parameter, mechanical loading angle and material gradient parameter on the possible fracture behavior.     

Dynamic crack tip fields and dynamic crack propagation characteristics of anisotropic material

Based on mechanics of anisotropic material, the dynamic crack propagation problem of I/II mixed mode crack in an infinite anisotropic body is investigated. Expressions of dynamic stress intensity factors for modes I and II crack are obtained. Components of dynamic stress and dynamic displacements around the crack tip are derived. The strain energy density theory is used to predict the dynamic crack extension angle. The critical strain energy density is determined by the strength parameters of anisotropic materials. The obtained dynamic crack tip fields are unified and applicable to the analysis of the crack tip fields of anisotropic material, orthotropic material and isotropic material under dynamic or static load. The obtained results show Crack propagation characteristics are represented by the mechanical properties of anisotropic material, i.e., crack propagation velocity M and fiber direction α. In particular, the fiber direction α and the crack propagation velocity M give greater influence on the variations of the stress fields and displacement fields. Fracture angle is found to depend not only on the crack propagation but also on the anisotropic character of the material.     

铁中Ⅰ型裂纹裂尖场与裂尖形变机制分析

应用有限元及连续介质力学计算了平面应交条件下体心立方铁中不同取向裂纹的裂尖应力场和弹性应变能密度分布,从宏微观相结合角度分析了裂尖场分布与裂尖微结构演化的相互关联.指出裂尖塑性变形的具体形式与裂尖滑移面上分切应力的大小密切相关,并从能量角度解释了裂尖相变产生的原因,最后结合相关分子动力学研究成果探讨了裂尖奇异性区域.     

Moving crack for functionally grated material in an infinite length strip under antiplane shear

Considered is a Yoffe crack in an infinite strip of functionally grated material (FGM) subjected to antiplane shear. The shear moduli in two directions of FGM are assumed to be of exponential form. The dynamic stress intensity factor and strain energy density factor at the crack tip are obtained by using integral transforms and dual-integral equations. The numerical results show that the decrease of the strain energy density factor varies with the shear moduli gradient, and the increase of the strain energy density factor varies with the increase of the moving crack speed. The ratio of shear moduli in material vertical orientation has a great influence on the strain energy density factor.     

Impact failure prediction of Mode III crack in orthotropic functionally graded strip

Mode III impact of a crack in an orthotropic functionally graded strip is investigated. The shear moduli in two directions of the material are assumed to vary proportionately with gradient. Laplace transform and Fourier cosine transform are used to reduce the problem to solving a Fredholm integral equation. The crack tip stress field is obtained by considering the asymptotic behavior of Bessel function. Energy density factor criterion is applied to obtain the maximum of minimum energy density and direction of crack initiation. Numerical results are given graphically. The effects of orthotropy, nonhomogeneity and height of the strip on the energy density factor are discussed.