Accumulative plastic deformation due to repeated loading is crucial to the lives of many mechanical components, such as gears, stamping dies, and rails in rail-wheel contacts. This paper presents a three-dimensional numerical model for simulating the repeated rolling or sliding contact of a rigid sphere over an elasto-plastic half-space. This model is a semi-analytical model based on the discrete convolution and fast Fourier transform algorithm. The half-space behaves either elastic-perfectly plastically or kinematic plastically. The analyses using this model result in histories of stress, strain, residual displacement, and plastic strain volume integral (PV) in the half-space. The model is examined through comparisons of the current results with those from the finite element method for a simple indentation test. The results of rolling contact obtained from four different hardening laws are presented when the load exceeds the theoretical shakedown limit. Shakedown and ratchetting behaviors are discussed in terms of the PV variation. The effect of friction coefficient on material responses to repeated sliding contacts is also investigated.

Interaction of nominally flat engineering surfaces that leads to a large contact area exists in many mechanical systems. Considering periodic similarity of surface geometry, a numerical three-dimensional elastoplastic contact model can be used to simulate the contact behaviors of two nominally flat surfaces with the assistance of the continuous convolution and Fourier transform (CC-FT) algorithm. This model utilizes the analytical frequency response functions (FRF) of elastic/plastic responses of materials and provides contact performance results, including the average surface gap, the contact area ratio, and the volume of plastically deformed material, which may be defined as performance variables. Following the digital filtration technology, rough surfaces can be numerically generated with specified autocorrelation length and the first four orders of statistical moments. A group of contact simulations are conducted with various working conditions. The effects of topographic and material properties on the contact behaviors are discussed. With a multivariables regression method, empirical formulas are developed for the performance variables as functions of surface statistical characteristics, material properties, a hardening parameter, and the applied load in terms of pressure.