On some aspects of numerical solutions of thin-film and mixed elastohydrodynamic lubrication

Abstract

Numerical solution of the thin-film and mixed elastohydrodynamic lubrication (EHL) is of great importance to the study of lubrication transition and breakdown, as well as surface failures caused by contact. However, in this region, converged and accurate numerical solutions become difficult, and effects of computational mesh density and differential schemes appear to be more significant. Also, there is currently a debate on how surface contact should be defined and analysed, and whether it is indeed possible to model contact through a grid-converged solution of the EHL equation system. This paper presents a set of analysed sample cases with different computational meshes, different cut-off values for handling contact, on different levels of central film thickness, from several hundred nanometres down to zero. Based on the results, some key issues are discussed (such as whether the contact can be modelled with the Reynolds equation system, how a contact should be defined and analysed, how computational meshes should be selected and managed, what the reasonable and satisfactory convergence accuracy should be, and so on). In addition, a progressive mesh densification (PMD) method is presented and compared with the multi-grid (MG) method. It is demonstrated that the PMD method may be able to significantly reduce computing time. On the other hand, the MG method may be good for thick-film cases with smooth surfaces or low frequency roughness, but is not desirable for the ultra-thin film and mixed EHL, especially when moderate and high frequency surface roughness is involved.

Keywords

elastohydrodynamic lubrication mixed lubrication surface contact numerical solution film thickness thin film behaviour

Get full access to this article

View all access options for this article.

References

Ranger

A. P.

Ettles

C. M. M.

Cameron

The solution of the point contact elastohydrodynamic problem

Proc. R. Soc. Lond., 1975, A346, 227–244.

Hamrock

B. J.

Dowson

Isothermal elastohydrodynamic lubrication of point contacts part 1 - Theoretical formulation

ASME J. Lubr. Technol., 1976, 98, 223–229.

K. P.

Rohde

S. M.

Numerical solution of the point contact problem using the finite element method

Int. J. Numer. Meth. Eng., 1977, 11, 1507–1518.

Evans

H. P.

Snidle

R. W.

Inverse solution of Reynolds equation of lubrication under point contact elastohydrodynamic conditions

ASME J. Lubr. Technol., 1981, 103, 539–546.

Zhu

Wen

S. Z.

A full numerical solution for the thermo-elastohydrodynamic problem in elliptical contacts

ASME J. Tribol., 1984, 106, 246–254.

Lubrecht

A. A.

The numerical solution of elastohydrodynamic lubricated line and point contact problems using multigrid techniques. PhD Thesis, University of Twente, The Netherlands, 1987.

Venner

C. H.

Multilevel solution of the EHL line and point contact problems. PhD Thesis, University of Twente, The Netherlands, 1991.

X. L.

Numerical analysis of elastohydrodynamically lubricated line and point contacts with rough surfaces by using semi-system and multigrid methods. PhD Thesis, Northwestern University, IL, USA, 1993.

Sadeghi

Thermal EHL analysis of circular contacts with measured surface roughness

ASME J. Tribol., 1996, 118, 473–483.

10.

Zhu

X. L.

Point contact EHL based on optically measured 3-Dimensional rough surfaces

ASME J. Tribol., 1997, 119, 375–384.

11.

Zhu

Y. Z.

The study of transition from full film elastohydrodynamic to mixed and boundary lubrication. The advancing frontier of engineering tribology. In Proceedings of the 1999 STLE/ASME H. S. Cheng

Tribology surveillance, October 1999, pp. 150–156 (STLE, Park Ridge, IL).

12.

Y. Z.

Zhu

A full numerical solution to the mixed lubrication in point contacts

ASME J. Tribol., 2000, 122, 1–9.

13.

Zhu

Y. Z.

A computer program package for prediction of EHL and mixed lubrication characteristics, friction, subsurface stresses and flash temperatures based on measured 3-D surface roughness

Tribol. Trans., 2001, 44, 383–390.

14.

Zhu

A design tool for selection and optimization of surface finish in mixed lubrication.

In Proceedings of the 29th Leeds-LyonSymposium on Tribology, 2002, pp. 703–711 (Elsevier B.V., Amsterdam).

15.

Morales-Espejel

G. E.

Dumont

M. L.

Lugt

P. M.

Olver

A. V.

A limiting solution for the dependence of film thickness on velocity in EHL contacts with very thin films

Tribol. Trans., 2005, 48, 317–324.

16.

Zhu

Open discussion in response to Morales-Espejel et al. [15]

Tribol. Trans., 2005, 48, 324–327.

17.

Liu

Y. C.

Wang

W. Z.

Y. Z.

Zhu

Effects of differential scheme and mesh density on EHL film thickness in point contacts. Presented at the STLE 2005 Annual Meeting, Las Vegas, May 2005, ASME J. Tribol., 2005, 128, 641–653.

18.

Venner

C. H.

EHL film thickness computations at low speeds: Risk of artificial trends as a result of poor accuracy and implications for mixed lubrication modeling

Proc. IMechE, Part J: J. Engineering Tribology, 2005, 219, 285–290.

19.

Venner

C. H.

Lubrecht

A. A.

Multilevel methods in lubrication, Elsevier Tribology Series, 2000, vol. 37 (Elsevier B.V., Amsterdam).

20.

Zhao

J. X.

Sadeghi

Hoeprich

M. H.

Analysis of EHL circular contact start up: Part I - Mixed contact model with pressure and film thickness results

ASME J. Tribol., 2001, 123, 67–74.

21.

Holmes

M. J. A.

Evans

H. P.

Snidle

R. W.

Comparison of transient EHL calculations with a start-up experiment.

In Proceedings of the 29th Leeds-Lyon Symposium on Tribology, 2002, pp. 79–89 (Elsevier B.V., Amsterdam).

22.

Holmes

M. J. A.

Evans

H. P.

Snidle

R. W.

Analysis of mixed lubrication effects in simulated gear tooth contacts

ASME J. Tribol., 2005, 127, 61–69.

23.

Popovici

Venner

C. H.

Lugt

P. M.

Effects of load system dynamics on the film thickness in EHL contacts during start up

J. Tribol., 2004, 126, 258–266.

24.

Zhu

Hartl

Krupka

A comparative study on the EHL film thickness results from both numerical solutions and experimental measurements. Presented at the STLE Annual Meeting, Toronto, Canada, 19 May 2004.

25.

Liu

Wang

Liu

A versatile method of discrete convolution and FFT (DC-FFT) for contact analyses

Wear, 2000, 243, 101–111.

26.

Liu

Wang

A Three-dimensional thermomechanical model of contact between non-conforming rough surfaces

ASME J. Tribol., 2001, 123, 17–26.

27.

Wang

W. Z.

Wang

Liu

Y. C.

Y. Z.

Zhu

A comparative study of the methods for calculation of surface elastic deformation

Proc. Instn Mech. Engrs, Part J: J. Engineering Tribology, 2003, 217, 145–153.

28.

Booser

E. R.

(Ed.) Tribology data handbook, 1997 (CRC Press LLC, New York).

29.

Evans

H. P.

Snidle

R. W.

The isothermal elastohydrodynamic lubrication of spheres

ASME J. Lubr. Technol., 1981, 103, 547–557.

30.

Chittenden

R. J.

Dowson

Dunn

J. F.

Taylor

C. M.

A theoretical analysis of the isothermal Elastohydrodynamic lubrication of concentrated contacts (Part 1 and Part 2)

Proc. R. Soc. Lond., 1985, A397, 245–294.

31.

Venner

C. H.

Lubrecht

A. A.

Numerical simulation of a transverse ridge in a circular EHL contact

ASME J. Tribol., 1994, 116, 751–761.

32.

Chevalier

Lubrecht

A. A.

Film thickness in starved EHL point contacts

ASME J. Tribol., 1998, 120, 126–133.

33.

Y. Z.

Barber

G. C.

Zhu

Numerical analysis for the elastic contact of real rough surfaces

Tribol. Trans., 1999, 42 (3), 443–452.

34.

35.

Xiong

S. W.

Wang

Liu

W. K.

Yang

Q. M.

Vaidyanathan

Zhu

Approaching mixed elastohydrodynamic lubrication of smooth journal bearing systems with low rotating speed

Tribol. Trans., 2006, 49, 598–610.