Picture of Open Access badges

Discover Open Access research at Strathprints

It's International Open Access Week, 24-30 October 2016. This year's theme is "Open in Action" and is all about taking meaningful steps towards opening up research and scholarship. The Strathprints institutional repository is a digital archive of University of Strathclyde research outputs. Explore recent world leading Open Access research content by University of Strathclyde researchers and see how Strathclyde researchers are committing to putting "Open in Action".


Image: h_pampel, CC-BY

The Hertz contact problem, coupled Volterra integral equations and a linear complementarity problem

Gauthier, A. and Knight, P.A. and McKee, S. (2007) The Hertz contact problem, coupled Volterra integral equations and a linear complementarity problem. Journal of Computational and Applied Mathematics, 206 (1). pp. 322-340. ISSN 0377-0427

Full text not available in this repository. (Request a copy from the Strathclyde author)


This paper is concerned with the indentation of an elastic half-space by an axisymmetric punch under a monotonically applied normal force and under the assumption of Coulomb friction with coefficient μ in the region of contact. Within an inner (unknown) circle the contact is adhesive, while in the surrounding annulus the surface moves inwards with increasing load. In this paper it is shown how this problem is equivalent to two coupled Abel's equations with an unknown free point, the inner circumference of the annulus. It is further shown that a product integration finite difference approximation of those integral equations leads to a mixed linear complementarity problem (mixed LCP). A method based on Newton's method for solving non-smooth nonlinear equations is demonstrated to converge under restrictive assumptions on the physical parameters defining the system; and numerical experimentation verifies that it has much wider applicability. The method is also validated against the approach of Spence. The advantage of the mixed LCP formulation is that it provides the radius of the inner adhesive circle directly using the physical parameters of the problem.