Next Article in Journal
A Time-Frequency Based Approach for Acoustic Emission Assessment of Sliding Wear
Previous Article in Journal
Methylenealkane-Based Low-Viscosity Ester Oils: Synthesis and Outlook
Previous Article in Special Issue
Misalignment-Induced Micro-Elastohydrodynamic Lubrication in Rotary Lip Seals
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Elastohydrodynamic Lubrication

by
James A. Greenwood
Department of Engineering, University of Cambridge, Cambridge CB2 IPZ, UK
Lubricants 2020, 8(5), 51; https://doi.org/10.3390/lubricants8050051
Submission received: 30 December 2019 / Revised: 19 April 2020 / Accepted: 26 April 2020 / Published: 6 May 2020
(This article belongs to the Special Issue Elastohydrodynamic (EHD) Lubrication)

Abstract

:
The development of EHL theory from its tentative beginnings is outlined, with an account of how Ertel explained its relation to Hertz contact theory. The problems caused by the failure of the early numerical analysts to understand that the film thickness depends on only two variables are emphasised, and answers of the form H = F ( P , S ) given. Early methods of measuring the film thickness are described, but these became archaic with the development of optical EHL. The behaviour of surface roughness as it passes through the high pressure region and suffers elastic deformation is described, and the implication for the traditional Λ -ratio noted. In contrast, the understanding of traction is far from satisfactory. The oil in the high pressure region must become non-Newtonian: the early explanation that the viscosity reduction is the effect of temperature proved inadequate. There must be some form of shear thinning (perhaps according to the Eyring theory), but also a limiting shear stress under which the lubricant shears as an elastic solid. It seems that detailed, and difficult, measurements of the high pressure, high shear-rate behaviour of individual oils are needed before traction curves can be predicted.

1. Introduction

Gear engineers have of course long known that when gears are removed from service after long use, the teeth, despite having been wiped (for of course only at the pitch point do they roll!) across each other for millions of passes, still show the original grinding marks. There must have been an oil film protecting the surface and preventing contact. However, classical lubrication theory when applied to very localised contacts, such as those between gear teeth, predicts film thicknesses way below the surface roughness. However, the classical theory satisfactorily explains the behaviour of thrust and journal bearings, predicting film thicknesses adequately greater than the surface roughness of the components. Why the difference?
Classical lubrication theory naturally assumes that a lubricating oil is a Newtonian fluid with a fixed viscosity, and that the bearings are elastic solids, which might bend under oil pressure or thermally, but could never deform locally. Surely that must be the case??

1.1. Classical Theory

Early investigations studied the simple two-dimensional problem: that is, where the lubricant flow is in one direction, and so should be governed by the one-dimensional Reynolds equation [1], in its simple integrated form d p d x = 6 η   u 1 h h h 3 , where p is the pressure, η the viscosity and u 1 the speed of the moving surface. h arises as an integration constant, but then is recognisable as the film thickness at the location of the pressure maximum. While in a thrust pad or a journal bearing normally only one of the two bearing surfaces will be moving, in gear contacts invariably both teeth are moving. Until we look at friction and heat production only the sum, or average, of the two speeds matters. The average (or “rolling”) speed u ¯ ( 1 / 2 ) ( u 1 + u 2 ) will be used throughout this paper: accordingly the Reynolds equation must be written d p d x = 12 η   u ¯ h h h 3 . We note that h then has an additional significance, for u ¯   h is the lubricant flow rate through the bearing.
Figure 1 (drawn for an inlet shape we shall meet later) shows the flow pattern in the inlet: the moving surfaces carry in more lubricant than can pass through the gap, and the surplus must be rejected and pushed back; and of course it is this surplus oil supply which causes the pressure build-up.
In a concentrated “contact” between a circular roller (radius R ) and a plane, it is easy to show that locally the gap is approximately h = h 0 + x 2 / 2 R (the parabolic approximation) [2]: the pressure is then readily found (analytically) to be as in Figure 2; the variations depending on the value of the ratio h / h 0 .
The two curves of interest are the one through the origin ( h * / h 0   = 4 / 3 ), which continues as the reflection of the curve from to 0 and becomes zero again at + , (known as the Sommerfeld solution), and the one becoming zero at the pressure minimum ( h * / h 0   = 1.2257 ), the Reynolds solution. The Sommerfeld solution is useless here, for the symmetry means that no load is carried: the Reynolds solution supports a load (per unit width) of W = 4.895   η   u ¯ R / h 0 (the “half-Sommerfeld” solution—ignoring the negative pressures—supports a load W = 4   η   U R / h 0 ; but violates flow continuity); in other words, the minimum film thickness h 0 will be 4.895   η   u ¯ R / W (the coefficient can be shown to be 6/H*). As an example of applying this to a gear tooth contact, take two steel teeth with relative radius of curvature R = 0.02   m with a tooth load W = 5 × 10 5   N / m . If this were a Hertzian line contact with no lubricant, the maximum Hertz pressure, given by p 0 = { W E / π   R } , with E = 115   GN / m 2 would be an acceptable 0.96   GN / m 2 . However, with a speed u ¯ = 2   m / s and oil of viscosity η = 0.07   Ns / m 2 , we find a film thickness of 0.03 µm, less than a realistic surface finish.

1.2. Digression: Factors of 2: Sums and Averages

We noted above that the Reynolds equation could be written using either the average speed u ¯ (or U ), or the sum speed u Σ , usually without warning. Either may be referred to as the “entraining“ speed, and so provide an opportunity of an error of a factor 2!
In the above Hertz calculation, we used the standard notation in Contact Mechanics [3] for the contact modulus 1 E 1 ν 1 2 E 1 + 1 ν 2 2 E 2 , derived from the fundamental Boussinesq result that a point load P on an elastic half space gives a displacement w ( r ) = P · ( 1 ν 2 ) π   E   r , so that the separation between two solids from the mutual effect of a force is the sum of two such terms. The combination E   /   ( 1 ν 2 ) is the plane strain modulus, conventionally denoted by E . A valuable consequence is that a solution of a contact problem between a rigid surface and an elastic solid, expressed in terms of E , can immediately be converted into the solution for two elastic solids by simply replacing E by E . Unfortunately, in EHL [4] the basic elastic quantity is defined by 1 E 1 2 [ 1 ν 1 2 E 1 + 1 ν 2 2 E 2 ] , a confusing double use of the symbol E . Thus, the maximum pressure in a Hertzian line contact with a rigid indenter will be p 0 = W E / π   R in contact mechanics, but p 0 = W E / 2 π   R in EHL!
We note with some surprise that although the local gap between two parallel rollers of radii R 1 and R 2 is z = x 2 2 R 1 + x 2 2 R 2 , giving h = h 0 + x 2 2 R 1 + x 2 2 R 2 in a concentrated contact with a minimum gap h 0 , no-one seems to use the average curvature 1 2 [ 1 R 1 + 1 R 2 ] ; happily, all define the effective radius of curvature by 1 R 1 R 1 + 1 R 2 and write h = h 0 + x 2 2 R . (However, be careful: for the contact between two crossed rollers of radii R 1 and R 2 , consult Johnson [3].

2. First Attempts

In the solution above, the lubricating oil is, of course, taken to be a Newtonian fluid with a fixed viscosity. We all know that the viscosity of oil decreases seriously with temperature, and this must be considered in the design of ordinary bearings: any effect of pressures, of the order of a few MN/m2 is (correctly) ignored. However, as noted above, in a dry gear contact there would be pressures of order 1 GN/m2; what would such a pressure due to the oil viscosity? There is good evidence that the Barus equation η = η 0 exp ( α   p ) is at least a reasonable approximation [4] (chapter 4), with a pressure–viscosity index (for a typical mineral oil at 30 °C, but decreasing as the temperature rises) of α = 20 × 10 9 m2/N, so at 1 GN/m2 the viscosity would increase by a factor of e 20 = 5 × 10 8 !
The Martin analysis [2] can easily be repeated using this viscosity law: the Reynolds equation becomes d p d x = 12 η 0 e α   p   U h h h 3 or e α   p   d p d x = 12 η 0 U h h h 3 , and the solution of the equation d q d x = 12 η 0 U h h h 3 is shown in Figure 1. Let e α   p   d p d x = d q d x , and choose q = 0 when p = 0 ; then q = 0 p e α   p d p = ( 1 / α ) · [ 1 e α   p ] . Thus, p = ( 1 / α ) ln ( 1 α   q ) , so that the q-plot is readily converted into a p-plot as in Figure 3.
The solid (“isoviscous”) curve is taken from Figure 2; the others are found from the q p equation. The pressures are substantially increased, and if the greatest pressure in the isoviscous solution ( q max = 2.1509 ·   η     U R h o 3 / 2   ) reaches 1/α, the real pressure becomes infinite. The minimum film thickness will then be h 0 = [ 2.1509 α η 0 U R ] 2 / 3 = 1.6663   ( α η 0 U ) 2 / 3 R 1 / 3 (Gatcombe [5]), while h = 2.04246   ( α η 0 U ) 2 / 3 R 1 / 3 , and this is in fact quite close to what we now know to be correct. However, the very high pressures are local, and although the load is increased from the isoviscous load, it is only by a factor of 2.3 (i.e., W = ( 2.3 × 4.985 )   η 0 U R / h 0 ), so this seems unable to explain high load behaviour (if only Blok, while integrating the logarithmic pressures, had taken them seriously, and wondered what the response of the material to an infinite pressure spike would be, we might have had a Dutch EHL theory in 1952) [6].

Kapitza’s Solution

It is convenient here to divert temporarily from the discussion of line contacts to include another rigid contact solution. The three-dimensional Reynolds equation for flow in the x-direction is x h 3 12 η p x + y h 3 12 η p y = U h x .
For a point contact the film thickness is approximately h = h 0 + x 2 / 2 R x + y 2 / 2 R y . Kapitza [7] noticed that p ( x , y ) = C x / h 2   is a solution of the Reynolds equation, provided that C = 4 η   U / ( 1 + 2 R x / 3 R y ) . Thus, for a sphere ( R x = R y ) we have C = ( 12 / 5 )   η U   , while for a roller ( R y = ) we get C = 4 η   U : and we recognise that p ( x ) = 4 η   U   x / h 2 is just the Sommerfeld solution for line contact that we noted before. The maximum pressure is of course on the centreline y = 0 , at x = ( 2 / 3 )   h 0 R x ; [ h = ( 4 / 3 ) h 0 ] and equals p max = η   U R x 1 + 2 R x / 3 R y 3 3 2 2 h 0 3 / 2 .
A convenient film thickness parameter is the ratio h   /   [ ( α   η 0   U ) 2 / 3 R x 1 / 3 ]   , which is denoted by H K commemorating Kapitza [7] (why Kapitza? For noting the simple hydrodynamic solution p = K   x / h 2 , but particularly for providing the first point contact solution). With this notation, the Gatcombe solution [5] becomes H K   = 2.042 .
It can be seen that the maximum pressure along the centreline of a circular contact is reduced from that of a band contact by the factor ( 1 + 2 R x 3 R y ) 1 , and this is obviously attributable to “side-leakage”; not all the fluid which enters the inlet continues through the contact. Of course the exact factor is specific to the Kapitza analysis, but the same factor is a useful general guide.

3. The “Grubin” Theory and Its Extension

It should be emphasised that the “Grubin” theory is not by Grubin [8], but by Ertel [9]. However, Ertel defected from the USSR, so in the political climate of the time neither his name nor his work could be mentioned. Grubin, the head of the department on which Ertel had done the work, correctly decided that the theory was too important not to be published.
Ertel saw that under the loads occurring in gear contacts, the metal parts must deform elastically to much the same way as in a Hertzian contact. In addition, indeed, at the pressure levels needed, the oil viscosity would be so enormous that the product η   · ( h h ) in the Reynolds equation could only remain of an acceptable magnitude if ( h h ) is vanishingly small: the oil must flow along a substantial parallel channel of width h h —the geometry of a Hertzian contact, but with the surfaces separated by a constant film thickness. It must be in the inlet to this parallel channel that the oil pressure builds up to Hertzian levels.
The shape of the gap at the edge of a Hertz contact is well known: for a band contact of width 2 b between two rollers with relative radius of curvature R , the gap is z = ( b 2 / 4 R )   a   ( x / b ) where the shape function is a   ( ξ ) 2 [   | ξ | ξ 2 1   cosh 1 | ξ |   ] . So we need to integrate the Reynolds equation from x = to x = b along a channel of width h = h +   ( b 2 / 4 R )   a   ( x / b ) . It is convenient to do this for a constant viscosity oil and then convert the reduced pressures q ( x ) so found into real pressures using the transformation p = ( 1 / α ) ln ( 1 α   q ) introduced above. Then, by equating the reduced pressure q ( b ) at the end of the inlet to 1 / α , we obtain a film thickness formula; the corresponding load is taken to be the Hertz load for this contact. Thus, for a Hertz load W we have b 2 = 8 R W π   E and 1 α = 12 η 0 U b h h h 3   d x = 12 η 0 U 1   ( b 2 / 4 R )   a ( ξ ) [ h + ( b 2 / 4 R )   a ( ξ ) ] 3   b d ξ .
Defining H 0 4 R h b 2 , we can write the equation as G b 3 8 α   η 0 U   R 2 = 24 1   a ( ξ ) [ H 0 + a ( ξ ) ] 3 d ξ f ( H 0 ) .
Ertel evaluated the integral numerically for a few values of H 0 (presumably laboriously by hand, so some inaccuracy is excusable) and made a log-log plot of G ( H 0 ) which he fitted by a power law. However, in fact the curve is not straight, the slope changing from −3/2 to −4/3 as H 0 increases, and it is better to consider H 0 · G 2 / 3 h ( α η o U ) 2 / 3 R 1 / 3 , (which is the Kapitza film thickness parameter H K ), and plot it against G (see Figure 4).
Also shown in Figure 4 is a solution to the problem by Crook [10], which avoids numerical integration by using the asymptotic approximation for the entrance to a Hertzian contact: h ( b 2 3 R ) ( 2 3 / 2 ) ( x / b 1 ) 3 / 2 . The integration can then be done analytically, giving H K = 2.084   G 1 / 12 . This gives reasonable answers for G > 1 , but is misleading for G < 1 where numerical integration of the exact entry shape shows H K becoming constant. Cameron suggested replacing Crook’s use of the asymptotic equation for the inlet gap by a curve-fit h ( b 2 / 2 R ) · ( 2.114   ( x / b 1 ) 1 / 55 (no extra work is involved). The agreement with Ertel is improved—but again only for G > 1 .
Ertel went on to consider the exit to the parallel region. To get the pressure back down from its high values, a negative pressure gradient is required, which from the Reynolds equation requires the film thickness to be less than the value h : there must be an exit constriction. Ertel argued that this implies a pressure spike (I believe this to be the first occurrence of the concept, predating Petrusevich’s encounter with it in his numerical solution), but this is explained more completely by Greenwood’s extension of the theory [11].

Extension to Ertel’s Theory

The Hertzian geometry provides no opportunity for the pressure to return to zero. However, suppose the Hertzian flat is moved forward [11] from its position on the centre-line of the rollers? The result is as shown in Figure 5.
A Hertzian flat centred on the original point of closest approach requires the standard Hertzian pressure distribution p ( x ) = ( E / 4 R ) b 2 x 2 . However, a flat moved forward a distance d , as shown in Figure 5 ( μ d / b ), can be obtained by applying a pressure distribution p ( x + d ) = E 4 R [ b 2 x 2 + d   ( b + x b x ) 1 / 2 ] .
Figure 6 shows examples of these pressures. When μ d / b is small, the pressure closely follows the Hertz curve, except for a narrow pressure spike at the end of the flat: then as μ increases the resemblance to Hertz disappears and the spike dominates (it should be noted that it is always a square root singularity).
Both the deformed shapes and the pressure distributions of Figure 5 and Figure 6 strongly resemble those found by the numerical solutions to be described later.
The inlet shape resulting from these pressures is known from contact mechanics, so just as for the Ertel analysis above, the Reynolds equation is integrated (for a constant viscosity), then interpreting the result as the reduced pressure q = ( 1 exp ( α p ) ) / α ) instead of p , and setting its terminal value (at the beginning of the flat) equal to 1 / α gives a film thickness equation, which now depends on both the contact width 2 b and the shift d.
To find d / b , a second integration over the exit constriction is performed, terminating at the point where the film thickness returns to h : then matching the two integrals gives d . Then just as for the Ertel solution we could plot G 2 / 3 H 0 h ( α η o U ) 2 / 3 R 1 / 3 , G b 3 8 α   η 0 U   R 2 against (the relation between contact width and Hertz pressure is modified; now, p 0 / E = ( b / 2 R )   ( 1 + 2 μ ) ).
The figure shows that the Ertel extension, while explaining and locating the pressure spike, does not greatly alter the predicted film thickness: indeed, for H K < 1.8 ( G > 1 ), it is hardly changed; and even for G < 0.01 , where there is a 10% reduction in film thickness, the prime feature of H K becoming almost a constant (i.e., the film thickness varying as the two-thirds power of the speed, h U 2 / 3 ) is unchanged.
Please note that now G has two (slightly) different meanings. The factor b is always the half-width of the flat, so in the extended theory it is no longer the Hertzian half-width b, and its relation to p 0 has changed from b = 4 R p 0 / E , to b = 4 R p 0 / E 1 + 2 μ for the half-width of the off-centre flat.
Of course, none of this is a real solution of the EHL problem. Figure 6 shows only the pressure distributions found from elastic theory: it does not show the hydrodynamic pressures, particularly not the inlet pressures with their own, separate, pressure spike at the end of the inlet. The offence may be reduced by modifying Ertel’s condition that the reduced pressure q 1 / α ; we may require it to tend to a lower value so that the hydrodynamic pressure tends to, for example 0.4 p 0 as used in Figure 7.
The outlet pressure distribution neatly continues the elastic spike, so offends less, and we get a plausible approximation of a complete pressure distribution (though now for a specific α p 0 , and with a slightly reduced value of G when α p 0 is small (as here)). However, both contributions must increase the surface deformation, and so invalidate the integration by which they were found!

4. Numerical Solutions and Non-Dimensional Groups

The need for a full solution using hydrodynamic pressures to produce the film shape needed to generate those pressures was clear: a simultaneous solution of the Reynolds equation and the film shape, including the elastic deflection, is needed. The equations (for the line contact problem and assuming the Barus law governs the viscosity are [4,12]:
d p d x = 12 η 0 e α p   u ¯ h h h 3 where h h   ( b ) and p ( b ) = 0 .
h   ( x ) = h 0 + x 2 2 R + u ( x ) u ( 0 ) (for the 2D problem, the deflection has an unknown datum)
u ( x ) = 4 π   E b p ( x ) ln ( | x x | )   d x w = b p ( x )   d x .
The first difficulty is the number of parameters involved in the problem: there are seven physical variables [ η 0 ,   α ,   u ,   R ,   h 0 ,   E ,   w ], so a solution h = F ( U , R , η   0 , w , E , α ) is needed. (Either the film thickness h 0 at x = 0 , or the load per unit width w will be known.)
Dimensional analysis helps, and leads to h / R = F ( α E   ,   η 0 u ¯ / E R ,   W / E R ) , where the three independent variables are the Dowson and Higginson parameters G D α E ,   U D η 0 u ¯ / E R and     W D w / E R . The first published numerical solutions all used this form, with unfortunate consequences, for the only way to report the answers was a grand curve-fitting exercise to obtain h / R = C ·   ( G D ) m (   U D ) n   ( W D ) p . Much of the early literature is concerned with researchers [13] comparing their values for the three indices. This was despite Blok and his colleagues pointing out that by “optimal similarity analysis” it can be shown that only two independent variables are needed, so that film thicknesses can be plotted on a map and curve-fitting avoided [14]. One wonders if this was so long ignored because of the intimidating name [15]! For of course dimensional analysis is a valuable tool for the experimenter but quite unnecessary for a theoretician when all the equations are available (and so of course have dimensional analysis built in!); instead, just scale the equations and see how many (few!) governing parameters are needed. (Appendix A shows the procedure.)
The final two independent variable formats arrived at by Blok and his colleagues were the Moes variables H M = F ( M , L ) where H M h E / η   0 U Σ   R , M w / ( η 0 U Σ E R )   , and L α E   ( η 0   U Σ / E   R ) 1 / 4 (using U Σ 2 U !). These groups may of course be combined in many different ways, though a very desirable measure of the load is the Hertz pressure p 0 ; for if the answers for line contacts and point contacts are to be compared (using the same point contact parameter W D W / E R x 2 for all elliptical point contacts with varying a/b, using a fixed W is highly misleading) we cannot use the load per unit width w for one but the load W itself for the other. This suggests the combination α   p 0 (for line contacts, α   p 0 M 1 / 2 L ). When the pressure distribution is roughly Hertzian, α   p 0 is a useful indicator of the viscosity enhancement, but for low loads (see Gatcombe’s solution!) it is merely a convenient load parameter (for line contacts P α   p 0 is found as α w   E / 2 π   R , for circular point contacts, α   p 0 must be found as ( α / π ) ( 3 W   E 2 / 2 R 2 ) 1 / 3 . For elliptical contacts, approximate answers can be found by replacing R by R 1 R 2 ).
The Dowson and Higginson groups have the great virtue of having speed and load as separate, independent variables, which suits the experimenter; the Blok/Moes group M is less suitable, so more appropriate partners are S α E   ( η 0 u ¯ E R ) 1 / 4 ( S 2 1 / 4 · L ) and P α   p 0 . Finally, the film thickness parameter should vary no more than necessary (charts using the parameter W   h / η U R (to avoid using E or α ) cover a range of 1 to 500), and it has become clear that the Kapitza combination h   /   ( α η 0 U ) 2 / 3 R 1 / 3 ( H K   H M /   L 2 / 3 ) is ideal. So the simplest form is H K = F ( P ,   S ) with H K H M / ( L / 2 ) 2 / 3 = 2 1 / 2 · H M / S 2 / 3 . The fact that the problem requires only two independent variables means that the indices in a three-group equation cannot be independent; writing h / R = C   ( α   E ) k ( η 0 U / E   R ) n ( W / E   R ) m , the indices must satisfy 2 n + m k / 2 = 1 . Thus, the original Dowson and Higginson equation [7] (chapter 7) with k = 0.6 ,   n = 0.7 ,   m = 0.13 has 2 n + m k / 2 = 0.97 and is impossible (uncorrected in the second edition of Elastohydrodynamic Lubrication (1977), sixteen years after the first edition). In the end, the results were reviewed, and the exponent k reduced to 0.54, to give (from ESDU (line contacts) 1985) h min / R = 2.65   G D 0.54 U D 0.70 W D 0.13 or in full h min = 2.65   α 0.54 ( η 0 U ) 0.70 R 0.43 W 0.13 E 0.03 (Note how slight is the dependence on E , and indeed, in the original equation the exponent is positive: E +   0.03 ). The corrected form can be written H K = 2.65   ( G D U D 1 / 4 ) 0.13 ( G D W D 1 / 2 ) 0.26 or substituting G D W D 1 / 2 = 2 π   P , as H K = 2.09     ( S ) 0.13 ( P ) 0.26 ; note the dependence is on the combination ( S / P 2 ) (There is a direct correspondence between the Dowson and Higginson load and speed variables and the preferred ones; thus, for line contacts P = ( G / 2 π )   W D 1 / 2 , S = G U D 1 / 4 ; then H K = ( G U D ) 2 / 3 H D ).
(Even later Pan and Hamrock [16] offered {for ( H D ) min } k = 0.568 ,   n = 0.694 ,   m = 0.128 with a checksum 0.976 instead of 1; fitting their data in a two-variable form leads to k = 0.562 ,   n = 0.695 ,   m = 0.153 , with no loss of accuracy.) The results of the Ertel analysis can also of course be written in terms of S and P . For the Ertel solution, G = 8 P 3 / S 4 or P / S 4 / 3 = 0.5   G 1 / 3 , while for the extended solution, P / S 4 / 3 = 0.5     G 1 / 3 · ( 1 + 2 μ ) 1 / 2 .
The plot is similar to Figure 6, now replacing the abscissa G by 0.5   G 1 / 3 · ( 1 + 2 μ ) 1 / 2 . For P / S 4 / 3 > 0.4 ( d / b < 0.2 ) the two models predict the same film thicknesses.
It is of interest to consider the limiting behaviour of the off-centre flat model as d / b , that is, as the length of the flat tends to zero. Then the deformation becomes very local, and the body effectively a rigid body for which the Gatcombe (like Kapitza, but with the correct Reynolds exit condition) analysis applies. This gives H K * = 2.0425 , and as Figure 8 shows, this is a plausible asymptote to the answers found.

4.1. Regimes of Lubrication

Attempts have been made in the past to identify different “regimes” [17], moving from the classical Rigid-Isoviscous (RI) to Rigid-Piezoviscous (RP) and Elastic-Isoviscous (EI) before reaching the full Elastic-Piezoviscous (EP) regime. (Piezoviscous being a convenient shorthand for “viscosity increasing with pressure”.) One method of determining the boundaries is to locate where the curve-fits obtained for the different regimes become equal; this becomes problematic if the curve-fits are inaccurate. (One chart found in this way located an EP region for which its curve-fit equation was based almost entirely on calculated film thicknesses at points which then proved to lie outside the region! (see [18]).) When maps giving numerical values of the film thickness for substantial ranges of load and speed were not available, such that power law curve-fits had to be employed, it was vital to identify where they were valid, and so to identify the regions noted above, and use the power law for that region. However, now consulting the ESDU charts [19] is easier and more reliable.

4.2. Line Contacts

“Once more detailed (and probably more accurate) values of the film thickness in line contacts became available (Venner [20], it is possible to investigate how accurate are the power law curve fits. Figure 9 compares the Dowson & Higginson power law discussed above with Venner’s calculated values.”
The two certain features in Figure 9 are that the film thickness decreases when the load increases, and that it increases when the speed increases. Clearly there is general agreement between the two sets of values in the area where the Dowson contours H K < 2.2 have been drawn; note in particular the Venner points near the Dowson H K = 1.4 contour. In contrast, the contours of both Ertel and extended Ertel solutions will be straight lines of slope 3/4; very different from the slope 2 of the Dowson contours. However, see Figure 10, showing the interpolated contours of the Venner values! The H K = 2.2 contour is already failing to match the increase of the Venner results, showing that at lighter loads or higher speeds, the film thicknesses are considerably above the H K 1   to   2 range of the right side. The Kapitza parameter is no longer the best film-thickness parameter, and we are approaching the rigid-isoviscous region where the classical (Martin) equation h 0 = 4.895   η   u ¯ R / W holds. Indeed, rewriting the Martin equation in the present terms it becomes H K = 1.56   S 4 / 3 / P 2 . The condition for this to be double the “universal” Kapitza prediction H K min = 1.5 is the line S 1.6 P 3 / 2 , and this is perhaps a better boundary to the EHL region than the Dowson contour. To the right of the H K = 2 contour, or better the line S 1.6 P 3 / 2 , the film thicknesses all more or less agree, but no more. The Dowson contours have the best overall slope, but this is just an average between the slightly negative slopes at low speeds and the positive ( 1 ) slopes on the upper right; there is certainly no suggestion that the points lie along straight lines. The Ertel 3/4 slopes become more plausible in the upper right-hand corner.
When contours are drawn from the Venner [20] data, they are very far from straight, and clearly a power law curve fit, which necessarily assumes straight contours, simply cannot reproduce the answers accurately. There is a strong resemblance to the contours given by the ESDU (point contact) item [19]. Note how the gentle spacing of the contours in the 1 < H K < 2 dramatically changes with the H K = 5 , 10 contours: the move into the rigid isoviscous region is clear.
More generally, we see that only such a plot enables the differences between sets of solutions to be seen: in contrast, a loglog plot of theoretical h / R against measured h / R , covering a range from 10−6 to 10−3 is dominated by the speed dependence, and indeed, one published example, although showing results differing by factors up to 2, “[can be taken as showing] that the equation can be used with some confidence”. Of course, often this may be good enough!

4.3. Circular Point Contacts

Numerical analysis of EHL point contacts inevitably lagged behind the solution for line contacts, but advances in computing power and (perhaps more importantly) in computing techniques mean that film thickness values, and indeed, film thickness maps, are now available. Not surprisingly, the dependence of film thickness on load and speed is similar, and the differences can to some extent be explained as being due to “side leakage”; not all the oil drawn into the inlet needs to be passed downstream through the “Hertz” region, though this leakage is impeded by the side constrictions which develop (and may house the point of minimum film thickness). Figure 11 provides a direct comparison with Figure 10, although it is for the (more confidently determined) central film thickness.
We see the same feature of a large region on the right where the film thickness varies gently from H K = 1 to H K = 2 , bordered on the upper left by the region where H K rises violently. This of course merely confirms the unsuitability of our film thickness parameter for this region. The actual film thickness behaves much as predicted by Kapitza’s initial isoviscous analysis h / R ( η 0 u ¯ R / W ) 2 , so very steeply using our variables.

4.4. Elliptical Point Contacts and Side-Leakage

The study of circular point contacts is very convenient both theoretically and experimentally, but contacts in ball bearings are often elliptical. A film thickness equation drawn from Hamrock and Dowson [21] is often used: H D = 2.69 R x U D 0.67 G D 0.53 W D 0.067 ( 1 0.61 e 0.73 ε ) .
Here the usual Dowson and Higginson variables (now W D W / E R x 2 , where R x is the radius of curvature in the entraining (the minor “b”) direction) are supplemented by an approximate ellipticity factor ε = 1.03 ( R y / R x ) 0.64 (i.e., ε a / b ). This is based on 25 values for k = 6 for varying U D ,   W D , supplemented by 10 values for varying k (from 1 to 8) at U D = 1.68 × 10 12 ,   W D = 0.1106 × 10 6 . Hamrock & Dowson (H&D) [21] report that the curve fit value for a circular contact is 9.4% low, but remarkably Lubrecht et al. [22] report good general agreement with their calculations.
However, of course W D = 0.1106 × 10 6 represents a very different Hertz pressure for a circular contact ( α   p 0 = 5.73 ) than it does for a highly elliptical ( k = 8 ,   α   p 0 = 1.55 —almost a line!) contact: and we know that film thickness decreases when the Hertz pressure is increased. This raises the question, how much of the film thickness reduction below the line-contact value is side leakage, and how much is the effect of “load”, i.e., of Hertz pressure?
To test this, we calculate the Hertz pressure for Hamrock’s 10 points for varying k , and for the 10 points for k = 6 for the same speed, and compared the film thicknesses. Figure 12 shows the result.
The k = 6 points are for varying load ( W D ), but would have negligible side leakage [( 1 0.61 exp ( 0.73 k ) ) = 0.992 ]; the k 6 are for a single W D . Clearly there is no side leakage (and some inaccurate computing!); this merely shows that the film thickness decreases when the Hertz pressure increases. Note the general agreement with a curve fit for all 34 points H K 2.62   · S 0.078 P 0.234 , which is simpler then the H&D equation, and a (slightly) better fit to their values, with a maximum deviation of 6% compared with their “less than 10%” (i.e., 9.7%!).
It is of course well-known that in point contacts, side-lobes develop within the Hertzian area, already for near-circular contacts, but very much more for slender contacts (ones with entrainment along the major axis). For these the minimum film thickness will be much lower that the minimum along the centre-line, although compressibility can sometimes offset this, and give a central film thickness less than the “minimum”. The side lobes will certainly impede side flow, but this is within the “flat”, and so acts after the film thickness has already been determined in the inlet. So it is not clear why the side flow factor ( 1 + 2 R x 3 R y ) 1 predicted by the Kapitza analysis seems to be absent.
The recognition that for “wide” contacts ( k a / b > 1 ) increasing k while keeping W D fixed will reduce the film thickness by increasing the Hertz pressure, and make it harder to identify any effect of side-leakage was late to appear. Perhaps the first statement was by Wheeler et al. [23] 2016, who wrote:
“Yet, Chittenden et al. varied the ellipticity ratio while keeping constant the dimensionless parameters (U, G, and W as defined by Hamrock and Dowson). Consequently, the actual load was varied in a large range together with the ellipticity, and the influence of the latter on film thickness could not be highlighted in itself.”
Yet, they then calculated film thicknesses for fixed speeds and loads (ue = 0.5 m/s, 2 m/s, 10 m/s and w = 120 N, 800 N, 2500 N.), maintaining fixed Hertz pressure by varying R x as k varied, so replacing the confusion between P and k by confusion between S and k. It is not clear how their (no doubt correct) conclusion that the greatest value of h c is when k 2 / 3 could be used.
A direct demonstration of the absence of side leakage (for k > 1 ) is provided by the calculated film thickness values reported in [22]. These are for R y / R x = 1, 2, 5, 10…100, with the loads being adjusted to keep the Hertz pressure constant. Two combinations were studied: [P = 7.9; S = 5.15] and [P = 15.8, S = 10.30]. For each set, the film thicknesses were almost constant, with perhaps slightly lower values for R y / R x = 1, 2 (see Table 1).
It seems that when the Hertz pressure is constant, the ellipticity has no effect on the film thickness. For wide contacts ( k > 1 ) there is no evidence of side-leakage.
Perhaps a more important contribution from Wheeler et al. [24] is their investigation of the accuracy of six different published curve-fit equations for the central film thickness of elliptical contacts. They calculated the film thickness for some dozen combinations, and compared the answers with their own. They conclude that for circular, and “wide” (k = 2.92) elliptical contacts, all the equations overestimate the central film thickness, although for circular contacts the best (Chittenden et al. [CDDT] [25]) are only 2% high. (Hamrock and Dowson (HD) [21] averages 9% high). For the k = 2.92 contacts they find CDDT averages 11%, HD 14%: the best being a recent (2015) set from Mostofi and Khonsari [26] at 5%.
For slender contacts (k = 0.34) the errors are larger. The only equation found (CDDT [27]) gives answers (for h c e n ) which are too low, by between 10% and 30%. (Perhaps more seriously, their answers for h min are too high by between 20% and 200%, but it will be clear to the reader that because of shear thinning the author distrusts all Newtonian h min calculations.)
It should, of course, be remembered that it may be the Wheeler et al. values that are inaccurate. In addition, in fairness to the authors whose equations are studied, note that Wheeler et al. use a different pressure-viscosity law (“an improved Yasutomi correlation”) and a different compressibility law (the Murnaghan equation), so there may be no errors: simply results for a different fluid (see discussion of Sperka et al. 2018 below).

5. Experimental

Direct confirmation of these film thickness predictions was made using disc machines, initially from oil flow measurements [28], then more comprehensive measurements were made by shining X-rays through the oil gap [29], or by measuring the capacitance between the discs [10], which gives an average film thickness, or monitoring the varying capacitance between evaporated strips on the discs, which, at least qualitatively, confirms the predicted shape with its exit constriction [30]. General satisfaction with the agreement between experimental values, numerical solutions, and Ertel’s theory was perhaps partly due to the use of a film thickness parameter with a range of 100, where as noted above, a factor 2 discrepancy seemed minor. Comparison of the exponents in power law curve-fits again showed only “general agreement”. However, the important point, then, was that the oil-film existed!
Extending the capacitance techniques to point contacts is unsatisfactory because of the stray capacitances (much more serious in three dimensions), but fortunately this technique became unnecessary with the introduction of optical interference methods—essentially the use of Newton’s rings.

5.1. Optical EHL

Newton’s rings depend on the path difference between light reflected at the glass/air interface and that reflected at the air/steel interface, but it is essential that the amplitudes of the two beams are comparable. For a glass/air interface, the two refractive indices are substantially different and there is a strong reflected beam, but oil and glass have similar refractive indices, giving a weak reflected beam and consequently poor fringes. To overcome this, a thin, semi-transparent reflecting layer must be deposited on the glass (or preferably synthetic sapphire, with a higher elastic modulus). Much of the development of the techniques was in the Imperial College Lubrication Laboratory [31], but optical EHL rigs are now common, almost invariably for contact between a ball and a rotating disc.
The fringes show a central area recognisable as the “Hertz” contact, but within this a central band which is almost flat, but terminating in an exit constriction with a reduced film thickness much as predicted for a line contact, but in addition there are side lobes with usually, but not invariably, a very much reduced film thickness. Advances in computing techniques and computer speeds have enabled matching contour maps to be found numerically.
We note a minor inconvenience of optical EHL: the contours are unlabelled! To some extent this is circumvented by using white light, for the contours are separated by an integer multiple of the wavelength, resulting in recognisable differences between the appearance of the first few contours; but the main technique is simply to increase the running speed gently, keeping careful count as the contours change. There is also the need to rely on the Lorenz-Lorentz relation between the refractive index and density to predict the refractive index at the relevant pressures. A third inconvenience may offend the practical engineer: one surface must be transparent, so steel/steel contacts cannot be studied. That the elastic modulus of sapphire discs is not too far below that of steel is some defence.

5.2. Pressure Spikes and Minimum Film Thickness

The visual impact of optical EHL is no reason to forget the discoveries of experiments on rollers using sensitive pressure and capacitance transducers bonded to one of the rollers. Hamilton and S.L. Moore [32] believed they had constructed a pressure transducer capable of resolving the predicted pressure spikes: but found only much smaller spikes (the first pressure spikes ever observed); their simultaneous film shape measurements found h min to be much closer to h than predicted. The exit constriction was correspondingly longer than calculated. Hirst & A.J.Moore [33] pointed out that even in pure rolling, the Reynolds equation implies shear stresses reaching ( h / 2 ) ( d p / d x ) ; they argued that the critical shear stress found in their traction measurements (typically 10 MN/m2) was below the shear stress required to generate the pressure spike (which needed 10 MN/m2 on the leading edge of the predicted pressure spike but 100 MN/m2 behind it!). They concluded that because of the modest pressures in the inlet where the central film thickness is determined (it is the lubricant flow rate ρ h which is determined in the inlet), the values found by analysis using a Newtonian fluid should be acceptable, but that to predict the minimum film thickness (or the size of the pressure spike) a non-Newtonian model must be used. Despite this, much effort, both numerical and analytical, has been put into studying the height (or finiteness) of the pressure spike.
Even more have attempts to predict the minimum film thickness using a Newtonian fluid. There is some agreement that for point contacts the correct way to find the minimum film thickness is through a detailed study of the ratio h min / h c e n (minimum to central film thickness). Sperka et al. [34] investigated the ratio for circular contacts for an extensive range of values of L and M, and offered a curve-fit to their calculated results h c e n / h min = 1 + 0.1 α f 0.128 M 0.38 M   . [ ( α f 0.2 ln ( L ) 3 ) / 22.7 ] 2 ] .
The term α f needs some explanation. When the Barus viscosity law η = η 0 exp ( α p ) is replaced by a more complicated law such as the Roelands law, then in effect the constant α becomes a variable. Blok long ago suggested that to make calculated results from different viscosity laws comparable, an equivalent pressure/viscosity index must be used, and he proposed a quantity α defined by ( α ) 1 = I ( ) 0 η 0 η ( p ) d p , and it is usually this quantity, not α 0 , which is used in G D = α E ,   P = α   p 0 etc. Bair et al. [35] argued that since the film thickness is determined by the pressures in the inlet, Blok’s upper limit is too high, and the integration should only be to 3 / α replacing Blok’s I ( ) by I ( 3 / α ) 0 3 / α η 0 η ( p ) d p . However, if the example of Blok’s definition α 1 / I ( ) were followed and α f defined as 1 / I ( 3 / α ) , then for a Barus fluid (where α = α ), we should have the unsatisfactory result α f = α / ( 1 e 3 ) . To avoid this, Bair defines α f ( 1 e 3 ) / I ( 3 / α ) . α f is proposed as the new universal pressure/viscosity index.
Thus, for the three fluids considered in [34], with α 0 = 11 ,   22 ,   33 GPa−1 the Roelands law used gives α f = 8.7 ,   20.6 ,   32.7 GPa−1, and it is these three numbers, with no units, which are to be inserted into the equation. More accurately, what is inserted into the equation is the ratio ( α f / I ) where I is the dimensioned quantity 1 GPa−1; units have been introduced into the non-dimensional equation. To put this differently, the ratio of film thicknesses depends on the length of a bar of platinum in a cellar in Paris! Now of course the authors’ method works: but what can the calculated film thickness ratio really depend on? An immediate possibility is on the ratio α r a t α f i l m / α 0 , and happily we find that ( α f ) 0.2 2 α r a t , so the equation can be rewritten h c e n / h min = 1 + 0.1     ( 2 α r a t ) 0.625 M 0.38 M   · [ ( 2 α r a r ) · ln ( L ) 3 ) / 22.7 ] 2 ] .
We find that this equation is just as good a fit to the data tabulated by the authors as the original. (The power in the second term has been adjusted from 5 × 0.128 to (slightly) improve the fit (actually slightly better: [2.86%, 3.52%, 3.75%] instead of [3.02%, 3.76%, 3.71%]).)
Figure 13 shows how the film thickness ratio varies with load and speed for the authors’ case (3); α 0 = 33   GPa - 1 . It differs from the authors’ figure only in the choice of axes, though that brings out that very extensive ranges of L and M do not cover what might be thought a useful (P, L) area.
Please note that the classical belief that the minimum is 20% to 30% less than the central value (contours 1.2 to 1.4) holds only in a very limited region. The white area is because the data in [34] is limited to M ≤ 1000, and the plot may already have been extrapolated too far: surely the film thickness ratio cannot increase further? The curved contours near the origin hint that they might return into the white area, so the ratio may well decrease there.
Sperka et al. [34] carried out optical film thickness experiments to accompany the numerical work. The agreement between the two is only moderate except for case (3) (that shown in the figure), for which 2.2% accuracy was found. However, undoubtedly film thickness ratios up to 2.6 are measured: so perhaps some lubricants may show Newtonian behaviour.

5.3. Effects of Sliding and Oil Properties

When the two surfaces move at different speeds, as of course gears do except at the pitch point, the shearing of the oil produces heat. This produces a temperature rise in the oil as well as raising the temperature of the bearing/gear surfaces, both of which reduce the oil viscosity. It should be clear from the above that the film thickness is determined by the behaviour in the inlet, so the direct effect on the film thickness is small, but temperatures in the parallel channel will be higher. Calculations suggest that the only effect on the central film thickness will be due to the resulting change in η 0 (and perhaps in α ), and that the film thickness should be calculated using the oil properties at the inlet temperature, estimated as best one can! (Experimentally, a thermocouple dragging on the surface helps.)
The oil is also heated to some extent in the inlet, even in a rolling contact, since the reduction in the gap height in the inlet is itself a form of shearing: but this effect is only serious at very high speeds. A good estimate of the reduction factor is given by the equation f = ( 1 + 0.46 X ) 0.45 , where X ( u ¯ 2 / K ) | d η 0 / d θ | , K being the thermal conductivity of the oil, itself dependent on temperature and pressure [36].
As noted above, lubricating oils are shear-thinning (we shall see later in the discussion of traction that the problem of non-Newtonian behaviour is much worse in the high pressure region). This does not usually affect the central film thickness, except that the viscosity increase in mineral oils due to polymer additives may disappear completely under EHL conditions.

6. Surface Roughness

The desire to know the film thickness was of course to ensure that it exceeded the surface roughness, and traditionally the “lambda ratio” Λ = h min / R q has been studied, where h min is the predicted minimum film thickness for smooth surfaces and R q is the centre-line average (cla) roughness. Values of Λ 3 were recommended, but were often found to be unachievable and not necessary, and lower values ( Λ ~ 1 2 ) were allowed, with warnings about reliability and to avoid substantial sliding speeds. However, of course Λ < 1 —(a film thickness less than the roughness!) is absurd—and yet successful operation is usually achieved. So what is wrong with the lambda ratio concept?
Some effort has been put into finding and solving a modified Reynolds equation applicable to rough surface lubrication, but always regarding the roughness as given. However, when computing techniques advanced to the point when it became possible to solve the EHL problem of sliding over a single, stationary, rough surface, surprising results were found [37]. (All deformation was assumed to be elastic, and heating effects were ignored.) The roughness was largely eliminated, to give almost a uniform gap with the predicted smooth surface film thickness h (and little sign of an exit constriction). Optical EHL experiments gave general support to the numerical findings [38]. In contrast, the pressures resembled a grossly magnified surface roughness profile, superimposed on the Hertz pressure distribution. Thus, the problem was transformed from worrying about scuffing failure to worrying whether the violent pressure fluctuations could lead to a fatigue failure.
To understand this behaviour, we must return to the Ertel argument that at the high pressures of an EHL contact, the viscosity is so high that the Poiseuille term in the Reynolds equation can be neglected, so the line contact Reynolds equation becomes simply h = h ; because the fluid flow rate is constant along the channel, the film thickness must also be constant, even if the surface is rough! There must be additional pressure fluctuations to provide the necessary elastic deformation.
To allow a simple treatment [39,40], we consider a sinusoidal roughness z ( x ) = z 1     sin ( 2 π   x / λ ) on the stationary surface: then the elastic deflection u ( x ) must equal the roughness z 1     sin ( 2 π   x / λ ) . To find the pressures needed to do this, we use the answer for the pressures needed to flatten an infinite wavy surface, where it is well-known that a pressure p 1     sin ( 2 π   x / λ ) gives a deformation u ( x ) = ( 2 λ p 1 / π E )     sin ( 2 π   x / λ ) , so immediately p 1     = ( π   E / 2 λ )     z 1 , and this is a good approximation to the answers obtained in full numerical analyses.
However, this is not the complete story. The Reynolds equation is not properly a statement about the volume flow rate but about the mass flow rate, and should be x ρ h 3 12 η p x + y ρ h 3 12 η p y = U ( ρ h ) x , and for EHL pressures we can no longer treat the density ρ as a constant. Then with Ertel’s observation that the Poiseuille terms can be neglected, the equation reduces to a simple U ( ρ h ) x = 0 ; it is the product ρ   h which is constant, not the film thickness itself.
If the density obeys the law usually assumed in EHL studies, ρ / ρ 0 = ( 1 + γ   p ) / ( 1 + β   p ) , then the bulk compressibility of the fluid B ρ ( d p / d ρ ) will be B p ¯ / C where C = 1 / ( 1 + β   p ¯ )       1 / ( 1 + γ   p ¯ ) . Then for small pressure variations the density variations can be taken as ρ / ρ ¯ 1 + C ( p p ¯ ) / p ¯ (with the usual values γ = 2.266 × 10 9 ,     β = 1.683 × 10 9 m 2 / N , we get C = 0.0665 when p ¯ = 1   G N / m 2 ).
Consider now a smooth geometry with a film thickness h ¯ and a uniform pressure p (which could be the mean Hertz pressure { π   p 0 / 4 in line contact}). We assume the additional pressures to be small (although we shall find they are not!), so that if the pressure fluctuations are also sinusoidal ( p p ¯ ) = p 1 sin ( k x ) , then so will be the density variations: ρ / ρ ¯ 1 + C ( p 1 sin k x ) / p ¯ .
Then ρ   h = ρ ¯   h ¯ gives ρ   ρ ¯ =   h ¯ h ; and h = h ¯ z + v so h ¯ h 1   +     z h ¯ v h ¯ . Thus, 1 + C   p 1 sin ( 2 π   x / λ ) / p ¯ = 1 + ( z 1 / h ¯ ) sin ( 2 π   x / λ ) ( 2 λ   p 1 / π   E   h ¯ ) sin ( 2 π   x / λ ) or C ( p 1 / p ¯ )   +   ( 2 λ   p 1 / π   E   h ¯ ) = ( z 1 / h ¯ ) .
Thus, p 1 p ¯ = 1 C + A z 1 h ¯ where A 2 λ p ¯ π   h ¯ E , relating the amplitude of the pressure fluctuations to the amplitude of the roughness ripples, and perhaps more impressively, h h ¯ z 1 = C C + A . Since C A ( 1 1 + β   p ¯ 1 1 + γ   p ¯ ) · π   E 2 p ¯ h ¯ λ = 24 h ¯ λ for p ¯ = 1   G N / m 2 , and decreases if p ¯ increases, we conclude that the compressibility contribution is usually small compared to the elastic contribution, except for very short wavelengths. Thus, compressibility does prevent the complete reduction of the roughness, but leaves the prediction of the amplitude of the pressures largely unchanged.
We note (applying the results far beyond the validity of the analysis!) that if z 1 = h ¯ when we naively expect solid contact at the roughness maxima, no problems would arise.
The analysis above can easily be extended to a roughness of any form, or to a bump or dent, by representing the roughness by a Fourier series and using superposition (for elastic theory is linear and we have linearised the compressibility). Good agreement is obtained for small amplitude features [40], but when the pressure variation is so large that pressures approaching (or below!) zero are predicted, the agreement becomes poor.
However, this proves to be only the beginning of the story. Hooke [41] derived more complete equations for the attenuated roughness, which for line contact become: h h ¯ z 1 = 1 i C Q 1 i Q i C Q , where C π   E h / ( 2 λ   B ) (and so equals C / A ); and Q = 12 π 2 η   u s h 1 E λ 2 h 2 , where u s is the sliding speed. (The complex form accounts for the phase change needed when reconstructing the result for a general roughness from its Fourier components.). If Q this reduces to h h ¯ z 1 = C 1 + C , which is just the earlier result [39]. In addition, for a Newtonian fluid, Q is indeed enormous. However, note the factor η   u s / h , which is the shear stress. With the enhanced viscosity η 0 exp ( α p ) this is usually orders of magnitude greater than the yield stress of steel: the oil cannot be a Newtonian fluid under these conditions, but must be shear-thinning. This will be discussed at length under “Friction”; here we note only that Hooke [42] showed that the shear stress could be replaced by a limiting value τ 0 , independent of speed, and representing either the “Eyring stress” or a limiting shear stress. Then Q becomes Q = 6 π 2 τ 0 E λ 2 h 2 and is finite, so the full form of the attenuation equation must be used. In addition, by some very impressive optical EHL experiments, Sperka et al. [43] found good agreement with the theory.

6.1. Moving Roughness

If the roughness is on a moving surface, the response becomes more complicated, and indeed, when the first numerical solution for sliding with a moving dent appeared, the results are very hard to believe [44,45]. For the (attenuated) dent and the associated pressure fluctuation became completely separated! The problem is no longer one of steady motion, and the relevant (line contact) Reynolds equation becomes x ρ h 3 12 η p x = u ¯ ( ρ h ) x + ( ρ h ) t or, again dropping the Poiseuille term, u ¯ ( ρ h ) x + ( ρ h ) t = 0 . This is a simple transport equation, with a solution ρ   h = f ( x u ¯ t ) , as well as the simple solution ρ   h = const . This is often called the “particular integral”, and is indeed the solution we just found for stationary roughness, except that now it (the pressures and the displacements) moves through the contact with the speed of the surface carrying the roughness. In effect, the moving roughness (speed u 1 ) is replaced by a moving pressure pulse (speed u 1 ). However, any solution of the form ρ   h = f ( x u ¯ t ) (the “complementary solution”) can be added [46].
Numerical solutions showing typical behaviour are shown in [44]. To understand this, we need to think about the inlet, and unfortunately, we know rather little about it. However, it is clear that each ridge of the surface roughness; however much attenuated, is going to deliver a varying amount of fluid to the parallel gap, and that this must be accommodated. A sine wave entering the inlet moving at speed u 1 must deliver fluid at a frequency u 1 / λ , but it must be carried down the channel at a speed u ¯ , so the pumped fluid must have a wavelength λ u ¯ / ( u 1 / λ ) = λ     ( u ¯ / u 1 ) . Some success in predicting its amplitude has been achieved [47], but currently it seems safer to rely on numerical solutions [48], which suggest the fluid flow can be represented by an attenuated roughness with amplitude given by z 2 / z 1 = ( 1 + 0.125 ¯ + 0.04 ¯ 2 ) 1 where ¯ ( λ / b ) ( ( 2 π ) 3 / 4 P 3 / 2 / L 2 u 1 / u ¯   ) (recall L = 2 1 / 4 S ). The surface amplitudes given by this Lubrecht-Venner equation are just the elastic deflections due to the varying pressures, so directly the amplitude of the pressure (it might be more fundamental to regard the action of the input as pumping a fluctuating volume of lubricant into the channel, where it must be accommodated just as for the removal of the actual roughness by a combination of compression and deformation. If we knew the pumped mass flow < ρ   h >       r sin ( 2 π ( x u 1 t ) / λ ) , the amplitude of the necessary pressure variations would then be p 2 p ¯ = 1 C + A r ρ ¯   h ¯ · ) is p 2 = z 2   [ 2 λ / π   E ] 1 , or p 2 p ¯ = 1 A z 2 h ¯ · where A is now based on the new wavelength λ λ     ( u ¯ / u 1 ) and so becomes A 2 p ¯ π   h ¯ E · λ   u ¯ u 1 .
In Figure 14, the left-hand pair represents the Particular Integral, which moves at the speed u 1 of the rough surface: below, the film thickness (much smaller than the roughness because of the elastic deformation) and above, the large pressure variations needed. The centre pair is the C.F., which moves at the mean speed u m = 1.5 u 1 ( u 2 = 2 u 1 ). Below is the attenuated initial roughness (attn × 0.3 ) and above, the relatively small pressures needed to produce it. The right hand pair shows the combination; both curves vary because of the differing speeds, but because of the relative magnitudes, the film thickness appears to move at u m (and be long wavelength), while the pressure variations move with the rough surface at u 1 .
This is relatively easy to understand. However, using Fourier analysis, a discrete feature such as a bump or a dent can be represented by a set of waves, again with the same two speeds. This can lead to the pressures triggered by the feature reaching the end of the channel and disappearing, while the feature, still in the channel, travels along with “no” effect (for both surfaces moving, but the rough one the slower, so u m > u 1 ). However, if the rough surface is faster, then u m < u 1 and the feature may exit the channel, leaving behind an “inexplicable” pressure pulse!
However, it is not clear how much of the above is real and not just “computer games”! For we shall see when we consider traction that we cannot assume the fluid in the “Hertzian” channel to be a Newtonian fluid, with the enormous viscosity we have been assuming to explain the film thickness, for then the friction would become impossibly high. Studies of waviness using one form of fluid behaviour, an Eyring fluid (see below), show that while the pressure predictions remain, the roughness is no longer removed. However experiments [45] using optical EHL (with a mineral “bright stock” oil) do indeed seem to support the theory (But, of course they cannot confirm the detachment of dent from pressure; only the film shape is available); later experiments [49,50] show agreement with a theory modified for a non-Newtonian fluid [42]. Corresponding experiments examining the “complementary function”, by choosing a feature and speeds to separate the CF from the PI [49] also show agreement with theory, but it is not clear which of many versions of the theory, especially of the modification for a point contact, should be preferred.

6.2. The Λ Ratio

There remains the vital question of how much of the roughness on a real gear tooth can be removed elastically. My belief (but is it just a hope?) is that most of the long wavelength components can be removed, and usually these contribute the majority of the “roughness amplitude” (the cla roughness). So a Λ-ratio less than 1 is not necessarily disastrous.
Perhaps the feature that comes out most clearly from investigations of the critical Λ-ratio is the lack of attention to the nature of surface roughness. It is absurd to associate contact with rms roughness; what possible connection with contact do the valleys have? Back in 1933 Abbott and Firestone knew better, when they invented the Abbott curve, and divided it into peak, medial and valley zones. Contact occurs at high peaks—wear occurs at high peaks—but the main contribution to rms roughness is from the medial roughness (divide the total height range into the top 10%, the middle 80% and the bottom 10%; then the surface area in the middle section is the medial roughness). Hansen et al. [51], by relocation profilometry, showed clearly where contact and wear has taken place, and will not take place again; the minimal change to the rms roughness is self-evident. Not surprisingly, they report successful operation with a Λ-ratio of 0.33.
Other assessments are perhaps more realistic. The conclusion reached by Cann et al. [52] “the Λ criterion [is] at low values inapt, leading, perhaps, to an overly conservative or widely optimistic view of the level of lubrication” after they decide that the film-thickness prediction h is untrustworthy because of non-Newtonian effects (especially with grease), and starvation, and that the roughness σ is irrelevant because the roughness is composed of a spectrum of amplitudes, all contributing to the cla (or rms), but not all present in the high pressure region. Jacobson [53] reinforces the worries about the effect of starvation on the film thickness: but also emphasises that surfaces will often run-in to a well-lubricated state, but that sometimes they do not!
Despite problems over whether in general the Λ-ratio does or does not indicate contact, attempts have been made to link contact with friction. Clarke et al. [54], using electrical contact resistance as an indicator of contact, plot a delightful pair of graphs, the first showing how there is negligible contact for Λ 1 , but a dramatic fall to full contact (actually to zero contact resistance) as the rolling speed is reduced and Λ falls to 0.5. The second shows the measured friction, falling slightly as Λ increases from 0.5, with no signs of any distinction between the values for Λ 1 (full film?) and those for Λ between 0.5 and 1 (partial contact?). Zapletal et al. [55] show the relative contact area (based on optical EHL measurements) rising from zero for Λ < 0.7 , but friction rising steadily from the full-film value for Λ < 3 . So it appears that roughness can increase friction without contact occurring! (Both studies are for a ball on a plane. We note wryly that while [54] uses the Chittenden et al. film thickness equation, based on point contact film thickness calculations for many geometries including circular contacts, [55] uses the earlier Hamrock and Dowson equation, based almost entirely on elliptical contact solutions ( a / b = 6 ) and on a single solution for a circular contact, and it differs by 4% from the curve fit. However, such a change is not important.)

7. Friction

It is not clear how much notice need be taken of film thicknesses, whether calculated or measured: engineers built and used gears very happily when lubrication theory predicted film thicknesses far below the known surface roughness; and while it is perhaps satisfactory to know that gear lubrication is possible, it is no great surprise. Of course, an approximate value is needed, for the shear stress is the ratio of the sliding velocity to the film thickness. It does not depend, as most authors seem to believe, on the slide/roll ratio. This may be unfair, for the film thickness depends approximately on the two-thirds power of the rolling speed, so | u 2 u 1 | h | u 2 u 1 | | u 2 + u 1 | , roughly!

7.1. Thermal Effects

More important is the prediction of friction, and now that we have EHL theory and an exponential viscosity increase, this is a problem. For we now know that we essentially have a long parallel channel containing oil with an enormous viscosity [ η = η 0 exp ( α   p 0 ) ], suggesting a friction force ( 2 b ) η   U s / h , where U s is the sliding speed U s | u 1 u 2 | , and this is impossibly large. However, this energy dissipation increases the oil temperature, and so decreases the viscosity. Crook [56] showed that although for thrust or journal bearings the heat is carried away by the hot oil leaving the bearing (and where of course the oil temperature rise is often a design parameter), under EHL conditions the heat flow is across the film to the bearing surfaces, and the heat is then convected away by the moving parts.
Crook [57] derived a simple implicit equation for the temperature rise across the film:
U s = θ 1 θ c 2 K η   ( θ ) d θ + θ 2 θ c 2 K η   ( θ ) d θ where θ 1 ,   θ 2 are the temperatures of the surfaces, and θ c is the maximum oil temperature and K the thermal conductivity of the oil, and derived a corresponding equation for the shear stress. Here we shall take the two surface temperatures to be the same, so Crook’s equations become U s = 2   { θ 1 θ c 2 K η   ( θ ) d θ } 1 / 2 and τ · h = 2 θ 1 θ c 2 K   { θ θ c 2 K η   ( θ ) d θ } 1 / 2 d θ .
If we assume the oil thermal conductivity to be constant and the viscosity to obey η = η 1 exp ( γ   θ ) , so that η 1 is the viscosity at the temperature of the surfaces, still increased by the pressure term exp ( α   p ) ], we obtain U s = 2 2 K γ   η 1   X   and τ · h = 2 2 K η 1 γ sinh 1 X 1 + X 2 where γ   ( θ c θ 1 ) = ln ( 1 + X 2 ) . Figure 15 shows the resulting traction.
It must be emphasised that this graph is very much a half-truth. As well as the simplification of assuming the two surfaces to be at equal temperatures, the temperatures have been taken as fixed. In practice they will not be equal, and will be higher when the sliding speeds and thus the heat generation increase: at high sliding speeds the curve will be misleadingly high. However, it does illustrate the important point that one must not refer to the oil film temperature: the markers give the difference between the maximum oil temperature (in this case mid-film) and the oil temperature at the surface. With a typical coefficient γ = 0.02 as assumed by Crook, there is already, at the curve maximum, a 50 °C temperature difference within (across!) the oil!
However, of course a disc machine does not measure shear stress, but only traction, and this shear stress prediction is for a particular combination of pressure and surface temperatures. To predict the traction, the shear stress must be integrated along the contact, taking account of the pressure variation (including the pressure spike!). However, the surface temperatures ( θ 1 and θ 2 ) must also be calculated by the usual methods of flash temperature calculation for sliding contacts. The general results are unchanged: the expected proportionality between sliding speed and shear stress holds only for low sliding speeds. The shear stress and the traction reach a maximum, and then fall as the sliding speed is increased further. For moderate pressures and rolling speeds the final predictions agree qualitatively with experiment [57,58].

7.2. Non-Newtonian Behaviour

For thinner films (higher pressures and lower rolling speeds) the predicted traction is too high, and the maximum occurs at lower sliding speeds than predicted: the oil must be showing some form of non-Newtonian behaviour!
One explanation is that the oil is an “Eyring fluid” with a shear-stress/shear-strain behaviour predicted by the Eyring thermal activation energy theory of fluid flow. According to this, the strain rate γ ˙ obeys a relation ( η     γ ˙ ) = τ 0 sinh ( τ / τ 0 ) , where η and τ 0 are related to molecular properties. Eyring’s theory predicts the exponential increase of η with pressure, and an exponential decrease with (1/T). It seems clear that in the physical chemistry world the predicted hyperbolic sine was not taken seriously; that τ / τ 0 was assumed to be small so that sinh ( τ / τ 0 ) τ / τ 0 , and this was an explanation of Newtonian behaviour. In contrast, in EHL the Eyring theory is treated seriously, and the full η     γ ˙ = τ 0 sinh ( τ / τ 0 ) is used, typically with τ 0 = 4   MN / m 2 , but probably increasing with pressure.
When Crook’s thermal analysis is repeated using the Eyring relation [59], the results are U s = 2   τ τ   { θ 1 θ c 2 K η   ( θ ) d θ } 1 / 2 and h · τ = 2   τ τ θ 1 θ c 2 K   { θ θ c 2 K η   ( θ ) d θ } 1 / 2 d θ , where τ τ 0 sinh ( τ / τ 0 ) ; these reduce to the equations for a Newtonian fluid when τ / τ 0 is small, so that τ τ . We write these as U s =   τ τ   F ( θ c ) and h · τ =   τ τ G ( θ c ) . To compare Eyring and Newtonian behaviour, we assume the same η   ( θ ) law for both fluids, and denote the shear stress and the corresponding sliding speed of the solution for a Newtonian fluid by τ N and V s , so that V s =     F ( θ c ) and h · τ N =   G ( θ c ) . Then U s / V s =   τ τ   and τ / τ N =   τ τ τ / τ 0 sinh ( τ / τ 0 ) . The last equation may be solved iteratively to obtain τ / τ 0 as a function of τ / τ N , and hence to get τ / τ N as a function of τ N / τ 0 . Thus, we can scale the Newtonian τ N ( V s ) to obtain the Eyring curve τ ( U s ) by shifting each point to the right by the factor τ / τ N , and downwards by the factor τ N / τ (Figure 16). Please note that while the Newtonian curve is linear almost to γ   Δ θ = 0.5 , the Eyring curve is non-linear well before this (the parameter h τ 0 / M used here is the reciprocal of the parameter χ of the original Johnson and Greenwood paper. Μ is not related to the Moes parameter mentioned earlier).
The reduced shear stresses due to the Eyring behaviour shown above (and much larger reductions can occur) raise the question of whether both methods of reducing the shear stresses must be considered?
An isothermal Eyring analysis for a parallel channel is easily done [60]: the equilibrium condition gives the shear stress across the film as τ = τ 1 + z · P ) where P is used for the pressure gradient d p / d x , and substituting in the Eyring equation gives u z = τ 0 η sinh ( τ 1 + z · P τ 0 ) . Integrating and substituting in the boundary conditions on speed leads to u 2 u 1 = τ 0 2 η P [ cosh { τ 1 + h P / 2 τ 0 } cosh { τ 1 h P / 2 τ 0 } ] = 2 τ 0 2 η P [ sinh { τ 1 τ 0 } sinh { h P 2 τ 0 } ] from which τ 1 = τ 0 · arcsinh { η   ( u 2 u 1 ) h   τ 0 · Φ 1 ( { h · | P | / 2 τ 0 } ) } where Φ 1 ( V ) V / sinh ( V ) .
The rapid increase of sinh(V) with V suggests the possibility of a greatly reduced effective shear rate, but in fact the pressure gradient d p / d x will be of the order p 0 / b , so the term will be of the order ( h / b ) ( p 0 / τ 0 ) . The second factor will be large, typically 500 (2 GN/m2/4 MN/m2), but the first factor will be small (0.2 mm/0.5 µm) giving V 1 and Φ 1 , so a reasonable approximation is simply τ 1 = τ 0 · arcsinh { η   ( u 2 u 1 ) h   τ 0 } (especially near the maximum pressure location, where it matters most). Thus, the shear stress steadily increases with sliding speed; the traction decrease found experimentally by Crook, with a maximum found to be too high and at too high a sliding speed in Johnson’s experiments, does not occur. To account for it, both shear-thinning and temperature rise must be allowed for.
Please note that the shear stress on the moving surfaces is not τ 1 : just as for a Newtonian fluid it is ( h / 2 ) ( d p / d x ) ± τ 1 , depending on which surface is the faster. However, with h constant and the pressure effectively zero at each end of the channel, so that b b ( d p / d x )   d x = 0 , only τ 1 contributes to the traction.
As noted above, the shear stress curve is merely the first step towards predicting the traction curve, for that requires the integration over the contact area, taking into account the pressure and temperature variations. To some extent an approximate traction curve can be obtained using the shear stress curve for the mean pressure and estimated surface temperatures: and this basis has led experimenters to believe that oil behaviour is even more complex, with in addition to the non-linear relation between shear stress and strain, also a limiting shear stress (and of course steels have a limiting shear stress, so why not oils?). Although earlier it was believed that these were alternatives, favoured by different groups, it now seems that both a shear-thinning mechanism and a limiting shear stress are needed. Even this is not all; under some conditions it seems the oil has an elastic modulus and so gives a viscoelastic response.
Thus, the prediction of the friction force requires many choices.
(1)
It is generally agreed that the Barus equation for the pressure dependence ( exp { α   p } ) of the oil viscosity is inadequate. The initial exponential increase is not continued, and this is described by the Roelands (2-parameter) equation [61]. This correctly describes the viscosity up to 0.4 GN/m2 and is adequate for (central) film thickness predictions. However, the region where the traction is developed has much higher pressures, and direct measurements up to 1.4 GN/m2 in a falling ball viscometer show an upturn in the log(viscosity) curve, so that the Roelands viscosities are then too low, and the Yasutomi equation [62] should be used.
(2)
An equation for the dependence of the viscosity on temperature is needed to allow for the surface temperature rise. The simple η exp ( γ   θ ) is unjustified at high pressures, and values from the full Roelands equation (which then becomes a four-parameter equation) may be better. However, it is claimed that experimental values require the use of the full Yasutomi equation (which then has eight parameters). A full thermal analysis of the whole configuration to find the surface temperatures is needed.
(3)
Unless the neglect of temperature variations within the oil film is justified, an acceptable equation for the variation of the oil thermal conductivity with pressure and temperature is needed.
(4)
Which non-Newtonian model should be used, and what are its parameters? Are the Eyring stress τ 0 (or the limiting shear stress) constants, or do they depend on temperature and pressure?
All the above must then be combined with a full numerical EHL analysis to determine the film thickness and the complete pressure distribution.
What are actually available are full isothermal analyses, using the Roelands equation for the viscosity and the Eyring shear stress model (Jacod [63]). A curve fit over all the manifold cases studied led to a beautifully simple master curve. Oddly, this is published in a form reminiscent of Hooke’s famous publication of “ut tension, sic vis” as an anagram, and it was 10 years later before one of the authors published the solution [64]: that the friction depends on the simple combination τ c 2   { η ( p 0 ) · u s / h     }     /     τ 0 ) , where u s is the sliding speed | ( u 2 u 1 ) | (I malign them: they use the product of the rolling speed and the slide/roll ratio, not the sliding speed), so that this is just (twice) the ratio of the Newtonian shear stress at the pressure maximum to the Eyring stress τ 0 —the combination already introduced above! Then the master curve is μ = ( τ 0 / p 0 )   sinh 1 ( τ c / 5 ) . This is shown to be a good fit to most of the numerical examples, except for slow sliding speeds, where it is suggested that an additional contribution comes from rolling friction, and may generally be neglected. More importantly, experimental measurements were found to lie on the master curve, although in fitting them, the estimates of τ 0 from the experimental traction curves (by the original authors) were used. We might note that the factor 5 is of little importance, for because of the factor η ( p 0 ) , the values of τ c are enormous, and sinh 1 ( τ c / 5 ) ~ ln ( 2 · τ c / 5 ) = ln ( τ c ) 0.9 ; the shift ( 0.9 ) is insignificant compared with perhaps ln ( 10 10 ) = 23 .
However, the master curve shows the traction steadily increasing with shear rate: there is no indication of the maximum reported by the experimenters (either a thermal effect or due to reaching a limiting shear stress), or of the subsequent decreasing branch of the curve (surely thermal?). A later publication presents a master curve for a limiting shear stress fluid [65]—this is 1 μ = p 0 τ 0 [ 1 + ( 4 / 3 )   sinh 1 ( τ c / 5 )   ] —but now τ c is calculated not using η ( p 0 ) as τ c , but from the viscosity at a modified (lower) pressure p . Experimental results are given, agreeing and showing the traction curve beginning to flatten out. However, in contrast, results from Habchi et al. [66] show a picture more resembling Figure 16, although their calculations include a different form of shear-thinning (see below) and a limiting shear stress as well as thermal effects. Both experiments and calculations (for a steel ball sliding on a steel disc at 2 m/s, p0 = 0.85 GN/m2) show the friction rising rapidly before levelling out at μ 0.07 (when u s 0.2 m / s ), then the experiment and the full solution show a gentle decrease to μ 0.06 at u s 4   m / s . They also present an isothermal solution, and this gives a continuing gentle rise to μ 0.075 . However, the authors emphasise the need to use the exact (independently measured) true thermal properties (thermal conductivity and specific heat) of the oil, varying with temperature and pressure, instead of constant values. They emphasise the importance of using real values of viscosity variation throughout the contact rather than using the Roelands equation to describe its dependence on pressure and temperature. They cast doubt on the Crook thermal analysis, arguing that heat convection along the channel is significant, and also assert that heat generated from fluid compression where the pressure is rising, and absorbed during decompression when the pressure falls, must be taken into account. Perhaps the most recent (2019) contribution [67] shows experimental and calculated traction curves for squalane, again using measured fluid properties. There is reasonable agreement, again with a maximum friction coefficient (0.05 at p0 = 1 MN/m2 or 0.07 at p0 = 2 MN/m2), followed by a gentle decrease.
Both the two last cited contributions [66,67] use yet another shear-thinning law, referred to as the “Carreau equation”. The Carreau equation η = η + ( η 0 η ) ( 1 + ( λ γ ˙ ) a ) ( n 1 ) / a (but in a general tensor form) is widely used to describe rheology in chemical engineering. Liu et al. [67] do indeed use this, in its basic form with a = 2 and omitting the term η   , so that no “second Newtonian regime” is included (but is believed to be needed with polymer-thickened oils). With n < 1 this does describe shear-thinning (they take n = 0.463). In contrast, the equation used by Habchi et al. [66] is η   ( τ ) = η ( 0 ) . [ 1 + ( τ / G ) a ] ( 1 1 / n ) / a (with, for Shell T9, G = 7 MN/m2, a = 5, n = 0.35.), with the shear stress τ replacing the shear rate γ ˙ as the independent variable, does not seem to be a Carreau equation, but with three fitting constants certainly provides a convenient framework for fitting experimental data.
For comparison, note that the Eyring equation may be written η   ( τ ) = η   ( 0 ) . [ ( τ / τ 0 )   /   sinh ( τ / τ 0 ) ] , the Eyring stress τ 0 replacing G : the two equations differ only qualitatively, and the solution procedures will be the same.
However, can it really be true that unless oil companies supply oil accompanied by a complete set of parameters, traction cannot be predicted? (Habchi et al. [66] detailed 18 oil parameters, presumably all oil-specific, plus 4 material parameters). Traction coefficients all seem to reach maxima in the band 0.05 to 0.1, after a linear range at low sliding speeds, and followed by a gentle decrease as the sliding speed increases. Can a simple guide not be derived? Habchi et al. [68] make perhaps the first step at establishing what must be included by eliminating in turn the various factors; in contrast Liu et al. [67] obtain acceptable predictions by simply assuming a Hertzian pressure distribution. Is there hope?

8. Conclusions

EHL theory is now able to predict film thicknesses in heavily loaded contacts, and to explain how these can happen. For this, the classical idea of a lubricating oil as Newtonian fluid, but with its viscosity strongly pressure dependent (to accompany the known temperature dependence) is all that is needed. There remains some doubt as to whether the famous, and embarrassing, pressure spikes are real, and correspondingly the predicted minimum film thicknesses may be somewhat pessimistic.
Surface roughness is seen to be a less immediate threat than traditionally believed, as its amplitude as it passes through the high pressure region will be greatly reduced by elastic deformation: but the associated pressure ripples may “assist” fatigue failure.
In the end, it seems that at present no useful guide to traction for the practising engineer is available. Claims in the literature that traction curves can only be predicted by using full detailed information obtained by direct, independent tests to obtain the detailed variation of oil properties must, one hopes, be pessimistic. Certainly, traction initially rises with sliding speed before levelling off: almost certainly it then falls as the sliding speed increases further. However, the details of this behaviour—the actual values and the pressures and speeds at which they occur—vary between individual oils, and at this point no general advice is available.

Funding

This research received no funding.

Conflicts of Interest

There is no conflict of interest.

Appendix A. Minimising the Number of Governing Parameters

The problem:
d p d x = 12 η 0 e α p   u ¯ h h h 3 where h h   ( b ) and p ( b ) = 0 .
h   ( x ) = h 0 + x 2 2 R + u ( x ) u ( 0 ) (For 2D problem, the deflection has an unknown datum) u ( x ) = 4 π   E b p ( x ) ln ( | x x | )   d x w = b p ( x )   d x .
We assume that h 0 is known, and the corresponding load is to be found.

Reduction

First refer all vertical distances to h 0 : thus H ( x ) = h ( x ) / h 0 and U ( x ) = u ( x ) / h 0 . Then H   ( x ) = 1 + x 2 2 R h 0 + U ( x ) U ( 0 ) , U ( x ) = 4 π   E h 0 b p ( x ) ln ( | x x | )   d x and the Reynolds equation becomes d p d x = 12 e α p   η 0   u ¯ h 0 2 H H H 3 .
A less obvious move is to choose our own units for the internal quantities x and p ; why should we use m and N/m2? So let x = A ξ and p = B q . (Of course, normally A = 1 m, B = 1 GN/m2: the computer only accepts numbers). Then H = 1 + A 2 ξ 2 2 R h 0 + U ( ξ ) U ( 0 ) ; U ( ξ ) U ( 0 ) = 4 A B π   E h 0 A b q ( ξ ) [ ln ( | ξ ξ | ) ln | ξ |   ]   d ξ and d q d ξ = 12 e α B p   A η 0   u ¯ B h 0 2 H H H 3 , so the problem is governed by the four parameters A 2 / R h 0 , A B / E h 0 , α   B and A η 0 u ¯ / B h 0 2 . However, A and B are at our disposal, so we may eliminate two of the four. An obvious choice is B = 1 / α , so leaving A 2 / R h 0 , A / α   E h 0 , and A α   η 0 u ¯ / h 0 2 . Any one of the three may be removed: for example, A = R h 0 . Then only two parameters govern the problem: Y ( α   E ) 1 R / h 0   and Z α   η 0 u ¯ R / h 0 3 / 2 . With this choice, the equations to be solved are H   ( ξ ) = 1 + 1 2 ξ 2 + U ( ξ ) U ( 0 ) ; U ( ξ ) = 4   Y π   b q ( ξ ) ln ( | ξ ξ | )   d ξ and e q d q d ξ = 12 Z   H H H 3 , the boundary conditions remaining effectively the same.
It may be reassuring to note that this is merely a particular case of the original problem: identical equations would have arisen if we had chosen R = 0.01 m, h 0 = 10 6 m, E = 100   GN / m 2 , α = 10 8 m 2 / N with Y = 0.1, and u ¯ = 1   m / s , η 0 = 0.1   Ns / m 2 with Z = 0.1: but of course now we have not only the answer for these values but for many other combinations.
It remains to find the load: our solution gives W = b q ( ξ )   d ξ , but we want w = b p ( x )   d x Then w = A B W = R h 0 α W ; so the load variable is W α   w R h 0 , and we have found W = F ( Y , Z ) . However, we can manipulate these variables: we have ( Z / Y 3 ) 1 / 4 = ( α   E ) ( η 0 u ¯ / E R ) 1 / 4 , which is our speed variable S , and ( W / Y ) 1 / 2 = ( α   ) ( w E / R ) 1 / 2 , which except for a factor 1 / 2 π   is our load variable P. Then since Z 2 / 3 = h 0 / ( α η 0 u ¯ ) 2 / 3 R 1 / 3 and the solution will have given the full film shape, with in particular the value of h h ( b ) , in the form H K = F 1 ( P ,     S ) .

References

  1. Williams, J.A. Engineering Tribology (Chapter 7); OUP: Oxford, UK, 1994. [Google Scholar]
  2. Martin, H.M. Lubrication of Gear Teeth. Eng. Lond. 1916, 102, 119–121. [Google Scholar]
  3. Johnson, K.L. Contact Mechanics (Chapter 4); CUP: Cambridge, UK, 1985. [Google Scholar]
  4. Dowson, D.; Higginson, G.R. Elasto-hydrodynamic Lubrication; Pergamon: Oxford, UK, 1977. [Google Scholar]
  5. Gatcombe, E.K. Lubrication characteristics of involute spur gears. Trans. ASME 1945, 67, 195. [Google Scholar]
  6. Blok, H. Discussion. The lubrication of gears. Gear Lubrication Symposium. J. Inst. Petrol. 1952, 38, 673. [Google Scholar]
  7. Kapitza, P.L. Hydrodynamic theory of lubrication during rolling. Zhurnal Tekhnicheskoi Fiziki 1955, 25, 747–762. [Google Scholar]
  8. Grubin, A.N. Investigation of the contact of machine components. Central Scientific Research Institute for Technology and Mechanical. Eng. Mosc. 1949, 30. [Google Scholar]
  9. Mohrenstein-Ertel, A. Berechnung der Hydrodynamischen Schmierung Gekrummter Oberflachen unter Hoher Belastung und Relativbewegung. In Fortschrittsberichte VDI.; Lang, O.R., Oster, P., Eds.; VDI Verlag: Dusseldorf, Germany, 1984. [Google Scholar]
  10. Crook, A.W. The lubrication of rollers - II. Film thickness with relation to viscosity and speed. Philos. Trans. R. Soc. Lond. Ser. A 1961, 254, 223–237. [Google Scholar]
  11. Greenwood, J.A. An extension of the Grubin theory of elastohydrodynamic lubrication. J. Phys. D Appl. Phys. 1972, 5, 2195–2211. [Google Scholar] [CrossRef]
  12. Dowson, D.; Higginson, G.R. A Numerical Solution to the Elasto-Hydrodynamic Problem. J. Mech. Eng. Sci. 1959, 1, 6. [Google Scholar] [CrossRef]
  13. Cheng, H.S. Isothermal Elastohydrodynamic Theory in Full Range of Pressure/Viscosity Coefficient. ASME J. Lubr. Technol. 1972, 94, 35–43. [Google Scholar] [CrossRef] [Green Version]
  14. Moes, H. Discussion of a paper by D. Dowson. Proc. Inst. Mech. Engrs. 1965, 180, 244–245. [Google Scholar]
  15. Moes, H. Optimum similarity analysis with applications to EHL. Wear 1992, 159, 57–66. [Google Scholar] [CrossRef] [Green Version]
  16. Pan, P.; Hamroch, B.J. Simple formulae for performance parameters used in elastohydrodynamic lubricated line contact. ASME J. Tribol. 1989, 111, 246–251. [Google Scholar] [CrossRef]
  17. Johnson, K.L. Regimes of elasto-hydrodynamic lubrication. J. Mech. Eng. Sci. 1970, 12, 9–16. [Google Scholar] [CrossRef]
  18. Greenwood, J.A. Film thicknesses in circular elastohydrodynamic contacts. Proc. I. Mech. Engrs. 1988, 202, 11–16. [Google Scholar] [CrossRef]
  19. ESDU 89045. Film Thicknesses in Lubricated Hertzian Contacts. Part 2: Point Contacts; IHS ESDU: London, UK, 1989; ISBN 978-0-85679-719-4. [Google Scholar]
  20. Venner, C.H. Multilevel Solution of the EHL Line and Point Contact Problems. Ph.D. Thesis, University of Twente, Enschede, The Netherlands, 1991. [Google Scholar]
  21. Hamrock, B.J.; Dowson, D. Isothermal elastohydrodynamic lubrication of point contacts: Part III—Fully flooded results. J. Lubr. Technol. 1977, 99, 264–275. [Google Scholar] [CrossRef]
  22. Lubrecht, A.A.; Venner, C.H.; Colin, F. Film thickness calculation in elasto-hydrodynamic lubricated line and elliptical contacts: The Dowson, Higginson, Hamrock contribution. Proc. I. Mech. Engrs. Part J: J. Eng. Tribol. 2009, 223, 511–515. [Google Scholar] [CrossRef]
  23. Wheeler, J.-D.; Fillot, N.N.; Vergne, P.; Philippon, D.; Morales Espejel, G.E. On the crucial role of ellipticity on elastohydrodynamic film thickness and friction. Proc. I. Mech. E Part J: J Eng. Tribol. 2016, 230, 1503–1515. [Google Scholar] [CrossRef]
  24. Wheeler, J.D.; Vergne, P.; Fillot, N.; Philippon, D. On the relevance of analytical film thickness EHD equations for isothermal point contacts: Qualitative or quantitative predictions? Friction 2016, 4, 369–379. [Google Scholar] [CrossRef] [Green Version]
  25. Chittenden, R.J.; Dowson, D.; Taylor, C.M. Elastohydrodynamic film thickness in concentrated contacts: Part 1: Experimental investigation for lubricant entrainment aligned with the major axis of the contact ellipse. Proc. R. Soc. Lond. 1985, 397, 245–269. [Google Scholar] [CrossRef]
  26. Masjedi, M.; Khonsari, M.M. On the effect of surface roughness in point-contact EHL: Formulas for film thickness and asperity load. Tribol. Intl. 2015, 82, 228–244. [Google Scholar] [CrossRef]
  27. Chittenden, R.J.; Dowson, D.; Dunn, J.F.; Taylor, C.M. Elastohydrodynamic lubrication of concentrated contacts—Part 2: General case, with lubricant entrainment along either principal axis of the hertzian contact ellipse. Proc. R. Soc. Lond. A 1985, 397, 271–294. [Google Scholar]
  28. Crook, A.W. The lubrication of rollers. Philos. Trans. R. Soc. Lond. Ser. A 1958, 250, 387. [Google Scholar]
  29. Sibley, L.B.; Orcutt, F.K. Elasto-hydrodynamic lubrication of rolling contact surfaces. Trans. Amer. Soc. Lub. Engrs. 1961, 4, 234. [Google Scholar] [CrossRef]
  30. Crook, A.W. Elasto-hydrodynamic lubrication of rollers. Nature 1961, 190, 1182. [Google Scholar] [CrossRef]
  31. Foord, C.A.; Wedeven, L.D.; Westlake, F.J.; Cameron, A. Optical Elastohydrodynamics (Part 1). Proc. I. Mech. Engrs. 1969, 184, 487–505. [Google Scholar] [CrossRef]
  32. Hamilton, G.M.; Moore, S.L. Deformation and pressure in an elastohydrodynamic contact. Proc. Roy. Soc. 1971, 322, 313–330. [Google Scholar]
  33. Hirst, W.; Moore, A.J. Non-Newtonian behaviour in elastohydrodynamic lubrication. Proc. Roy. Soc. 1974, 337, 101–121. [Google Scholar]
  34. Sperka, P.; Krupka, I.; Hartl, M. Analytical Formula for the Ratio of Central to Minimum Film Thickness in a Circular EHL Contact. Lubricants 2018, 6, 80. [Google Scholar] [CrossRef] [Green Version]
  35. Bair, S.; Liu, Y.; Wang, Q.J. The pressure–viscosity coefficient for Newtonian EHL film thickness with general piezoviscous response. J. Tribol. 2006, 128, 624–631. [Google Scholar] [CrossRef]
  36. Dyson, A. Discussion of Inlet shear heating in elastohydrodynamic lubrication by Greenwood J.A. & Kauzlarich J.J. J. Lub. Tech. 1973, 95, 417–423. [Google Scholar]
  37. Kweh, C.C.; Evans, H.P.; Snidle, R.P. Micro-Elastohydrodynamic Lubrication of an Elliptical Contact With Transverse and Three-Dimensional Sinusoidal Roughness. J. Tribol. 1989, 111, 577–584. [Google Scholar] [CrossRef]
  38. Kaneta, M.; Sakai, T.; Nishikawa, H. Effects of Surface Roughness on Point Contact EHL. Tribol. Trans 1993, 36, 605–612. [Google Scholar] [CrossRef]
  39. Greenwood, J.A.; Johnson, K.L. The Behaviour of Transverse Roughness in a Sliding Elastohydrodynamically Lubricated Contact. Wear 1992, 153, 107–117. [Google Scholar] [CrossRef]
  40. Greenwood, J.A.; Morales-Espejel, G.E. The behaviour of real transverse roughness in a sliding EHL contact. In 19th Leeds-Lyon Symposium: Thin Films in Tribology; Dowson, D., Taylor, C.M., Childs, T.H.C., Godet, M., Dalmaz, G., Eds.; Elsevier Science Publishers B.V.: Amsterdam, The Netherlands, 1993. [Google Scholar]
  41. Hooke, C.J. Surface Roughness Modification in Elastohydrodynamic Line Contacts operating in the Elastic Piezoviscous Regime. Proc. Inst. Mech. Engrs: Part J 1998, 212, 145–162. [Google Scholar] [CrossRef]
  42. Hooke, C.J. The Behaviour of Low-Amplitude Surface Roughness under Line Contacts: Non-Newtonian Fluids. Proc. Inst. Mech. Engrs. Part J: J. Eng. Tribol. 2000, 214, 253–265. [Google Scholar] [CrossRef]
  43. Šperka, P.; Kǔpka, I.; Hartl, M. Experimental study of roughness effect in a rollingsliding EHL contact. Part I: Roughness deformation. Tribol. Trans. 2016, 59, 267–276. [Google Scholar] [CrossRef]
  44. Venner, C.H.; Lubrecht, A.A. Transient Analysis of Surface Features in an EHL Line Contact in the Case of Sliding. J. Tribol. 1994, 116, 186–193. [Google Scholar] [CrossRef]
  45. Kaneta, M.; Kanada, T.; Nishikawa, H. Optical Interferometric Observations of the Effects of a Moving Dent on Point Contact EHL. In Elastohydrodynamics 96, Proceedings of the 23rd Leeds-Lyon Symposium on Tribology; Elsevier Tribology Series 32; Elsevier: Amsterdam, The Netherlands, 1997. [Google Scholar]
  46. Greenwood, J.A.; Morales-Espejel, G.E. The behaviour of transverse roughness in EHL. Proc. Instn. Mech. Engrs (Part J): J. Eng. Tribol. 1994, 208, 121–132. [Google Scholar] [CrossRef]
  47. Hooke, C.J.; Li, K.Y.; Morales-Espejel, G. Rapid Calculation of the Pressures and Clearances in Rough, Rolling-Sliding Elastohydrodynamically Lubricated Contacts. Part 1: Low-Amplitude, Sinusoidal Roughness. Proc. I. Mech. Engrs. Part C: J. Mech. Eng. Sci. 2007, 221, 535–550. [Google Scholar] [CrossRef]
  48. Lubrecht, A.A.; Graille, D.; Venner, C.H.; Greenwood, J.A. Waviness Amplitude Reduction in EHL Line Contacts Under Rolling-Sliding. J. Tribol. 1998, 120, 705–709. [Google Scholar] [CrossRef]
  49. Šperka, P.; Kǔpka, I.; Hartl, M. Experimental study of roughness effect in a rolling-sliding EHL contact. Part II: Complementary effects. Tribol. Trans. 2016, 59, 277–285. [Google Scholar] [CrossRef]
  50. Morales-Espejel, G.E.; Quiñonez, A.F. On the complementary function amplitude for film thickness in Micro-EHL. Tribol. Int. 2019, 131, 631–636. [Google Scholar] [CrossRef]
  51. Hansen, J.; Björling, M.; Larsson, R. Topography transformation due to running-in of rolling-sliding non-conformal surfaces. Tribol. Int. 2020, 144, 106126. [Google Scholar] [CrossRef]
  52. Cann, P.; Ionnides, E.; Jacobsen, B.; Lubrecht, A.A. The lambda ratio—A critical re-examination. Wear 1994, 175, 177–188. [Google Scholar] [CrossRef]
  53. Jacobson, B. Thin film lubrication of real surfaces. Tribol. Int. 2000, 33, 205–210. [Google Scholar] [CrossRef]
  54. Clarke, A.; Weeks, I.J.; Evans, H.P.; Snidle, R.W. An investigation into mixed lubrication conditions using electrical contact resistance techniques. Tribol. Int. 2016, 93, 709–716. [Google Scholar] [CrossRef]
  55. Zapletal, T.; Šperka, P.; Kǔpka, I.; Hartl, M. The effect of surface roughness on friction and film thickness in transition from EHL to mixed friction. Tribol. Int. 2018, 128, 356–364. [Google Scholar] [CrossRef]
  56. Crook, A.W. The lubrication of rollers Pt. III: A theoretical discussion of friction and the temperatures in the oil film. Philos.Trans. R. Soc. Lond. Ser. A 1961, 254, 237–258. [Google Scholar]
  57. Crook, A.W. The lubrication of rollers Pt. IV: Measurements of friction and effective viscosity. Philos. Trans. R. Soc. Lond. Ser. A 1963, 255, 281–312. [Google Scholar]
  58. Johnson, K.L.; Cameron, R. Shear behaviour of elastohydrodynamic oil films at high rolling contact pressure. Proc. Inst. Mech. Engrs 1967, 182, 12–24. [Google Scholar] [CrossRef]
  59. Johnson, K.L.; Greenwood, J.A. Thermal analysis of an Eyring fluid in elastohydrodynamic traction. Wear 1980, 61, 353–374. [Google Scholar] [CrossRef]
  60. Conry, T.F.; Wang, S.; Cusano, C. A Reynolds-Eyring Equation for Elastohydrodynamic Lubrication in Line Contacts. ASME J. Tribol. 1987, 109, 648–654. [Google Scholar] [CrossRef]
  61. Roelands, C.J.A.; Vlugter, J.C.; Waterman, H.I. The viscosity-pressure-temperature relationships of lubricating oils. Trans. ASME J. Basic Eng. 1963, 85, 601–607. [Google Scholar] [CrossRef]
  62. Yasutomi, S.; Bair, S.; Winer, W. An Application of a Free Volume Model to Lubricant Rheology. ASME J. Tribol. 1984, 106, 291–303. [Google Scholar] [CrossRef]
  63. Jacod, B.; Venner, C.H.; Lugt, P.M. A Generalized Traction Curve for EHL Contacts. ASME J. Tribol. 2001, 123, 248–253. [Google Scholar] [CrossRef] [Green Version]
  64. Lugt, P.M.; Morales-Espejel, G.E. A review of elastohydrodynamic lubrication theory. Tribol. Trans. 2011, 54, 470–496. [Google Scholar] [CrossRef]
  65. Jacod, B.; Venner, C.H.; Lugt, P.M. Extension of the Friction Mastercurve to Limiting Shear Stress Models. ASME J. Tribol. 2003, 125, 739–746. [Google Scholar] [CrossRef] [Green Version]
  66. Habchi, W.; Vergne, P.; Bair, S.; Andersson, O.; Eyheramendy, D.; Morales-Espejel, G.E. Influence of pressure and temperature dependence of thermal properties of a lubricant on the behaviour of circular TEHD contacts. Tribol. Int. 2010, 43, 1842–1850. [Google Scholar] [CrossRef]
  67. Liu, H.C.; Zhang, B.B.; Bader, N.; Li, X.M.; Terwey, T.; Poll, G. Fast Traction Prediction In Rolling/Sliding Ehl Contacts. In Proceedings of the 46th Leeds-Lyon Symposium on Tribology, Lyon, France, 2–4 September 2019. [Google Scholar]
  68. Habchi, W.; Bair, S.; Vergne, P. On friction regimes in quantitative elastohydrodynamics. Tribol. Int. 2013, 58, 107–117. [Google Scholar] [CrossRef]
Figure 1. Inlet flow pattern for pure rolling. (N.B. schematic: the stagnation point is where h = 3 h .
Figure 1. Inlet flow pattern for pure rolling. (N.B. schematic: the stagnation point is where h = 3 h .
Lubricants 08 00051 g001
Figure 2. Classical pressure distributions for roller lubrication.
Figure 2. Classical pressure distributions for roller lubrication.
Lubricants 08 00051 g002
Figure 3. Reynolds exit condition: effect of exponential viscosity dependence.
Figure 3. Reynolds exit condition: effect of exponential viscosity dependence.
Lubricants 08 00051 g003
Figure 4. Ertel and extended Ertel solutions, and analytical approximations.
Figure 4. Ertel and extended Ertel solutions, and analytical approximations.
Lubricants 08 00051 g004
Figure 5. An off-centre flat implies an exit constriction, with a minimum film thickness 0.8 h .
Figure 5. An off-centre flat implies an exit constriction, with a minimum film thickness 0.8 h .
Lubricants 08 00051 g005
Figure 6. Elastic pressure distributions for an off-centre flat.
Figure 6. Elastic pressure distributions for an off-centre flat.
Lubricants 08 00051 g006
Figure 7. Complete pressure distribution, showing the ignored hydrodynamic pressures.
Figure 7. Complete pressure distribution, showing the ignored hydrodynamic pressures.
Lubricants 08 00051 g007
Figure 8. Film thickness h from Ertel and extended Ertel theories.
Figure 8. Film thickness h from Ertel and extended Ertel theories.
Lubricants 08 00051 g008
Figure 9. Numerical solutions for minimum film thickness in line contacts. points: Venner [20]; lines: Dowson’s equation.
Figure 9. Numerical solutions for minimum film thickness in line contacts. points: Venner [20]; lines: Dowson’s equation.
Lubricants 08 00051 g009
Figure 10. Contours interpolated from the Venner minimum film thicknesses.
Figure 10. Contours interpolated from the Venner minimum film thicknesses.
Lubricants 08 00051 g010
Figure 11. Central film thickness for a circular point contact based on Venner’s data.
Figure 11. Central film thickness for a circular point contact based on Venner’s data.
Lubricants 08 00051 g011
Figure 12. Variation of film thickness with Hertz pressure.
Figure 12. Variation of film thickness with Hertz pressure.
Lubricants 08 00051 g012
Figure 13. Film thickness ratio h c e n / h min for circular contacts.
Figure 13. Film thickness ratio h c e n / h min for circular contacts.
Lubricants 08 00051 g013
Figure 14. Pressure and film thickness: smooth surface faster.
Figure 14. Pressure and film thickness: smooth surface faster.
Lubricants 08 00051 g014
Figure 15. Variation of shear stress with sliding speed (after Crook [54]).
Figure 15. Variation of shear stress with sliding speed (after Crook [54]).
Lubricants 08 00051 g015
Figure 16. Comparison of shear stresses for Newtonian and Eyring oils.
Figure 16. Comparison of shear stresses for Newtonian and Eyring oils.
Lubricants 08 00051 g016
Table 1. Hk calculated by LVC [22] for varying Ry/Rx but fixed speed and Hertz pressure.
Table 1. Hk calculated by LVC [22] for varying Ry/Rx but fixed speed and Hertz pressure.
Ry/Rx125102050100Mean
P = 7.9; S = 5.151.6311.6901.7111.6641.6621.6591.6571.668
P = 15.8; S = 10.31.4711.5361.5611.5531.5531.5451.5401.537

Share and Cite

MDPI and ACS Style

Greenwood, J.A. Elastohydrodynamic Lubrication. Lubricants 2020, 8, 51. https://doi.org/10.3390/lubricants8050051

AMA Style

Greenwood JA. Elastohydrodynamic Lubrication. Lubricants. 2020; 8(5):51. https://doi.org/10.3390/lubricants8050051

Chicago/Turabian Style

Greenwood, James A. 2020. "Elastohydrodynamic Lubrication" Lubricants 8, no. 5: 51. https://doi.org/10.3390/lubricants8050051

APA Style

Greenwood, J. A. (2020). Elastohydrodynamic Lubrication. Lubricants, 8(5), 51. https://doi.org/10.3390/lubricants8050051

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop