This is a repository copy of Zinc dialkyl dithiophosphate antiwear tribofilm and its effect on the topography evolution of surfaces: A numerical and experimental study.
White Rose Research Online URL for this paper: http://eprints.whiterose.ac.uk/100856/
Version: Accepted Version
Article: Ghanbarzadeh, A orcid.org/0000-0001-5058-4540, Piras, E, Nedelcu, I et al. (5 more authors) (2016) Zinc dialkyl dithiophosphate antiwear tribofilm and its effect on the topography evolution of surfaces: A numerical and experimental study. Wear, 362-36. pp. 186-198. ISSN 0043-1648 https://doi.org/10.1016/j.wear.2016.06.004
© 2016, Elsevier. Licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International http://creativecommons.org/licenses/by-nc-nd/4.0/
Reuse Unless indicated otherwise, fulltext items are protected by copyright with all rights reserved. The copyright exception in section 29 of the Copyright, Designs and Patents Act 1988 allows the making of a single copy solely for the purpose of non-commercial research or private study within the limits of fair dealing. The publisher or other rights-holder may allow further reproduction and re-use of this version - refer to the White Rose Research Online record for this item. Where records identify the publisher as the copyright holder, users can verify any specific terms of use on the publisher’s website.
Takedown If you consider content in White Rose Research Online to be in breach of UK law, please notify us by emailing [email protected] including the URL of the record and the reason for the withdrawal request.
[email protected] https://eprints.whiterose.ac.uk/ 1 Zinc Dialkyl Dithiophosphate Antiwear Tribofilm and its Effect on the Topography 2 Evolution of Surfaces: A Numerical and Experimental Study
3 Ali Ghanbarzadeh1, Elio Piras1,2, Ileana Nedelcu2, Victor Brizmer2, Mark C. T. Wilson1, 4 Ardian Morina1, Duncan Dowson1, Anne Neville1
5 1 Institute of Functional Surfaces, School of Mechanical Engineering, University of Leeds, 6 Leeds, UK 7 2 SKF Engineering and Research Centre, the Netherlands
8 Corresponding authors: 9 E-mail address: *[email protected] (Ali Ghanbarzadeh) 10 *[email protected] (Elio Piras) 11 Abstract 12 A modelling framework has recently been developed which considers tribochemistry in 13 deterministic contact mechanics simulations in boundary lubrication. One of the capabilities of 14 the model is predicting the evolution of surface roughness with respect to the effect of 15 tribochemistry. The surface roughness affects the behaviour of tribologically loaded contacts 16 and is therefore of great importance for designers of machine elements in order to predict 17 various surface damage modes (e.g. micropitting or scuffing) and to dessign more efficient 18 tribosystems. The contact model considers plastic deformation of the surfaces and employs a 19 modified localized version of Archard’s wear equation at the asperity scale that accounts for 20 the thickness of the tribofilm. The evolution of surface topography was calculated based on the 21 model for a rolling/sliding contact and the predictions were validated against experimental 22 results. The experiments were carried out using a Micropitting Rig (MPR) and the topography 23 measurements were conducted using White Light Interferometry. Numerically, it is shown that 24 growth of the ZDDP tribofilm on the contacting asperities affects the topography evolution of 25 the surfaces. Scanning Electron Microscopy (SEM) and X-ray Photoelectron Spectroscopy 26 (XPS) have been employed to confirm experimentally the presence of the tribofilm and its 27 chemistry. The effects of the contact load and surface hardnesses on the evolution of surface 28 topography have also been examined in the present work.
29 Key words: Surface Roughness, Boundary Lubrication, Tribochemistry, Wear, ZDDP
1
30 1. Introduction 31 Running-in is an important term used in the field of tribology. Due to the complexity and 32 diversity of the phenomena occurring in this period various definitions of the term can be found 33 in the literature (1). As described by Blau (2), running-in is a combination of processes that 34 occur prior to the steady-state when two surfaces are brought together under load and with 35 relative motion, and this period is characterized by changes in friction, wear and physical and 36 chemical properties of surfaces. During the running-in period, surface micro topography is 37 subjected to various changes. In boundary and mixed lubrication conditions, the height of the 38 asperities of rough surfaces normally decrease (3-6). However, in the case of very smooth 39 surfaces, an increase in the roughness value is observed (7, 8). During the process of change in 40 the roughness of the surfaces, load carrying capacity is increased due to the gradual 41 development of asperity-level conformity. The increase in the conformity of the surfaces is 42 significant as peaks and valleys of the surfaces correspond to each other, and the overall 43 performance of the system is improved (9-11).
44 Predicting changes in the topography of contacting bodies is important for designers to be able 45 to predict the mechanical and chemical behaviour of surfaces in loaded tribological systems. 46 The roughness of operating surfaces influences the efficiency of mechanical parts. In the 47 design of machine elements and selection of materials, the film thickness parameter, known as 48 the ratio, which is a representation of the severity of the contact, is important (12). Its value 49 is inversely proportional to the composite roughness of the two surfaces in contact. It is also 50 widely reported that fatigue life of bearing components is dependent on their functional 51 surfaces’ characteristics, such as roughness. Optimization of the surface roughness can help in 52 increasing the lifetime of bearings based on their applications (13, 14). Surface roughness can 53 enhance stress concentrations that can lead to surface-initiated rolling contact fatigue (15). 54 Therefore, it is important to be able to predict the surface topography changes in a 55 tribologically-loaded system.
56 The changes in topography can be either due to plastic deformation of the surface asperities or 57 due to removal, loss or damage to the material, which is known as wear. Evaluating wear in 58 boundary lubrication has been the subject of many studies. There are almost 300 equations for 59 wear/friction in the literature which are for different conditions and material pairs but none of 60 them can fully represent the physics of the problem and offer a universal prediction (16, 17). 61 Some examples of these models are the Suh delamination theory of wear (18), the Rabinowicz 62 model for abrasive wear (19) and the Archard wear equation (20, 21). Wear occurs by different
2
63 interfacial mechanisms and all these mechanisms can contribute to changes in the topography. 64 It has been widely reported that third body abrasive particles play an important role in changes 65 in the topography of surfaces. There are several parameters that govern the wear behaviour in 66 this situation such as wear debris particle size or shape, configuration of the contact and contact 67 severity etc. (22-24). It was reported by Godet (25) that a comprehensive mechanical view of 68 wear should consider the third body abrasive particles and their effect on wear and topography 69 changes. A study of abrasive wear under three-body conditions was carried out by Rabinowicz 70 et al (26). They proposed a simple mathematical model for third body abrasive wear rate and 71 showed that the wear rate in this situation is about ten times less than two-body abrasive wear. 72 It was reported by Williams et al. (27) that lubricant is used to drag the wear debris inside the 73 interface and the abrasive wear action then depends on the particle size, its shape and the 74 hardness of the materials. They reported that a critical ratio of particle size and film thickness 75 can define the mode of surface damage. Despite the importance of a three-body abrasive wear 76 mechanism there is no comprehensive mechanistic model to describe such a complicated 77 mechanism. In the mild wear regime in lubricated contacts the effect of third body abrasive is 78 often assumed to be insignificant.
79 Most of the work in the literature is based on using the well-known Archard wear equation to 80 evaluate wear in both dry and lubricated contacts. Olofsson (28-30) used Archard’s wear 81 equation to evaluate wear in bearing applications and observed the same behaviour between 82 model and experiments. Flodin (31) showed that Archard’s wear equation is good enough to 83 predict wear in spur helical gears application. Andersson et al (32) tested and reviewed different 84 wear models and reported that Archard’s wear model can predict wear of lubricated and 85 unlubricated contacts and is able to predict the surface topography both in macro and micro- 86 scales. They tested their generalized Archard’s wear model for random rough surface contact 87 (33). The Archard wear equation was widely used in numerical studies in order to predict the 88 wear and topography at different scales (34-47).
89 Hegadekatte et al. (39) developed a multi-time-scale model for wear prediction. They used 90 commercial codes to determine the contact pressure and deformations and then used Archard’s 91 wear equation to calculate wear. Andersson et al. (47) have employed the Archard wear 92 equation to predict wear in a reciprocating ball-on-disc experiment. They used a wear model 93 and implemented Fast Fourier Transforms (FFT) based contact mechanics simulations to 94 calculate contact pressure and deformations. However, in all these implementations of 95 Archard’s wear equation in numerical models, which resulted in reasonably good agreement
3
96 with experimental results, the effect of lubrication and lubricant properties was neglected. 97 Recently, there have been some attempts to consider the lubrication effects in boundary 98 lubrication modelling that could affect modifications to Archard’s wear equation.
99 Bosman & Schipper (48) proposed a numerical model for mild wear prediction in boundary 100 lubricated systems. They assumed that the main mechanisms that protect the boundary 101 lubricated system are the chemically-reacted layers and when these layers are worn off, the 102 system will restore the balance and the substrate will react with the oil to re-establish the 103 tribofilm. They also proposed a transition from mild wear to more severe wear by making a 104 complete wear map. In another recent work by Andersson et al. (49), contact mechanics of 105 rough surfaces was used to develop a chemo-mechanical model for boundary lubrication. They 106 used an Arrhenius-type thermodynamic equation to develop a mathematical model for 107 formation of the tribofilm on the contacting asperities. They have also employed the 108 mechanical properties of the antiwear tribofilm and used Archard’s wear equation to predict 109 wear of the surfaces. The coefficient of wear was assumed to be the same for the areas where 110 the tribofilm is formed with the areas without the tribofilm.
111 Recent work by Morales Espejel et al. (50) used a mixed lubrication model to predict the 112 surface roughness evolution of contacting bodies by using a local form of Archard’s wear 113 equation, and the model results show good agreement with experimental data. They used a 114 spatially and time-dependent coefficient of wear that accounts for lubricated and unlubricated 115 parts of the contact. The same modelling framework was used in other works of those authors 116 to predict wear and micropitting (51).
117 A range of experimental work has investigated changes in surface roughness during 118 tribological contacts. Karpinska (7) studied the evolution of surface roughness over time for 119 both base oil and base oil with ZDDP. She also studied the wear of surfaces at different instants 120 during running-in. It was suggested that a ZDDP tribofilm significantly affects the 121 topographical changes of surfaces during running-in. Blau et al. (1) stated that friction and wear 122 in running-in are time-dependent and related to the nature of energy dissipation in the contacts; 123 they are governed by a combination of different mechanical and chemical processes. They 124 showed that roughness evolution of contacting surfaces might have different patterns for both 125 surfaces, depending on several parameters.
126
4
127 Despite the importance and the attempts in the literature to monitor and predict the roughness 128 evolution of surfaces, there is no reported work that addresses the effect of tribochemistry. 129 However, a modelling framework has recently been developed by the authors (52),(53) that is 130 capable of predicting changes in surface topography under boundary lubrication conditions, 131 taking into account the simultaneous dynamics of an anti-wear tribofilm. The present paper 132 therefore seeks to test and exploit this model to explore the effect of a ZDDP tribofilm on the 133 evolution of surface topography. The topography evolution of both contacting surfaces is 134 predicted, taking into account not only the effects of plastic deformation and mild wear but 135 also the coupled development and influence of a ZDDP antiwear tribofilm.
136 Experimental results from a Micropitting Rig (identical to the one used in Ref (50)) are used 137 to validate the model in terms of general prediction of topography changes and growth of the 138 ZDDP tribofilm on the contacting asperities. The numerical model is described briefly in 139 Section 2, while the experimental set-up that the numerical model is adapted to is explained in 140 Section 3. The numerical results based on the model are then reported and discussed in Section 141 4, where special attention is given to the different parameters in the model that affect the surface 142 topography evolution. The importance of the growth of ZDDP tribofilm in changing the 143 topography of surfaces in the model is shown in that section. The experimental results from the 144 MPR and surface roughness measurements are reported in Section 5, where the thickness of 145 the tribofilm and the evolution of surface roughness are compared with the numerical results 146 of Section 4. Using the validated model, the effects of two important physical parameters – the 147 hardness and the load – are studied numerically in Sections 6 and 7.
148 2. Numerical model 149 The numerical model used in this work is the one reported by the authors in Ref (52). The 150 model is adaptable to tribosystems with different configurations which makes it possible to 151 investigate different problems. The model consists of three important parts:
152 (i) A contact mechanics code for rough surfaces assuming an elastic-perfectly plastic 153 material response; 154 (ii) A semi-analytical tribofilm growth model which includes both tribofilm formation 155 and partial removal; and 156 (iii) A modified Archard’s wear equation which accounts for the local thickness of the 157 ZDDP tribofilm.
5
158 A brief description of the model is given below; further details on parts (i) and (ii) can be found 159 in Ref (52) and an expanded discussion on part (iii) and the wear model validation is reported 160 in (53). The whole numerical model algorithm is shown schematically in Figure 1.
Initial roughness patterns, Contact inputs, End of simulation Material properties
Yes
Apply contact mechanics and No Is time = End? calculate contact pressures
(time increment)
Calculate elastic and plastic Calculate the local wear and deformations modify the surfaces
Apply plastic deformation Modify the Archard wear on surfaces equation with respect to (permanent surface the tribofilm (equation 3) modifications)
Contact temperature calculation, Modify the surfaces with Calculate local tribofilm tribofilm on top growth on surfaces
161 Figure 1. Flow chart for the numerical procedure.
162 Digitized surfaces are important inputs for the deterministic contact mechanics simulations. 163 Surfaces are generated in this work based on the model developed by Tonder (54) which is 164 based on the digital filtering of a Gaussian input sequence of numbers. This method generates 165 surfaces with desired roughness and asperity lateral size. Using artificial surface topographies 166 instead of actual surface scans does have some drawbacks, but also several advantages. For
6
167 practical reasons, it is good to use the measured surface topography to study different cases for 168 bearings, gears and other machine element parts. In addition, some important real 169 characteristics of the engineering surface might be missed if using artificial generated surfaces. 170 However, artificial surfaces can be used for scientific studies in different lubrication regimes, 171 and they are particularly good for performing parametric studies independently of specific 172 experimental datasets. It is possible to see the effect of different topographic properties on the 173 tribofilm formation, plastic deformation, wear and topography evolution. The method is simple 174 to manipulate and, since this work seeks to evaluate the performance of the model as a general 175 tool for exploring topography evolution, is appropriate.
176 The contact mechanics model is based on the complementary potential energy concept using 177 the model of Tian et al. (55). An elastic-perfectly plastic approach is then employed with the 178 hardness of the material to be the criterion for the plastic flow. This assumption neglects the 179 work hardening behaviour of the asperities as they deform. This method gives realistic contact 180 pressure estimation but the deformations due to the pressure can be different from real values 181 as the influence coefficients of an elastic material are used, neglecting the non-linear behaviour 182 of the yielding subsurface material. Despite its shortcomings, this method has been widely used 183 in the literature for modelling pressures and surface deformations in the elastic-plastic contacts. 184 The model is based on the contact model by Sahlin et al (56). The surfaces move relative to 185 each other and the slide-to-roll ratio determines the speed of the movement of surfaces. 186 Therefore an asperity of one surface can contact with a number of asperities on its way. The 187 movement of surfaces is periodic and is carried out by shifting the elements of matrices 188 containing the values for the asperity heights.
189 The thermal model used to calculate the flash temperature is based on the Blok theory (57). It 190 is obtained from calculation of the frictional heating. It should be noted that only the maximum 191 temperature is important in this work and using Blok’s theory seems reasonable. The 192 formulation used in this work is based on the equations reported in Kennedy and Tian (58, 59).
193 A tribochemical model was developed in the previous work of the authors that can capture the 194 growth of the tribofilm on the asperity level. The tribofilm growth is taken to be a combination 195 of the formation and partial removal at the same time (52). The formation of the tribofilm is 196 assumed to be due to tribochemical reactions and follows the reaction kinetics based on the 197 non-equilibrium thermodynamics of interfaces. The formation model is combined with a
7
198 phenomenological term that accounts for the simultaneous partial removal of the tribofilm, 199 giving the net development of the tribofilm thickness as a function of time:
200 201 where and are the Boltzmann and the Planck constants respectively, is the flash 202 temperature and and are constants accounting for the continuous partial removal of the 203 tribofilm. Hence the tribofilm is modelled as a dynamic system such that even when an 204 equilibrium thickness is achieved, formation and partial removal continue but balance each 205 other to maintain that thickness. Note that t in equation (1) refers to a local time for each point 206 in the domain, and starts increasing at a given point once an asperity contact occurs there. The 207 term captures the effect of rubbing in inducing the tribochemical reactions. As reported 208 recently (60), mechanical activation plays an important role in the growth behaviour of the 209 tribofilm on a single asperity. Gosvami et al. (60) showed that the rate of tribochemical reaction 210 is highly dependent on the pressure applied on a single asperity. This is in line with the model 211 presented above, where the role of mechanical activation is represented by the term . A 212 detailed discussion can be found in the previous work (52).
213 At this point it is important to distinguish between the wear of the tribofilm, which is captured 214 by the ‘partial removal’ term in Equation (1), and the wear of the substrate itself. Henceforth 215 in this paper the term ‘wear’ will be used to refer to the (mild) wear of the substrate underneath 216 the tribofilm – i.e. the wear of the substrate in the presence (and, initially, absence) of the 217 tribofilm. Studies of ZDDP tribofilms on steel show that the tribofilm contains substrate atoms 218 at a concentration that decreases towards the top of the tribofilm. Hence material from the 219 substrate is consumed in forming (and maintaining) the tribofilm, and therefore if part of the 220 tribofilm is removed due to the contact, this corresponds to an effective removal of material 221 from the substrate. This principle is the basis of the mild wear model of Bosman & Schipper 222 (Ref (48)), who linked the rate of substrate wear to the rate of tribofilm removal by considering 223 the volumetric percentage of iron as a function of depth in the tribofilm. The functional form 224 was determined from X-ray photoelectron spectroscopy analysis.
225 In the present work, the link between the substrate wear and the tribofilm takes a different form. 226 The details are explained elsewhere (53) but, in brief, the wear model is a modified version of 227 Archard’s wear equation in which the wear coefficient is related to the local tribofilm thickness. 228 The local wear depth of each point at the surface is given by:
8
229 (2) 230 in which , , , , and are the material hardness, dimensionless Archard’s wear 231 coefficient, local contact pressure, sliding speed, and time step respectively. All the parameters 232 in Equation (2) except are calculated in the contact mechanics simulation. The decrease 233 in the concentration of substrate atoms within the tribofilm as the distance from the 234 substrate/tribofilm interface increases supports the fact that less wear of the substrate occurs if 235 a thicker tribofilm exists. Hence it is assumed that the coefficient of wear is at its maximum for 236 steel-steel contact (i.e. when no tribofilm is present) and at its minimum when the tribofilm has 237 its maximum thickness. Assuming, in addition, a linear variation with tribofilm thickness h, 238 the coefficient of wear is given by:
239 (3) 240 where is the coefficient of wear for a substrate covered by a tribofilm with thickness .
241 and are the coefficients of wear for steel and for the maximum ZDDP tribofilm 242 thickness respectively, and is the maximum tribofilm thickness. The values of and 243 are determined from a single calibration experiment as described in Section 4. Figure 1 244 shows a flow chart of the numerical model.
245 3. Experimental setup 246 The tests were carried out using a Micropitting rig (MPR) which is shown schematically in 247 Figure 2.
248
249 Figure 2. The load unit of the Micropitting Rig (MPR)
250 The principle part of the Micropitting rig is the load unit which includes a spherical roller 251 12mm in diameter (taken from a spherical roller bearing), which comes into contact with three
9
252 larger counterbodies (which are in fact inner rings of a cylindrical roller bearing). The roller 253 and the rings can be ground and/or honed and finished to any desired roughness on the surface. 254 With this configuration, the roller accumulates loading cycles 13.5 times faster than any of the 255 three rings during a test. The oil is carried into the contact by the two bottom rings, creating 256 the conditions of fully lubricated contact – depending, of course, on the desired lubrication 257 regime.
258 The maximum load that can be applied on the machine is 1250 N, corresponding to a maximum 259 Hertzian pressure of 2.75 GPa, which is high enough for bearing studies. The temperature can 260 be controlled up to 135ºC and the maximum tangential speed is 4 m/s.
261 The roller and rings are driven by two independent motors, which means that a controlled slide- 262 to-roll ratio of ±200% can be reached. The size of the Hertzian contact varies with the load and 263 the transverse radius of the roller, but the typical values are around 0.244 1.016 mm (in the 264 rolling and transverse direction, respectively) corresponding to the maximum Hertzian contact 265 pressure of 1.5 GPa. The operating conditions used in the present work are listed in Table 1. 266 The oil temperature, load, speed and slide-to roll ratio were maintained constant during the 267 experiments. The roughness was changed (on the rings only), to simulate different lubrication 268 conditions.
269 Over a given experiment duration, the number of rolling cycles experienced is different for the 270 roller and the rings; in fact each cycle of the rings corresponds to 13.5 cycles of the roller.
271 The lubricant used was a synthetic model oil (poly-alpha-olephine, PAO) mixed with 1% 272 weight of primary zinc dialkyl dithiophosphate (ZDDP). The properties of the oil and additive 273 are listed in Table 2.
274
275
276
277
278
279
10
280 Table 1. List of experimental test conditions for the Micropitting rig Temperature 90ºC
Entrainment speed ( ) 1 m/s