Numerical Solution Program for Elastohydrodynamic Lubrication in Journal Bearings: Pressure, Film Thickness, and More

Numerical Solution Methods for Elastohydrodynamic Lubrication Problems in Journal Bearings

Elastohydrodynamic Lubrication (EHL) serves as the fundamental theory for studying lubricant film behavior in mechanical components like sliding bearings under high-pressure conditions. The core challenge involves coupled solutions of the following physical fields:

Reynolds Equation This partial differential equation describes the pressure distribution within the lubricant film, requiring consideration of viscosity-pressure relationships (e.g., Barus formula) and density variations with pressure. For steady-state conditions, it simplifies to a two-dimensional form, while transient analysis necessitates incorporating time-dependent terms. In code implementation, this typically involves finite difference discretization with appropriate boundary conditions (e.g., ambient pressure at edges, cavitation boundaries).

Elastic Deformation Equation Surface elastic deformation under high-pressure film loading is commonly calculated using Boussinesq integrals or Fast Fourier Transform (FFT) methods. Young's modulus and Poisson's ratio of metallic materials are critical input parameters. Programmatically, FFT-based approaches significantly accelerate deformation computations through convolution theorem applications, often implemented with zero-padding techniques to avoid cyclic errors.

Film Thickness Equation This equation defines the geometric lubricant film thickness, incorporating initial bearing surface geometry (such as arc radius), elastic deformation amounts, and potential thermal expansion components. Minimum film thickness directly determines bearing wear resistance. Algorithmically, this requires precise coordinate mapping between deformed surfaces, often handled through influence coefficient matrices in numerical schemes.

Numerical Solution Strategy Multigrid methods are frequently employed to accelerate convergence, with pressure and film thickness fields solved through iterative coupling. Relaxation factor selection is crucial for stability, following a typical workflow: assume initial film thickness → solve Reynolds equation → update deformations → recalculate film thickness → check convergence criteria. Code implementations often feature adaptive relaxation factors based on residual norms, with convergence thresholds typically set between 10^-4 to 10^-6 for relative errors.

Industrial software implementations typically include material databases (lubricant viscosity-temperature characteristic curves), parallel computing modules (for handling 3D contact problems), and post-processing visualization capabilities (pressure contour plots, film thickness isosurfaces). Key outputs include pressure peak locations, secondary pressure peak phenomena, and minimum film thickness variations with load/speed parameters. Modern implementations often integrate Python/Matlab scripting interfaces for parametric studies and C++/Fortran kernels for computationally intensive simulations.