Comparative analysis of fatigue life predictions in central notched specimens of Al–Mg–Si alloys

The fatigue behaviour of an Al–Mg–Si alloy was studied using notched specimens. Fatigue tests were conducted at two stress ratios R= 0 and R= 0.4 on thin plates with a central hole. Constant and block variable loading amplitudes were applied to the specimens using a servo‐hydraulic machine. The applicability of the local strain approach method to the prediction of the fatigue life was investigated for this type of discontinuity. Two methods, the equivalent strain energy density approach and a modified stress–strain intensity field approach, were used to predict the fatigue strength. For the second one an elastic–plastic finite element analysis was carried out in order to obtain the local strain and stress distributions near the notch root. Based on Miner's rule an equivalent stress was used to correlate the fatigue lives for the variable amplitude histories. The experimental results were compared with the predicted results obtained by the two methods investigated and better agreement was found with the stress–strain field intensity approach, while the strain energy approach gave more conservative results. Miner's rule gives a good correlation between the variable amplitude and constant amplitude results.     

Multiaxial notch fatigue life prediction based on pseudo stress correction and finite element analysis under variable amplitude loading

A new computational methodology is proposed for fatigue life prediction of notched components subjected to variable amplitude multiaxial loading. In the proposed methodology, an estimation method of non‐proportionality factor (F) proposed by authors in the case of constant amplitude multiaxial loading is extended and applied to variable amplitude multiaxial loading by using Wang‐Brown's reversal counting approach. The pseudo stress correction method integrated with linear elastic finite element analysis is utilized to calculate the local elastic‐plastic stress and strain responses at the notch root. For whole local strain history, the plane with weight‐averaged maximum shear strain range is defined as the critical plane in this study. Based on the defined critical plane, a multiaxial fatigue damage model combined with Miner's linear cumulative damage law is used to predict fatigue life. The experimentally obtained fatigue data for 7050‐T7451 aluminium alloy notched shaft specimens under constant and variable amplitude multiaxial loadings are used to verify the proposed methodology and equivalent strain‐based methodology. The results show that the proposed methodology is superior to equivalent strain‐based methodology.     

Anisotropic Gurson‐Tvergaard‐Needleman plasticity and damage model for finite element analysis of elastic‐plastic problems

Implementation and analysis of the anisotropic version of the Gurson‐Tvergaard‐Needleman (GTN) isotropic damage criterion are performed on the basis of Hill's quadratic anisotropic yield theory with the definition of an effective anisotropic coefficient to represent the elastic‐plastic behavior of ductile metals. This study aims to analyze the extension of the GTN model suitable for anisotropic porous metals and to investigate the GTN model extension. An anisotropic damage model is implemented using the user material subroutine in ABAQUS/standard finite element code. The implementation is verified and applied to simulate a uniaxial tensile test on a commercially produced aluminum sheet material for three‐dimensional and plane stress test cases. Spherical and ellipsoidal micro voids are considered in the matrix material, and their effects on the uniaxial stress‐strain response of the material are analyzed. Hill's quadratic anisotropic yield theory predicts substantially large damage evolution and a low stress‐strain curve compared with those predicted by the isotropic model. An approximate model for anisotropic materials is proposed to avoid increased damage evolution. In this approximate model, Hill's anisotropic constants are replaced with an effective anisotropy coefficient. All model‐generated stress‐strain predictions are compared with the experimental stress‐strain curve of AA6016‐T4 alloy.     

State‐of‐the‐art review on extended stress intensity factor concepts

The stress intensity factor concept for describing the stress field at pointed crack or slit tips is well known from fracture mechanics. It has been substantially extended since Williams' basic contribution (1952) on stress fields at angular corners. One extension refers to pointed V‐notches with stress intensities depending on the notch opening angle. The loading‐mode‐related simple notch stress intensity factors K₁, K₂ and K₃ are introduced. Another extension refers to rounded notches with crack shape or V‐notch shape in two variants: parabolic, elliptic or hyperbolic notches (‘blunt notches’) on the one hand and root hole notches (‘keyholes’ when considering crack shapes) on the other hand. Here, the loading‐mode‐related generalised notch stress intensity factors K_1ρ, K_2ρ and K_3ρ are defined. The concepts of elastic stress intensity factor, notch stress intensity factor and generalised notch stress intensity factor are extended into the range of elastic–plastic (work‐hardening) or perfectly plastic notch tip or notch root behaviour. Here, the plastic notch stress intensity factors K_1p, K_2p and K_3p are of relevance. The elastic notch stress intensity factors are used to describe the fatigue strength of fillet‐welded attachment joints. The fracture toughness of brittle materials may also be evaluated on this basis. The plastic notch stress intensity factors characterise the stress and strain field at pointed V‐notch tips. A new version of the Neuber rule accounting for the influence of the notch opening angle is presented.     

Generalized stress intensity factors due to steady and transient thermal loads with applications to welded joints

This paper deals with the elastic and plastic stress fields induced by thermal loads in the vicinity of sharp V‐notch tips in plates. Under the hypothesis of steady‐state heat transfer and plane‐strain conditions, the thermal and mechanical problem requires the numerical solution of an ordinary differential equation (ODE) system, obtained by extending the ‘stress function approach’. The intensity of the stress distributions ahead of V‐notch tips can be expressed in terms of thermal notch stress intensity factors (thermal NSIFs), as for external loads. The problem becomes much more demanding in the presence of transient thermal loads. The residual asymptotic stress distribution arising from the solidification of a fusion zone during an arc welding process is obtained by considering different boundary conditions. An aluminium butt‐welded joint is analysed after having modelled the weld toe region as a sharp V‐notch. A finite element (FE) simulation of the welding process is carried out by means of SYSWELD code (version 2004.1) modelling the arc welding torch by means of Goldak's source. Near the weld toe, the intensity of the residual stress field is given in terms of elastic or elastic—plastic generalized NSIFs.     

The Equivalent Strain Energy Density approach re-formulated and applied to sharp V-shaped notches under localized and generalized plasticity

In the case of a rounded notch, the stress and strain at the notch tip can be determined by the traditional Neuber rule or by the Equivalent Strain Energy Density (ESED) approach, as formulated by Glinka and Molski. In the latter case the elastoplastic strain energy density at the notch tip is thought of as coincident with that determined under purely elastic conditions. For sharply V‐shaped notches this approach is not directly applicable, since the strain energy density at the notch tip tends toward infinity both for a material obeying an elastic law and a material obeying a power hardening law. By using the notch stress intensity factors, the present paper suggests a re‐formulation of the ESED approach which is applied no longer at the notch tip but to a finite size circular sector surrounding the notch tip. In particular we have adopted the hypothesis that, under plane strain conditions, the value of the energy concentration due to the notch is constant and independent of the two constitutive laws. When small scale yielding conditions are present, such a hypothesis immediately results in the constancy of the strain energy averaged over the process volume. As a consequence, plastic notch stress intensity factors valid for sharp V‐shaped notches can be predicted on the basis of the linear elastic stress distributions alone.     

Two‐dimensional stress distributions in tensioned orthotropic plates weakened by blunt V‐shaped notches

In this work, a comprehensive analytical and numerical study on the two‐dimensional stress fields in orthotropic plates with blunt notches is carried out. Initially, a quick review of the available solutions for lateral notches is presented. Later on, using Lekhnitskii's approach, a new analytical solution is derived which explicitly accounts for the local geometry of the notch, such as the notch root radius and the notch opening angle, as well as the elastic properties of the material. Theoretical predictions based on the new solution are compared with the results from a bulk of 2D finite element analyses carried out on tensioned plates weakened by lateral hyperbolic and blunt V‐shaped notches, showing a very satisfactory agreement.     

A low‐order,hexahedral finite element for modelling shells

A thin, eight‐node, tri‐linear displacement, hexahedral finite element is the starting point for the derivation of a constant membrane stress resultant, constant bending stress resultant shell finite element. The derivation begins by introducing a Taylor series expansion for the stress distribution in the isoparametric co‐ordinates of the element. The effect of the Taylor series expansion for the stress distribution is to explicitly identify those strain modes of the element that are conjugate to the mean or average stress and the linear variation in stress. The constant membrane stress resultants are identified with the mean stress components, and the constant bending stress resultants are identified with the linear variation in stress through the thickness along with in‐plane linear variations of selected components of the transverse shear stress. Further, a plane‐stress constitutive assumption is introduced, and an explicit treatment of the finite element's thickness is introduced. A number of elastic simulations show the useful results that can be obtained (tip‐loaded twisted beam, point‐loaded hemisphere, point‐loaded sphere, tip‐loaded Raasch hook, and a beam bent into a ring). All of the gradient/divergence operators are evaluated in closed form providing unequivocal evaluations of membrane and bending strain rates along with the appropriate divergence calculations involving the membrane stress and bending stress resultants. The fact that a hexahedral shell finite element has two distinct surfaces aids sliding interface algorithms when a shell folds back on itself when subjected to large deformations. Published in 2004 by John Wiley & Sons, Ltd.     

Sharp V‐notches in viscoplastic solids: Strain energy rate density rule and fracture toughness

Different from Neuber's rule or Glinka's energy method which are always adopted to characterize the notch tip field under elastoplastic condition, in this paper, the strain energy rate density (SERD) rule is used for viscoplastic materials. In particular, based on the definition of generalized notch stress intensity factor (G‐NSIF) for sharp V‐notch in viscoplastic solids, the concept of SERD for sharp V‐notch in viscoplastic solids is presented. Subsequently, by taking as a starting point the SERD, the averaged strain energy density (SED) for sharp V‐notch in viscoplastic solids is derived with integration of time. The fracture toughness relation between sharp V‐notch specimens and crack specimen in viscoplastic materials is given based on the transformation of SERD. A numerical approach is presented to compute the SERD and SED based on finite element method. Some crucial comments on the G‐NSIF have been discussed. Some typical solutions for SERD and SED for sharp V‐notched specimens are investigated.     

Formulation and numerical treatment of incompressibility constraints in large strain elastic–plastic analysis

The present paper is concerned with an efficient framework for a nonlinear finite element procedure for the rate‐independent finite strain analysis of solids undergoing large elastic‐isochoric plastic deformations. The formulation relies on the introduction of a mixed‐variant metric deformation tensor which will be multiplicatively decomposed into a plastic and an elastic part. This leads to the definition of an appropriate logarithmic strain measure which can be additively decomposed into the exact isochoric (deviatoric) and volumetric (spheric) strain measures. This fact may be seen as the basic idea in the formulation of appropriate mixed finite elements which guarantee the accurate computation of isochoric strains. The mixed‐variant logarithmic elastic strain tensor provides a basis for the definition of a local isotropic hyperelastic stress response whereas the plastic material behavior is assumed to be governed by a generalized J₂ yield criterion and rate‐independent isochoric plastic strain rates are computed using an associated flow rule. On the numerical side, the computation of the logarithmic strain tensors is based on higher‐order Padé approximations. To be able to take into account the plastic incompressibility constraint a modified mixed variational principle is considered which leads to a quasi‐displacement finite element procedure. Finally, the numerical solution of finite strain elastic‐plastic problems is presented to demonstrate the efficiency and the accuracy of the algorithm. Copyright © 1999 John Wiley & Sons, Ltd.     

A new expression for crack opening stress determined based on maximum crack opening displacement under tension–compression cyclic loading

The maximum crack opening displacement is introduced to investigate the effect of compressive loads on crack opening stress in tension–compression loading cycles. Based on elastic–plastic finite element analysis of centre cracked finite plate and accounting for the effects of crack geometry size, Young's modulus, yield stress and strain hardening, the explicit expression of crack opening stress versus maximum crack opening displacement is presented. This model considers the effect of compressive loads on crack opening stress and avoids adopting fracture parameters around crack tip. Besides, it could be applied in a wide range of materials and load conditions. Further studies show that experimental results of da/dN ? ΔK curves with negative stress ratios could be condensed to a single curve using this crack opening stress model.     

Driving forces for localized corrosion‐to‐fatigue crack transition in Al–Zn–Mg–Cu

Research on fatigue crack formation from a corroded 7075‐T651 surface provides insight into the governing mechanical driving forces at microstructure‐scale lengths that are intermediate between safe life and damage tolerant feature sizes. Crack surface marker‐bands accurately quantify cycles (N_i) to form a 10–20 μm fatigue crack emanating from both an isolated pit perimeter and EXCO corroded surface. The N_i decreases with increasing‐applied stress. Fatigue crack formation involves a complex interaction of elastic stress concentration due to three‐dimensional pit macro‐topography coupled with local micro‐topographic plastic strain concentration, further enhanced by microstructure (particularly sub‐surface constituents). These driving force interactions lead to high variability in cycles to form a fatigue crack, but from an engineering perspective, a broadly corroded surface should contain an extreme group of features that are likely to drive the portion of life to form a crack to near 0. At low‐applied stresses, crack formation can constitute a significant portion of life, which is predicted by coupling macro‐pit and micro‐feature elastic–plastic stress/strain concentrations from finite element analysis with empirical low‐cycle fatigue life models. The presented experimental results provide a foundation to validate next‐generation crack formation models and prognosis methods.     

Creep fatigue of a directionally solidified Ni‐base superalloy – smooth and cylindrically notched specimens

Turbine blade life modelling is complicated by the presence of notches, dwells, high temperatures, thermal cycles and temperature gradients. Furthermore, directionally solidified (DS) Ni‐base superalloys are highly anisotropic. This work seeks to characterize the response of the DS Ni‐base superalloy CM247LC subjected to isothermal low cycle fatigue at either 750 or 950 °C. This study considers the effects of strain rate, dwells at the maximum temperature, and stress concentrations. Experiments were conducted under uniaxial loading on smooth and cylindrically notched round‐bar specimens in both longitudinal and transverse orientations. The location of the creep‐fatigue crack is at the maximum Hill's effective stress in the notched specimens. In addition, the notch behaviour is discussed in light of finite element analysis using an anisotropic elastic‐crystal viscoplastic material model.     

Effects of approximations in analyses of beams of open thin‐walled cross‐section—part II: 3‐D non‐linear behaviour

In a companion paper, the effects of approximations in the flexural‐torsional stability analysis of beams was studied, and it was shown that a second‐order rotation matrix was sufficiently accurate for a flexural‐torsional stability analysis. However, the second‐order rotation matrix is not necessarily accurate in formulating finite element model for a 3‐D non‐linear analysis of thin‐walled beams of open cross‐section. The approximations in the second‐order rotation matrix may introduce ‘self‐straining’ due to superimposed rigid‐body motions, which may lead to physically incorrect predictions of the 3‐D non‐linear behaviour of beams. In a 3‐D non‐linear elastic–plastic analysis, numerical integration over the cross‐section is usually used to check the yield criterion and to calculate the stress increments, the stress resultants, the elastic–plastic stress–strain matrix and the tangent modulus matrix. A scheme of the arrangement of sampling points over the cross‐section that is not consistent with the strain distributions may lead to incorrect predictions of the 3‐D non‐linear elastic–plastic behaviour of beams. This paper investigates the effects of approximations on the 3‐D non‐linear analysis of beams. It is found that a finite element model for 3‐D non‐linear analysis based on the second‐order rotation matrix leads to over‐stiff predictions of the flexural‐torsional buckling and postbuckling response and to an overestimate of the maximum load‐carrying capacities of beams in some cases. To perform a correct 3‐D non‐linear analysis of beams, an accurate model of the rotations must be used. A scheme of the arrangement of sampling points over the cross‐section that is consistent with both the longitudinal normal and shear strain distributions is needed to predict the correct 3‐D non‐linear elastic–plastic behaviour of beams. Copyright © 2001 John Wiley & Sons, Ltd.     

Strain rate concentration and dynamic stress concentration for double‐edge‐notched specimens subjected to high‐speed tensile loads

Engineering plastics provide superior performance to ordinary plastics for wide range of the use. For polymer materials, dynamic stress and strain rate may be major factors to be considered when the strength is evaluated. Recently, high‐speed tensile test is being recognized as a standard testing method to confirm the strength under dynamic loads. In this study, therefore, high‐speed tensile test is analysed by the finite element method; then, the maximum dynamic stress and strain rate are discussed with varying the tensile speed and maximum forced displacement. The maximum strain rate increases with increasing the tensile speed u/t, but the strain rate concentration factor is found to be constant independent of tensile speed, which is defined as the maximum strain rate appearing at the notch root over the average nominal strain rate at the minimum section . It is found that the strain rate at the notch root depends on the dynamic stress rate at the notch root and independent of the notch root radius ρ. It is found that the difference between the static and dynamic maximum stress concentration (σ_yA,max ? σ_yA,st) at the notch root is proportional to the tensile speed when u/t = 5000 mm/s. Strain rate concentration factors are also discussed with varying the notch depth and specimen length. Based on the elastic strain rate concentration factor, the master curve is obtained useful for understanding the impact fracture of polycarbonate for the wide range of temperature and impact speed.     

Notch fatigue life prediction considering nonproportionality of local loading path under multiaxial cyclic loading

An innovative numerical methodology is presented for fatigue lifetime estimation of notched bodies experiencing multiaxial cyclic loadings. In the presented methodology, an evaluation approach of the local nonproportionality factor F for notched specimens, which defines F as the ratio of the pseudoshear strain range at 45° to the maximum shear plane and the maximum shear strain range, is proposed and discussed deeply. The proposed evaluation method is incorporated into the material cyclic stress‐strain equation for purpose of describing the nonproportional hardening behavior for some material. The comparison between multiaxial elastic‐plastic finite element analysis (FEA) and experimentally measured strains for S460N steel notched specimens shows that the proposed nonproportionality factor estimation method is effective. Subsequently, the notch stresses and strains calculated utilizing multiaxial elastic‐plastic FEA are used as input data to the critical plane‐based fatigue life prediction methodology. The prediction results are satisfactory for the 7050‐T7451 aluminum alloy and GH4169 superalloy notched specimens under multiaxial cyclic loading.     

A rational elasto‐plastic spatially curved thin‐walled beam element

Torsion is one of the primary actions in members curved in space, and so an accurate spatially curved‐beam element needs to be able to predict the elasto‐plastic torsional behaviour of such members correctly. However, there are two major difficulties in most existing finite thin‐walled beam elements, such as in ABAQUS and ANSYS, which may lead to incorrect predictions of the elasto‐plastic behaviour of members curved in space. Firstly, the integration sample point scheme cannot capture the shear strain and stress information resulting from uniform torsion. Secondly, the higher‐order twists are ignored which leads to loss of the significant effects of Wagner moments on the large twist torsional behaviour. In addition, the initial geometric imperfections and residual stresses are significant for the elasto‐plastic behaviour of members curved in space. Many existing finite thin‐walled beam element models do not provide facilities to deal with initial geometric imperfections. Although ABAQUS and ANSYS have facilities for the input of residual stresses as initial stresses, they cannot describe the complicated distribution patterns of residual stresses in thin‐walled members. Furthermore, external loads and elastic restraints may be applied remote from shear centres or centroids. The effects of the load (and restraint) positions are important, but are not considered in many beam elements. This paper presents an elasto‐plastic spatially curved element with arbitrary thin‐walled cross‐sections that can correctly capture the uniform shear strain and stress information for integration, and includes initial geometric imperfections, residual stresses and the effects of the load and restraint positions. The element also includes elastic restraints and supports, which have to be modelled separately as spring elements in some other finite thin‐walled beam elements. Comparisons with existing experimental and analytical results show that the elasto‐plastic spatially curved‐beam element is accurate and efficient. Copyright © 2006 John Wiley & Sons, Ltd.     

A superposition approach to the stress fields for various blunt V‐notches

This paper deals with the problems of blunt V‐notch with various notch shapes. The purpose is to develop a new method capable of obtaining more accurate solutions for the stress fields around a blunt V‐notch tip under opening and sliding modes. The key method is to use the principle of superposition for linear elastic materials. On the basis of the superposition method and the conventional stress fields for a sharp V‐notch, the stress fields useful for any shapes of blunt V‐notch is proposed. The notch stress intensity factors are estimated by the numerical analysis with finite element analysis, and then the effectiveness and validation of the proposed superposition approach are discussed by comparison with the results from the literature.     

Vergleichende Untersuchungen zur Anrisslebensdauervorhersage gekerbter Proben mit dem Örtlichen Konzept und dem Nennspannungskonzept am Beispiel des Stahls Cm15

Comparative Investigations on Service Life Assessment of Notched Specimens Based on the Local Strain and the Nominal Stress Approach to Fatigue for a Steel SAE 1017 It is still unclear whether the strain based approach to fatigue or the stress based approach to fatigue should be preferred for service life assessment of notched components. In order to clarify the similarities and differences between these concepts stress and strain controlled fatigue experiments have been performed with notched specimens. It has been found, that stress and strain controlled fatigue testing results in the same number of cycles until failure. Essential for this correlation is that the cyclic stable strain amplitude at the notch root is taken for the entry into the strain‐life diagram in both cases. Starting from an elastic‐plastic analysis of the material behaviour at the notch root it is shown, how the strain‐life curve can be converted into a stress‐life curve. Based on that result service‐life is calculated from both approaches mentioned above. The calculation gives nearly the same service‐lives for both cases, but overestimates the measured data. It becomes obvious, that a S‐N curve determined under one‐level loading doesn’t provide a proper basis for service life assessment. While strain or stress‐life curves always contain crack initiation phase as well as crack propagation phase, the fatigue process under irregular loads is mainly governed by crack propagation. As a consequence, the damage per cycle is underestimated for loads near the fatigue limit, if Miner’s rule is used.