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.     

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.     

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.     

Numerical modelling of crack path in the lubricated rolling–sliding contact problems

A two-dimensional computational model for simulation of contact fatigue of lubricated rolling–sliding contact problems is presented. In the model it is assumed that the initial crack of length 0.02 mm is initiated at the surface due to previous mechanical or heat treatment of the material as well as a consequence of running process in an early stage of exploitation. The discretised model with the initial crack is then subjected to the normal contact pressure and tangential loading due to friction between the contacting surfaces. The model also considers the moving contact of the contacting surfaces and fluid trapped in the crack. The crack propagation path is predicted with the MTS and modified MTS criterion, which takes into account the influence of the stress intensity factor K_I and K_II, T-stress, stress on the crack surface caused by lubricant pressure inside the crack and critical distance ahead the crack tip. The numerical results correspond well with available experimental data.     

Stress intensity factor K and the elastic T-stress for corner cracks

The stress intensity factor K and the elastic T-stress for corner cracks have been determined using domain integral and interaction integral techniques. Both quarter-circular and tunnelled corner cracks have been considered. The results show that the stress intensity factor K maintains a minimum value at the mid-plane where the T-stress reaches its maximum, though negative, value in all cases. For quarter-circular corner cracks, the K solution agrees very well with Pickard's (1986) solution. Rapid loss of crack-front constraint near the free surfaces seems to be more evident as the crack grows deeper, although variation of the T-stress at the mid-plane remains small. Both K and T solutions are very sensitive to the crack front shape and crack tunnelling can substantially modify the K and T solutions. Values of the stress intensity factor K are raised along the crack front due to crack tunnelling, particularly for deep cracks. On the other hand, the difference in the T-stress near the free surfaces and at the mid-plane increases significantly with the increase of crack tunnelling. These results seem to be able to explain the well-observed experimental phenomena, such as the discrepancies of fatigue crack growth rate between CN (corner notch) and CT (compact tension) test pieces, and crack tunnelling in CN specimens under predominantly sustained load.     

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.     

Cracked asphalt pavement under traffic loading – A 3D finite element analysis

An asphalt pavement containing a transverse top-down crack is investigated under traffic loading using 3D finite element analysis. The stress intensity factors (SIFs) and the T-stress are calculated for different distances between the crack and the vehicle wheels. It is found that all the three Modes (I, II and III) are present in the crack deformation. The signs and magnitudes of K_I, K_II, K_III and T are significantly dependent on the location of the vehicle wheels with respect to the crack plane. The magnitude of T-stress is considerable, if compared to the stress intensity factors, when one of the wheels is very close to the crack plane.     

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.     

Quasi-static crack tip fields in rate-sensitive FCC single crystals

In this work, the effects of loading rate, material rate sensitivity and constraint level on quasi-static crack tip fields in a FCC single crystal are studied. Finite element simulations are performed within a mode I, plane strain modified boundary layer framework by prescribing the two term (K − T) elastic crack tip field as remote boundary conditions. The material is assumed to obey a rate-dependent crystal plasticity theory. The orientation of the single crystal is chosen so that the crack surface coincides with the crystallographic (010) plane and the crack front lies along [10[`1]][10\overline 1] direction. Solutions corresponding to different stress intensity rates [(K)\dot]\dot{{K}}, T-stress values and strain rate exponents m are obtained. The results show that the stress levels ahead of the crack tip increase with [(K)\dot]\dot{{K}} which is accompanied by gradual shrinking of the plastic zone size. However, the nature of the shear band patterns around the crack tip is not affected by the loading rate. Further, it is found that while positive T-stress enhances the opening and hydrostatic stress levels ahead of crack tip, they are considerably reduced with imposition of negative T-stress. Also, negative T-stress promotes formation of shear bands in the forward sector ahead of the crack tip and suppresses them behind the tip.     

A fracture mechanics model for the analysis of micro-pitting in regard to lubricated rolling–sliding contact problems

A computational model is presented for the analysis of micro-pitting in regard to lubricated rolling–sliding contact problems. This model assumes the appearance of an initial microcrack on the contact surface due to the mechanical or thermal treatment of the material, and as a consequence of an on-going process in early the stage of exploitation. The discretised model of the contacting mechanical elements is subjected to normal loading (Hertzian contact pressure), tangential loading (friction between contacting surfaces) and internal pressure to the crack surfaces. Crack propagation is predicted as follows: (1) using modified maximum tangential stress criterion, which takes into account the influence of stress intensity factors K_I and K_II, T-stress, stress on the crack’s surface caused by lubricant pressure inside the crack, and the critical distance ahead of the crack tip and (2) the classical maximum tangential stress criterion, which only takes into account the influence of the stress intensity factors K_I and K_II. The stress intensity factor based on these two criteria is then used in a short crack growth theory to determine the fatigue life of an initial crack to extent up to micro-pit. The developed model is applied to a real spur gear pair.     

A thermomechanical fracture modeling and simulation for functionally graded solids using a residual-strain formulation

This paper addresses mixed-mode crack growth in two-dimensional functionally graded solids under thermomechanical loads, and investigates the effect of mechanical and thermal loads as well as the T-stress on their crack growth behavior. A novel residual strain-based formulation in the interaction integral method is developed and used for the accurate evaluation of mixed-mode stress intensity factors and/or the T-stress. Simulation of mixed-mode crack propagation in functionally graded materials including solid oxide fuel cells under thermomechanical loads is performed by means of the finite element method and the generalized interaction integrals in conjunction with a remeshing algorithm. An iterative procedure is used for crack growth simulation including the calculation of mixed-mode stress intensity factors and/or the T-stress by means of the generalized interaction integral method, determination of crack growth direction and crack initiation condition based on selected fracture criteria, and local automatic remeshing along the crack path. The present approach employs a user-defined crack increment at the beginning of the simulation. Crack trajectories and fracture parameters obtained by the present simulation for thermomechanical loads are assessed for some numerical examples in comparison with those for mechanical loads.     

Numerical simulation of constraint effects in fatigue crack growth

The computational analysis of constraint effects on fatigue crack growth is discussed. An irreversible cohesive zone model is used in the computations to describe the processes of material separation under cyclic loading. This approach is promising for the investigation of fatigue crack growth under constraint as the energy dissipation due to the formation of new crack surface and cyclic plastic deformation is accounted for independently. Fatigue crack growth in multi-layer structures under consideration of different levels of T-stress are conducted with a modified boundary layer model. Fatigue crack growth is computed as a function of layer thickness and T-stress for constant and variable amplitude loading cases.     

On fracture of kinked interface cracks – The role of T-stress

This paper presents a modified maximum tangential stress criterion (MMTS) for prediction of the fracture initiation conditions in kinked bi-material cracks. The criterion takes into account the effect of T-stress as well as the stress intensity factors (K_I and K_II) to predict the mixed mode fracture toughness of interface cracked specimens. First the fracture criterion is developed and the effect of sign and magnitude of T-stress on mixed mode fracture toughness is studied analytically. Then, the suggested criterion is evaluated using the experimental data reported for some epoxy/Aluminum Brazil-nut-sandwich specimens in the literature. The MMTS criterion is also compared with the conventional maximum tangential stress (MTS) criterion and hence, significantly improved estimates were achieved for mixed mode fracture toughness of the tested specimens.     

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 in-plane and out-of-plane stress constraint factors and K-T-Tz description of stress fields near the border of a quarter-elliptical corner crack

The elastic T-stress and stress intensity factor K for quarter-elliptical corner cracks have been investigated in elastic plates by detailed three-dimensional finite-element calculations. The distributions of normalized K and T-stress have been obtained along the crack front with aspect ratios (a/c) of 0.2, 0.3, 0.4, 0.5, 0.6, 0.8 and 1.0, and far-field tension and the effect of Poisson's ratio have also been considered. The normalized K increases and the normalized T-stress decreases with the increase of Poisson's ratio v. For v= 0.3, the normalized K gradually increases in the range of crack-face angle φ≥ 22.5° and decreases in the range of φ≤ 22.5° with the increase of a/c. The normalized T-stress increases in the beginning and then decreases with increasing φ except for a/c= 0.2 and a/c= 0.3. By fitting the numerical results with the least squares method, empirical formulae have been given for the convenience of engineering applications. Combining with the corresponding out-of-plane constraint factor T_z, the three-parameter K-T-T_z approach has been provided, which can accurately describe the stress field around the crack front.     

Numerical investigation of crack tip strain localization under cyclic loading in FCC single crystals

In this work, the crack tip strain localization in a face centered cubic single crystal subject to both monotonic and cyclic loading was investigated. The effect of constraint was implemented using T-stress and strain accumulation was studied for both isotropic and anisotropic elastic cases with the appropriate application of remote displacement fields in plane strain. Modified boundary layer simulations were performed using the crystal plasticity finite element framework. The consideration of elastic anisotropy amplified the effect of constraint level on stress and plastic strain fields near the crack tip indicating the importance of its use in fracture simulations. In addition, to understand the cyclic stress and strain behavior in the vicinity of the crack tip, combined isotropic and kinematic hardening laws were incorporated, and their effect on the evolution of yield curves and plastic strain accumulation were investigated. With zero-tension cyclic load, the evolution of plastic strain and Kirchhoff stress components showed differences in magnitudes between isotropic and anisotropic elastic cases. Furthermore, under cyclic loading, ratcheting was observed along the localized slip bands, which was shown to be affected by T-stress as well as elastic anisotropy. Negative T-stress increased the accumulation of plastic strain with number of cycles, which was further amplified in the case of elastic anisotropy. Finally, in all the cyclic loading simulations, the plastic strain accumulation was higher near the \(55^0 \) slip band.     

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.     

Analysis of the effects of compressive stresses on fatigue crack propagation rate

In this study, the effects of compressive stresses on the crack tip parameters and its implication on fatigue crack growth have been studied. Elastic–plastic finite element analysis has been used to analyse the change of crack tip parameters with the increase of the applied compressive stress level.The near crack tip opening displacements and the reverse plastic zone size around the crack tip have been obtained. The finite element analysis shows that when unloading from peak tensile applied stress to zero applied stress, the crack tip is still kept open and the crack tip opening displacement gradually decreases further with the applied compressive stress. It has been found that for a tension–compression stress cycle these crack tip parameters are determined mainly by two loading parameters, the maximum stress intensity K_max in the tension part of the stress cycle and the maximum compressive stress σ_maxcom in the compression part of the stress cycle.Based on the two parameters, K_max, and σ_maxcom, a fatigue crack propagation model for negative R ratios only has been developed to include the compressive stress effect on the fatigue crack propagation rate.Experimental fatigue crack propagation data sets were used for the verification of this model, good agreements have been obtained.     

Variation of elastic T-stresses along slightly wavy 3D crack fronts

The non-singular terms in the series expansion of the elastic crack-tip stress field, commonly referred to as the elastic T-stresses, play an important role in fracture mechanics in areas such as the stability of a crack path and the two-parameter characterization of elastic-plastic crack-tip deformation. In this paper, a first order perturbation analysis is performed to study some basic properties of the T-stress variation along a slightly wavy 3D crack front. The analysis employs important properties of Bueckner-Rice 3D weight function fields. Using the Boussinesq-Papkovitch potential representation for the mode I weight function field, it is shown that, for coplanar cracks in an infinite isotropic and homogeneous linear elastic body, the mean T-stress along an arbitrary crack front is independent of the shape and size of the crack. Further, a universal relation is discovered between the mean T-stress and the stress field at the same crack front location under the same loading but in the absence of a crack. First-order-accurate solutions are given for the T-stress variation along a slightly wavy crack front with nearly circular or straight confifurations. Specifically, cosine wave functions are adopted to describe smooth polygonal and slightly undulating planar crack shapes. The results indicate that T ₁₁, the 2D T-stress component acting normal to the crack front, increases with the curvature of the crack front as it bows out but T ₃₃, acting parallel to the crack front, decreases with the crack front curvature.     

Application of a conservation integral to an interface crack interacting with singularities

A mutual integral, which has the conservation property is applied to the problem of an interfacial crack. The stress intensity factors K ₁, K ₂, K ₃ and T-stress for the problem in an infinite medium are easily obtained by using the mutual integral without solving the boundary value problem. In doing so the auxiliary solutions are required and they are taken from the known asymptotic solutions. This method is amenable to numerical evaluation of the stress intensity factors and T-stress if the interfacial crack in a finite medium is considered.