Research Papers: Elastohydrodynamic Lubrication

Coupling Strategies for Finite Element Modeling of Thermal Elastohydrodynamic Lubrication Problems

[+] Author and Article Information
W. Habchi

Department of Industrial and
Mechanical Engineering,
Lebanese American University,
Byblos, Lebanon
e-mail: wassim.habchi@lau.edu.lb

Contributed by the Tribology Division of ASME for publication in the JOURNAL OF TRIBOLOGY. Manuscript received May 3, 2016; final manuscript received September 24, 2016; published online April 4, 2017. Assoc. Editor: Xiaolan Ai.

J. Tribol 139(4), 041501 (Apr 04, 2017) (12 pages) Paper No: TRIB-16-1149; doi: 10.1115/1.4034956 History: Received May 03, 2016; Revised September 24, 2016

This paper investigates coupling strategies for finite element modeling (FEM) of thermal elastohydrodynamic lubrication (TEHL) problems. The TEHL problem involves a strong coupling between several physics: solid mechanics, fluid mechanics, and heat transfer. Customarily, this problem is split into two parts (elastohydrodynamic (EHD) and thermal) and the two problems are solved separately while an iterative procedure is established between their respective solutions. This weak coupling strategy involves a loss of information, as each problem is not made intimately aware of the evolution of the other problem's solution during the resolution procedure. This typically leads to slow convergence rates. The current work offers a full coupling strategy for the TEHL problem, i.e., both the EHD and thermal parts are solved simultaneously in a monolithic system. The system of equations is generated from a finite element discretization of the governing field variables: hydrodynamic pressure, solids elastic deformation, and temperature. The full coupling strategy prevents any loss of information during the resolution procedure leading to very fast convergence rates (solution is attained within a few iterations only). The performance of the full coupling strategy is compared to that of different weak coupling strategies. Out of simplicity, only steady-state line contacts are considered in this work. Nevertheless, the proposed methodology, results, and findings are of a general nature and may be extrapolated to circular or elliptical contacts under steady-state or transient conditions.

Copyright © 2017 by ASME
Your Session has timed out. Please sign back in to continue.



Grahic Jump Location
Fig. 1

Geometry of a line contact

Grahic Jump Location
Fig. 2

Computational domain of the EHL part

Grahic Jump Location
Fig. 3

Computational domain of the thermal part

Grahic Jump Location
Fig. 4

“Extra coarse” mesh case for the EHL domain (a) and the thermal domain (b)

Grahic Jump Location
Fig. 5

Film thickness and temperature convergence with respect to mesh size. Left: moderately loaded case (ph = 0.7 GPa) and right: highly loaded case (ph = 2.2 GPa).

Grahic Jump Location
Fig. 6

Dimensionless pressure and film thickness distribution (left), temperature rise over the plane's surface (Z = 0), the midlayer of the lubricant film (Z = 0.5), and the cylinder's surface (Z = 1) (middle) and dimensionless lubricant shear stress distribution on the plane's surface over the contact width (right) for the typical test case

Grahic Jump Location
Fig. 7

Flowcharts of weak coupling strategies




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In