Derivation of famous analytical solution of EHL equations is presented here. In 1949, Grubin obtained a solution for so called elasto-hydrodynamic lubrication (EHL) line contact problem with simplifications known as Grubin approximation. He was the first to combine both elastic deformation and lubricant hydrodynamic flow, including pressure-viscosity dependence. Although his solution did not satisfy strictly hydrodynamic equations of EHL, his analysis was recognized as very useful under highly loaded conditions, not only in stationary case, but also in transient conditions. Corresponding MATLAB code can be found here.

Here, the solution approach is described for the case of stationary line contact (1D Reynolds equation). The base for the validity of Grubin’s approximation is the fact that at highly loaded contacts, the hydrodynamic pressure is almost the same as Hertz contact pressure obtained for dry contacts (cylinder on plane). On top of it, the film thickness is almost constant within Hertz contact area, as shown in Fig. 1.

Following system of equations is considered:

(1)

where are hydrodynamic film thickness, pressure, viscosity, and pressure-viscosity coefficient correspondingly. Integration of Reynolds equation gives:

where C is an integration constant. Since the pressure drops to zero in the infinity, there exist a point , so that . Hence,, where . Following expression can be obtained from the integrated Reynolds equation:

(2)

Further, a following variable is introduced: . Then, the equation is true:

(3)

When pressure is very high, and .

From this it follows that and film is constant within the Hertz contact area. It is clear from Fig. 1 that outside the Hertz area, the hydrodynamic pressures are small compared to those within the contact and therefore it can be assumed that all the elastic defection is due to Herztian pressure. Therefore, the film thickness can be written in the form:

(4)

where is the combined elastic deflection of the surfaces (given by Hertz solution), is the applied load per unit length of the cylinder, is the combined elastic modulus and is given by the following expression:

(5)

From the expression of the film thickness and hence (equation 2):

(6)

By introducing dimensionless variables and , the equation reduces to the following form:

(7)

Integration of this equation from to gives . This equation was integrated by Grubin numerically for various values of and further fitted with the following simple formula:

(8)

Now, for the line contact, with being equivalent radius of the bodies. At in dimensional coordinates . Hence,

(9)

Substituting the expression for dimensionless pressure one can obtain following:

(10)

And finally, the film thickness can be expressed as

(11)