An FFT-based method for rough surface contact

Elastic contact between a rigid plane and a halfspace whose surface height is described by a bandwidth-limited Fourier series is considered. The surface normal displacements and contact pressures are found by a numerical technique that exploits the structure of the Fast Fourier Transform (FFT) and an exact result in linear elasticity. The multiscale nature of rough surface contact is implicit to the method, and features such as contact agglomeration and asperity interaction-a source of difficulty for asperity-based models-evolve naturally. Both two-dimensional (2-D) and three-dimensional (3-D) contact ape handled with equal ease. Finally, the implementation is simple, compact, and fast.