An optimization-based algorithm for Coulomb's frictional contact
ESAIM. Proceedings, Tome 27 (2009), pp. 54-69
Cet article a éte moissonné depuis la source EDP Sciences
The main goal of this paper is to propose a stable algorithm to compute friction forces governed by Coulomb's law in the course of the simulation of a nonsmooth Lagrangian dynamical system. The problem appears in computational mechanics, to simulate the dynamics of granular materials, robots, etc . Using a classical impulse-velocity formulation of Coulomb's law to model friction, and a semi-implicit time discretization scheme, we get a set of linear, non-linear and complementarity equations which has to be solved at each timestep. Two mutually dual parametric convex optimization problems coupled with a fixed point equation appear naturally. By solving (one of) these optimization problems iteratively within a damped nonsmooth-Newton algorithm, it is possible to decrease some infeasibility criterion and hopefully converge to a solution of the system. Numerical results are provided, which show that the number of iterations needed by the algorithm is very small in general and that the method is stable.
@article{EP_2009_27_a3,
author = {Florent Cadoux},
title = {An optimization-based algorithm for {Coulomb's} frictional contact},
journal = {ESAIM. Proceedings},
pages = {54--69},
year = {2009},
volume = {27},
doi = {10.1051/proc/2009019},
language = {en},
url = {http://geodesic.mathdoc.fr/articles/10.1051/proc/2009019/}
}
Florent Cadoux. An optimization-based algorithm for Coulomb's frictional contact. ESAIM. Proceedings, Tome 27 (2009), pp. 54-69. doi: 10.1051/proc/2009019
Cité par Sources :