Material constraint parameters for the assessment of shallow defects in structural components. Part II: Constraint-based assessment of shallow defects

This Part II paper describes a series of constraint-based R6 assessments of shallow-cracked specimens under equibiaxial loading using material constraint parameters obtained from lookup tables presented in the companion Part I paper. Parameters are derived for an A533B-1 steel forging from a knowledge of the yield and flow behaviour of the material and a calibration of the Beremin model parameter m. Resulting assessments, which use both T- and Q-stress to quantify constraint, are found to be conservative with respect to the available experimental data. The results are also used to predict the influence of T-stress on the Master Curve reference temperature. Available data are broadly in agreement with the predicted trend curve. Finally, the results demonstrate that out-of-plane biaxial loading increases constraint to the extent that the inherent conservatism of the elastic T-stress at high L_r is eroded. Out-of-plane constraint effects are only properly accounted for by the hydrostatic Q-stress and for this reason, the use of T to assess biaxial loading situations should be undertaken with a degree of caution.     

Experimental T33-stress formulation of test specimen thickness effect on fracture toughness in the transition temperature region

This paper describes a study of the test specimen thickness effect on fracture toughness of a material, in the transition temperature region, for CT specimens. In addition we studied the specimen thickness effect on the T₃₃-stress (the out-of-plane non-singular term in the series of elastic crack-tip stress fields), expecting that T₃₃-stress affected the crack-tip triaxiality and thus constraint in the out-of-plane direction. Finally, an experimental expression for the thickness effect on the fracture toughness using T₃₃-stress is proposed for 0.55% carbon steel S55C. In addition to the fact that T₃₃ (which was negative) seemed to show an upper bound for large B/W, these results indicate the possibility of improving the existing methods for correlating fracture toughness obtained by test specimen with the toughness of actual cracks found in the structure, using T₃₃-stress.     

The T-stress solutions for through-wall circumferential cracks in cylinders subjected to general loading conditions

The elastic T-stress is a parameter used to define the level of constraint at a crack tip. It is important to provide T-stress solutions for practical geometries to apply the constraint-based fracture mechanics methodology. In the present work, T-stress solutions are provided for circumferential through-wall cracks in thin-walled cylinders. First, cylinders with a circumferential through-wall crack were analyzed using the finite element method. Three cylinder geometries were considered; defined by the mean radius of the cylinder (R) to wall thickness (t) ratios: R/t = 5, 10, and 20. The T-stress was obtained at eight crack lengths (θ/π = 0.0625, 0.1250, 0.1875, 0.2500, 0.3125, 0.3750, 0.4375, and 0.5000, θ is the crack half angle). Both crack face loading and remote loading conditions were considered including constant, linear, parabolic and cubic crack face pressures and remote tension and bending. The results for constant and linear crack face pressure were used to derive weight functions for T-stress for the corresponding cracked geometries. The weight functions were validated against several linear and non-linear stress distributions. The derived weight functions are suitable for T-stress calculations for circumferential cracks in cylinders under complex stress fields.     

The effective T-stress estimation and crack paths emanating from U-notches

The concept of the T-stress as a local constraint factor has been extended to U-notch tip stress distribution as the effective T-stress. The effective T-stress has been estimated as the average value of the T-stress in the region corresponding to the effective (characteristic) distance ahead of the notch tip. The T-stress is evaluated by finite element method using the experimental load for crack initiation and computing the difference between principal stresses along ligament. A large range of critical effective T-stress values is investigated for different specimen configurations and notch aspect ratios. Crack stabilisation and crack bifurcation for fracture emanating from notches according to the critical effective T-stress is discussed. A model involving the influence of the critical effective T-stress on void growth for ductile failure in the vicinity of the notch tip has been proposed.     

Estimations of the T-stress for small cracks at notches

The T-stress is increasingly being recognized as an important additional stress field characterizing parameter in the analyses of cracked bodies. Using T-stress as the constraint parameter, the framework of failure assessments including the constraint effect has been established; and the effect of T-stress on fatigue crack propagation rate has been investigated by several researchers. In this paper, a simple method for determining the T-stress for small notch-emanating cracks is presented. First, the background on the T-stress calculation using the superposition principle and the similarities between the elastic notch-tip stress fields described by two parameters: the stress concentration factor (K_t) and the notch-tip radius (ρ), are summarized. Then, the method of estimating T-stress for both short and long cracks at the notches is presented. The method is used to predict T-stress solutions for cracks emanating from an internal hole in a wide plate, and cracks emanating from an U-shaped edge notch in a finite thickness plate. The results are compared to the T-stress results in the literature, and the T-stresses solutions obtained from finite element analysis. Excellent agreements have been achieved for small cracks. The method presented here can be used for a variety of notch crack geometries and loading conditions.     

T-stress of cracks loaded by near-tip tractions

T-stress solutions were derived for tractions acting on the crack-faces near a crack tip. Such solutions are of interest for the determination of the leading term of a weight function representation of T-stresses and the computation of an “intrinsic” T-stress for cracks growing in a material with a rising crack growth resistance. First, the type of the Green’s function for T-stresses is theoretically established. Then, results of finite element computations are reported for edge-cracked bars, DCB and CT specimens, which are suited for the determination of the first series term. As an application of the Green’s functions, the T-stresses caused by bridging interactions very close to the crack tip are computed.     

Quasicontinuum simulation of crack propagation in bcc-Fe

We have investigated fracture in bcc-Fe through multiscale simulations. The quasicontinuum (QC) method with an embedded atom method (EAM) interatomic potential is applied. The analyses have been carried out assuming different crystallographic orientations and different T-stress under Mode I loading. Both anisotropic and isotropic formulations of the modified boundary layer (MBL) approach has here been investigated and compared. The results show that the mechanisms at the crack tip and the critical stress intensity factor K_Ic are sensitive to both the crystallographic orientation and whether or not the formulation of the boundary conditions are isotropic or anisotropic. Mechanisms such as cleavage crack propagation, twinning, and dislocation emission are observed in the analyses. A short literature review on atomistic and multiscale simulations of fracture in bcc-Fe has been performed and evaluated, and also compared with the current results.     

Cracked Brazilian disc specimen subjected to mode II deformation

The centrally cracked Brazilian disc specimen has been used by many researchers to study mode I and mode II brittle fracture in different materials. However, the experimental results obtained in the past from this specimen indicate that the fracture toughness ratio (K_IIc/K_Ic) is always significantly higher than the theoretical predictions. It is shown in this paper that the increase in the ratio K_IIc/K_Ic can be predicted if a modified maximum tangential stress (MTS) criterion is used. The modified criterion takes into account the effect of T-stress in addition to the conventional singular stresses. The fracture toughness ratio K_IIc/K_Ic is calculated for two brittle materials using the modified criterion and is compared with the relevant published experimental results obtained from fracture tests on the cracked Brazilian disc specimen. A very good agreement is shown to exist between the theoretical predictions and the experimental results.     

Elastic T-stress solutions for semi-elliptical surface cracks in finite thickness plates

Three-dimensional finite element analyses have been conducted to calculate the elastic T-stress for semi-elliptical surface cracks in finite thickness plates. Far-field tension and bending loads were considered. The analysis procedures and results were verified using both exact solutions and approximate solutions. The T-stress solutions are presented along the crack front for cracks with a/t values of 0.2, 0.4, 0.6 or 0.8 and a/c values of 0.2, 0.4, 0.6 or 1.0. Based on the present finite element calculations for T-stress, empirical equations for the T-stress at three locations: the deepest, the surface and the middle points of the crack front under tension or bending are presented. The numerical results are approximated by empirical formulae fitted with an accuracy of 1% or better. They are valid for 0.2?a/c?1 and 0?a/t?0.8. These T-stress results together with the corresponding K or J values for surface cracks are suitable for the analysis of constraint effects for surface cracked components.     

Elastic T-stress solutions for flat elliptical cracks under tension and bending

Exact solutions for elastic T-stress of a flat elliptical crack in an infinite body under tension and bending are obtained in this paper. Many papers have been devoted to the problems for elliptical cracks in an elastic medium, but all their attention has been concentrated on the determination of stress intensity factors. In the current paper, elastic T-stress solutions are derived by means of the potential method and a specific collection of harmonic functions. The formulas of the elastic T-stress for a penny-shaped crack [Wang X. Elastic T-stress solutions for penny-shaped cracks under tension and bending. Engng Fract Mech 2004;71:2283-98] follow from the present results as a special case. It is obtained that under tension loading, the elastic T-stress is always compressive along the elliptical crack front. In both tension and bending cases, T-stress essentially depends on the Poisson’s ratio of the material, a parametric angle and semi-axes of the ellipse.     

Elastic T-stress for cracks in test specimens subjected to non-uniform stress distributions

Finite element analyses have been conducted to calculate elastic T-stress solutions for cracked test specimens. The T-stress solutions are presented for single edge cracked plates, double edge crack plates and centre cracked plates. Uniform, linear, parabolic or cubic stress distributions were applied to the crack face. The results for uniform and linear stress distributions were used to derive weight functions for T-stress for the corresponding specimens. The weight functions for T-stress are then verified against several linear and non-linear stress distributions. The derived weight functions are suitable for the T-stress calculation for cracked specimens under any given stress field.     

T-stress effects on steady crack growth in a thin, ductile plate under small-scale yielding conditions: Three-dimensional modeling

The non-singular T-stress provides a first-order estimate of geometry and loading mode, e.g. tension vs. bending, effects on elastic–plastic, crack-front fields under mode I conditions. The T-stress has a pronounced effect on measured crack growth resistance curves for ductile metals – trends most computational models confirm using a two-dimensional setting. This work examines T-stress effects on three-dimensional (3D), elastic–plastic fields surrounding a steadily advancing crack for a moderately hardening material in the framework of a 3D, small-scale yielding boundary-layer model. A flat, straight crack front advances at a constant quasi-static rate under near invariant local and global mode I loading. The boundary-layer model has thickness B that defines the only geometric length-scale. The material flow properties and (local) toughness combine to limit the in-plane plastic-zone size during steady growth to at most a few multiples of the thickness (conditions obtainable, for example, in large, thin aluminum components). The computational model requires no crack growth criterion; rather, the crack front extends steadily at constant values of the plane-stress displacements imposed on the remote boundary for the specified far-field stress intensity factor and T-stress. The specific numerical results presented demonstrate similarity scaling of the 3D near-front stresses in terms of two non-dimensional loading parameters. The analyses reveal a strong effect of T-stress on key stress and strain quantities for low loading levels and less effect for higher loading levels, where much of the plastic zone experiences plane-stress conditions. To understand the combined effects of T-stress on stresses and plastic strain levels, normalized values from a simple void-growth model, computed over the crack plane for low loading, clearly reveal the tendency for crack-front tunneling, shear-lip formation near the outside surfaces, and a minimum steady-state fracture toughness for T = 0 loading.     

Quantifying crack tip displacement fields with DIC

Crack paths under both fatigue and fracture conditions are governed by the crack tip displacement field and the material deformation characteristics, including those influenced by metallurgical anisotropy. Experimental techniques such as thermoelasticity and photoelasticity have been successfully used to characterise the elastic stress fields around cracks but they do not take into account either plasticity or anisotropy. Considerable work has been carried out to characterise crack tip stress fields from displacement measurements. The current method of choice for obtaining displacement field data is digital image correlation (DIC) which has undergone significant advances in the recent years. The ease of use and capabilities of the technique for full field displacements has led to improved methods for characterising crack tip displacement fields based on data obtained from DIC. This paper gives an overview of some of the applications of DIC for crack tip characterisation such as K, T-stress and crack tip opening angle (CTOA) measurements as well as data obtained from 3D measurements of a propagating crack.     

Two-parameter characterization of elastic-plastic crack front fields: Surface cracked plates under tensile loading

In this paper the J-Q two-parameter characterization of elastic-plastic crack front fields is examined for surface cracked plates under uniaxial and biaxial tensile loadings. Extensive three-dimensional elastic-plastic finite element analyses were performed for semi-elliptical surface cracks in a finite thickness plate, under remote uniaxial and biaxial tension loading conditions. Surface cracks with aspect ratios a/c = 0.2, 1.0 and relative depths a/t = 0.2, 0.6 were investigated. The loading levels cover from small-scale to large-scale yielding. In topological planes perpendicular to the crack fronts, the crack stress fields were obtained. In order to facilitate the determination of Q-factors, modified boundary layer analyses were also conducted. The J-Q two-parameter approach was then used in characterizing the elastic-plastic crack front stress fields along these 3D crack fronts. Complete distributions of the J-integral and Q-factors for a wide range of loading conditions were obtained. It is found that the J-Q characterization provides good estimate for the constraint loss for crack front stress fields. It is also shown that for medium load levels, reasonable agreements are achieved between the T-stress based Q-factors and the Q-factors obtained from finite element analysis. These results are suitable for elastic-plastic fracture mechanics analysis of surface cracked plates.     

T-Stress solutions for a semi-elliptical axial surface crack in a cylinder subjected to mode-I non-uniform stress distributions

This paper presents the T-stress solutions (T₁₁ and T₃₃) for semi-elliptical axial surface cracks in a cylinder subjected to mode-I non-uniform stress on the crack surface. Two cylindrical geometries with inner radius (R_i) to wall thickness (t) ratios R_i/t = 5 and 10 were considered. The T-stresses were applied along the crack front for normalized crack depth values a/t of 0.2, 0.4 and 0.5 and aspect ratios a/c of 0.2, 0.4, 0.6 and 1.0. Three stress distribution; uniform, linear and parabolic were applied to the crack face. In addition to these solutions, concrete formulation of the superposition principle is given for the T₃₃-stress, which is known as an elastic parameter that describes the out-of-plane crack tip constraint effect. Then, the validity of the formulation was shown through application of our T-stress solutions to the problem of an axial semi-elliptical surface crack in a cylinder subjected to internal pressure, and checking that the principle of superposition holds for the problem.     

Methods for including constraint effects within the SINTAP procedures

To enable constraint-based approaches to be applied to the practical assessment of defective components, development work within SINTAP has addressed a number of issues: constraint parameters for surface defects, constraint parameters for secondary stresses, simplified weight function methods for calculation of T-stress and relevant fracture toughness data. This paper describes these developments and indicates how the methods have been included in the overall SINTAP procedures.     

A finite element analysis of stable crack growth in an aluminium alloy

Extensive stable cracking has been observed in large test pieces of 25 mm thick weldable AlMgZn alloy which is used in the construction of a portable bridge. Standard fracture specimens produced valid K_IC values, with short cracks exhibiting unstable fracture. Finite element analysis of the large specimens determined a valid J_R-curve that can increase the effective K_C by several times the K_IC value. The R-curve has an unusual ‘concave’ shape that is associated with the change from initially flat fracture to fully slant fracture. The early stages of the R-curve are affected by in-plane constraint that can be indexed by the T-stress. The R-curve can be used to explain the stability of long cracks in full-scale tests on a bridge prototype, compared with the instability of short cracks in small, standard test pieces.     

Mode I fatigue crack growth under biaxial loading

This paper is centred on the role of the T-stress during mode I fatigue crack growth. The effect of a T-stress is studied through its effect on plastic blunting at crack tip. As a matter of fact, fatigue crack growth is characterized by the presence of striations on the fracture surface, which implies that the crack grows by a mechanism of plastic blunting and re-sharpening (Laird C. The influence of metallurgical structure on the mechanisms of fatigue crack propagation. In: Fatigue crack propagation, STP 415. Philadelphia: ASTM; 1967. p. 131–68 [8]). In the present study, plastic blunting at crack tip is a global variable ρ, which is calculated using the finite element method. ρ is defined as the average value of the permanent displacement of the crack faces over the whole K-dominance area. The presence of a T-stress modifies significantly the evolution of plastic deformation within the crack tip plastic zone as a consequence of plastic blunting at crack tip. A yield stress intensity factor K_Y is defined for the cracked structure, as the stress intensity factor for which plastic blunting at crack tip exceeds a given value. The variation of the yield stress intensity factor was studied as a function of the T-stress. It is found that the T-stress modifies significantly the yield point of the cracked structure and that the yield surface in a (T, K_I) plane is independent of the crack length. Finally, a yield criterion is proposed for the cracked structure. This criterion is an extent of the Von-Mises yield criterion to the problem of the cracked structure. The proposed criterion matches almost perfectly the results obtained from the FEM. The evolution of the yield surface of the cracked structure in a (T, K_I) plane was also studied for a few loading schemes. These results should develop a plasticity model for the cracked structure taking into account the effect of the T-stress.     

Mode I cracks subjected to large T-stresses

There are several criteria for predicting brittle fracture in mode I and mixed mode loading. In this paper, the modified maximum tangential stress criterion originally proposed for mixed mode loading, is employed to study theoretically brittle fracture for mode I cracks. In particular, the effect of the non-singular term of stress, often known as the T-stress, on the angle of initiation of fracture and the onset of crack growth is explored. The T-stress component of the tangential stress vanishes along the crack line. Therefore, it is often postulated for linear elastic materials that the effect of T-stress on mode I brittle fracture can be ignored. However, it is shown here that the maximum tangential stress is no longer along the line of initial crack when the T-stress exceeds a critical value. Thus, a deviation in the angle of initiation of fracture can be expected for specimens having a large T-stress. It is shown that the deviation angle increases for larger values of T-stress. Theoretical results show that the apparent fracture toughness decreases significantly when a deviation in angle occurs. Earlier experimental results are used to corroborate the findings. The effect of large T-stresses is also explored for a crack specimen undergoing moderate scale yielding. The elastic-plastic investigation is conducted using finite element analysis. The finite element results reveal a similar deviation in the angle of maximum tangential stress for small to moderate scale yielding.