Study on crack-inclusion interaction using digital gradient sensing method

The interaction between matrix crack and a round inclusion was studied by the method of digital gradient sensing. First, the stress fields at the matrix crack tip in the neighbor of a round inclusion were derived based on transformation toughening theory and Eshelby inclusion method, and the effect of the inclusion on the stress intensity factor of the matrix crack was analyzed. Then, the non-contact optical measurement system of digital gradient sensing was built up, and a three-point bending test was carried out using a single-edge cracked specimen. The mode I stress intensity factor was extracted from the angular deflection of the light rays. Finally, the effect of the inclusion on the angular deflection fields and the stress intensity factor at the crack tip was analyzed experimentally. These results will play an important role for evaluating the fracture mechanism of crack-inclusion interaction in composites.     

Failure initiation at a blunt V-notch tip under mixed mode loading

A criterion to predict crack onset at a sharp V-notch tip in homogeneous brittle materials under a mixed-mode loading was presented and validated by experimental observations in a previous paper by the authors. This criterion slightly underestimates the experimental loads causing failure which is attributed to a small notch tip radius that blunts the sharp corner. This discrepancy is rigorously analyzed mathematically in this paper by means of matched asymptotics involving 2 small parameters: a micro-crack increment length and the notch tip radius. A correction is brought to the initial prediction and a better agreement is obtained with experiments on PMMA notched specimens.     

K‐dominance of static crack tip in functionally gradient materials

K‐dominance of static crack tip in functionally gradient materials (FGMs) with a crack oriented along the direction of the elastic gradient is studied through coherent gradient sensing (CGS), digital speckle correlation method (DSCM) and finite element method (FEM). In the direction of crack propagation, the shear modulus has a linear variation with constant mass density and Poisson's ratio. First, the CGS and DSCM governing equations related to the measurements and the elastic solutions at mode I crack in FGMs are obtained in terms of the stress intensity factor, material constants and graded index. Secondly, two kinds of FGMs specimens and one homogenous specimen are prepared to observe the influences of the property variation on the K‐dominance. Then, CGS and DSCM experiments using three‐point‐bending of FGMs and homogenous beams are performed. Thirdly, based on the results of the experiments, the stress intensity factors of three kinds of specimens are calculated by CGS and DSCM. Meanwhile, the stress intensity factors are obtained by FEM. Finally, comparing the results from CGS, DSCM and FEM, the K‐dominance of mode‐I static crack tip in FGMs is discussed in detail.     

Determination of Stress Intensity Factor for Cracks in Orthotropic Composite Materials using Digital Image Correlation

Abstract: This paper focuses on the application of the digital image correlation (DIC) technique to determine the stress intensity factor (SIF) for cracks in orthotropic composites. DIC is a full‐field technique for measuring the surface displacements of a deforming object and can be applied to any type of material. To determine the SIF from full‐field displacement data, the asymptotic expansion of the crack‐tip displacement field is required. In this paper the expansion of the crack tip displacement field is derived from an existing solution for strain fields. Unidirectional fibre composite panels with an edge crack aligned along the fibre were tested under remote tensile loading and the displacements were recorded using DIC. The SIF was calculated from the experimental data by fitting the theoretical displacement field using the least squares method. The SIF thus determined was in good agreement with theoretical results and therefore demonstrates the applicability of the derived displacement field and DIC technique for studying fracture in composites.     

Evaluating crack tip stress field in a thin glass plate under thermal load

The stress field around a propagating crack tip in a quenched thin glass plate is discussed through experimental and theoretical analyses. Instantaneous phase-stepping photoelasticity using a CCD camera equipped with a pixelated micro-retarder array is used for measuring the crack tip stress field. From the successive phase maps of principal direction, the position and the velocity of the crack tip are evaluated. On the other hand, the fracture parameters, that is, the stress intensity factors and the T-stress are determined from the phase maps of the retardation. Experimental results obtained for a straight crack show good agreement with those obtained by theory of elasticity. The results also indicate that the direction of the crack propagation arising in the quenching process is not determined by the direction of the maximum principal stress. Furthermore, the results show that the T-stress criterion is inappropriate to evaluate the crack path instability in a quenched thin glass plate.     

Factors affecting the intensity distribution in high stress gradient zones in digital photoelasticity

In digital photoelastic analysis it is observed that in high stress gradient zones the fringes become brighter. It leads to erroneous results, if not accounted for, in those techniques which use quantitative values of intensity to evaluate the experimental parameters. In this note, the reasons for the above phenomenon is investigated and suitable recommendations are given.     

Stress intensity factor analysis for an interface crack between dissimilar isotropic materials under thermal stress

Thermal stresses, one of the main causes of interfacial failure between dissimilar materials, arise from different coefficients of linear thermal expansion. Two efficient numerical procedures in conjunction with the finite element method (FEM) for the stress intensity factor (SIF) analysis of interface cracks under thermal stresses are presented. The virtual crack extension method and the crack closure integral method are modified using the superposition method. The SIF analyses of some interface crack problems under mechanical and thermal loads are demonstrated. Very accurate mode separated SIFs are obtained using these methods.     

Stress Intensities at Spot Welds

The stress intensities (notch stress, stress intensity factors and J-integral) at spot welds under typical loads of tensile-shear, cross-tension and coach-peel are derived as a number of simple formulas on the basis of an analytic solution where the stress intensities at spot welds are generally determined by the stresses around the spot welds and of some analytic solutions to circular rigid inclusions in plates with the inclusions simulating the weld nuggets. The derived formulas show consistently the trends in the stress intensities with the design parameters for spot welds such as nugget diameter and sheet thickness and additionally with spacing of force for cross-tension spot welds and load eccentricity for coach-peel spot welds. The stress intensities at spot welds under general loading conditions are estimated in terms of the forces and moments transferred by the spot welds based on the derivations. The theoretical predictions from the formulas are compared favorably with the finite element results. As an application example, some fatigue test data for spot welds in the form of load range versus life to failure are transferred into the form of stress intensities range versus life to failure with the scatterband of the fatigue test data being substantially reduced. This revised version was published online in August 2006 with corrections to the Cover Date.     

聚合物梯度材料黏弹性断裂的双控制参数

针对组分材料体积分数任意分布的聚合物功能梯度材料,研究其在蠕变加载条件下Ⅰ型裂纹应力强度因子（SIFs）和应变能释放率的时间相依特征。由Mori-Tanaka方法预测等效松弛模量,在Laplace变换域中采用梯度有限元法和虚拟裂纹闭合方法计算断裂参数,由数值逆变换得到物理空间的对应量。分析边裂纹平行于梯度方向的聚合物功能梯度板条,分别考虑均匀拉伸和三点弯曲蠕变加载。结果表明,聚合物梯度材料应变能释放率随时间增加,其增大的程度与黏弹性组分材料体积分数正相关;材料的非均匀黏弹性性质产生应力重新分布,导致裂纹尖端应力场强度随时间变化,当裂纹位于黏弹性材料含量较低的一边时,应力强度因子随时间增加,反之,随时间减小。而且,材料的应力强度因子与时间相依的变化范围和体积分数分布以及加载方式有关,当体积分数接近线性分布时,变化最明显,三点弯曲比均匀拉伸的变化大。SIFs随时间的延长增加或减小、加剧或减轻裂纹尖端部位的“衰坏”,表明黏弹性功能梯度裂纹体的延迟失稳需要联合采用应力强度因子与应变能释放率作为双控制参数。     

Stress intensity factor computation using the method of fundamental solutions: mixed‐mode problems

The method of fundamental solutions is applied to the computation of stress intensity factors in linear elastic fracture mechanics. The displacements are approximated by linear combinations of the fundamental solutions of the Cauchy–Navier equations of elasticity and the leading terms for the displacement near the crack tip. Two algorithms are developed, one using a single domain and one using domain decomposition. Numerical results are given. Copyright © 2006 John Wiley & Sons, Ltd.     

NiCr/ZrO2功能梯度复合材料中混合型准静态裂纹启裂的数字散斑相关方法实验研究

采用粉末冶金方法制备出NiCr/ZrO₂功能梯度材料FGMs。通过2种断裂试件研究了材料梯度对混合型断裂行为的影响(FGM-A试件,裂纹位于试件的弹性模量较大一侧;FGM-B试件,裂纹位于试件的弹性模量较小一侧)。对2种断裂试件在非对称载荷下进行准静态断裂实验,并利用数字散斑相关方法测得Ⅰ、Ⅱ型应力强度因子。结果表明：FGM-A的裂纹的开裂角小于FGM-B的开裂角;FGM-A的弹性梯度对静态裂纹有保护作用;弹性模量的梯度变化和裂尖局部材料的断裂韧性会影响混合型裂纹的启裂。     

Automatic Evaluation of Mixed-mode Stress Intensity Factors Utilizing Digital Image Correlation

Abstract:  A method for determining mixed-mode stress intensity factors from displacement fields obtained by digital image correlation is studied. To perform automatic evaluation, not only stress intensity factors but crack-tip location, higher-order terms in the series expansion of displacement fields and rigid-body displacement components are also determined simultaneously using nonlinear least squares based on the Newton–Raphson method. Experimental results show that the mixed-mode stress intensity factors are evaluated accurately from either radial or circumferential displacement components. As mixed-mode stress intensity factors can be evaluated easily, simply and automatically by the technique with digital image correlation and nonlinear least squares, it is expected that the proposed method can be applied to solve various fracture problems.     

Stress intensity factor solution for a semi-elliptical crack in an arbitrarily distributed stress field

An equation for the stress intensity factor (SIF) for semi-elliptical crack has been developed. It is based on the Newman-Raju's solution for the crack in a plate under bending or tension. The equation can be applied when a stress distribution is described by a power function. Using the approach outlined, the SIF for a surface crack in a T-butt welded connection has been estimated. The results obtained can be used in a fracture-mechanics-based fatigue analysis.     

Arbitrarily oriented microcracks near the tip of an interface macrocrack in bonded dissimilar anisotropic materials

It is well known that microcracking in brittle materials results in a reduction of the stress intensity factor (SIF) and energy release rate (ERR). The reduced SIF or ERR represents crack tip shielding which is of significant interest to micromechanics and material science researchers. However, the effect of microcracking on the SIF and ERR is a complicated subject even for isotropic homogeneous materials, and becomes much more formidable in case of interface cracks in bonded dissimilar solids. To unravel the micromechanics of interface crack tip shielding in bonded dissimilar anisotropic solids, an interface crack interacting with arbitrarily oriented subinterface microcracks in bonded dissimilar anisotropic materials is studied. After deducing the fundamental solutions for a subinterface crack under concentrated normal and tangential tractions, the present interaction problem is reduced to a system of integral equations which is then solved numerically. A J‐integral analysis is then performed with special attention focused on the J₂‐integral in a local coordinate system attached to the microcracks. Theoretical and numerical results reassert the conservation law of the J‐integral derived for isotropic materials 1 , 2 also to be valid for bonded dissimilar anisotropic materials. It is further concluded that there is a wastage when the remote J‐integral transmits across the microcracking zone from infinity to the interface macrocrack tip. In order to highlight the influence of microstructure on the interfacial crack tip stress field, the crack tip SIF and ERR in several typical cases are presented. It is interesting to note that the Mode I SIF at the interface crack tip is quite different from the ERR in bonded dissimilar anisotropic materials.     

高温合金CT试样裂纹前沿应力分布及第二相影响的模拟计算

利用有限元方法分析了IN718合金CT试样裂纹区域的应力应变分布、裂纹尖端应力强度因子及第二相状态对它的影响，计算结果与实测数据基本一致，表明采用有限元方法研究含裂纹体构件的蠕变行为是一种较好的研究手段。通过计算认为，IN718合金在时效初期虽然有α-Cr相的析出而不至于促进蠕变裂纹扩展，是由于α—Cr于δ相周围γ″贫化区(微塑性区)内析出而处于微塑性区内使得裂纹尖端钝化，不利于裂纹扩展。     

Stress intensity factors for circumferential through‐wall cracks in thin‐walled cylindrical shells subjected to tension and torsion

In order to assess the structural integrity of tubular members or pipes containing circumferential through‐wall cracks, their stress intensity factor solutions are required. While stress intensity factors for tension and bending are available, few solutions exist for the case of torsion, even though these components may also be subjected to torque. In this paper, the finite element method is used to compute the stress intensity factors for this geometry under tension and torsion. Shell elements are employed to compute the results for thin shells by the means of the displacement extrapolation technique. The computed results indicate that the available analytical solution for torsional loading, which is based on shallow shell theory, is nonconservative for long cracks in thin shells. Shallow shell theory is in general not applicable to long cracks, and the present work is therefore able to provide solutions for a wider range of crack lengths than what is currently available.     

An accurate scheme for mixed‐mode fracture analysis of functionally graded materials using the interaction integral and micromechanics models

The interaction integral is a conservation integral that relies on two admissible mechanical states for evaluating mixed‐mode stress intensity factors (SIFs). The present paper extends this integral to functionally graded materials in which the material properties are determined by means of either continuum functions (e.g. exponentially graded materials) or micromechanics models (e.g. self‐consistent, Mori–Tanaka, or three‐phase model). In the latter case, there is no closed‐form expression for the material‐property variation, and thus several quantities, such as the explicit derivative of the strain energy density, need to be evaluated numerically (this leads to several implications in the numerical implementation). The SIFs are determined using conservation integrals involving known auxiliary solutions. The choice of such auxiliary fields and their implications on the solution procedure are discussed in detail. The computational implementation is done using the finite element method and thus the interaction energy contour integral is converted to an equivalent domain integral over a finite region surrounding the crack tip. Several examples are given which show that the proposed method is convenient, accurate, and computationally efficient. Copyright © 2003 John Wiley & Sons, Ltd.