A numerical model of friction between rough surfaces

This paper describes a computational method to calculate the friction force between two rough surfaces. In the model used, friction results from forces developed during elastic deformation and shear resistance of adhesive junctions at the contact areas. Contacts occur between asperities and have arbitrary orientations with respect to the surfaces. The size and slope of each contact area depend on external loads, mechanical properties and topographies of surfaces. Contact force distribution is computed by iterating the relationship between contact parameters, external loads, and surface topographies until the sum of normal components of contact forces equals the normal load. The corresponding sum of tangential components of contact forces constitutes the friction force. To calculate elastic deformation in three dimensions, we use the method of influence coefficients and its adaptation to shear forces to account for sliding friction. Analysis presented in Appendix A gives approximate limits within which influence coefficients developed for flat elastic half-space can apply to rough surfaces. Use of the method of residual correction and a successive grid refinement helped rectify the periodicity error introduced by the FFT technique that was used to solve for asperity pressures. The proposed method, when applied to the classical problem of a sphere on a half-space as a benchmark, showed good agreement with previous results. Calculations show how friction changes with surface roughness and also demonstrate the method's efficiency.