Microsoft Word - numero_33_art_49 F. Castro et alii, Frattura ed Integrità Strutturale, 33 (2015) 444-450; DOI: 10.3221/IGF-ESIS.33.49 444 Focussed on multiaxial fatigue Estimation of fretting fatigue life using a multiaxial stress-based critical distance methodology F. C. Castro, J. A. Araújo, M. S. T. Pires Department of Engineering Mechanics, University of Brasilia fabiocastro@unb.br L. Susmel Department of Civil and Structural Engineering, The University of Sheffield l.susmel@sheffield.ac.uk ABSTRACT. This work presents a methodology for life estimation of mechanical couplings subjected to fretting fatigue. In this approach, a stress-based multiaxial fatigue parameter is evaluated at a critical distance below the contact surface. The fatigue parameter is based on an improved formulation of the Modified Wöhler Curve Method, in which the shear stress amplitude is measured via the Maximum Rectangular Hull method. To apply the Theory of Critical Distances in the medium-cycle fatigue regime, the critical distance is assumed to depend on the number of cycles to failure. Available fretting fatigue data, conducted on a cylinder-plane contact configuration made of Al alloy 4% Cu, were used to assess the methodology. Most of the life estimates were within an error band given by a factor of 2. KEYWORDS. Fretting fatigue; Multiaxial fatigue; Life estimation; Theory of Critical Distances. INTRODUCTION retting fatigue refers to the conjoint action of a small oscillatory motion between contacting bodies and a cyclic remote loading. The oscillatory motion often leads to surface damage phenomena that may speed up the formation of micro-cracks. Due to the remote loading, these cracks may propagate until catastrophic failure occurs. Many engineering assemblies are prone to fretting fatigue problems as, for instance, the blade/disk interfaces of gas turbine engines, wires of overhead conductors and riveted joints. Fretting fatigue usually displays high stress gradients that affect observed lives [1, 2]. Among the different formulations that account for the stress gradient in fretting, non-local approaches have received considerable attention in the literature [3]. The appeal of such approaches derives from their relative simplicity and from satisfactory correlations with experimental data [2, 4-7]. The Theory of Critical Distances (TCD) [8-10] has been one of the most used non-local approaches in the last ten years. In this approach, the size of the fatigue process zone is related to fatigue thresholds of cracked or sharply notched specimens. The success of the TCD in correlating fatigue thresholds of notched members has been widely reported in the literature (see [9, 10] and references therein). Due to the similarities between notch and fretting fatigue, an attempt to estimate fretting fatigue thresholds using the TCD has been carried out by Araújo et al. [6]. In that paper, the Modified F F. Castro et alii, Frattura ed Integrità Strutturale, 33 (2015) 444-450; DOI: 10.3221/IGF-ESIS.33.49 445 Wöhler Curve Method (MWCM) [11] was applied at a point located at a critical distance below the trailing edge of the contact. Estimates of fretting fatigue thresholds fell within an error interval of ±20% when compared to experimental data. The methodology proposed by Araújo et al. [6] is only applicable to design situations involving threshold conditions. In this paper, an extension of this approach to the medium-cycle fatigue regime is presented. A first attempt to assess the new methodology is carried out based on available fretting fatigue tests [1]. MULTIAXIAL FATIGUE LIFE ESTIMATION he MWCM [10-12] is a multiaxial stress-based critical plane approach where the driving parameters for crack nucleation are the maximum shear stress amplitude, a, and the maximum normal stress acting on the maximum shear stress plane, n,max. Once the values of a and n,max have been evaluated, the stress ratio  is defined as n,max a     (1) It is noteworthy that  is sensitive not only to mean stresses, but also to the degree of multiaxiality and non- proportionality of the stress path [10]. Also, it is worth recalling that for an unnotched specimen under fully reversed uniaxial loading the  ratio is equal to unity, whereas under fully reversed torsional loading  equals zero [10]. The MWCM is based on a modified Wöhler diagram where a is plotted against the number of cycles to failure, Nf (Fig. 1). This diagram is made of different fatigue curves, each one corresponding to a certain value of  ratio and being unambiguously described by its negative inverse slope  and by a reference shear stress amplitude, A,Ref, corresponding to an appropriate number of cycles to failure, NA. To obtain the diagram, one should properly define the  vs.  and A,Ref vs.  relationships, and correctly calibrate them by running appropriate experiments. The following linear relationships have been found by Lazzarin and Susmel [12] to correlate a wide range of experimental data: ( ) ba    (2) A,Ref ( ) c d    (3) where a, b, c and d are material constants. When these constants are obtained from fully-reversed uniaxial and torsional tests on plain specimens, Eqs. (2) and (3) are stated as Figure 1: Modified Wöhler diagram. ( ) [ ( ) (1 0)] ( )0             (4) 0 A,Ref 0 0( ) 2            (5) where (=1) and 0 are, respectively, the inverse slope of the modified Wöhler curve and the fatigue limit under uniaxial loading condition, whereas (=0) and 0 are the corresponding quantities for torsional loading. It is noteworthy that, for T F. Castro et alii, Frattura ed Integrità Strutturale, 33 (2015) 444-450; DOI: 10.3221/IGF-ESIS.33.49 446 materials that do not exhibit a fatigue limit, 0 and 0 must be defined as endurance limits corresponding to an appropriate number of cycles to failure. After a proper calibration of Eqs. (2) and (3), any curve of the modified Wöhler diagram can be obtained. Hence, the number of cycles to failure can be estimated as A,Ref ( ) f,e A a ( ) N N            (6) Usually, the value of a in the MWCM is determined via the Minimum Circumscribed Circle (MCC) method [13]. Here, the Maximum Rectangular Hull (MRH) approach is adopted, since fatigue estimates based on the MWCM are improved when a is measured by the MRH rather than by the MCC [14]. The MRH method is schematically represented in Fig. 2. The halves of the sides of the rectangular hull with orientation  are calculated as 1 max , m( ) ( in , , 1( 2 2 ) ) ,i i i tt t ia t           (7) where i(,t) (i=1,2) are the components of the shear stress vector (t) with respect to the -oriented frame. The amplitude of each -oriented rectangular hull can be evaluated as 2 2 a 1 2( ) )) (( a a    (8) The shear stress amplitude is then defined as the maximum value of Eq. (8) among all  -oriented rectangular hulls: 2 2 a MRH 1 2 90 º0 max ( ) ( )a a        (9) Figure 2: Amplitudes of the -oriented rectangular hull bounding the shear stress path. CRITICAL DISTANCE APPROACH he stress gradient effect is a crucial aspect in the design of stress raisers such as notches and mechanical contacts. Among the formulations that account for this effect, the TCD [8-10] has been recognized as one of the most attractive due to its simplicity and good results for a number of notch configurations. The central idea of the TCD is the definition of an effective stress, eff, based on an averaging procedure over a volume surrounding the stress raiser. Fatigue failure is expected to occur if eff exceeds a reference material strength, ref. Simplified methods can also be formulated by considering averages over an area or a line (Area and Line Methods, respectively) or the effective stress of the point located at a critical distance, L, from the stress raiser (Point Method). In this paper attention is focused on the Point Method, as schematically illustrated in Fig. 3. Taylor [8] has developed fitting procedures to determine the critical distance, which are based on the fatigue thresholds of cracked or sharply notches specimens. In particular, the critical distance for the Point Method has been found to be expressed as 2 th 0 1 2 L K          (10) T F. Castro et alii, Frattura ed Integrità Strutturale, 33 (2015) 444-450; DOI: 10.3221/IGF-ESIS.33.49 447 where Kth is the threshold stress intensity factor range and 0 is the uniaxial plain fatigue limit range. Further work [10, 15] has investigated the extension of the TCD to stress raisers under multiaxial loadings, concluding that the combination of the MWCM with the Point Method requires the same critical distance given by Eq. (10). Figure 3: Schematic representation of the Point Method. An extension of the TCD to estimate fatigue life of stress raisers has been developed in Refs. [16, 17]. To explain the formulation, it should be recalled that the appropriate critical distance to estimate static failure of notched members is given as [9, 18]: I 2 S c r 1 2 K L          (11) where KIc is the plane strain material fracture toughness and r is a reference material constant which can be equal or larger than the ultimate tensile strength, UTS [9]. As the values of the critical distances at the threshold and static conditions are usually different, it can be assumed that the critical distance at the medium-cycle fatigue regime, LM, depends on the number of cycles to failure, Nf. In particular, a power law relationship between LM and Nf has been proposed by Susmel and Taylor as follows: M f f( ) BL N AN (12) where A and B are material constants. Two fitting procedures have been developed to obtain these constants: one based on critical distances determined at threshold and static conditions (Eqs. (10) and (11), respectively), and another based on fatigue curves of plain and sharply notched specimens. Although the latter procedure has proved to provide more accurate notch life estimates when compared to experimental data, the former is simpler to work with as the material constants required to determine A and B are usually available or can be extracted from empirical correlations. APPLICATION TO FRETTING FATIGUE everal attempts to address the fretting fatigue problem using notch methodologies have been investigated in the literature [2-7] due to the similarities between notch and fretting fatigue. Indeed, the stress fields in both problems are characterized by stress gradients and multiaxial stresses. In the fretting fatigue problem, however, there is also a wear process due to the relative motion between the contacting surfaces. In this setting, notch methodologies could be applied to fretting fatigue if the surface damage could be regarded as negligible. This approximation is considered in this paper, at least for the partial slip case where the amount of wear debris is usually small [19]. The use of the methodology for a typical fretting fatigue problem is shown in Fig. 4. The contact configuration involves a normal force P, a cyclic tangential loading Q(t), and a cyclic remote stress (t). As a first step, the crack initiation point must be determined. This task can be accomplished, for instance, by searching the point where a given fatigue parameter achieves its maximum value. Normally, the crack initiation point occurs at the trailing edge of the contact, as illustrated in Fig. 4a. Subsequent analysis is carried out on the straight line that emanates from the crack initiation point and is perpendicular to the contact surface. In order to obtain the number of cycles to failure, the critical distance corresponding to a trial number of cycles to failure, N, is determined as follows: M ( ) BL N AN (13) S F. Castro et alii, Frattura ed Integrità Strutturale, 33 (2015) 444-450; DOI: 10.3221/IGF-ESIS.33.49 448 The stress quantities a, n,max, and  are then evaluated at this critical distance, and the number of cycles to failure, Nf,e, is calculated using Eq. (6). If the calculated Nf,e is different from the trial value N, a new iteration is started considering N = Nf,e. This process is repeated until convergence between the trial and calculated fatigue lives has been attained. Figure 4: Procedure to estimate fretting fatigue life: (a) location of the critical distance and (b) flowchart of the iterative procedure for life estimation. COMPARISON WITH EXPERIMENTAL DATA vailable fretting fatigue data [1] were used to assess the methodology. Such tests were carried out with a pair of cylindrical pads pressed against a flat dog-done specimen, both made of an Al alloy 4% Cu. In each series of tests, the tests were run (each one with a different pad radius) at constant values of peak contact pressure, tangential force and remote stress. Hence, the magnitude of the stress field of each test was the same, but different stress gradients were induced by the contact. Constants for the methodology are given in Tab. 1. The constants A and B in the LM versus Nf relationship, Eq. (12), were determined using L values calculated at the threshold and static conditions, Eqs. (10) and (11), respectively. The values of  and A,Ref were extracted from fully reversed uniaxial and torsional fatigue curves. These curves were estimated by fitting fatigue data for 103 cycles and 107 cycles, using empirical relationships [20] to obtain the fatigue strengths at 103 cycles from the ultimate tensile strength. UTS 0 KIC Kth (=0) A,Ref (=0) (=1) A,Ref (=1) (MPa) (MPa) (MPam0.5) (MPam0.5) (MPa) (MPa) 500 124 34 4.4 12.8 161 12.8 115 Table 1: Constants for the methodology for Al alloy 4% Cu. Due to the geometry of the experimental setup and the applied loadings, analytical techniques [21] were employed to solve the elastic contact problem. The surface tractions are characterized by a Hertzian contact pressure distribution, and by shear tractions that are similar to the Mindlin-Cattaneo one except that the stick zone is not symmetrical with respect to the center of the contact zone but shifted due to the presence of an alternating remote stress. Once the surface tractions have been determined, subsurface stresses can be obtained by using a Muskhelishvili potential. The time varying elastic stress field in any material point in the specimen can be finally calculated by superposing the effects of contact pressure, shear traction and remote stress. A F. Castro et alii, Frattura ed Integrità Strutturale, 33 (2015) 444-450; DOI: 10.3221/IGF-ESIS.33.49 449 Estimated and observed number of cycles to failure is shown in Fig. 5. The solid diagonal line corresponds to a perfect correlation between estimated and observed fatigue lives, and the two dashed lines define the factor of 2 bandwidth. As can be clearly seen in this figure, most of the life estimates fall within a error band given by a factor of 2. Considering the well known scatter that characterizes the fatigue phenomena this can be considered a very good correlation. Figure 5: Observed and estimated fatigue lives for fretting fatigue tests conducted by Nowell [1]. CONCLUSIONS n engineering methodology for fatigue life estimation of mechanical couplings subjected to fretting fatigue was presented in this paper. Constants for the methodology are relatively simple to determine, as only experimental data from conventional fatigue tests on plain and sharply notched (or cracked) specimens are required. It can also be easily incorporated in a Finite Element Method based software for fast assessment of fretting fatigue life in more realistic type of mechanical couplings. Fatigue life estimates correlated well with the experimental results [1], falling in most cases within a factor of two bandwidth. Although some of the required material constants were obtained either from literature or from empirical relations, the methodology was still capable of producing satisfactory life estimates. However, further assessment considering different contact configurations and other materials is still required to corroborate the proposed life methodology. ACKNOWLEDGMENTS he supports provided by CNPq (contracts 310845/2013-0 and 309748/2013-5) and by FINATEC are gratefully acknowledged. REFERENCES [1] Nowell, D., An analysis of fretting fatigue. Ph.D. thesis, Oxford University, (1988). [2] Amargier, R., Fouvry, S., Chambon, L., Schwob, C., Poupon, C., Stress gradient effect on crack initiation in fretting using a multiaxial fatigue framework, Int J Fatigue, 32 (2010) 1904–1912. [3] Nowell, D., Dini, D., Hills, D.A., Recent developments in the understanding of fretting fatigue, Eng Fract Mech, 73 (2006) 207–222. A T F. Castro et alii, Frattura ed Integrità Strutturale, 33 (2015) 444-450; DOI: 10.3221/IGF-ESIS.33.49 450 [4] Araújo, J.A., Nowell, D., The effect of rapidly varying contact stress fields on fretting fatigue, Int J Fatigue, 24 (2002) 763–776. [5] Araújo, J.A., Castro, F.C., A comparative analysis between multiaxial stress and ΔK-based short crack arrest models in fretting fatigue Engineering Fracture Mechanics, 93 (2012) 34–47. [6] Araújo, J.A., Susmel, L., Taylor, D., Ferro J.C.T., Mamiya, E.N., On the use of the theory of critical distances and the Modified Wöhler curve method to estimate fretting fatigue strength of cylindrical contacts, Int J Fatigue, 29 (2007) 95–107. [7] Navarro, C., Muñoz, S., Domínguez, J., On the use of multiaxial fatigue criteria for fretting fatigue life assessment, Int J Fatigue, 30 (2008) 32–44. [8] Taylor, D., Geometrical effects in fatigue: a unifying theoretical model. Int J Fatigue, 21 (1999) 413–20. [9] Taylor, D., The Theory of Critical Distances: a new perspective in fracture mechanics. Elsevier BV, (2007). [10] Susmel, L., Multiaxial Notch Fatigue: from nominal to local stress-strain quantities. Woodhead & CRC, Cambridge, UK, (2009). [11] Susmel, L., Lazzarin, P., A bi-parametric Wöhler curve for high cycle multiaxial fatigue assessment, Fatigue Fract Eng Mater Struct, 25 (2002) 63–78. [12] Lazzarin, P., Susmel, L., A stress-based method to predict lifetime under multiaxial fatigue loading, Fatigue Fract Eng Mater Struct, 26 (2003) 1171–1187. [13] Castro, F.C., Araújo, J.A., Zouain, N. On the application of multiaxial high-cycle fatigue criteria using the theory of critical distances. Engng Fract Mech. 2009; 76:512–24. [14] Araújo, J.A., Dantas A.P., Castro F.C., Mamiya E.N., Ferreira J.L.A., On the characterization of the critical plane with a simple and fast alternative measure of the shear stress amplitude in multiaxial fatigue, Int J Fatigue, 33 (2011) 1092– 1100. [15] Susmel, L., Taylor, D., Two methods for predicting the multiaxial fatigue limits of sharp notches, Fatigue Fract Engng Mater Struct, 26 (2003) 821–833. [16] Susmel, L., Taylor, D., A novel formulation of the theory of critical distances to estimate lifetime of notched components in the medium-cycle fatigue regime, Fatigue Fract Engng Mater Struct, 30 (2007) 567–581. [17] Susmel, L., Taylor, D., The Modified Wöhler Curve Method applied along with the Theory of Critical Distances to estimate finite life of notched components subjected to complex multiaxial loading paths, Fatigue Fract Engng Mater Struct, 31 (2008) 1047–1064. [18] Susmel L, Taylor D., On the use of the theory of critical distances to predict static failures in ductile metallic materials containing different geometrical features, Eng Fract Mech, 75 (2008) 4410–4421. [19] Vingsbo, O., Söderberg, S., On fretting maps, Wear, 126 (1988) 131-147. [20] Susmel, L., Fatigue Design: Summary, Class notes, (2012). [21] Hills, D.A., Nowell, D., Mechanics of Fretting Fatigue. Dordrecht: Kluwer Academic Publishers, (1994).