TY - JOUR
T1 - Aspects of modeling and numerical simulation of dry point contacts between viscoelastic solids
AU - Zhao, Yan
AU - Morales-Espejel, G.E.
AU - Venner, C.H.
N1 - Funding Information:
The authors would like to thank Mr. Bernie van Leeuwen, SKF Research and Technology Development Director, for his kind permission to publish this article. The authors would like to thank Dr. Ehret Pascal for sharing his knowledge on this topic and thank Dr. Binbin Zhang and Dr. Haichao Liu for proofreading. The authors would thank Dr. C.W.Visser for his helpful discussion and comments. The first author acknowledges the China Scholarship Council (CSC) for providing the PhD scholarship.
Publisher Copyright:
© 2021 The Authors
PY - 2022/1
Y1 - 2022/1
N2 - The virtual absence of pressure driven flow inside an elastohydrodynamic lubrication (EHL) contact leads to remarkable similarity between lubricated film thickness and separation by a viscoelastic layer as is suggested by experimental and theoretical results, see [1]. This offers interesting opportunities for mixed lubrication modeling without needing the flow continuity-based Reynolds equation. In this paper numerical simulation of viscoelastic concentrated contact modeling is revisited, first considering the problem of a viscoelastic half-space and a rigid sphere, both in static (squeeze) as well as (frictionless) rolling/sliding conditions. In particular modeling aspects are discussed so that the efficiency of computation of the viscoelastic deformation remains close to the efficiency of computation of an elastic deformation. This is achieved by rewriting Hunter's [2] deformation equation from an integral form into a differential form by applying the Leibnitz integral theorem, even for complex materials with multiple relaxation times to represent complex viscoelastic properties of materials. The approach can be straightforwardly implemented in any (EHL) contact solver regardless of the numerical method. Detailed results are presented, such as the effect of the Deborah number, the relation between the extreme cases of the rolling problem and the static problem, and discussed in detail. In addition to the rigid sphere against a viscoelastic half-space, also the case of two viscoelastic bodies is considered.
AB - The virtual absence of pressure driven flow inside an elastohydrodynamic lubrication (EHL) contact leads to remarkable similarity between lubricated film thickness and separation by a viscoelastic layer as is suggested by experimental and theoretical results, see [1]. This offers interesting opportunities for mixed lubrication modeling without needing the flow continuity-based Reynolds equation. In this paper numerical simulation of viscoelastic concentrated contact modeling is revisited, first considering the problem of a viscoelastic half-space and a rigid sphere, both in static (squeeze) as well as (frictionless) rolling/sliding conditions. In particular modeling aspects are discussed so that the efficiency of computation of the viscoelastic deformation remains close to the efficiency of computation of an elastic deformation. This is achieved by rewriting Hunter's [2] deformation equation from an integral form into a differential form by applying the Leibnitz integral theorem, even for complex materials with multiple relaxation times to represent complex viscoelastic properties of materials. The approach can be straightforwardly implemented in any (EHL) contact solver regardless of the numerical method. Detailed results are presented, such as the effect of the Deborah number, the relation between the extreme cases of the rolling problem and the static problem, and discussed in detail. In addition to the rigid sphere against a viscoelastic half-space, also the case of two viscoelastic bodies is considered.
KW - Contact mechanics
KW - Lubrication modelling
KW - Multi-level multi-integration method (MLMI)
KW - Viscoelasticity
KW - UT-Hybrid-D
UR - http://www.scopus.com/inward/record.url?scp=85114737446&partnerID=8YFLogxK
U2 - 10.1016/j.triboint.2021.107245
DO - 10.1016/j.triboint.2021.107245
M3 - Article
AN - SCOPUS:85114737446
SN - 0301-679X
VL - 165
JO - Tribology international
JF - Tribology international
M1 - 107245
ER -