Iteration framework for solving mixed lubrication computation problems

Abstract

The general discrete scheme of time-varying Reynolds equation loses the information of the previous step, which makes it unreasonable. A discretization formula of the Reynolds equation, which is based on the Crank–Nicolson method, is proposed considering the physical message of the previous step. Gauss–Seidel relaxation and distribution relaxation are adopted for the linear operators of pressure during the numerical solution procedure. In addition to the convergent criteria of pressure distribution and load, an estimation framework is developed to investigate the relative error of the most important term in the Reynolds equation. Smooth surface with full contacts and mixed elastohydrodynamic lubrication is tested for validation. The asperity contact and sinusoidal wavy surface are examined by the proposed discrete scheme. Results show the precipitous decline in the boundary of the contact area. The relative error suggests that the pressure distribution is reliable and reflects the accuracy and effectiveness of the developed method.

Graphical abstract

Keywords

mixed lubrication / discretization formula / relative error / Reynolds equation / asperity

References

[1] Morales-Espejel G E, Rycerz P, Kadiric A. Prediction of micropitting damage in gear teeth contacts considering the concurrent effects of surface fatigue and mild wear. Wear, 2018, 398–399: 99–115

[2] Brizmer V, Matta C, Nedelcu I, . The influence of tribolayer formation on tribological performance of rolling/sliding contacts. Tribology Letters, 2017, 65(2): 57 doi:10.1007/s11249-017-0839-3

[3] Brizmer V, Gabelli A, Vieillard C, . An experimental and theoretical study of hybrid bearing micropitting performance under reduced lubrication. Tribology Transactions, 2015, 58(5): 829–835

[4] Patir N, Cheng H S. An average flow model for determining effects of three-dimensional roughness on partial hydrodynamic lubrication. Journal of Tribology, 1978, 100(1): 12–17

[5] Geng K H, Geng A N, Wang X, . Frictional characteristics of the Vane–Chute pair in a rolling piston compressor based on the second-order motion. Tribology International, 2019, 133: 111–125

[6] Liu C, Lu Y, Wang P, . Numerical analysis of the effects of compression ring wear and cylinder liner deformation on the thermal mixed lubrication performance of ring-liner system. Mechanics & Industry, 2018, 19(2): 203 doi:10.1051/meca/2018015

[7] Cui S H, Gu L, Fillon M, . The effects of surface roughness on the transient characteristics of hydrodynamic cylindrical bearings during startup. Tribology International, 2018, 128: 421–428

[8] Zhang S W, Zhang C H, Hu Y Z, . Numerical simulation of mixed lubrication considering surface forces. Tribology International, 2019, 140: 105878

[9] Azam A, Ghanbarzadeh A, Neville A, . Modelling tribochemistry in the mixed lubrication regime. Tribology International, 2019, 132: 265–274 doi:10.1016/j.triboint.2018.12.024

[10] Jiang X, Hua D Y, Cheng H S, . A mixed elastohydrodynamic lubrication model with asperity contact. Journal of Tribology, 1999, 121(3): 481–491

[11] Hu Y Z, Zhu D. A full numerical solution to the mixed lubrication in point contacts. Journal of Tribology, 2000, 122(1): 1–9 doi:10.1115/1.555322

[12] Hu Y Z, Wang H, Wang W Z, . A computer model of mixed lubrication in point contacts. Tribology International, 2001, 34(1): 65–73

[13] Wang X P, Liu Y C, Zhu D. Numerical solution of mixed thermal elastohydrodynamic lubrication in point contacts with three-dimensional surface roughness. Journal of Tribology, 2017, 139(1): 011501

[14] He T, Zhu D, Yu C J, . Mixed elastohydrodynamic lubrication model for finite roller-coated half space interfaces. Tribology International, 2019, 134: 178–189

[15] Zhu D, Wang J, Ren N, . Mixed elastohydrodynamic lubrication in finite roller contacts involving realistic geometry and surface roughness. Journal of Tribology, 2012, 134(1): 011504

[16] Xie Z J, Xue Q H, Wu J Q, . Mixed-lubrication analysis of planetary roller screw. Tribology International, 2019, 140: 105883

[17] Lu X Q, Dong Q B, Zhou K, . Numerical analysis of transient elastohydrodynamic lubrication during startup and shutdown processes. Journal of Tribology, 2018, 140(4): 041504

[18] Xiao S B. Numerical analysis of transient thermal mixed lubrication in point contact. Thesis for Master’s Degree. Wuhan: Wuhan University of Science and Technology, 2019 (in Chinese)

[19] He Y C. Numerical analysis for the thermal EHL in point contacts of cam-tappet pair. Thesis for Master’s Degree. Qingdao: Qingdao University of Technology, 2019 (in Chinese)

[20] Marti J, Ryzhakov P. An explicit/implicit Runge–Kutta-based PFEM model for the simulation of thermally coupled incompressible flows. Computational Particle Mechanics, 2020, 7(1): 57–69

[21] Ezzatneshan E, Hejranfar K. Simulation of three-dimensional incompressible flows in generalized curvilinear coordinates using a high-order compact finite-difference lattice Boltzmann method. International Journal for Numerical Methods in Fluids, 2019, 89(7): 235–255

[22] Huang P. Lubrication Numerical Calculation Methods. Beijing: Higher Education Press, 2012

[23] Hamid Y, Usman A, Afaq S K, . Numeric based low viscosity adiabatic thermo-tribological performance analysis of piston-skirt liner system lubrication at high engine speed. Tribology International, 2018, 126: 166–176

[24] Wen S Z, Huang P. Principles of Tribology. 2nd ed. Beijing: Tsinghua University Press, 2002

[25] Zhu D. On some aspects of numerical solutions of thin-film and mixed elastohydrodynamic lubrication. Proceedings of the Institution of Mechanical Engineers. Part J, Journal of Engineering Tribology, 2007, 221(5): 561–579 doi:10.1243/13506501JET259

[26] Anderson J D. Computational Fluid Dynamics: The Basics with Applications. New York: McGraw-Hill, 1995

[27] Liao S K, Zhang Y, Chen D. Runge–Kutta finite element method based on the characteristic for the incompressible Navier–Stokes equations. Advances in Applied Mathematics and Mechanics, 2019, 11(6): 1415–1435 doi:10.4208/aamm.OA-2018-0150

[28] Fischels M V, Rajagopalan R G. Family of Runge–Kutta-based algorithms for unsteady incompressible flows. AIAA Journal, 2017, 55(8): 2630–2644

[29] Zhang X N, Glovnea R, Morales-Espejel G E, . The effect of working parameters upon elastohydrodynamic film thickness under periodic load variation. Tribology Letters, 2020, 68(2): 62

[30] Liu S B, Wang Q, Liu G. A versatile method of discrete convolution and FFT (DC-FFT) for contact analyses. Wear, 2000, 243(1–2): 101–111 doi:10.1016/S0043-1648(00)00427-0

[31] Bair S, Liu Y, Wang Q J. The pressure-viscosity coefficient for Newtonian EHL film thickness with general piezoviscous response. Journal of Tribology, 2006, 128(3): 624–631

[32] Venner C H, Lubrecht A A. Multilevel Methods in Lubrication. Amsterdam: Elsevier, 2000

[33] Lv F R, Jiao C X, Ta N, . Mixed-lubrication analysis of misaligned bearing considering turbulence. Tribology International, 2018, 119: 19–26

[34] Xiang G, Han Y F, Wang J X, . Coupling transient mixed lubrication and wear for journal bearing modeling. Tribology International, 2019, 138: 1–15

Acknowledgements

This study was financially supported by the National Natural Science Foundation of China (Grant No. U1637206), Shanghai Academy of Spaceflight Technology Projects (Grant Nos. SAST2017-079 and USCAST2019-25), and the State Key Laboratory of Mechanical System and Vibration (Grant No. MSVZD201912).

RIGHTS & PERMISSIONS

2021 Higher Education Press