Fracture mechanics analysis of a crack in a residual stress field

The standard definition of the J integral leads to a path dependent value in the presence of a residual stress field, and this gives rise to numerical difficulties in numerical modelling of fracture problems when residual stresses are significant. In this work, a path independent J definition for a crack in a residual stress field is obtained. A number of crack geometries containing residual stresses have been analysed using the finite element method and the results demonstrate that the modified J shows good path-independence which is maintained under a combination of residual stress and mechanical loading. It is also shown that the modified J is equivalent to the stress intensity factor, K, under small scale yielding conditions and provides the intensity of the near crack tip stresses under elastic-plastic conditions. The paper also discusses two issues linked to the numerical modelling of residual stress crack problems-the introduction of a residual stress field into a finite element model and the introduction of a crack into a residual stress field.     

A fracture parameter for welded structures with residual stresses

In this paper, a finite element procedure to determine a fracture parameter J _ld is presented for welded structures with consideration of residual stresses. The method is based on the energy difference of two cracked solids with slightly different crack sizes. Our computational results show that J _ld and the J integral agree well for a cracked plate without consideration of residual stresses. When the residual stresses are considered, the values of J _ld for different contours close to the crack tip in the cracked plate subject to remote tensile stresses are in good agreement. The computational results also indicate that for the given residual stress distribution, the values of J _ld with consideration of residual stresses are lower than those without consideration of residual stresses for the cracked plate subject to large remote tensile stresses.     

Analytical stress intensity factor solutions for resistance and friction stir spot welds in lap-shear specimens of different materials and thicknesses

In this paper, analytical stress intensity factor and J integral solutions for resistance and friction stir spot welds without and with gap and bend in lap-shear specimens of different materials and thicknesses are developed. The J integral and stress intensity factor solutions for spot welds are first presented in terms of the structural stresses for a strip model. Analytical structural stress solutions for spot welds without and with gap and bend in lap-shear specimens are then developed based on the closed-form structural stress solutions for a rigid inclusion in a finite thin plate subjected to various loading conditions. With the available structural stress solutions, the analytical J integral and stress intensity factor solutions can be obtained as functions of the applied load, the elastic material property parameters, and the geometric parameters of the weld and specimen. The analytical stress intensity factor solutions are selectively validated by the results of three-dimensional finite element analyses for a spot weld with ideal geometry and for a friction stir spot weld with complex geometry, gap and bend. The stress intensity factor and J integral solutions at the critical locations of spot welds in lap-shear specimens of dissimilar magnesium, aluminum and steel sheets with equal and different thicknesses are then presented in the normalized forms as functions of the ratio of the specimen width to the weld diameter. Finally, general trends and simple estimation methods of the stress intensity factor and J integral solutions at the critical locations of spot welds in lap-shear specimens of different materials and thicknesses are given for convenient engineering applications.     

Residual stresses and fracture mechanics analysis of a crack in welds of high strength steels

This study presented the characteristics of residual stresses in welds of high strength steels (POSTEN60, POSTEN80) whose tensile strengths were 600 MPa and 800 MPa, respectively. Three-dimensional thermal elastic-plastic analyses were conducted to investigate the characteristics of welding residual stresses in welds of high strength steels through the thermal and mechanical properties at high temperatures obtained from the elevated temperature tensile tests. A finite element analysis method which can calculate the J-integral for a crack in a residual stress field was developed to evaluate the J-integral for a centre crack when mechanical stresses were applied in conjunction with residual stresses.The results show that the volumetric changes associated with the austenite to martensite phase transformation during rapid cooling after welding of high strength steels significantly influence on the development of residual stresses in the weld fusion zone and heat-affected zone. For a centre crack in welds of high strength steels where only residual stresses are present, increased tensile strength of the steel, increased the J-integral values. The values of the J-integral for the case when mechanical stresses are applied in conjunction with residual stresses are larger than those for the case when only residual stresses are present.     

Approximate J estimates for tension-loaded plates with semi-elliptical surface cracks

This paper provides a simplified engineering J estimation method for semi-elliptical surface cracked plates in tension, based on the reference stress approach. Note that the essential element of the reference stress approach is the plastic limit load in the definition of the reference stress. However, for surface cracks, the definition of the limit load is ambiguous (“local” or “global” limit load), and thus the most relevant limit load (and thus reference stress) for the J estimation should be determined. In the present work, such limit load solution is found by comparing reference stress based J results with those from extensive 3-D finite element (FE) analyses. Based on the present FE results, the global limit load solution proposed by Goodall for surface cracked plates in combined bending and tension was modified, in the case of tension loading only, to account for a weak dependence on w/c and was defined as the reference normalizing load. Validation of the proposed equation against FE J results based on actual experimental tensile data of a 304 stainless steel shows excellent agreements not only for the J values at the deepest point but also for those at an arbitrary point along the crack front, including at the surface point. Thus the present results provide a good engineering tool for elastic-plastic fracture analyses of surface cracked plates in tension.     

Effects of weld geometry and sheet thickness on stress intensity factor solutions for spot and spot friction welds in lap-shear specimens of similar and dissimilar materials

In this paper, three-dimensional finite element analyses for spot welds with ideal geometry in lap-shear specimens of different materials and thicknesses were first conducted. The computational results indicate that the stress intensity factor and J integral solutions based on the finite element analyses agree well with the analytical solutions and that the analytical solutions can be used with a reasonable accuracy. Three-dimensional finite element analyses based on the micrographs of an aluminum 6111 resistance spot weld, an aluminum 5754 spot friction weld, and a dissimilar Al/Fe spot friction weld were also conducted. The computational results indicate that the stress intensity factor and J integral solutions based on the finite element analyses for the aluminum 6111 resistance spot weld and aluminum 5754 spot friction weld with complex geometry are in good agreement with the analytical solutions for the equivalent spot welds with ideal geometry. However, the stress intensity factor and J integral solutions based on the finite element analysis for the Al/Fe spot friction weld with complex geometry are completely different from the analytical solutions for the equivalent spot weld with ideal geometry. Different three-dimensional finite element analyses based on the meshes that represent different features of the complex geometry of the Al/Fe spot friction weld were then conducted. The computational results indicate that the stress intensity factor and J integral solutions for the Al/Fe spot friction weld based on the finite element analysis agree reasonably well with the analytical solutions for the equivalent spot weld with consideration of gap and bend. The computational and analytical results suggest that the stress intensity factor and J integral solutions based on the finite element analysis and the analytical solutions with consideration of gap and bend may be used to correlate with the fatigue crack growth patterns of Al/Fe spot friction welds observed in experiments.     

Creep-fatigue crack propagation in lead-free solder under cyclic loading with various waveforms

Crack propagation tests of lead-free solder were conducted using center-notched plate specimens under cyclic tension-compression of three load waveforms: pp waveform having fast loading and unloading, cp-h waveform having a hold time under tension, and cc-h waveform having a hold time under tension and compression. In the case of fatigue loading, i.e. pp waveform, the path of crack propagation was macroscopically straight and perpendicular to the maximum principal stress direction, showing tensile-mode crack propagation. The introduction of the creep components by hold time in cc-h and cp-h waveforms promoted shear-mode crack propagation. For fatigue loading of pp wave, the crack propagation rate was expressed as a power function of the fatigue J integral and the relation was identical for load-controlled and displacement-controlled conditions. The creep component due to the hold time greatly accelerates the crack propagation rate when compared at the same values of the fatigue J integral or the total J integral (the sum of fatigue J and creep J integrals). The creep crack propagation rate was expressed as a power function of the creep J integral for each case of cp-h and cc-h waveforms. The crack propagation rate for cp-h waveform is higher than that for cc-h waveform. The predominant feature of fracture surfaces was striations for pp waveform and grain boundary fracture for cp-h waveform. Grain fragmentation was abundantly observed on the fracture surface made under cc-h waveform.     

The energy dissipation rate approach to tearing instability

Stable and unstable tearing in metals is currently analysed by J integral theory, or by the G_R curve approach. This paper explains an alternative analysis route based on energy dissipation rate, D. It is shown that the implication of increasing toughness with crack growth in G_R and J_R curves is misleading. Even in small scale yielding (SSY), it is possible to have stable tearing under increasing G or J whilst at the same time D is constant (or even reducing) with crack growth. New terms: C for crack driving force, D^* for geometry normalised D, D_ssy for D in SSY, and crack stability index are explained. A D based fracture analysis diagram is introduced. Comparisons are made between energy dissipation rate, J integral, and G_R curve instability prediction methods. It is shown that, in most instances, these different approaches are compatible; but that the use of J_R curves derived in fully yielded test pieces to predict failure in SSY has the potential to lead to an unconservative instability prediction. The practical advantage of the energy dissipation rate approach is that it can be applied to all product thicknesses at any extent of crack growth. The major advantage compared to the G_R approach is that toughness measurements can be made on much smaller specimens.     

DIC‐based J‐integral evaluation of laser repaired cracks with micro/nanomaterial addition

Cracks generated by a wire electric discharge machining on compact tension specimens are first repaired by high powered laser beam with different weight fractions of nano‐WC added at the crack tip. Digital image correlation method combined with J‐integral theory is used to measure and calculate J‐integrals of the repaired specimens. Fracture properties of the repaired specimens are compared and studied. The residual stress of specimens after being repaired by laser was also studied to analyze the feasibility of J‐integral to evaluate the fracture properties of the repaired parts. It is found that the influence of residual stress can be neglected when calculating J‐integral in a certain region in case of when the residual stress is small. The J‐integral obtained by the digital image correlation measurement is accurate and can effectively evaluate repairing effects. The paper provides guidance for laser repaired cracks with nanomaterial addition and an effective method for the evaluation of repairing effect.     

Advances in J‐integral estimation of circumferentially surface cracked pipes

This paper describes enhanced J‐integral estimation schemes for pipes with circumferential semi‐elliptical cracks subjected to tensile loading, global bending and internal pressure. These schemes are given in two different forms to cover the wide ranges of geometries and material parameters; the modified GE/EPRI method and the modified reference stress method. In the former method, new plastic influence functions for fully plastic J‐integral estimation are developed based on extensive three‐dimensional finite element calculations. In the latter method, new optimized reference loads are suggested and utilized to predict the J values. To verify the feasibility of these two schemes, J‐integral values obtained from further detailed FE analyses are compared to those from the proposed schemes. Because the estimated J‐integrals agree fairly well with the detailed FE analysis results, the new solutions can be applied for accurate structural integrity assessment of different size pipes with a circumferential surface crack.     

Approximate elastic-plastic J estimates of cylinders with off-centred circumferential through-wall cracks

This paper provides approximate J estimates for off-centred, circumferential through-wall cracks in cylinders under bending and under combined tension and bending. The proposed method is based on the reference stress approach, where the dependence of elastic and plastic influence functions of J on the cylinder/crack geometry, the off-centred angle and strain hardening is minimised through the use of a proper normalising load. Based on published limited FE results for off-centred, circumferential through-wall cracks under bending, such normalising load is found, based on which the reference stress based J estimates are proposed for more general cases, such as for a different cylinder geometry and for combined loading. Comparison of the estimated J with extensive FE J results shows overall good agreements for different crack/cylinder geometries and for combined tension and bending, which provides sufficient confidence in the use of the proposed method for fracture mechanics analyses of off-centred circumferential cracks. Furthermore, the proposed method is simple to use, giving significant merits in practice.     

Three-dimensional stress intensity factor calibration for a stiffened cracked plate

Three-dimensional finite element analyses are used in this paper to calibrate the stress intensity factor in a cracked stiffened plate subjected to remote uniform traction. An accurate numerical determination of the stress field and stress intensity variation through the thickness of a central cracked plate was first carried out in order to evaluate three-dimensional effects. A stiffened cracked plate was then analysed, taking into account the results and the conclusions obtained in the previous study. Such a structure was chosen due to the growing interest for large integral metallic structures for aircraft applications, following the continuous need for low cost and the emergence of new technologies. The J-Integral technique was used to calculate the values of the stress intensity factor along the plate thickness. The plane strain behaviour near the crack front and the variation of the opening stress are discussed.     

Mode I–III Decomposition of the J‐integral from DIC Displacement Data

This paper presents the implementation of the decomposition method on digital image correlation (DIC) obtained displacement fields to obtain J‐integral results (J) and respective stress intensity factors (SIFs). DIC is increasingly used with the J‐integral approach in experimental mechanics to obtain J estimates from complex fracture processes. In this approach, the decomposition method is applied to DIC displacement fields for the first time. Here, displacement fields are separated before stresses and strains are computed, so that subsequent computation of separate J or SIF components may follow the classical full‐field J‐integral approach. The sensitivity of the decomposition method to experimental errors is investigated using synthetically generated errors imposed on crack tip displacement fields (Williams' series), from which improvements to the procedure are proposed. The method is experimentally tested on PMMA Arcan specimens under mode I, II, and III, and mixed‐mode I–III loading. Test results were compared to fracture toughness values obtained from ASTM tests and literature with close agreement.     

Evaluation of polymer fracture parameters by the boundary element method

The boundary element method (BEM) for two-dimensional linear viscoelasticity is applied to polymer fracture. The time-dependence of stress intensity factors is assessed for various viscoelastic models as well as loading and support conditions. Various representations of the energy release rate under isothermal conditions are adopted. Additional boundary integral equations for the displacement gradient in the domain are linked to algorithms for the evaluation of path-independent J-integrals. The consistency of BEM predictions and their good agreement with other analytical results confirms BEM as a valid modelling tool for viscoelastic fracture characterisation and failure assessment under complex geometric and loading conditions.     

An extended equivalent domain integral method for mixed mode fracture problems by the p-version of FEM

A new path-independent contour integral formula is presented to estimate the crack-tip integral parameter, J-value, for two-dimensional cracked elastic bodies which may quantify the severity of the crack-tip stress fields. The conventional J-integral method based on a line integral has been converted to an equivalent area or domain integral (EDI) by the divergence theorem. It is noted that the EDI method is very attractive because all the quantities necessary for computation of domain integrals are readily available in a finite element analysis. The details and its implementation are extended to the p-version FE model with hierarchic elements using integrals of Legendre polynomials. By decomposing the displacement field obtained from the p-version finite element analysis into symmetric and antisymmetric displacement fields with respect to the crack line, the Mode-I and Mode-II non-dimensional stress intensity factors can be determined by using the decomposition method. The example problems for validating the proposed techniques are centrally oblique cracked plates under tensile loading. The numerical results associated with the variation of oblique angles show very good agreement with the existing solutions. Also, the selective distribution of polynomial orders and the corner elements for automatic mesh generation are applied to improve the numerical solution in this paper. © 1998 John Wiley & Sons, Ltd.     

Approximate J estimates for circumferential cracked pipes under primary and secondary stresses

This paper presents finite element solutions for elastic-plastic J for circumferentially cracked pipes under combined mechanical and thermal loads in terms of the V/V_o factor used within the failure assessment diagram approach. Systematic analyses suggest that the two major variables affecting V/V_o are the relative magnitude of the secondary stress and the primary load magnitude. It has been found that the variations with these parameters obtained from the FE results agree well with current R6 predictions for modest thermal loads and low primary loads. For larger thermal stresses, the R6 predictions are overly conservative and more accurate predictions are suggested.     

Crack growth resistance behavior of a functionally graded material: computational studies

This paper describes crack growth resistance simulation in a ceramic/metal functionally graded material (FGM) using a cohesive zone ahead of the crack front. The plasticity in the background (bulk) material follows J₂ flow theory with the flow properties determined by a volume fraction based, elastic-plastic model (extension of the original Tamura-Tomota-Ozawa model). A phenomenological, cohesive zone model with six material-dependent parameters (the cohesive energy densities and the peak cohesive tractions of the ceramic and metal phases, respectively, and two cohesive gradation parameters) describes the constitutive response of the cohesive zone. Crack growth occurs when the complete separation of the cohesive surfaces takes place. The crack growth resistance of the FGM is characterized by a rising J-integral with crack extension (averaged over the specimen thickness) computed using a domain integral (DI) formulation. The 3-D analyses are performed using WARP3D, a fracture mechanics research finite element code, which incorporates solid elements with graded elastic and plastic properties and interface-cohesive elements coupled with the functionally graded cohesive zone model. The paper describes applications of the cohesive zone model and the DI method to compute the J resistance curves for both single-edge notch bend, SE(B), and single-edge notch tension, SE(T), specimens having properties of a TiB/Ti FGM. The numerical results show that the TiB/Ti FGM exhibits significant crack growth resistance behavior when the crack grows from the ceramic-rich region into the metal-rich region. Under these conditions, the J-integral is generally higher than the cohesive energy density at the crack tip even when the background material response remains linearly elastic, which contrasts with the case for homogeneous materials wherein the J-integral equals the cohesive energy density for a quasi-statically growing crack.     

J-integral evaluation for cases involving non-proportional stressing

Calculation of J for cases where the proportional stressing condition cannot be satisfied is investigated. A modified J definition is derived and implemented into an ABAQUS post-processing program for both 2-D and axisymmetric problems. The modified J-integral is path independent for cases of proportional and non-proportional stressing. For cases with proportional stressing, the modified integral gives the same value as does the standard ABAQUS J function. It is also found that the modified J is equivalent to the stress intensity factor for a linear elastic material and provides a measure of the intensity of the crack-tip fields for non-linear elastic and elastic-plastic materials. The modified J formulation is applied to the case of a cylinder with an external circumferential crack under various load conditions.     

Interaction of residual stress with mechanical loading in a ferritic steel

A well-defined residual stress field was introduced into modified single edge notched bend, SEN(B), specimens by the ‘in-plane compression’ procedure in order to investigate the interaction between residual stress and applied mechanical loading. Numerical predictions of the residual stress field arising from the in-plane compression procedure are given along with details of the numerical fracture modelling and experimental fracture test results made on A533B ferritic steel specimens in the lower transition region at −150 °C. Use was made of a recently developed finite element post-processor capable of determining path-independent J-integral values in the presence of residual stress fields. The paper compares the experimental results to predictions made using a probabilistic ‘global approach’ based on the conventional crack-tip parameters K and J and predictions made using a well-known structural integrity assessment code, R6 (Revision 4). It is shown that obtaining more accurate estimates of the crack driving force created by residual stresses leads to better correlation between experiments and predictions, and less conservatism in the assessment code.     

Analytical and numerical study of cohesive zones for a crack in an adhesive layer between identical isotropic materials

A crack in a thin adhesive elastic-perfectly plastic layer between two identical isotropic elastic half-spaces is considered. Uniformly distributed normal stress is applied to the substrates at infinity. First, stress distribution in the cohesive zones and the J-integral values are defined numerically by the finite element method (FEM). Further, a mathematical formulation of the problem is given and its analytical solution is proposed. It is assumed that, at the crack continuations, there exist cohesive zones. The interlayer thickness is neglected since it is much smaller than the crack length. The distribution of the normal stress, which was obtained by means of the FEM, is now approximated by a piecewise-constant function and assumed to be applied at the faces of the cohesive zones. The formulated problem is solved analytically and an equation for determination of the cohesive zone lengths is derived. Also, closed expressions for the crack tip opening displacement and for the J-integral are obtained in an analytical form. These parameters are found with respect to the values of the normal stress applied at infinity. Finally, a universal approximating function, which describes the stress distribution in the cohesive zones, is constructed. This function depends on the ratio between the interlayer thickness and the crack length and on the ratio between the normal stress applied at infinity and the yield limit of the interlayer’s material. Once again, the problem is solved analytically, but this time for the stress distribution prescribed by the universal approximating function. The cohesive zone lengths, the values of the crack tip opening displacement and of the J-integral are calculated. A comparative analysis of the obtained results is carried out. A good agreement of the J-integral values calculated by means of the developed analytical models and by the associated finite element analysis is demonstrated.