Fracture assessment of U-notches under mixed mode loading: two procedures based on the ‘equivalent local mode I’ concept

Two fracture criteria are proposed and applied to blunt-notched components made of brittle materials loaded under mixed mode; the former is based on the averaged strain energy density over a given control volume, the latter on the cohesive crack zone model. In both instances use of the equivalent local mode I hypothesis is made. Only two material properties are needed: the ultimate tensile strength and the fracture toughness. Numerical predictions of rupture loads from the two criteria are compared with experimental measurements from more than 160 static tests with notched beams. The samples are made of PMMA and tested at − 60°C to assure a bulk behaviour almost linear elastic up to rupture. Notch root radii range from 0.2 to 4.0 mm and load mixicity varies from pure mode I to a prevailing mode II. The good agreement between theory and experimental results adds further confidence to the proposed fracture criteria.     

A criterion for brittle fracture in U-notched components under mixed mode loading

A failure criterion is proposed for brittle fracture in U-notched components under mixed-mode static loading. The criterion, called UMTS, is developed based on the maximum tangential stress criterion and also a criterion proposed in the past for mode I failure of rounded V-shaped notches [Gomez FJ, Elices M. A fracture criterion for blunted V-notched samples. Int J Fracture 2004;127:239-64]. Using the UMTS criterion, a set of fracture curves are derived in terms of the notch stress intensity factors. These curves can be used to predict the mixed mode fracture toughness and the crack initiation angle at the notch tip. An expression is also obtained from this criterion for predicting fracture toughness of U-notched components in pure mode II loading. It is shown that there is a good agreement between the results of UMTS criterion and the experimental data obtained by other authors from three-point bend specimens.     

Strain energy-based evaluations of plastic notch stress intensity factors at pointed V-notches under tension

An analytical study is carried out on the existing link between elastic and plastic notch stress intensity factors at pointed V-notches in plates under tension.The frame is developed on the basis of the elastic and plastic energy concentration factors of the notch defined here as the ratio between the local and the nominal strain energy densities. The link varies under plane stress and plane strain conditions. The local strain energy density is evaluated over a control volume drawn by the energy contour lines ahead of the notch and allows plastic notch stress intensity factors to be predicted on the basis of an ideally linear elastic analysis, both under small and large scale yielding.     

Fracture of U-notched specimens under mixed mode: Experimental results and numerical predictions

The first part of the paper gives an account of 153 fracture tests on blunted notched specimens (with notches of root radius ranging from 0.3 to 4.0 mm), loaded under mixed mode (ranging from almost pure mode I to mode II, and beyond). Maximum loads and initial crack angles were measured as a function of notch root radius and loading mixity. Such results can help in evaluating numerical models of the fracture of notched components. The second part of the paper deals with the suitability of the cohesive crack concept for predicting fracture loads under mixed mode. Use of local mode I was considered for numerical computations. Comparison of experimental results with numerical predictions was significantly accurate. Diagrams of fracture loci for notched components loaded under mixed mode are discussed.     

Acoustic emission imaging of shear failure in large reinforced concrete structures

The averaged value of the strain energy density over a well-defined volume is used to predict the static strength of U-notched specimens under mixed-mode conditions due to combined bending and shear loads. The volume is centered in relation to the maximum principal stress present on the notch edge, by rigidly rotating the crescent-shaped volume already used in the literature to analyse U- and V-shaped notches subject to mode I loading. The volume size depends on the ultimate tensile strength σ _u and the fracture toughness K _IC of the material. In parallel, an experimental programme was performed. All specimens are made of polymethyl-metacrylate (PMMA), a material which exhibits quasi-brittle behaviour at -60°C. Good agreement is found between experimental data for the critical loads to failure and theoretical predictions based on the constancy of the mean strain energy density over the control volume.     

Plastic notch stress intensity factors for pointed V-notches under antiplane shear loading

The paper deals with a work-hardening, elastic–plastic, stress analysis of pointed V-notches under antiplane shear deformation loading both under small and large scale yielding. Stress and strain field intensities are expressed in terms of plastic Notch Stress Intensity Factors, which are analytically linked to the corresponding linear elastic ones under small scale yielding. The near tip stress and strain fields are then used to give closed-form expressions for the Strain Energy Density over a circular sector surrounding the notch tip, and for the J-integral parameter, both as a function of the relevant plastic NSIFs, these expressions being valid both under small and large scale yielding.     

Influence of effective stress intensity factor range on mixed mode fatigue crack propagation

ABSTRACT The behaviour of fatigue crack propagation of rectangular spheroidal graphite cast iron plates, each consisting of an inclined semi‐elliptical crack, subjected to axial loading was investigated both experimentally and theoretically. The inclined angle of the crack with respect to the axis of loading varied between 0° and 90°. In the present investigation, the growth of the fatigue crack was monitored using the AC potential drop technique, and a series of modification factors, which allow accurate sizing of such defects, is recommended. The rate of fatigue crack propagation db/dN is postulated to be a function of the effective strain energy density factor range, ΔS_eff. Subsequently, this concept is applied to predict crack growth due to fatigue loads. The mixed mode crack growth criterion is discussed by comparing the experimental results with those obtained using the maximum stress and minimum strain energy density criteria. The threshold condition for nongrowth of the initial crack is established based on the experimental data.     

Elastic notch stress intensity factors for sharply V-notched rounded bars under torsion

The paper deals with the determination of analytical expressions for the mode III notch stress intensity factors for circumferentially-sharply-notched rounded bars under torsion loading, starting from the theoretical stress concentration factors of the corresponding notch problem.An exact, closed-form solution for the NSIFs is obtained for deep notches; subsequently the solution is extended also to finite notched components taking advantage of a shape function determined by a numerical best fitting procedure.     

Review of the fatigue strength of welded joints based on the notch stress intensity factor and SED approaches

Several approaches exist for the fatigue strength assessment of welded joints. In addition to the traditional nominal stress approach, various approaches were developed using a local stress as fatigue parameter. In recent times, the N-SIF based approaches using the notch stress intensity at the weld toe or root have been developed. Based on this, the more practical strain energy density (SED) and the Peak Stress approaches were proposed. This paper reviews the proposed design S–N curves of the N-SIF and SED approaches questioning in particular the consideration of misalignment effects, which should be included on the load side of local approaches in order to consider them individually in different types of welded joints. A re-analysis of fatigue tests evaluated for the effective notch stress approach leads to slight changes of the design S–N curves and of the radius of the control volume used for averaging the SED at the notches. Further, on purpose fatigue tests of artificially notched specimens show that the fatigue assessment using a single-point fatigue parameter might be problematic because the crack propagation phase, being part of the fatigue life, is strongly affected by the stress distribution along the crack path that may vary considerably between different geometries and loading cases.     

Rapid evaluation of notch stress intensity factors using the peak stress method: Comparison of commercial finite element codes for a range of mesh patterns

The peak stress method (PSM) is an engineering, finite element (FE)‐oriented method to rapidly estimate the notch stress intensity factors by using the singular linear elastic peak stresses calculated from coarse FE analyses. The average element size adopted to generate the mesh pattern can be chosen arbitrarily within a given range. Originally, the PSM has been calibrated under pure mode I and pure mode II loadings by means of Ansys FE software. In the present contribution, a round robin between 10 Italian universities has been carried out to calibrate the PSM with 7 different commercial FE codes. To this aim, several two‐dimensional mode I and mode II problems have been analysed independently by the participants. The obtained results have been used to calibrate the PSM for given stress analysis conditions in (i) FE software, (ii) element type and element formulation, (iii) mesh pattern, and (iv) criteria for stress extrapolation and principal stress analysis at FE nodes.     

Multiaxial fatigue of V‐notched steel specimens: a non‐conventional application of the local energy method

The paper deals with the multi‐axial fatigue strength of notched specimens made of 39NiCrMo3 hardened and tempered steel. Circumferentially V‐notched specimens were subjected to combined tension and torsion loading, both in‐phase and out‐of‐phase, under two nominal load ratios, R=?1 and R= 0, also taking into account the influence of the biaxiality ratio, λ=τ_a/σ_a. The notch geometry of all axi‐symmetric specimens was a notch tip radius of 0.1 mm, a notch depth of 4 mm, an included V‐notch angle of 90° and a net section diameter of 12 mm. The results from multi‐axial tests are discussed together with those obtained under pure tension and pure torsion loading on plain and notched specimens. Furthermore the fracture surfaces are examined and the size of non‐propagating cracks measured from some run‐out specimens at 5 million cycles. Finally, all results are presented in terms of the local strain energy density averaged in a given control volume close to the V‐notch tip. The control volume is found to be dependent on the loading mode.     

Optimum location of a three strain gauge rosette for measuring mixed mode stress intensity factors

This paper analyzes the errors inherent to the determination of mixed mode stress intensity factors from data obtained by using a three strain gauge rosette. The analysis shows that the errors are mainly due the third characteristic value (3/2) and its corresponding coefficients. It is also shown that the errors do not depend on the orientation angle of the rosette, the angle between the strain gauges and the material properties. The error mainly depends on its location (radius, angle), being linear in the radius. For pure mode I, an angle of 90° will completely eliminate the error due to the angle, while for pure mode II, a 0° angle will minimize it. The normalized variation of the errors with the angle at any radius is shown for different ratios of the corresponding coefficients of the third characteristic value. The analytical results are applied to a numerical example of an edge crack subjected to mixed mode loading. From the numerical example, it is recommended to use two strain gauge rosettes at the same angle, and linearly extrapolate their results, if errors less than 15% for a mixed mode field are desired.     

Weight functions and stress intensity factors for pin‐loaded single‐edge notch bend specimen

A recently developed pin‐loaded single‐edge notch bend specimen provides an alternative to the single‐edge notch tension specimen commonly used for small‐crack growth testing. In this paper, weight functions for pin‐loaded single‐edge notch bend specimen are derived by using two methods, the classical analytical weight function method and the newly developed numerical weight function complex variable Taylor series expansion method. Excellent agreement between the two methods is achieved. Based on these weight functions, accurate stress intensity factors for two load cases, that is, pin‐loading and Dugdale loading, which is required for plasticity‐induced crack‐closure analysis based on the strip‐yield model, are determined.     

A phantom node formulation with mixed mode cohesive law for splitting in laminates

A phantom node method with mixed mode cohesive law is proposed for the simulation of splitting in laminates. With this method, a discontinuity in the displacement field can be modeled at arbitrary locations. The micromechanical phenomenon that splitting cracks grow parallel to the fiber, is incorporated on the mesolevel, i.e., in the homogenized ply, by setting the direction of the crack propagation equal to the fiber direction. A new mixed mode cohesive law is introduced for increased robustness of the incremental-iterative solution procedure. The model is validated with mixed mode bending tests, and its utility is illustrated with examples for a single ply and for a laminate.     

Analytical study of the elastic–plastic stress fields ahead of parabolic notches under antiplane shear loading

An analytical study is carried out on the elastic–plastic stress and strain distributions and on the shape of the plastic zone ahead of parabolic notches under antiplane shear loading and small scale yielding. The material is thought of as obeying an elastic-perfectly-plastic or a strain hardening law. When the notch root radius becomes zero, the analytical frame matches the solutions for the crack case due to Hult–McClintock (elastic-perfectly-plastic material) and Rice (strain hardening material). The analytical frame provides an explicit link between the plastic stress and the elastic stress at the notch tip. Neuber’solution for blunt notches under antiplane shear is also obtained and the conditions under which such a solution is valid are discussed in detail by using elastic and plastic notch stress intensity factors. Finally, revisiting Glinka and Molski’s equivalent strain energy density (ESED), these factors are used also to give, under antiplane shear loading, the increment of the strain energy at the notch tip with respect to the linear elastic case.     

Finite element evaluation of mixed mode stress intensity factors in functionally graded materials

This paper is directed towards finite element computation of fracture parameters in functionally graded material (FGM) assemblages of arbitrary geometry with stationary cracks. Graded finite elements are developed where the elastic moduli are smooth functions of spatial co‐ordinates which are integrated into the element stiffness matrix. In particular, stress intensity factors for mode I and mixed‐mode two‐dimensional problems are evaluated and compared through three different approaches tailored for FGMs: path‐independent J^*_k‐integral, modified crack‐closure integral method, and displacement correlation technique. The accuracy of these methods is discussed based on comparison with available theoretical, experimental or numerical solutions. Copyright © 2001 John Wiley & Sons, Ltd.     

Variations of stress intensity factor of a semi-elliptical surface crack subjected to mixed mode loading

Maximum stress intensity factors of a surface crack usually appear at the deepest point of the crack, or a certain point along crack front near the free surface depending on the aspect ratio of the crack. However, generally it has been difficult to obtain smooth distributions of stress intensity factors along the crack front accurately due to the effect of corner point singularity. It is known that the stress singularity at a corner point where the front of 3 D cracks intersect free surface is depend on Poisson's ratio and different from the one of ordinary crack. In this paper, a singular integral equation method is applied to calculate the stress intensity factor along crack front of a 3-D semi-elliptical surface crack in a semi-infinite body under mixed mode loading. The body force method is used to formulate the problem as a system of singular integral equations with singularities of the form r ⁻³ using the stress field induced by a force doublet in a semi-infinite body as fundamental solution. In the numerical calculation, unknown body force densities are approximated by using fundamental density functions and polynomials. The results show that the present method yields smooth variations of mixed modes stress intensity factors along the crack front accurately. Distributions of stress intensity factors are indicated in tables and figures with varying the elliptical shape and Poisson's ratio.     

A precise computation of stress intensity factor on the front of a convex planar crack

In this paper, we consider a general integral expression for mode I stress intensity factor along the fronts of convex planar cracks. For this integral approximation, we develop a simple numerical quadrature formula on every convex set Ω which allows a precise estimation of the error. This permits the use of extrapolation techniques for the accurate computation of the integral. Copyright © 2002 John Wiley & Sons, Ltd.     

Wide range data for crack tip parameters in two disc-type specimens under mixed mode loading

Disc-type specimens are among favorite test samples for determining mode I and mixed mode fracture toughness in brittle materials like rocks, brittle polymers, ceramics, etc. In this research, the finite element method is used to analyze two disc-type specimens: a semi-circular disc specimen containing an edge crack and subjected to three-point-bend loading (SCB specimen), and a centrally cracked circular disc subjected to diametral compressive loading, often called the Brazilian disc specimen. The crack parameters K_I, K_II and T are calculated for different mode mixities from pure mode I to pure mode II. Although the stress intensity factors K_I and K_II are presented mainly for validation of the analyses, they are also used for determining the crack angle corresponding to pure mode II for each specimen. It is shown that in general the T-stress increases for larger crack angles. While the T-stress in the Brazilian disc specimen is always negative for any combinations of mode I and mode II, the sign of T-stress in the SCB specimen depends on the mode mixity. A very good agreement is shown to exist between the calculated results for T and those very limited data presented in other papers.     

A non-local stress and strain energy release rate mixed mode fracture initiation and propagation criteria

A non-local stress condition for crack initiation and propagation in brittle materials is presented. This condition is expressed in terms of normal and tangential traction components acting on a physical plane segment (damage zone) of specified length. Next, a non-local strain energy release rate criterion is proposed. This condition is based on the assumption that initiation or propagation of cracking occurs when the maximal value of the function of opening and sliding energy release rates reaches a critical value. The value of energy release rates is determined for finite elementary crack growth. Mixed mode conditions are considered, for which both the critical load value and the crack orientation are predicted from the non-local stress and energy criteria, which are applicable to both regular and singular stress concentrations. The effect of non-singular second order term (T_σ-stress) on the crack propagation is discussed.