Logo eprints

Optimization algorithms for the solution of the frictionless normal contact between rough surfaces

Bemporad, Alberto and Paggi, Marco Optimization algorithms for the solution of the frictionless normal contact between rough surfaces. International Journal of Solids and Structures, 69-70. pp. 94-105. ISSN 0020-7683 (2015)

Download (658kB) | Preview


This paper revisits the fundamental equations for the solution of the frictionless unilateral normal contact problem between a rough rigid surface and a linear elastic half-plane using the boundary element method (BEM). After recasting the resulting Linear Complementarity Problem (LCP) as a convex quadratic program (QP) with nonnegative constraints, different optimization algorithms are compared for its solution: (i) a Greedy method, based on different solvers for the unconstrained linear system (Conjugate Gradient CG, Gauss–Seidel, Cholesky factorization), (ii) a constrained CG algorithm, (iii) the Alternating Direction Method of Multipliers (ADMM), and (iv) the Non-Negative Least Squares (NNLS) algorithm, possibly warm-started by accelerated gradient projection steps or taking advantage of a loading history. The latter method is two orders of magnitude faster than the Greedy CG method and one order of magnitude faster than the constrained CG algorithm. Finally, we propose another type of warm start based on a refined criterion for the identification of the initial trial contact domain that can be used in conjunction with all the previous optimization algorithms. This method, called cascade multi-resolution (CMR), takes advantage of physical considerations regarding the scaling of the contact predictions by changing the surface resolution. The method is very efficient and accurate when applied to real or numerically generated rough surfaces, provided that their power spectral density function is of power-law type, as in case of self-affine fractal surfaces.

Item Type: Article
Identification Number: 10.1016/j.ijsolstr.2015.06.005
Projects: ERC StG CA2PVM
Funders: European Research Council
Uncontrolled Keywords: Unilateral contact problem; Frictionless normal contact; Quadratic programming; Optimization algorithms; Boundary element method; Roughness
Subjects: T Technology > TJ Mechanical engineering and machinery
Research Area: Computer Science and Applications
Depositing User: Prof Marco Paggi
Date Deposited: 08 Oct 2015 08:03
Last Modified: 08 Oct 2015 08:04
URI: http://eprints.imtlucca.it/id/eprint/2763

Actions (login required)

Edit Item Edit Item