Microsoft Word - numero_33_art_47 J.A Araujo et al, Frattura ed Integrità Strutturale, 33 (2015) 427-433; DOI: 10.3221/IGF-ESIS.33.47 427 Focussed on multiaxial fatigue Equivalent configurations for notch and fretting fatigue J. A. Araújo, F. C. Castro University of Brasilia, UNB, Departamento de Engenharia Mecanica alex07@unb.br S. Pommier, J. Bellecave LMT-Cachan, ENS Cachan / CNRS / Université Paris Saclay sylvie.pommier@universite-paris-saclay.fr jbellecave@gmail.com J. Mériaux SNECMA, Group SAFRAN jean.meriaux@snecma.fr ABSTRACT. Under the typical partial slip conditions under which fretting fatigue takes place, the amount of superficial damage is small. Therefore, the substantial reduction in fatigue life caused by fretting, when compared to plain fatigue, may well be more associated with the stress concentration and the stress gradient phenomena generated by the contact problem than to the superficial loss of material. In this setting, notch stress-based methodologies could, in principle, be applied to fretting in the medium/high cycle fatigue regime. The aim of this work was to investigate whether it is possible to design fretting and notch fatigue configurations, which are nominally identical in terms of damage measured by a multiaxial fatigue model. The methodology adopted to carry out this search considered a cylindrical on flat contact and a V-notch. Load and geometry dimensions of both configurations were adjusted in order to try to obtain the “same” decay of the Multiaxial Fatigue Index from the hot spot up to a critical distance. Positive results of such simulations can lead us to design an experimental program that can bring more firm conclusions on the use of pure stress-based approaches, which do not include the wear damage, in the modeling of fretting fatigue. KEYWORDS. Fretting fatigue; Notch fatigue; Stress gradient; Multiaxial fatigue. INTRODUCTION he term fretting denotes a small oscillatory movement between two solid surfaces in contact. In many industrial designs, one or both components of the assembly may also be subjected to a bulk fatigue load, the phenomenon is then called fretting fatigue. Fretting failure of components such as splines, the dovetail fixing between blade and disc in fans of aeroengines [1] and riveted skins of the aircraft fuselage [2] have become a major design concern. Indeed, experimental evidence has shown that the conjoint action of fretting and fatigue may produce strength reduction factors T J.A Araujo et al, Frattura ed Integrità Strutturale, 33 (2015) 427-433; DOI: 10.3221/IGF-ESIS.33.47 428 varying from 2 up to 10 compared to plain fatigue [3]. The most severe reduction of life is observed in partial slip conditions [4]. Indeed gross sliding can fastly remove small nucleated cracks before they can further propagate. Under fretting conditions, the loads involved generate a time varying non-proportional multiaxial stress field under the contact [5]. There is usually a high stress concentration on the contact interface but it is limited to a small region and hence decays fastly as one moves from the surface to the interior of the component [6]. In this setting, the use of non- local approaches, initially developed to estimate the fatigue endurance of notched specimens, have been extended to the fretting fatigue problem [7]–[9]. Araújo et al. [8] used the Theory of Critical Distances (TCD) in conjunction with the Modified Whöler Curve Method (MWCM) and showed that it was capable of estimating the results of fretting fatigue experiments showing a contact size effect with a good degree of accuracy (±20% error band). It is worthy noticing that notched components and mechanical assemblies under time varying loads present quite similar characteristics. In both configurations there are usually severe stress gradients and the state of stress becomes complex (multiaxial) as the analysis moves inside the material. The general aim of this work is to investigate whether it is possible to design fretting and notch fatigue configurations, which are nominally identical in terms of damage measured by a multiaxial fatigue index. Further, these designs must be such that the specimens can be machined within a good level of accuracy and with tight tolerances in both configurations (notch and fretting), turning the future experimental campaign reproducible and reliable. The advantage of proving that such experimental campaign can be generated is that it consists in a way to somehow isolate the role of the surface fretting wear in diminishing the material fatigue resistance of mechanical assemblies under partial slip. As both fatigue problems are equivalent in terms of fatigue loading over a material process zone, their resistance should be essentially similar, unless the influence of the small surface damage caused by the fretting wear, and which does not exist in the notch configuration, is greater than the authors expect it will have. Another quite important aspect that such experimental campaign could clear out is the fact that one can use the same fatigue modeling approach to design either industrial components containing geometrical discontinuities or mechanical couplings. This would avoid the need for lengthy and costly experimental programs considering complex geometries and specific test rigs to calibrate the fatigue material constants. MATERIAL he material considered in this study was a 7050-T7451 aluminium alloy. Material properties were taken from Araújo et al. [10] and are presented in Tab. 1. Fatigue strength for fully reversed (stress ratio, R=-1) and repeated loading (R=0) are quoted for 107 cycles. Young’s Modulus Poisson’s ratio Yield strength 1 0 , 1th RK  73.4 GPa 0.33 454 MPa 161 MPa 120 MPa 4.5 MPam Table 1: Material properties of 7050-T7451 Al alloy. STRESS TENSOR FOR ELASTIC CONTACT OF CYLINDERS UNDER PARTIAL SLIP AND IN THE PRESENCE OF A BULK FATIGUE LOAD he first step towards a solution for the subsurface stress field is to solve the contact problem itself, i.e., to find the magnitude and distribution of the surface tractions. In the present problem, the pad radius, R , the normal load per unit length, P , and the specimen thickness were defined so that each solid could be considered as an elastic half-space and the solution for the pressure distribution was Hertzian [11]. For pure fretting, the time varying tangential force Q (t) generates a shear traction described firstly by Cattaneo [12] and independently by Mindlin [13]. In partial slip condition, the contact is characterized by a central stick zone bordered by two sliding zones. In fretting fatigue, the effect of the bulk tension is to offset this stick zone. By superposing the normal and tangential contributions, it is possible to directly evaluate the resulting stress tensor using Muskhelishvili’s potential theory [14]. For synchronous and in phase load combinations, at the instants of maximum and minimum bulk/shear loads, the stress tensor turns out to be: T T J.A Araujo et al, Frattura ed Integrità Strutturale, 33 (2015) 427-433; DOI: 10.3221/IGF-ESIS.33.47 429 0 0 0 0 0 , , , , , ( ) n t t B y y y yx x x x e t tca a a a a a c c f f p p fp a fp p                                                                 (1) For any other time instant t, during loading or unloading, the appropriate analytical solution is stated as: 0 0 0 0 0 0 ( ) ,, , , , , ( ) ( )( ) ( ) 2 tn t t B x e t yy y y yx x x x e t c t c tc t tca a a a a a c c f f f p p fp a fp a fp p                                                                             (2) In the above equations, 0p is the peak pressure, c and e are the stick zone half width and its offset from the center of the contact at the instant of maximum or minimum shear load. At any other time instant, 'c and 'e correspond to the stick zone half width and its offset from the center of the contact. The superscripts n and t stand for the stress components due to the normal and the tangential load, respectively. Finally, B is the stress tensor associated with the bulk fatigue load. Plane strain conditions are assumed. Explicit expressions to compute c , e , 'c , 'e , n and t are given in a convenient form by Hills and Nowell [15]. STRESS TENSOR AROUND A V-NOTCH he elastic solution for the stress field of plane V-notched specimens can be derived using complex potential [14] or bi-harmonic potential functions [16]. Recently, the Williams’ crack solution [17] was re-addressed by Lazzarin and Tovo [18] and Filippi et al. [19]. Such approach was used in this work to evaluate the stresses along the bisector of a notch loaded in mode I. MULTIAXIAL CRITERION usmel and Lazzarin [20] observed that the multiaxial High-Cycle Fatigue behaviour of metallic materials could successfully be estimated by using a simple a vs , /n max a  relationship. The model has been described by Susmel and co-workers in a series of articles [20-22]. Briefly, the so-called Modified Wöhler Curve Method (MWCM) can be formalized as follows: ,( , ) ( , )n maxc c c ca a            (3) where a is the equivalent shear stress amplitude in the critical plane ( , ) c c  , ,n max is the maximum stress perpendicular to this plane, and the parameters  and  are material constants that can be obtained from two fatigue strengths generated under different loading conditions. For instance, these constants can be evaluted using the fatigue limits 1 and 0 , generated respectively under fully-reversed (R =-1) and under repeated (R = 0) uniaxial loading, as follows [21]: 1 0 2       01 2     (4) which gives 20.8 MPa  and 101.5MPa  for the 7050-T7451 Al alloy. A key aspect in using such model for non- proportional loadings, as it is the case for fretting fatigue, is the computation of a relative to a material plane. We propose here to use the maximum rectangular hull concept developed by Araújo et al. [23-25]. It consists of rotating T S J.A Araujo et al, Frattura ed Integrità Strutturale, 33 (2015) 427-433; DOI: 10.3221/IGF-ESIS.33.47 430 rectangular hulls engulfing the shear stress vector path in a material plane, and computing the square root of the sum of the squares of each rectangle’ half-sides. The greatest of these values is then defined as the shear stress amplitude, a . IDENTIFICATION OF THE CRITICAL DISTANCE he use of the MWCM in conjunction with the Theory of Critical Distances (TCD) is based on the assumption that all the physical processes leading to crack initiation are confined within the so-called structural volume. The size of this volume is assumed not to be dependent on either the stress concentration feature weakening the component or the complexity of the stress field damaging the fatigue process zone [26]. In the TCD fatigue endurance is assumed to occur when: 1 ( ) V F dV c V   (5) where  is stress tensor within a volume of material V around the stress raiser, c is a material parameter associated with its fatigue resistance, and  .F denotes an effective stress that appropriately characterizes the fatigue loading. Eq. (5) can be simplified by substituting the material volume by a line (LM) or by a single point (PM) at a certain distance from the hot spot. To obtain this length parameter, the critical distance, a third material parameter is needed, the threshold stress intensity factor range under R=-1, thK . Considering a material point at a distance l from the crack tip and on the bisector line of a sharp crack it is possible to compute the stress tensor at this position at any time instant of the loading history. If now one computes the a and ,n max variables of the multiaxial model from this state of stress, it is possible to find equations relating l, thK and 1 for the point and the line methods: 1 2 1 2 t M h P K l          (6) 1 2 2 th LM K l          (7) Using material properties reported in Tab. 1, one can use Eq. (6) to calculate 0.031PMl  mm. An alternative manner to interpret the critical distance is to associate it with the length of non-propagating cracks in crack like notched specimens under fatigue limit conditions. Therefore, from a phenomenological point of view, some authors such as Susmel [27] and Taylor [28] assume the critical distance as a material property. However, it is clear that it depends on the effective stress used to model the multiaxial problem. METHODOLOGY he first step to try to establish equivalent notch/fretting fatigue experiments is to define a reference configuration. In this case, a fretting fatigue data tested in the apparatus of the University of Brasília was considered as the starting point. This test was carried out using two cylindrical fretting pads, which were loaded against a flat dogbone specimen. Tab. 2 lists the loads, pad radius, the coefficient of friction, f , and the life for complete fracture of the especimen for this test. The cross section of the fretting specimen is schematically shown in Fig. 1 (13 x 13 mm). The characteristics of the fretting apparatus are well detailed by Martins et al. [29], hence it will not be described here. Now one can plot the multiaxial fatigue index eqS , defined in Eq. (8), from the trailing edge of the contact (hot spot for the fatigue index) along the distance y in the interior of the fretting specimen. The computation is carried up to a distance equivalent to the characteristic length provided by the line method, LMl . ,n maxeq a a S       (8) T T J.A Araujo et al, Frattura ed Integrità Strutturale, 33 (2015) 427-433; DOI: 10.3221/IGF-ESIS.33.47 431 Fig. 2 (a) shows the gradient of eqS with respect to y for this test and the limit for crack initiation defined by the threshold parameter 101.5MPa  for this alloy. Notice that, considering the point method as the critical distance, the index is higher than  (at the distance PMl ), and is then expected the test would break before 107 cycles. Indeed the test broke around 1.2 x 106 cycles. Figure 1: Scheme of the methodology proposed to find equivalent notch/fretting fatigue configurations. To carry on with the analysis the following step is to search for the notch configuration that will provide an equivalent decay for the multiaxial fatigue index over LMl . To do so, a minimization technique based on the difference between the indexes ,eq notchS and ,eq frettingS along the distance LMl can be applied (function J in Fig. 1). Notch radius, notch opening angle and the remote fatigue load are the variables of the minimization. The range of variation of such variables must be such that one can guarantee that the notch specimens can be machined and mainly its notch radius lies within a tight tolerance. In this setting, notch radius could vary between 0.1 and 1 mm (in steps of 0.05 mm), its opening angle between 30° and 90° (15° increments) and the notch depth was fixed in 5 mm. The limitation for the choice of remote bulk load range was defined by the operational capability of the servohydraulic machine. It can accurately apply and control load amplitudes varying from 1 to 100 kN (a 1kN increment was defined). Pad radius Peak pressure a ( 1)R   /Q fP f Frequency Life (cycles) 70mm 300 MPa 55 MPa 0.33 0.54 5Hz 1119774 Table 2: Parameters and observed life of the fretting fatigue test used to design an equivalent notch fatigue test. Minimization of the difference between the two gradients was achieved by using a 0.5mm notch radius specimen, with a notch opening angle of 60° and a fully reverse fatigue stress amplitude of 42.5 MPa. Now Fig. 2 (a) is again invoked to depict the result of the optimized notch configuration. It can be seen that the variation of the multiaxial fatigue index against distance from the notch root matches quite well the curve obtained for the fretting fatigue configuration. Fig. 2(b) depicts the points representing the shear stress amplitude, a , against the normalized maximum stress, ,n max a  , for both configurations. These points do not lay one over the other but they keep the same distance from the continuous threshold line dividing safe (under the line) and unsafe domains. In Fig. 2(c) an interesting behavior is observed. The shear stress amplitude and the maximum normal stress are traced against the distance for the notch and for the fretting configuration. While the amplitude of shear stress on the critical plane for the fretting problem is always higher than for the notch one along the distance, the opposite happens to the maximum normal stress. Therefore, although the decay of the multiaxial fatigue index is essentially the same for both configurations, the same can not be said about the individual variables that are present in the index computation. J.A Araujo et al, Frattura ed Integrità Strutturale, 33 (2015) 427-433; DOI: 10.3221/IGF-ESIS.33.47 432 Figure 2: (a) Variation of the multiaxial fatigue index against distance in the specimen (from the hot spot) for the fretting fatigue test and for the optimized notch configuration. (b) Points representing the shear stress amplitude, a , against the normalized maximum stress, , /n max a  , for both configurations. (c) Variation of a and , /n max a  against distance for fretting and notch problems. CONCLUSION t was shown that it is possible to design a series of fretting and notch fatigue tests for a 7050-T7451 Al alloy, which are nominally equivalent in terms of multiaxial fatigue conditions along a process zone. The notch geometry and fretting configuration found can be machined and reproduced within a given tolerance, and closed form solutions for the cyclic stress field under plane strain are available for both. It was pointed out however that equivalent fatigue index parameters does not mean equivalent shear stress amplitude and maximum normal stress. The design of this test campaign can constitute a key initiative to improve our knowledge on the influence of the surface damage on the life of components under fretting fatigue conditions. Further, should the results of these tests confirm the equivalence hypothesis between this two apparently distinct fatigue problems, a single unified stress based approach could be used to design them. Future work involves the conduction of the experimental tests. ACKNOWLEDGMENTS he supports provided by Safran/Snecna, by CNPq (contracts 310845/2013-0 and 309748/2013-5) and by FINATEC are gratefully acknowledged. REFERENCES [1] Ruiz, C., Chen, K. C., Life assessment of dovetail joints between blades and discs in aero-engines, In: Proc. Int. Conf. Fatigue, (1986). [2] Harish G., Farris, T. N., Shell Modeling of Fretting in Riveted Lap Joints, AIAA J., 36(6) (1998) 1087-1093. [3] Lindley, T. C., Fretting fatigue in engineering alloys, Int. J. Fatigue, 19(93) (1997) 39-49. [4] Madge, J. J., Leen, S. B., Shipway, P. H., The critical role of fretting wear in the analysis of fretting fatigue, Wear, 263(1-6) (2007) 542-551. [5] Navarro, C., Muñoz, S., Domínguez, J., On the use of multiaxial fatigue criteria for fretting fatigue life assessment, Int. J. Fatigue, 30(1) (2008) 32-44. [6] Nowell, D., Dini, D., Hills, D. A., Recent developments in the understanding of fretting fatigue, Eng. Fract. Mech., 73(2) (2006) 207-222. [7] Fouvry, S., Kapsa, P., Vincent, L. Multiaxial fatigue analysis of fretting contact taking into account the size effect, in ASTM Special Technical Publication, 1367 (2000) 167-182. I T J.A Araujo et al, Frattura ed Integrità Strutturale, 33 (2015) 427-433; DOI: 10.3221/IGF-ESIS.33.47 433 [8] 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(1) (2007) 95-107. [9] Proudhon, H., Fouvry, S., Yantio, G. R., Determination and prediction of the fretting crack initiation: introduction of the (P, Q, N) representation and definition of a variable process volume, Int. J. Fatigue, 28(7) (2006)707-713. [10] Araújo, J. A., Castro, F. C., A comparative analysis between multiaxial stress and ΔK-based short crack arrest models in fretting fatigue, Eng. Fract. Mech., 93 (2012) 34-47. [11] Hertz, R., Über die Berührung fester elastischer Körper, J. für die reine und Angew. Math., 92 (1981) 156-171. [12] Cattaneo, C., Sul contatto di due corpi elastici: distribuzione locale degli sforzi. Rendiconti dell’Accademia nazionale dei Lincei, Rend. dell’Accademia Naz. dei Lincei, 27 (1938) 342-348. [13] Mindlin, R. D., Compliance of elastic bodies in contact, J. Appl. Mech., 16 (1949) 259-268. [14] Muskhelishvili, Some Basic Problems of Mathematical Theory of Elasticity. Noordhoff, Groningen, (1953). [15] Hills, D. A., Mechanics of fretting fatigue, Wear, 175(1-2) (1994) 107-113. [16] Neuber, H., Theory of Notch Stresses, Kerbspannungslehre, (1958). [17] Williams, M. L., On the stress distribution at the base of a stationary crack, J. Appl. Mech., 24 (1957) 109-114. [18] Lazzarin, P., Tovo, R., A unified approach to the evaluation of linear elastic stress fields in the neighborhood of cracks and notches, Int. J. Fract., 78 (1996) 3-19. [19] Filippi, S., Lazzarin, P., Tovo, R., Developments of some explicit formulas useful to describe elastic stress fields ahead of notches in plates, Int. J. Solids Struct., 39(17) (2002) 4543-4565. [20] Susmel, L., Lazzarin, P., A bi-parametric Wöhler curve for high cycle multiaxial fatigue assessment, Fatigue & Fract. Eng. Mater. & Struct., 25(1) (2002) 63-78. [21] Susmel, L., Tovo, R., Lazzarin, P., The mean stress effect on the high-cycle fatigue strength from a multiaxial fatigue point of view, Int. J. Fatigue, 27(8) (2005) 928-943. [22] Fatemi, A., Socie, D. F., Critical plane approach to multiaxial fatigue damage including out-of-phase loading, Fatigue Fract. Eng. Mater. Struct., 11(3) (1988) 149-165. [23] 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(8) (2011) 1092-1100. [24] Castro, F. C., Araújo, J. A., Mamiya, E. N., Zouain, N., Remarks on multiaxial fatigue limit criteria based on prismatic hulls and ellipsoids, Int. J. Fatigue, 31(11-12) (2009) 1875-1881. [25] Mamiya, E. N., Araújo, J. A., Castro, F. C., Prismatic hull: A new measure of shear stress amplitude in multiaxial high cycle fatigue, Int. J. Fatigue, 31(7) (2009) 1144-1153. [26] Susmel, L., A unifying approach to estimate the high-cycle fatigue strength of notched components subjected to both uniaxial and multiaxial cyclic loadings, Fatigue Fract. Eng. Mater. Struct., 27(5) (2004) 391-411. [27] Susmel, L., The theory of critical distances: a review of its applications in fatigue, Eng. Fract. Mech., 75(7) (2008) 1706-1724. [28] Taylor, D., A mechanistic approach to critical-distance methods in notch fatigue, Fatigue Fract. Eng. Mater. Struct., 24(4) (2001) 215-224. [29] Martins, L. H., Rossino, L. S., Bose Filho, W. W., Araújo, J. A., Detailed design of fretting fatigue apparatus and tests on 7050-T7451~Al alloy, Tribol. - Mater. Surfaces Interfaces, 2(1) (2008) 39-49.