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.     

Stress intensity factors for surface cracks in single‐edge notch bend specimen by a three‐dimensional weight function method

Stress intensity factors for half‐elliptical surface cracks at a semi‐circular notch in a recently developed single‐edge notch bend specimen are determined for a wide range of geometrical parameters using a three‐dimensional weight function method. Two load cases of pin loading and uniform remote tension are considered. The results are in good agreement with abaqus/franc3d finite element analysis. It is found that the Ziegler–Newman engineering similitude approach (programmed into the Fatigue Crack Growth Structural Analysis life‐prediction code) produces good results for a wide range in a/c ratios. Expressions by multi‐variable curve fitting to the weight function results are presented for easy engineering applications.     

A weight function method for mixed modes hole‐edge cracks

A weight function approach is proposed to calculate the stress intensity factor and crack opening displacement for cracks emanating from a circular hole in an infinite sheet subjected to mixed modes load. The weight function for a pure mode II hole‐edge crack is given in this paper. The stress intensity factors for a mixed modes hole‐edge crack are obtained by using the present mode II weight function and existing mode I Green (weight) function for a hole‐edge crack. Without complex derivation, the weight functions for a single hole‐edge crack and a centre crack in infinite sheets are used to study 2 unequal‐length hole‐edge cracks. The stress intensity factor and crack opening displacement obtained from the present weight function method are compared well with available results from literature and finite element analysis. Compared with the alternative methods, the present weight function approach is simple, accurate, efficient, and versatile in calculating the stress intensity factor and crack opening displacement.     

Estimation of three‐dimensional stress intensity factor for structural ‘T’ details

Three‐dimensional (3D) opening mode stress intensity factors (SIFs) for structural steel‐welded ‘T’ details were investigated by the finite element method. A 3D shape‐dependent correction factor is proposed for semi‐elliptical surface cracks. The aspect ratio (a/c) of a semi‐elliptical crack plays a key role in the approximation of 3D‐SIF values, and in the present study, it was estimated for a 3D crack analysis. The estimated 3D‐SIF was determined through a correlation between the a/c ratio and the two‐dimensional SIF for semi‐elliptical cracks in the thickness direction adjacent to the web‐flange junction of a welded ‘T’. The resulting equation can be used to estimate the 3D‐SIF values from the two‐dimensional SIF without much ambiguity.     

Automatic 3-D crack propagation calculations: a pure hexahedral element approach versus a combined element approach

This article presents an evaluation of two different crack prediction approaches based on a comparison of the stress intensity factor distribution for three example problems. A single edge notch specimen and a quarter circular corner crack specimen subjected to shear displacements and a three point bend specimen with a crack inclined to the mid-plane are examined. The stress intensity factors are determined from the singular stress field close to the crack front. Two different fracture criteria are adopted for the calculation of an equivalent stress intensity factor and crack deflection angle. The stress intensity factor distributions for both numerical methods agree well to available reference solutions. Deviations are recorded at crack front locations near the free surface probably due to global contraction effects and the twisting behaviour of the crack front. Crack propagation calculations for the three point bending specimen give results that satisfy intuitive expectations. The outcome of the study encourages further pursuit of a crack propagation tool based on a combination of elements.     

Experimental investigation on low cycle fatigue and fracture behaviour of a notched Ni‐based superalloy at elevated temperature

In order to investigate the effects of stress concentration on low cycle fatigue properties and fracture behaviour of a nickel‐based powder metallurgy superalloy, FGH97, at elevated temperature, the low cycle fatigue tests have been conducted with semi‐circular and semi‐elliptical single‐edge notched plate specimens at 550 and 700 °C. The results show that the fatigue life of the notched specimen decreases with the increase of stress concentration factor and the fatigue crack initiation life evidently decreases because of the defect located in the stress concentration zone. Moreover, the plastic deformation induced by notch stress concentration affects the initial crack occurrence zone. The angle α of the crack occurrence zone is within ±10° of notch bisector for semi‐circular notched specimens and ±20° for semi‐elliptical notched specimens. The crack propagation rate decreases to a minimum at a certain length, D, and then increases with the growth of the crack. The crack propagation rate of the semi‐elliptical notched specimen decelerates at a faster rate than that of the semi‐circular notched specimen because of the increase of the notch plasticity gradient. The crack length, D, is affected by both the applied load and the notch plasticity gradient. In addition, the fracture mechanism is shown to transition from transgranular to intergranular as temperature increases from 550 to 700 °C, which would accelerate crack propagation and reduce the fatigue life.     

Stress intensity factors for large aspect ratio surface and corner cracks at a semi-circular notch in a tension specimen

Stress intensity factor solutions for semi-elliptic surface and quarter-elliptic corner cracks emanating from a semi-circular notch in a tension specimen are presented. A threedimensional finite-element analysis in conjunction with the equivalent domain integral was used to calculate stress intensity factors (SIF). SIF solutions for surface or corner crack (crack length to depth ratio of 2) at a notch are presented for a wide range of crack sizes and notch radii. Results showed that the SIF are larger for larger crack lengths and for larger notch radii. The SIF are nearly constant all along the crack front for deep surface cracks and for all corner cracks analysed.     

Approximate weight function method for elliptical arc through cracks in mechanical joints

Mechanical joints such as bolted, riveted or pinned joints are widely used to join the constituent parts of structural components. Reliable stress intensity factor analysis of arbitrary cracks in mechanical joints is required for the safety evaluation or fracture mechanics design. It has been reported that cracks in mechanical joints usually nucleate as the corner crack and grow as the elliptical arc through crack. The weight function method is a useful technique to calculate the stress intensity factor using the appropriate weight function for a cracked body and the stress field of an uncracked body. In this paper, the weight function method for the two surface points of elliptical arc through cracks in mechanical joints is developed to analyze the mixed-mode stress intensity factors. Unknown coefficients included in the weight function are determined using the reference stress intensity factors obtained from finite element analysis.     

STRESS INTENSITY FACTORS FOR CORNER CRACKS AT A SEMI-CIRCULAR NOTCH UNDER STRESS GRADIENTS

Abstract— A weight function method, recently developed by the authors, is applied to calculate stress intensity factors for corner cracks emanating from a semi-circular notch under crack face polynomial pressure loading. A wide range of configuration parameters are considered. These results, combined with superposition principle, allow determination of stress intensity factors under general loading conditions. The approach is demonstrated by obtaining stress intensity factors for the load cases of remote tension and shot-peening residual stresses at the notch.     

Vorhersagemodell für die Wachstumsraten kurzer Risse auf der Basis von Kmax‐konstant‐Versuchen an M(T)‐Proben

Prediction model for the growth rates of short cracks based on K_max‐constant tests with M(T) specimens The fatigue crack growth behaviour of short corner cracks in the Aluminium alloys Al 6013‐T6 and Al 2524‐T351 was investigated. The aim was to determine the crack growth rates of small corner cracks at stress ratios of R = 0.1, R = 0.7 and R = 0.8 and to develop a method to predict these crack growth rates from fatigue crack growth curves determined for long cracks. Corner cracks were introduced into short crack specimens, similar to M(T)‐specimens, at one side of a hole (Ø = 4.8 mm) by cyclic compression (R = 20). The pre‐cracks were smaller than 100 μm (notch + precrack). A completely new method was used to cut very small notches (10–50 μm) into the specimens with a Focussed Ion Beam. The results of the fatigue crack growth tests with short corner cracks were compared with long fatigue crack growth test data. The short cracks grew at ΔK‐values below the threshold for long cracks at the same stress ratio. They also grew faster than long cracks at the same ΔK‐values and the same stress ratios. A model was developed on the basis of K_max‐constant tests with long cracks that gives a good and conservative prediction of the short crack growth rates.     

Weight functions for the determination of stress intensity factor and T‐stress for semi‐elliptical cracks in finite thickness plate

This paper presents the application of weight function method for the calculation of stress intensity factors (K) and T‐stress for surface semi‐elliptical crack in finite thickness plates subjected to arbitrary two‐dimensional stress fields. New general mathematical forms of point load weight functions for K and T have been formulated by taking advantage of the knowledge of a few specific weight functions for two‐dimensional planar cracks available in the literature and certain properties of weight function in general. The existence of the generalised forms of the weight functions simplifies the determination of specific weight functions for specific crack configurations. The determination of a specific weight function is reduced to the determination of the parameters of the generalised weight function expression. These unknown parameters can be determined from reference stress intensity factor and T‐stress solutions. This method is used to derive the weight functions for both K and T for semi‐elliptical surface cracks in finite thickness plates, covering a wide range of crack aspect ratio (a/c) and relative depth (a/t) at any point along the crack front. The derived weight functions are then validated against stress intensity factor and T‐stress solutions for several linear and nonlinear two‐dimensional stress distributions. These derived weight functions are particularly useful for the development of two‐parameter fracture and fatigue models for surface cracks subjected to fluctuating nonlinear stress fields, such as these resulting from surface treatment (shot peening), stress concentration or welding (residual stress).     

Modelling the residual strength of fatigue damage at a single semicircular edge notch: Semielliptical crack and through‐the‐thickness crack

In the present paper, computational framework for fatigue performance analysis of a semicircular edge notch with a through‐the‐thickness crack or a semielliptical crack is discussed. The failure behaviour of such configurations is theoretically examined through the stress‐intensity analysis and residual life estimation. The stress field of a damaged notch configuration is herein investigated by employing analytical and numerical approaches. Further, a fracture mechanics–based methodology, developed for fatigue life assessment, is taking into account the crack growth model proposed by Huang and Moan in which the stress ratio is involved. The efficiency of the obtained fatigue damage assessments, related to the edge notch configurations, is verified through appropriate experimental observations.     

Stress intensity factor analysis of elliptical corner cracks in mechanical joints by weight function method

Mixed-mode stress intensity factors at the surface and deepest points of quarter elliptical corner cracks in mechanical joints such as bolted or riveted joints are analyzed by weight function method. The weight function method is an efficient technique to calculate the stress intensity factors using uncracked stress field. The extended form of the weight function method for 2D mixed-mode problems to 3D mixed-mode is presented and the accuracy due to the number of terms included in the weight function is examined. The effects of the amount of clearance between the hole and the bolt or rivet on the stress intensity factors are investigated, and the critical angle causing the mode I stress intensity factor to be maximized is determined by analyzing the variation of the stress intensity factors along incline angle of crack.     

Rectangular Tensile Sheet with Double Edge Notch Cracks

This paper deals with the rectangular tensile sheet with symmetric double edge notch cracks. Such a crack problem is called an edge notch crack problem for short. By using a hybrid displacement discontinuity method (a boundary element method), two edge notch models are analyzed in detail. By changing the geometrical forms and parameters of the edge notch, and by comparing the stress intensity factors (SIFs) of the edge notch crack problem with those of the double edge cracked plate tension specimen (DECT), which is a model frequently used in fracture mechanics, the effect of the geometrical forms and parameters of the edge notch on the SIFs of the DECT specimen, is revealed. Some geometric characterestic parameters are introduced here, which are used to formulate the notch length and the branch crack length, which are to be determined in mechanical machining of the DECT specimen So we can say that the geometric characterestic parameters and the formulae used to determine the notch length and the branch crack length presented in this paper perhaps have some guidance role for mechanical machining of the DECT specimen.     

A crack propagation criterion based on ΔCTOD measured with 2D‐digital image correlation technique

The fatigue cracks growth rate of a forged HSLA steel (AISI 4130) was investigated using thin single edge notch tensile specimen to simulate the crack development on a diesel train crankshafts. The effect of load ratio, R, was investigated at room temperature. Fatigue fracture surfaces were examined by scanning electron microscopy. An approach based on the crack tip opening displacement range (ΔCTOD) was proposed as fatigue crack propagation criterion. ΔCTOD measurements were carried out using 2D‐digital image correlation techniques. J‐integral values were estimated using ΔCTOD. Under test conditions investigated, it was found that the use of ΔCTOD as a fatigue crack growth driving force parameter is relevant and could describe the crack propagation behaviour, under different load ratio R.     

Weight functions and stress intensity factors for quarter-elliptical corner cracks in fastener holes

Approximate weight functions for a quarter‐elliptical crack in a fastener hole were derived from a general weight function form and two reference stress intensity factors. Closed‐form expressions were obtained for the coefficients of the weight functions. The derived weight functions were validated against numerical data by comparison of stress intensity factors calculated for several nonlinear stress fields. Good agreements were achieved. These derived weight functions are valid for the geometric range of 0.5 ≤a/c≤ 1.5 and 0 ≤a/t≤ 0.8 and R/t= 0.5; and are given in forms suitable for computer numerical integration. The weight functions appear to be particularly suitable for fatigue crack growth prediction of corner cracks in fastener holes and fracture analysis of such cracks in complex stress fields.     

Calculation of mixed mode stress intensity factors for an elliptical subsurface crack under arbitrary normal loading

Normal loading causes mixed fracture modes in an elliptical subsurface crack because of the nonsymmetrical geometry with respect to the crack face. In this paper, mixed mode weight functions (MMWFs) for elliptical subsurface cracks in an elastic semi‐infinite space under normal loading are derived. Reference mixed mode stress intensity factors (MMSIFs), calculated by finite element analysis, under uniform normal loading are used to derive MMWFs. The cracks have aspect ratios and crack depth to crack length ratios of 0.2–1.0 and 0.05 to infinity, respectively. MMWFs are used to calculate MMSIFs for any point of the crack front under linear and nonlinear two‐dimensional (2D) loadings. So, in order to evaluate the fatigue crack growth phenomenon under complicated 2D stress distributions, MMWFs can be easily used. The comparison between the MMSIFs obtained from the MMWFs and finite element analysis indicates high accuracy.     

Crack growth predictions using a strip‐yield model for variable‐amplitude and spectrum loading

Fatigue‐crack‐growth tests were conducted on compact, C(T), specimens made of D16Cz aluminum alloy. Constant‐amplitude tests were conducted over a range of stress ratios (R = P_min/P_max = 0.1 to 0.75). Comparisons were made between test data from middle‐crack tension, M(T), specimens from the literature and C(T) specimens. A crack‐closure analysis was used to collapse the rate data from both specimen types into a fairly narrow band over many orders of magnitude in rates using proper constraint factors. Constraint factors were established from single‐spike overload and constant‐amplitude tests. The life‐prediction code, FASTRAN, which is based on the strip‐yield‐model concept, was used to calculate the crack‐length‐against‐cycles under constant‐amplitude (CA) loading and the single‐spike overload (OL) tests; and to predict crack growth under variable‐amplitude (VA) loading on M(T) specimens and simulated aircraft loading spectrum tests on both specimen types. The calculated crack‐growth lives under CA and the OL tests were generally within ±20 % of the test results, the predicted crack‐growth lives for the VA and Mini‐Falstaff tests on the M(T) specimens were short by 30 to 45 %, while the Mini‐Falstaff+ results on the C(T) specimens were within 10 %. Issues on the crack‐starter notch effects under spectrum loading are discussed, and recommendations are suggested on avoiding these notch effects.     

ANALYSES OF FOUR SYMMETRIC CORNER CRACKS IN ELASTIC PLATES

Abstract— With the aid of the three-dimensional weight function method stress intensity factors and crack opening displacements are determined for the problem of four symmetric quarter-elliptical corner cracks in a rectangular plate or a bar, which is subjected to remote tension. A wide range of configuration parameters is investigated. When the relative crack size is small, the results agree very well with those of a single corner crack determined by the finite element method. The significant difference between four symmetric corner cracks and a single corner crack is also illustrated when the relative crack size is large.     

A new variable‐order singular boundary element for two‐dimensional stress analysis

A new variable‐order singular boundary element for two‐dimensional stress analysis is developed. This element is an extension of the basic three‐node quadratic boundary element with the shape functions enriched with variable‐order singular displacement and traction fields which are obtained from an asymptotic singularity analysis. Both the variable order of the singularity and the polar profile of the singular fields are incorporated into the singular element to enhance its accuracy. The enriched shape functions are also formulated such that the stress intensity factors appear as nodal unknowns at the singular node thereby enabling direct calculation instead of through indirect extrapolation or contour‐integral methods. Numerical examples involving crack, notch and corner problems in homogeneous materials and bimaterial systems show the singular element's great versatility and accuracy in solving a wide range of problems with various orders of singularities. The stress intensity factors which are obtained agree very well with those reported in the literature. Copyright © 2002 John Wiley & Sons, Ltd.