Closed-form structural stress and stress intensity factor solutions for spot welds in commonly used specimens

Closed-form new structural stress and stress intensity factor solutions for spot welds in lap-shear, square-cup, U-shape, cross-tension and coach-peel specimens are obtained based on elasticity theories and fracture mechanics. The loading conditions for spot welds in the central parts of the five types of specimens are first examined. The resultant loads on the weld nugget and the self-balanced resultant loads on the lateral surface of the central parts of the specimens are then decomposed into various types of symmetric and anti-symmetric parts. Closed-form structural stress and stress intensity factor solutions for spot welds under various types of loading conditions are then adopted from the recent work of Lin and Pan to derive new closed-form structural stress and stress intensity factor solutions for spot welds in the five types of specimens. The selection of a geometric factor for square-cup specimens and the decompositions of the loads on the central parts of the U-shape, cross-tension and coach-peel specimens are based on the corresponding three-dimensional finite element analyses of these specimens. The new closed-form solutions are expressed as functions of the spot weld diameter, the sheet thickness, the width and the length of the five types of specimens. The closed-form solutions are also expressed as functions of the angular location along the nugget circumference of spot welds in the five types of specimens in contrast to the limited available solutions at the critical locations in the literature. The new closed-form solutions at the critical locations of spot welds in the five types of specimens are listed or can be easily obtained from the general closed-form solutions for fatigue life predictions.     

Geometric functions of stress intensity factor solutions for spot welds in cross-tension specimens

Stress intensity factor solutions for spot welds in cross-tension specimens are investigated by finite element analyses. Three-dimensional finite element models are developed to obtain accurate solutions. Various ratios of sheet thickness, half specimen width and half effective specimen length to nugget radius are considered. The computational results confirm the functional dependence on the nugget radius and sheet thickness of Zhang’s analytical solutions. The results also provide three geometric functions in terms of normalized half specimen width and normalized half effective specimen length to Zhang’s analytical solutions. Based on the analytical and computational results, the dimensions of cross-tension specimens and the corresponding approximate stress intensity factor solutions are suggested.     

Analytical stress intensity factor solutions for resistance and friction stir spot welds in lap-shear specimens of different materials and thicknesses

In this paper, analytical stress intensity factor and J integral solutions for resistance and friction stir spot welds without and with gap and bend in lap-shear specimens of different materials and thicknesses are developed. The J integral and stress intensity factor solutions for spot welds are first presented in terms of the structural stresses for a strip model. Analytical structural stress solutions for spot welds without and with gap and bend in lap-shear specimens are then developed based on the closed-form structural stress solutions for a rigid inclusion in a finite thin plate subjected to various loading conditions. With the available structural stress solutions, the analytical J integral and stress intensity factor solutions can be obtained as functions of the applied load, the elastic material property parameters, and the geometric parameters of the weld and specimen. The analytical stress intensity factor solutions are selectively validated by the results of three-dimensional finite element analyses for a spot weld with ideal geometry and for a friction stir spot weld with complex geometry, gap and bend. The stress intensity factor and J integral solutions at the critical locations of spot welds in lap-shear specimens of dissimilar magnesium, aluminum and steel sheets with equal and different thicknesses are then presented in the normalized forms as functions of the ratio of the specimen width to the weld diameter. Finally, general trends and simple estimation methods of the stress intensity factor and J integral solutions at the critical locations of spot welds in lap-shear specimens of different materials and thicknesses are given for convenient engineering applications.     

Modeling of failure mode of laser welds in lap-shear specimens of HSLA steel sheets

Failure mode of laser welds in lap-shear specimens of high strength low alloy (HSLA) steel sheets is investigated in this paper. The experiments for laser welds in lap-shear specimens under quasi-static loading conditions are briefly reviewed first. The experimental results showed that the laser welds failed in a ductile necking/shear failure mode and the ductile failure was initiated at a distance away from the crack tip near the boundary of the base metal and heat affected zone. In order to understand the failure mode of these welds, finite element analyses under plane strain conditions were conducted to identify the effects of the different plastic behaviors of the base metal, heat affected zone, and weld zone as well as the weld geometry on the ductile failure. The results of the reference finite element analysis based on the homogenous material model show that the failure mode is most likely to be a middle surface shear failure mode in the weld. The results of the finite element analysis based on the multi-zone non-homogeneous material models show that the higher effective stress–plastic strain curves of the weld and heat affected zones and the geometry of the weld protrusion result in the necking/shear failure mode in the load carrying sheet. The results of another finite element analysis based on the non-homogeneous material model and the Gurson yield function for porous materials indicate that the consideration of void nucleation and growth is necessary to identify the ductile failure initiation site that matches well with the experimental observations. Finally, the results of this investigation indicate that the failure mode of the welds should be examined carefully and the necking/shear failure mode needs to be considered for development of failure or separation criteria for welds under more complex loading conditions.     

Effects of weld geometry and sheet thickness on stress intensity factor solutions for spot and spot friction welds in lap-shear specimens of similar and dissimilar materials

In this paper, three-dimensional finite element analyses for spot welds with ideal geometry in lap-shear specimens of different materials and thicknesses were first conducted. The computational results indicate that the stress intensity factor and J integral solutions based on the finite element analyses agree well with the analytical solutions and that the analytical solutions can be used with a reasonable accuracy. Three-dimensional finite element analyses based on the micrographs of an aluminum 6111 resistance spot weld, an aluminum 5754 spot friction weld, and a dissimilar Al/Fe spot friction weld were also conducted. The computational results indicate that the stress intensity factor and J integral solutions based on the finite element analyses for the aluminum 6111 resistance spot weld and aluminum 5754 spot friction weld with complex geometry are in good agreement with the analytical solutions for the equivalent spot welds with ideal geometry. However, the stress intensity factor and J integral solutions based on the finite element analysis for the Al/Fe spot friction weld with complex geometry are completely different from the analytical solutions for the equivalent spot weld with ideal geometry. Different three-dimensional finite element analyses based on the meshes that represent different features of the complex geometry of the Al/Fe spot friction weld were then conducted. The computational results indicate that the stress intensity factor and J integral solutions for the Al/Fe spot friction weld based on the finite element analysis agree reasonably well with the analytical solutions for the equivalent spot weld with consideration of gap and bend. The computational and analytical results suggest that the stress intensity factor and J integral solutions based on the finite element analysis and the analytical solutions with consideration of gap and bend may be used to correlate with the fatigue crack growth patterns of Al/Fe spot friction welds observed in experiments.     

Modeling of failure near spot welds in lap-shear specimens based on a plane stress rigid inclusion analysis

In this paper, the failure mechanism of resistance spot welds in dual-phase steel lap-shear specimens is investigated based on experimental observations, two-dimensional elasticity theories and two-dimensional finite element analyses. Optical micrographs of the cross sections of spot welds in lap-shear specimens of a dual-phase steel before and after failure are first examined to understand the failure mechanism. The experimental results suggest that under lap-shear loading conditions, a necking failure is initiated near the middle of the nugget circumference in the base metal and then the failure propagates along the nugget circumference in the sheet to final fracture. Based on the stress function approach of the elasticity theory, an analytic solution for an infinite plate containing a rigid inclusion subjected to a resultant shear force is developed and used to investigate the stress and strain distributions near the nugget in lap-shear specimens. The results of the elastic analytic solution and those of a two-dimensional elastic finite element analysis indicate that the initial yielding starts on the two side edges of the inclusion in the sheet. However, the results of a two-dimensional elastic-plastic finite element analysis indicate that as the applied displacement increases, the maximum equivalent plastic strain shifts from the two side edges of the inclusion to the middle of the inclusion along the inclusion circumference in the sheet. The computational results suggest that the location of the initial necking failure should occur near the middle of the nugget circumference in the sheet as observed in experiments based on the forming limit diagram (FLD) for ductile sheet metals.     

A plasticity-corrected stress intensity factor for fatigue crack growth in ductile materials under cyclic compression

For prediction of the fatigue crack growth (FCG) behavior under cyclic compression, a plasticity-corrected stress intensity factor (PC-SIF) range ΔK_pc is proposed on the basis of plastic zone toughening theory. The FCG behaviors in cyclic compression, and the effects of load ratio, preloading and mean load, are well predicted by this new mechanical driving force parameter. Comparisons with experimental data showed that the proposed PC-SIF range ΔK_pc is an effective single mechanical parameter capable of describing the FCG behavior under different cyclic compressive loading conditions.     

Comparison of the simulation and experimental fatigue crack behaviors in the nanoseconds laser shocked aluminum alloy

This investigation was performed to compare the simulation and experimental results of the fatigue crack growth rates and behaviors of the 7050-T7451 aluminum alloy by nanoseconds laser shock processing (LSP). Forman–Newman–deKoning (FNK) model embedded in the Franc2D/L software was utilized to predict fatigue crack growth rate, which was conducted to weigh the stress intensity factor (SIF) changing on the surface cracks. LSP induced high compressive residual stresses that served to enhance fatigue properties by improving the resistance against fatigue crack initiation and propagation. The circulating times of crack growth obtained from the simulation and experimental values indicated a slower fatigue crack growth rates after LSP. The relationships between the elastic–plastic materials crack growth rates and the SIF changing after LSP are resolved.     

Effects of the transverse stress on biaxial fatigue crack growth predicted by plasticity-corrected stress intensity factor

The transverse stress has an important effect on the biaxial fatigue crack behavior. However, the experimental evidence has provided conflicting indications: it is sometimes considered to increase, decrease or have no effect. These complex phenomena cannot be rationally explained by the existing mechanical models. The effect of the transverse stress on the fatigue crack growth behavior is still one of the most puzzling questions in biaxial fatigue. Physically, this effect is a transverse stress induced plasticity phenomenon. In this paper, a plasticity-corrected stress intensity factor (PC-SIF) is proposed to describe the effect of transverse stress on biaxial fatigue. By use of this new crack driving force some important phenomena associated with transverse stress are predicted. Comparisons with experimental results showed that the PC-SIF as an effective mechanical parameter is capable of predicting the effects of the crack length, the stress level, cyclic stress ratio, biaxial stress ratio and phase difference on the biaxial fatigue crack growth. Consequently, the alleged conflicting experimental results have been rationally explained by the PC-SIF.     

Stress intensity factor solutions for cracks in railway axles

The aim of this paper is a collection of stress intensity factor solutions for cracks in railway axle geometries which the authors of the present special issue developed and/or used for damage tolerance analyses. These solutions comprise closed form analytical as well as tabled geometry functions and they refer to solid as well as hollow axles and various crack sites such as the T- and V-notch and the axle body.     

Stress intensity factor computation using the method of fundamental solutions

In this paper, we study the application of the method of fundamental solutions to the computation of stress intensity factors in linear elastic fracture mechanics. The displacements are approximated by linear combinations of the fundamental solutions of the Cauchy–Navier equations of elasticity and the leading terms for the displacement near the crack tip. The applicability of two formulations of the method is demonstrated on two mode I crack problems, where it is shown that accurate approximations for the stress intensity factors can be obtained with relatively few degrees of freedom. Parts of this work were undertaken while the first author was a Visiting Professor in the Department of Mathematical and Computer Sciences, Colorado School of Mines, Golden, Colorado 80401, U.S.A.     

Stress intensity factors for three‐dimensional weld toe cracks using weld toe magnification factors

In welded components, particularly those with complex geometrical shapes, evaluating stress intensity factors is a difficult task. To effectively calculate the stress intensity factors, a weld toe magnification factor is introduced that can be derived from data obtained in a parametric study performed by finite element method (FEM). Although solutions for the weld toe magnification factor have been presented, these are applicable only to non‐load‐carrying cruciform or T‐butt joints, due possibly to the requirement of very complicated calculations. In the majority of cases for various welded joints, the currently used weld toe magnification factors do not adequately describe the behaviour of weld toe cracks. In this study, the weld toe magnification factor solutions for the three types of welded joints such as cruciform, cover plate and longitudinal stiffener joints were provided through a parametric study using three‐dimensional finite elements. The solutions were formed with exponents and fractions that have polynomial functions in terms of a/c and a/t – that is, crack depths normalised by corresponding half crack lengths and specimen thickness. The proposed weld toe magnification factors were applied to evaluate the fatigue crack propagation life considering the propagation mechanisms of multiple‐surface cracks for all welded joints. It showed good agreement within a deviation factor of two between the experimental and calculated results for the fatigue crack propagation life.     

Stress intensity factor predictions: Comparison and round-off error

Practical steps required to obtain robust finite element triangular meshes for crack path and stress intensity calculation purposes are evaluated, and techniques to use their predictions to calculate fatigue lives, including load interaction effects, are discussed. These steps address: (a) how to simulate efficiently 2D crack paths under bi-axial loading using automatic remeshing schemes; (b) how to choose the best method to calculate stress intensity factors along the crack path; and (c) how the numerical problems associated with excessive FE mesh refinement along the crack path may affect predictions. Various modeling strategies are compared using different crack geometries and mesh refinements to quantify their performance, particularly when the elements around the crack tip are very small compared with the elements far from it. It is shown that, contrary to many other stress analysis applications, excessive mesh refinement may significantly degrade the calculation accuracy in crack problems. A limit for the elements size ratio is clearly established.     

An experimental investigation on fatigue crack growth of AL6XN stainless steel

The crack growth behavior of AL6XN stainless steel was experimentally investigated using round compact tension (CT) specimens. The influences of the R-ratio (the ratio of the minimum load over the maximum applied load in a cycle), the tensile and compressive overloads, and the loading sequence on crack growth were studied in detail. The results from the constant-amplitude experiments show a sensitivity of the crack growth rate to the R-ratio. The application of a tensile overload has a profound effect on crack growth, resulting in a significant retardation in the crack propagation rate. A compressive overload (underload) leads to a short-lived acceleration in crack growth. Results from the two-step high-low loading reveal a period of crack growth retardation at the beginning of the lower amplitude step, an effect similar to that of a single overload. A crack driving force parameter together with a modified Wheeler model is found to correlate the crack growth experiments well.     

On R ratio dependence of threshold stress intensity factor range for fatigue crack growth in type 316(N) stainless steel weld

Abstract

The influence of R ratio in the range 0·05–0·4 on the ambient temperature fatigue crack growth behaviour of an austenitic stainless steel weld, SS 316(N), has been studied. Results indicate that the cyclic threshold stress intensity factor ΔK_th increases with decreasing R ratio. The data are compared with those for SS 316, SS 316L and SS 316L(N) base materials from the literature, and various approaches dealing with the R ratio effects are examined. Zhang’s model considering the contribution of the crack tip plasticity to the fundamental fatigue crack propagation process does provide a consistent interpretation for the data.     

Stress intensity factor estimation for unbalanced rotating cracked shafts by artificial neural networks

The knowledge of the stress intensity factor (SIF) values along a crack front is essential to calculate the crack growth rate and the remaining life of a mechanical component. In the case of a rotating shaft, usually it presents disalignments, which modify the SIF data with regard to a balanced one. This paper presents the use of an artificial neural network (ANN) for estimating the SIF at the crack front in an unbalanced shaft under rotating bending, previously, a quasi‐static numerical (finite element) model, which simulates a rotating shaft, has been developed to create the training cases for the ANN. The obtained results allow to study the influence of the unbalance of rotating shafts in the crack breathing mechanism and will allow to predict the influence of this behaviour on the values of the SIF and in the propagation of cracks.     

An improved model for predicting the crack size and plasticity dependence of fatigue crack propagation

The angled crack problem has been given special attention in the recent years by fracture mechanics investigators due to its close proximity to realistic conditions in engineering structures. In this paper, an investigation of fatigue crack propagation in rectangular steel plates containing an inclined surface crack is presented. The inclined angle of the crack with respect to the axis of loading varied between 0° and 90°. During the fatigue tests, the growth of the fatigue crack was monitored using the AC potential drop technique. A series of modification factors, which allow accurate sizing of such defects, is recommended. Paris power law is normalized and adopted for data analysis. Subsequently, this concept is applied to predict crack growth due to fatigue loads. The results obtained are compared with those obtained using the commonly employed fracture criterion and the experimental data.     

Threshold stress intensity factor and crack growth rate prediction under mixed-mode loading

A new mixed-mode threshold stress intensity factor is developed using a critical plane-based multiaxial fatigue theory and the Kitagawa diagram. The proposed method is a nominal approach since the fatigue damage is evaluated using remote stresses acting on a cracked component rather than stresses near the crack tip. An equivalent stress intensity factor defined on the critical plane is proposed to predict the fatigue crack growth rate under mixed-mode loading. A major advantage is the applicability of the proposed model to many different materials, which experience either shear or tensile dominated crack growth. The proposed model is also capable to nonproportional fatigue loading since the critical plane explicitly considers the influence of the load path. The predictions of the proposed fatigue crack growth model under constant amplitude loading are compared with a wide range of fatigue results in the literature. Excellent agreements between experimental data and model predictions are observed.     

Abaqus implementation of extended finite element method using a level set representation for three-dimensional fatigue crack growth and life predictions

A three-dimensional extended finite element method (X-FEM) coupled with a narrow band fast marching method (FMM) is developed and implemented in the Abaqus finite element package for curvilinear fatigue crack growth and life prediction analysis of metallic structures. Given the level set representation of arbitrary crack geometry, the narrow band FMM provides an efficient way to update the level set values of its evolving crack front. In order to capture the plasticity induced crack closure effect, an element partition and state recovery algorithm for dynamically allocated Gauss points is adopted for efficient integration of historical state variables in the near-tip plastic zone. An element-based penalty approach is also developed to model crack closure and friction. The proposed technique allows arbitrary insertion of initial cracks, independent of a base 3D model, and allows non-self-similar crack growth pattern without conforming to the existing mesh or local remeshing. Several validation examples are presented to demonstrate the extraction of accurate stress intensity factors for both static and growing cracks. Fatigue life prediction of a flawed helicopter lift frame under the ASTERIX spectrum load is presented to demonstrate the analysis procedure and capabilities of the method.     

Stress intensity factor computation using the method of fundamental solutions: mixed‐mode problems

The method of fundamental solutions is applied to the computation of stress intensity factors in linear elastic fracture mechanics. The displacements are approximated by linear combinations of the fundamental solutions of the Cauchy–Navier equations of elasticity and the leading terms for the displacement near the crack tip. Two algorithms are developed, one using a single domain and one using domain decomposition. Numerical results are given. Copyright © 2006 John Wiley & Sons, Ltd.