Stress intensity factor error index for finite element analysis with singular elements

An error index for the stress intensity factor (SIF) obtained from the finite element analysis results using singular elements is proposed. The index was developed by considering the facts that the analytical function shape of the crack tip displacement is known and that the SIF can be evaluated from the displacements only. The advantage of the error index is that it has the dimension of the SIF and converges to zero when the actual error of the SIF by displacement correlation technique converges to zero. Numerical examples for some typical crack problems, including a mixed mode crack, whose analytical solutions are known, indicated the validity of the index. The degree of actual SIF error seems to be approximated by the value of the proposed index.     

Determination of SIFs, crack-tip coordinates and crack angle of anisotropic materials

This paper investigates the accuracy of a numerical method incorporating an image-correlation experiment for finding the crack-tip coordinates, the axes of the measurement and the mixed-mode stress intensity factors (SIFs) of anisotropic materials. Finite element and experimental simulations in this paper indicate that the proposed linear search method incorporating Powell's method is considerably stable for finding the unknown crack-tip coordinates and measurement axes, even though the initial input is far away from the actual data. Moreover, the displacement field very near the crack tip is not necessary. In the experimental validation, the maximum SIF error of the image-correlation experiment is about 9%, which should be acceptable for the mixed-mode fracture problems.     

Application of the weight function method to solving multiparametric three-dimensional fracture mechanics problems

A system of equations of the weight function method is obtained for the calculation of the stress intensity factors (SIF) along the crack front under arbitrary normal loading, as well as the crack opening displacement field (CODF) for cracks of arbitrary shape. The equations, in particular, make it possible to estimate the error of the known empirical formulas for the SIF, as well as to calculate the variation of the SIF values with the crack shape. It is shown that the SIF values for the problems with the mixed boundary conditions can be calculated making use of the weight functions obtained for the problems with the boundary conditions in stresses.     

Fatigue crack growth in CFRP-strengthened steel plates

In this paper fatigue crack growth in steel plates reinforced by using carbon fiber reinforced (CFRP) strips is investigated from the experimental, numerical and analytical point of view. Single edge notched tension (SENT) specimens were strengthened with different reinforcement configurations and tested at a stress ratio R of 0.4. Different initial damage levels were considered and the experimental results showed that the reinforcement application can effectively reduce the crack growth rate and significantly extend the fatigue life. Numerical models (finite elements) were also developed to evaluate the stress intensity factor (SIF) and the crack opening displacement (COD) profile. Based on the numerical results, an analytical model was proposed to predict the fatigue crack growth rate and the fatigue crack growth curves. The analytical results are validated by comparing the fatigue crack growth curves to the experimental ones.     

Approximate weight function technique using single-layer potential for a cracked circular disc

An efficient weight function technique using the indirect boundary integral method was presented for cracked circular discs. The crack opening displacement field was presented by a single layer whose kernel was a modified form of the fundamental solution in elastostatics. The application of a single-layer potential to the weight function method leads to a unique closed-form SIF (stress intensity factor) solution. The solution can be applied to a cracked circular discs with or without an internal hole or opening. For these crack geometries over a wide range of crack ratios, the SIF solution can be applied without any modification.

The calculation procedure of SIFs for the various cracked circular discs using only one analytical solution is very simple and straightforward. The information necessary in the analysis includes only two or three reference load cases. In most cases the SIF solution using two reference SIFs gives reasonably accurate results while the SIF solution with three reference load cases may be used to improve the solution accuracy of the crack configurations, with an internal opening or hole, compared with the solutions of the available literature.     

New a posteriori error estimator and adaptive mesh refinement strategy for 2-D crack problems

A posteriori error estimation and adaptive refinement technique for fracture analysis of 2-D/3-D crack problems is the state-of-the-art. The objective of the present paper is to propose a new a posteriori error estimator based on strain energy release rate (SERR) or stress intensity factor (SIF) at the crack tip region and to use this along with the stress based error estimator available in the literature for the region away from the crack tip. The proposed a posteriori error estimator is called the K-S error estimator. Further, an adaptive mesh refinement (h-) strategy which can be used with K-S error estimator has been proposed for fracture analysis of 2-D crack problems. The performance of the proposed a posteriori error estimator and the h-adaptive refinement strategy have been demonstrated by employing the 4-noded, 8-noded and 9-noded plane stress finite elements. The proposed error estimator together with the h-adaptive refinement strategy will facilitate automation of fracture analysis process to provide reliable solutions.     

Mode II stress intensity factors for central slant cracks with frictional surfaces in uniaxially compressed plates

The effect of crack surface friction on mode II stress intensity factor (SIF) of a central slant crack in a plate uniformly loaded in uniaxial compression is quantified. A previously developed two-dimensional finite element analysis was utilised after its modification to accommodate the friction between the crack surfaces. The plane strain state was assumed. A new numerical technique was devised to avoid the iteration procedures, which had to be employed due to the existence of frictional forces.

The crack inclination angle varied between zero and 75° measured from the horizontal direction. The coefficient of friction of the crack surfaces changed from zero to 1. In case of relatively sliding crack surfaces, mode II SIF existed. As is well known, the resulting mode II SIF decreased with increasing the coefficient of friction of the crack surfaces. Further, mode II SIF increased with increasing crack line inclination angle and then decreased after reaching a maximum value. The angle corresponding to that maximum SIF increased as the coefficient of friction of the crack surfaces increased.     

SGBEM analysis of crack-particle(s) interactions due to elastic constants mismatch

The effects of elastic constants mismatch on the interaction between a propagating crack and single or multiple inclusions in brittle matrix materials are investigated using numerical simulations. The simulations employ a quasi-static crack-growth prediction tool based upon the symmetric-Galerkin boundary element method (SGBEM) for multiregions, a modified quarter-point crack-tip element, the displacement correlation technique for evaluating stress intensity factors (SIFs), and the maximum principal stress criterion for crack-growth direction. It is shown that, even with this simple method for calculating SIF, the crack-growth prediction tool is both highly accurate and computationally effective. This is evidenced by results for the case of a single inclusion in an infinite plate, where the SGBEM results for the SIFs show excellent agreement with known analytical solutions. The simulation results for crack growth and stress intensity behaviors in particulate media are very stable. The crack-tip shielding and amplification behaviors, as seen in similar studies using other numerical approaches, can be clearly observed.     

A new analytical method for stress intensity factors based on insitu measurement of crack deformation under biaxial tension

A new approach for the calculation of stress intensity factors (SIF) for isotropic and orthotropic materials under biaxial tension loading was proposed in this paper. In order to determine SIF from the full-field displacement data, an asymptotic expansion of the crack tip displacement field was performed. The deforming shape and surface residual stress of the crack tip was obtained at the early extended stage of the loading process by using optical microscope and X-ray diffraction measurement. During this stage, a modified Dugdale Model, which takes into account the coupled effect at the crack tip, was proposed for the open displacement of the crack tip. In this paper, the SIFs of two types of silicon steel sheet with isotropic and orthotropic properties were calculated using the modified Dugdale Model based on the biaxial tension experimental data. From the results, it was found that analysis using the modified Dugdale Model is an effective way to evaluate SIF under biaxial stress.     

Near‐crack contour behaviour and extraction of log‐singular stress terms of the self‐regular traction boundary integral equation

This work contains an analytical study of the asymptotic near‐crack contour behaviour of stresses obtained from the self‐regular traction‐boundary integral equation (BIE), both in two and in three dimensions, and for various crack displacement modes. The flat crack case is chosen for detailed analysis of the singular stress for points approaching the crack contour. By imposing a condition of bounded stresses on the crack surface, the work shows that the boundary stresses on the crack are in fact zero for an unloaded crack, and the interior stresses reproduce the known inverse square root behaviour when the distance from the interior point to the crack contour approaches zero. The correct order of the stress singularity is obtained after the integrals for the self‐regular traction‐BIE formulation are evaluated analytically for the assumed displacement discontinuity model. Based on the analytic results, a new near‐crack contour self‐regular traction‐BIE is proposed for collocation points near the crack contour. In this new formulation, the asymptotic log‐singular stresses are identified and extracted from the BIE. Log‐singular stress terms are revealed for the free integrals written as contour integrals and for the self‐regularized integral with the integration region divided into sub‐regions. These terms are shown to cancel each other exactly when combined and can therefore be eliminated from the final BIE formulation. This work separates mathematical and physical singularities in a unique manner. Mathematical singularities are identified, and the singular information is all contained in the region near the crack contour. Copyright © 2003 John Wiley & Sons, Ltd.     

梯度复合材料应力强度因子计算的梯度扩展单元法

推导了一种适用于梯度复合材料断裂特性分析的梯度扩展单元, 采用细观力学方法描述材料变化的物理属性, 通过线性插值位移场给出了4节点梯度扩展元随空间位置变化的刚度矩阵, 并建立了结构的连续梯度有限元模型。通过将梯度单元的计算结果与均匀单元以及已有文献结果进行对比, 证明了梯度扩展有限元(XFEM)的优越性, 并进一步讨论了材料参数对裂纹尖端应力强度因子(SIF)的影响规律。研究结果表明: 随着网格密度的增加, 梯度单元的计算结果能够迅速收敛于准确解, 均匀单元的计算误差不会随着网格细化而消失, 且随着裂纹长度和属性梯度的增大而增大; 属性梯度和涂层基体厚度比的增大导致涂覆型梯度材料的SIF增大; 裂纹长度的增加和连接层基体厚度比的减小均导致连接型梯度材料的SIF增大。     

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.     

On the through-thickness crack with a curve front in center-cracked tension specimens

Three-dimensional finite element analyses are performed on through-thickness cracks with slightly wavy front in center-cracked plates. Considering there is an inherent relationship between the crack shape and the corresponding stress intensity factor (SIF) distribution of a crack, the curved configuration of the crack is determined using a heuristically derived iterative procedure if the SIF distribution function is known. Several simple SIF distribution functions, for instance the constant SIF distribution along the crack front, are assumed to determine the crack shape. Under the assumption that the rate of fatigue crack growth depends on the SIF range or the effective SIF range, possible effects of plate thickness, crack length and crack closure level gradient on the behaviour of crack tunneling are investigated. The stability of the curved shape of a through-thickness crack in fatigue is also discussed, i.e. whether a crack can maintain its shape satisfying the conditions of constant SIF distribution or other distribution along the crack front during fatigue growth. This study will be useful for a better understanding of the behaviour of crack tunneling and help to evaluate the validity of the two-dimensional linear elastic fracture mechanics in cracked plates.     

Fracture Analysis of Concrete Structural Components Accounting for Tension Softening Effect

This paper presents methodologies for fracture analysis of concrete structural components with and without considering tension softening effect. Stress intensity factor (SIF) is computed by using analytical approach and finite element analysis. In the analytical approach, SIF accounting for tension softening effect has been obtained as the difference of SIF obtained using linear elastic fracture mechanics (LEFM) principles and SIF due to closing pressure. Superposition principle has been used by accounting for non-linearity in incremental form. SIF due to crack closing force applied on the effective crack face inside the process zone has been computed using Green's function approach. In finite element analysis, the domain integral method has been used for computation of SIF. The domain integral method is used to calculate the strain energy release rate and SIF when a crack grows. Numerical studies have been conducted on notched 3-point bending concrete specimen with and without considering the cohesive stresses. It is observed from the studies that SIF obtained from the finite element analysis with and without considering the cohesive stresses is in good agreement with the corresponding analytical value. The effect of cohesive stress on SIF decreases with increase of crack length. Further, studies have been conducted on geometrically similar structures and observed that (i) the effect of cohesive stress on SIF is significant with increase of load for a particular crack length and (iii) SIF values decreases with increase of tensile strength for a particular crack length and load.     

Crack closure in fibre metal laminates

GLARE is a fibre metal laminate (FML) built up of alternating layers of S2-glass/FM94 prepreg and aluminium 2024-T3. The excellent fatigue behaviour of GLARE can be described with a recently published analytical prediction model. This model is based on linear elastic fracture mechanics and the assumption that a similar stress state in the aluminium layers of GLARE and monolithic aluminium result in the same crack growth behaviour. It therefore describes the crack growth with an effective stress intensity factor (SIF) range at the crack tip in the aluminium layers, including the effect of internal residual stress as result of curing and the stiffness differences between the individual layers. In that model, an empirical relation is used to calculate the effective SIF range, which had been determined without sufficiently investigating the effect of crack closure. This paper presents the research performed on crack closure in GLARE. It is assumed that crack closure in FMLs is determined by the actual stress cycles in the metal layers and that it can be described with the available relations for monolithic aluminium published in the literature. Fatigue crack growth experiments have been performed on GLARE specimens in which crack growth rates and crack opening stresses have been recorded. The prediction model incorporating the crack closure relation for aluminium 2024-T3 obtained from the literature has been validated with the test results. It is concluded that crack growth in GLARE can be correlated with the effective SIF range at the crack tip in the aluminium layers, if it is determined with the crack closure relation for aluminium 2024-T3 based on actual stresses in the aluminium layers.     

裂尖具线性分布约束应力的运动裂纹模型及其解析解

以恒定速度运动的Griffith裂纹解析解为著名的Yoffe解。静止裂纹的条状屈服模型即Dugdale模型,将其推广到运动裂纹模型时发现,当裂纹运动速度跨越Rayliegh波速时,裂纹张开位移COD趋于(∞,且表现为间断。通过在裂尖引入一个约束应力区及两个速度效应函数,假设约束应力为线性分布,采用复变函数方法,求得动态应力强度因子SIF与裂纹张开位移COD的解析解。新的结果,在Rayleigh波速下裂纹张开位移连续且为有限值。给出裂纹张开位移的一些数值结果,获得了一些有意义的结论。     

Reduction of Noise‐Induced Bias in Displacement Estimation by Linear Off‐Pixel Digital Image Correlation

Abstract: A linear digital image correlation algorithm is proposed to eliminate noise‐induced bias in one‐dimensional translation estimation using noisy images. The algorithm uses linear interpolation for both initial and current images at off‐pixel positions and solves directly the displacement parameter by minimizing a sum‐of‐squared‐differences coefficient. Both analytical results and numerical simulations using synthetic image sets show that there is indeed no noise‐induced bias in the displacement estimation using the proposed algorithm if the off‐pixel positions in both images are chosen properly according to the relative displacement between two images. When the displacement is only known initially within a range of ±0.5 pixels from the actual displacement, an iterative procedure using the algorithm is able to obtain the displacement estimation with a residual bias that converges to the noiseless subpixel approximation bias. A further refinement of the off‐pixel analysis algorithm will be needed so the remaining residual bias due to subpixel approximation can also be significantly reduced.     

Stress intensity factor analysis of friction sliding at discontinuity interfaces and junctions

A stress intensity factor (SIF) analysis for two-dimensional fractures with frictional contact (crack friction) is presented. This analysis is carried out using the symmetric-Galerkin boundary element method, and a modified quarter-point crack tip element. As in case of non-contact fracture, it is shown that highly accurate SIFs can be obtained, even with the simple Displacement Correlation SIF technique. Moreover, with the modified crack tip element, the mesh on the crack does not need to be excessively refined in order to achieve high accuracy. This meshing advantage is especially important in the context of the nonlinear frictional contact problem, as the computing time for the iterative process strongly depends on the number of elements used. Several numerical examples are presented and the SIF results are compared with available analytical or reference solutions. This research was supported in part by the University of South Alabama Research Council, and by the Applied Mathematical Sciences Research Program of the Office of Mathematical, Information, and Computational Sciences, U.S. Department of Energy under contract DE-AC05-00OR22725 with UT-Battelle, LLC.     

A new weight function approach using indirect boundary integral method

A new weight function approach to determine SIF (stress intensity factor) using the indirect boundary integral method has been presented. The crack opening displacement field was represented by one boundary integral in the form of a single-layer potential whose kernel was modified from the fundamental solution. The proposed method enables the calculation of SIF using only one SIF formula without any modification of the crack geometries symmetric in a two-dimensional plane, e.g. a center crack in a plate with or without an internal hole, double edge cracks, circumferential crack or radial cracks in a pipe. The application procedure for this variety of crack geometries is very simple and straghtforward with only one SIF formula. The necessary information in the analysis is two reference SIFs. The analysis results, using several examples, verified that the present closed-form solution was in good agreement with those of the literature and applicable to various crack geometries.     

A notch stress intensity approach applied to fatigue life predictions of welded joints with different local toe geometry

In the Notch Stress Intensity Factor (N‐SIF) approach the weld toe region is modelled as a sharp V‐shaped corner and local stress distributions in planar problems can be expressed in closed form on the basis of the relevant mode I and mode II N‐SIFs. Initially thought of as parameters suitable for quantifying only the crack initiation life, N‐SIFs were shown able to predict also the total fatigue life, at least when a large part of the life is spent as in the propagation of small cracks in the highly stressed region close to the notch tip. While the assumption of a welded toe radius equal to zero seems to be reasonable in many cases of practical interest, it is well known that some welding procedures are able to assure the presence of a mean value of the weld toe radius substantially different from zero. Under such conditions any N‐SIF‐based prediction is expected to underestimate the fatigue life. In order to investigate the degree of conservatism, a total of 128 fillet welded specimens are re‐analysed in the present work by using an energy‐based N‐SIF approach. The local weld toe geometry, characterised by its angle and radius, has been measured with accuracy for the actual test series. The aim of the work is to determine if the N‐SIF‐based model is capable of taking into account the large variability of the toe angle, and to quantify the inaccuracy in the predictions due to the simplification of setting the toe radius equal to zero.