Effect of residual stresses on crack behaviour in single edge bending specimens

Residual stresses due to manufacturing processes, such as welding, change the load bearing capacity of cracked components. The effects of residual stresses on crack behaviour in single edge bending specimens were investigated using Finite element analyses. Three parameters (J, Q and R) were used to study the crack behaviour. The J‐integral predicts the size scale over which large stresses and strains exist, the constraint parameter Q describes the crack‐tip constraint as a result of geometry, loading mode and crack depth and the constraint parameter R is used to describe the constraint resulting from residual stresses. To carry out a systematic investigation on the effect of residual stresses on the J‐integral and crack‐tip constraints, models under different combinations of residual stresses and external loads with different crack depths were analysed. It has been shown that the crack‐tip constraint R increased by tensile residual stresses around the crack‐tip. On the other hand, the constraint parameter R decreased and tended to zero at high external load levels.     

Crack-tip constraint in mode II deformation

Studies of cracked specimens loaded in mode I have shown that the stresses near the crack tip depend significantly on the level of constraint. The stresses can be determined near the crack tip using the HRR solution, but only for high constraint specimens. For other levels of constraint, O'Dowd and Shih's Q parameter may be used to adjust the stresses derived from the HRR solution. Only limited research has been carried out to study the effect of constraint in mode II. In this paper a mode II boundary layer formulation is used to study the effect of far field elastic stresses on the size and shape of the plastic zone around the crack tip and on the stresses inside the plastic zone. It is shown that in mode II, both positive and negative values of remote T-stress influence the tangential stress along the direction of maximum tangential stress. In the spirit of O'Dowd and Shih, a dimensionless parameter Q _II is introduced to quantify the constraint for mode II specimens failing by brittle fracture. The relation between Q _II and T/₀ is determined for different values of the strain hardening coefficient n. To investigate the range of validity of the Q–T diagram for real specimens, the constraint parameter Q _II is calculated directly from finite element analysis for three mode II specimens and compared with the evaluation using the Q–T diagram.     

Prediction of fatigue crack growth under constant amplitude loading and a single overload based on elasto-plastic crack tip stresses and strains

It is generally accepted that the fatigue crack growth (FCG) depends mainly on the stress intensity factor range (ΔK) and the maximum stress intensity factor (K_max). The two parameters are usually combined into one expression called often as the driving force and many various driving forces have been proposed up to date. The driving force can be successful as long as the stress intensity factors are appropriately correlated with the actual elasto-plastic crack tip stress-strain field. However, the correlation between the stress intensity factors and the crack tip stress-strain field is often influenced by residual stresses induced in due course.A two-parameter (ΔK_tot, K_max,tot) driving force based on the elasto-plastic crack tip stress-strain history has been proposed. The applied stress intensity factors (ΔK_appl, K_max,appl) were modified to the total stress intensity factors (ΔK_tot, K_max,tot) in order to account for the effect of the local crack tip stresses and strains on fatigue crack growth. The FCG was predicted by simulating the stress-strain response in the material volume adjacent to the crack tip and estimating the accumulated fatigue damage. The fatigue crack growth was regarded as a process of successive crack re-initiations in the crack tip region. The model was developed to predict the effect of the mean and residual stresses induced by the cyclic loading. The effect of variable amplitude loadings on FCG can be also quantified on the basis of the proposed model. A two-parameter driving force in the form of: was derived based on the local stresses and strains at the crack tip and the Smith-Watson-Topper (SWT) fatigue damage parameter: D = σ_maxΔε/2. The effect of the internal (residual) stress induced by the reversed cyclic plasticity manifested itself in the change of the resultant (total) stress intensity factors controlling the fatigue crack growth.The model was verified using experimental fatigue crack growth data for aluminum alloy 7075-T6 obtained under constant amplitude loading and a single overload.     

J–R curves corrected by load-independent constraint parameter in ductile crack growth

The constraint effect on J–resistance curves of ductile crack growth is considered under the condition of two-parameter J–Q^* controlled crack growth, where Q^* is a modified parameter of Q in the J–Q theory. Both J and Q^* are used to characterize the J–R curves with J as the loading level and Q^* as a constraint parameter. It is shown that Q^* is independent of applied loading under large-scale yielding or fully plastic deformation, and so Q^* is a proper constraint parameter during crack growth. An approach to correct constraint effects on the J–R curve is developed, and a procedure of transferring the J–R curves determined from standard ASTM procedure to nonstandard specimens or real cracked structures is outlined.The test data of fracture toughness, J_IC, and tearing modulus, T_R, by Joyce and Link (Engng. Fract. Mech. 57(4) (1997) 431) for a single-edge notched bend specimen with various depth cracks are employed to demonstrate the efficiency of the present approach. The variation of J_IC and T_R with the constraint parameter Q^* is obtained, and then a constraint-corrected J–R curve is constructed for the test material of HY80 steel. Comparisons show that the predicted J–R curves can match well with the experimental data for both deep and shallow cracked specimens over a reasonably large amount of crack extension.Finally, the present approach is applied to predict the J–R curves of ductile crack growth for five conventional fracture specimens. The results show that the effect of specimen geometry on the J–R curves is generally much larger than the effect of specimen sizes, and larger specimens tend to have lower crack growth resistance curves.     

RESIDUAL STRESS FIELDS DURING FATIGUE CRACK GROWTH

This study outlines the distinction between (1) residual stresses at an ideal crack tip, undergoing reversed deformation in the absence of crack closure, and (2) additional residual stresses generated due to plasticity induced closure upon fatigue crack growth. Residual stresses resulting from reversed deformation in plane strain were higher compared to the plane stress case, while residual stresses generated behind the crack tip were more significant in plane stress compared to plane strain. The origin of these residual stresses was studied for two specimen geometries over a wide range of loading conditions. We define a new crack tip parameter, S_tt as the applied stress level that corresponds to the development of tensile stresses immediately ahead of crack tips. The S_tt levels were significantly higher for a fatigue crack than for an ideal crack. We attribute the difference in S_tt levels between these two cases to plasticity induced closure. The results demonstrate the importance of the S_tt parameter, since the stresses ahead of crack tips could remain compressive even when the crack surfaces are open. Moreover, the study emphasizes the need, when describing fatigue crack growth, to incorporate both the closure concept and residual stress field ahead of a crack tip.     

Experimental and numerical analysis of in-plane and out-of-plane crack tip constraint characterization by secondary fracture parameters

In this study, several two-parameter- concepts are analyzed experimentally and numerically with respect to their capability of characterizing in-plane and out-of-plane crack tip constraint effects. Different approaches utilizing the second term T _stress of the linear-elastic crack tip stress field, a higher term A ₂ of the power-law hardening crack tip stress field, a hydrostatic correction term Q for a reference stress field or the local triaxiality parameter h are compared. Experimental results for a pressure vessel steel 22NiMoCr3-7 are investigated by means of the different approaches regarding their capability of constraint characterization for enhanced transferability. Theoretical aspects are investigated in a modified boundary layer analysis and in three-dimensional nonlinear elastic-plastic finite element analyses of the specimens. It is found that, with respect to their capability of quantifying combined in-plane and out-of-plane constraint effects, the investigated concepts differ significantly.     

A constraint based parameter for quantifying the crack tip stress fields in welded joints

By means of finite element analyses of plane strain crack tip stress fields from homogeneous and heterogeneous modified boundary layer formulations, as well as homogeneous and mismatched full field solutions, a new constraint parameter β_m has been established for overmatched welded joints, allowing the material mismatching effect on the crack tip stress fields to be quantified. In the case of complete specimens, both geometry and material mismatching affect the crack tip stress fields, and a total constraint parameter β_T can be defined. This approach allows to quantify the stress fields directly from the values of the remote applied load.     

A study of stable crack growth using experimental methods,finite elements and fractography

Systematic analysis of the in-plane constraint influence on J-resistance curves is presented. J_R curves were also recorded and analyzed beyond the limits of crack extension inside which the stress field can be assumed to be dominated by J-integral. Three steels and four types of specimen: SEN(B), SEN(T), CCT and DENT have been tested. Along with the J_R curves the fracture mechanisms have been analyzed with the help of scanning microscopy. The numerical, finite element analysis has been adopted to compute the Q-stresses, as a measure of the in-plane constraint prior to the onset of crack growth. The analysis of the stress field in front of the crack has been performed to check whether the state of stress prior to the crack growth can predetermine the way the crack will grow. It turns out that characteristic features in the J_R curves runs can be predicted qualitatively from the Q(a/W) and Q(J) curves. However, there is a good correlation between Q-stress and voids diameters on fractured surfaces. Several patterns in J_R curves runs have been observed for tested specimens; e.g. no influence of specimen thickness on J_R curves runs was observed for side-grooved specimens. Strong influence of specimen thickness on J_R curve shape was observed for non-side-grooved specimens. J_R curve run higher for thinner specimens unless they are dominated by plane stress. For bent specimens J_R curves run higher for shorter cracks but they run lower for specimens in tension.     

Quantitative characterization of creep constraint induced by crack depths in compact tension specimens

In this paper, the C(t) and C∗ integrals, stress redistribution time t_red and creep crack-tip stress distributions in the CT specimens with various crack depths have been calculated by the finite element method (FEM), and the creep constraint induced by crack depths are quantitatively investigated in detail. The results show that the creep constraint could be characterized by the new constraint parameter R. The constraint effect induced by crack depths at non-steady-state creep is more pronounced than that at steady-state creep. The effects of the crack depths, load levels (C∗) and distances from the crack tips on the creep constraint parameter R are analyzed.     

Effect of loading rate on dynamic fracture initiation toughness of brittle materials

Numerical simulation is carried out to investigate the effect of loading rate on dynamic fracture initiation toughness including the crack-tip constraint. Finite element analyses are performed for a single edge cracked plate whose crack surface is subjected to uniform pressure with various loading rate. The first three terms in the Williams’ asymptotic series solution is utilized to characterize the crack-tip stress field under dynamic loads. The coefficient of the third term in Williams’ solution, A ₃, was utilized as a crack tip constraint parameter. Numerical results demonstrate that (a) the dynamic crack tip opening stress field is well represented by the three term solution at various loading rate, (b) the loading rate can be reflected by the constraint, and (c) the constraint A ₃ decreases with increasing loading rate. To predict the dynamic fracture initiation toughness, a failure criterion based on the attainment of a critical opening stress at a critical distance ahead of the crack tip is assumed. Using this failure criterion with the constraint parameter, A ₃, fracture initiation toughness is determined and in agreement with available experimental data for Homalite-100 material at various loading rate.     

Modellierung des Ermüdungsrisswachstums in Schweißverbindungen unter Berücksichtigung von Schweißeigenspannungen

The present paper contains research results determined within the framework of a project called IBESS (?Integrale Bruchmechanische Ermittlung der Schwingfestigkeit von Schweißverbindungen“) by the Materials Mechanics Group of the Technische Universität Darmstadt [1]. Aim is to calculate the fatigue life of welded joints by taking into account the effect of residual stresses and the influence of the weld toe geometry. Here, the fatigue life is regarded as period of short fatigue crack growth. Two and three dimensional finite element models, with cracks as initial defects, are constructed for this purpose. Fatigue crack growth analyses are performed by using the node release technique together with the finite element program ABAQUS. The welding residual stresses as well as the plasticity induced crack closure effects are considered. Structural calculations are performed in order to introduce residual stress fields in finite element models. The calculated compressive residual stress field matches the measured one especially in the weld notch area. The effective cyclic J‐integral (ΔJ_eff) is used as crack tip parameter in a relation similar to the Paris equation for the calculation of the fatigue life. For this purpose, a Python code was written for the determination of ΔJ_eff at every crack length phase. The calculated fatigue lives were compared with experimental data and a good accordance between both results was achieved. The impact of welding residual stresses on ΔJ_eff as well as on the fatigue life during short crack growth was investigated. As expected, results revealed that at lower stress amplitude, a compressive residual stress field is favorable to the fatigue life, whilst a tensile residual stress field is unfavorable. The influence of residual stresses can be neglected only for large load amplitudes.     

A new constraint based fracture criterion for surface cracks

A new methodology for predicting the location of maximum crack extension along a surface crack front in ductile materials is presented. Three-dimensional elastic-plastic finite element analyses were used to determine the variations of a constraint parameter (α_h) based on the average opening stress in the crack tip plastic zone and the J-integral distributions along the crack front for many surface crack configurations. Monotonic tension and bending loads are considered. The crack front constraint parameter is combined with the J-integral to characterize fracture, the critical fracture location being the location for which the product Jα_h is a maximum. The criterion is verified with test results from surface cracked specimens.     

J-T characterized stress fields of surface-cracked metallic liners bonded to a structural backing - I. Uniaxial tension

Surface crack-tip stress fields in a tensile loaded metallic liner bonded to a structural backing are developed using a two-parameter J-T characterization and elastic-plastic modified boundary layer (MBL) finite element solutions. The Ramberg-Osgood power law hardening material model with deformation plasticity theory is implemented for the metallic liner. In addition to an elastic plate backed surface crack liner model, elastic-plastic homogeneous surface crack models of various thicknesses were tested. The constraint effects that arise from the elastic backing on the thin metallic liner and the extent to which J-T two parameter solutions characterize the crack-tip fields are explored in detail. The increased elastic constraint imposed by the backing on the liner results in an enhanced range of validity of J-T characterization. The higher accuracy of MBL solutions in predicting the surface crack-tip fields in the bonded model is partially attributed to an increase in crack-tip triaxiality and a consequent increase in the effective liner thickness from a fracture standpoint. After isolating the effects of thickness, the constraint imposed by the continued elastic linearity of the backing significantly enhanced stress field characterization. In fact, J and T along with MBL solutions predicted stresses with remarkable accuracy for loads beyond full yielding. The effects of backing stiffness variation were also investigated and results indicate that the backing to liner modulus ratio does not significantly influence the crack tip constraint. Indeed, the most significant effect of the backing is its ability to impose an elastic constraint on the liner. Results from this study will facilitate the implementation of geometric limits in testing standards for surface cracked tension specimens bonded to a structural backing.     

Surface flaws behavior under tension,bending and biaxial cyclic loading

Fatigue surface crack growth and in-plane and out-of-plane constraint effects are studied through experiments and computations for the aluminum alloy D16T. A tension/bending central notched plate and cruciform specimens under different biaxial loadings with external semi-elliptical surface cracks are studied. The variation of the fatigue crack growth rate and surface crack paths is studied under cyclic tension, bending and biaxial tension–compression loading. For the experimental surface crack paths in the tested specimens, the T-stress, out-of-plane T_z factor, local triaxiality parameter h and the governing parameter for the 3D-fields of the stresses and strains at the crack tip in the form of the I_n-integral are calculated as a function of the aspect ratio by finite element analysis to characterize the constraint effects along the semi-elliptical crack front. The plastic stress intensity factor approach is applied to the fatigue crack growth on the free surface, as well as at the deepest point of the semi-elliptical surface crack front, of the tested tension/bending plate and cruciform specimens. From the results, characteristics of the fatigue surface crack growth rate as a function of the loading conditions are established.     

Constraint effects at brittle fracture initiation in a cast ferritic steel

The fracture resistance of a cast low carbon manganese ferritic steel intended for containers for spent nuclear fuel has been analysed by combining several approaches. Based on data from three-point bend specimens with shallow and deep cracks the effect of crack tip constraint at brittle fracture initiation has been followed. Q-parameter was used for the constraint quantification. The crack length effect on the fracture toughness–temperature diagram has been analysed and peculiarities of fracture behaviour in the lower shelf region have been explained. The role of cleavage fracture stress in brittle fracture initiation under the influence of crack tip constraint has been analysed.     

Analysis of fatigue crack growth behavior of SS 316(N) welds using the unified approach

Fatigue crack growth (FCG) behavior of SS 316(N) weld has been evaluated at different R‐ratios at room temperature and compared with that of the base metal. The FCG resistance of weld is better than that of the base material and is due to the residual stresses developed during the welding. The data were analyzed using the unified approach that considers the two‐parametric (ΔK and K_max) nature of fatigue. The R‐ratio effects in both the base and weld metals are accounted for without invoking the extrinsic parameters, such as plasticity‐induced crack closure. Since the residual stresses are of the monotonic type, they affect the crack growth via the K_max‐parameter. The crack growth trajectory plots were developed, and they show how the two crack tip driving forces, ΔK and K_max, change to overcome the FCG resistance of the weld in relation to that of the base metal. The results also show that the effects from the compressive residual stresses are more dominant at low R‐values and occur via the K_max parameter.     

A parameter for quantitative analysis of plasticity induced crack closure

Numerical models have been successfully developed to predict plasticity induced crack closure (PICC). However, despite the large research effort a full understanding of the links between physical parameters, residual plastic wake and PICC has not been achieved yet. The plastic extension of material behind crack tip, Δy_p, obtained by the integration of vertical plastic deformation perpendicularly to crack flank, is proposed here to quantify the residual plastic field. The values of Δy_p and PICC were obtained numerically in a M(T) specimen using the finite element method. An excellent correlation was found between PICC and Δy_p which indicates that this parameter controls the phenomenon, and can be used to quantify the effect of physical parameters. An empirical model was developed to predict PICC assuming that the residual plastic field is a set of vertical plastic wedges, that the linear superposition principle applies and that the influence of a particular wedge exponentially decreases with distance to crack tip. The model was applied successfully to predict PICC for different residual plastic fields which provided an additional validation of Δy_p as the parameter controlling PICC.     

Two-Parameter (J-M) Description of Crack Tip Stress-Fields for an idealized weldment in small scale yielding

The crack tip stress-field in a trimaterial finite element model has been examined. The model represents an idealised steel weldment with a crack located at the fusion line. The model was loaded with a K _I displacement field to simulate small scale yielding conditions. The effect of changing the weld metal plastic properties and the HAZ layer thickness on the crack tip stress-field was studied, keeping the material properties of the HAZ and base metal constant. The results show that the calculated J-integral remains path independent in the trimaterial model. It is confirmed that the crack tip stress-fields can be normalised by the J-integral. The mismatch constraint can be characterised by a difference field, which is independent of the normalised distance from the crack tip. The results show that changes of HAZ thickness only have a small effect on the stress-fields close to the crack tip. The hardenability of the weld metal influences on the slope of the crack tip stress distribution, but for small changes in hardenability, this effect can be neglected. The results indicate that the difference fields show some radial dependence when a homogeneous reference field is used, but the radial dependence was removed by introducing an inhomogeneous reference field. The effect of changes in the weld metal yield strength has been described with a two parameter (J-M) formulation using the inhomogeneous reference field.     

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.