Three-dimensional rolling/sliding contact on a viscoelastic layered half-space

Abstract

This paper proposes an efficient contact model for a viscoelastic layered half-space where coating and substrate have different creep functions (i.e. different viscoelastic behaviours). The problem is formulated in three dimensions for an imposed pressure and shear field. The influence coefficients are calculated with the Papkovich-Neuber potentials and computation performed using Discrete Convolution and Fast Fourier Transform (DC-FFT) algorithms. From there, an elastic/viscoelastic correspondence is applied to move the elastic solution to a viscoelastic one. The Generalized Maxwell viscoelastic model is used and presented in the form of Prony series. A validation is performed using not only results within the literature but also by comparison with a finite element analysis. Next, a parametric study is done to analyse the influence of different parameters such as the elastic modulus ratio, the relaxation time ratio, the coating thickness or the rolling velocity. The model allows the analysis of the transient regime for which most of the known models struggle. It allows also to calculate 3D stresses and so to extract the stresses at the interface between the coating and the substrate so one can estimate the risk of failure such as delamination.