We focus on the challenging problem of simulating thin elastic rods in contact, in the presence of friction. Most previous approaches in computer graphics rely on a linear complementarity formulation for handling contact in a stable way, and approximate Coulombs's friction law for making the problem tractable. In contrast, following the seminal work by Alart and Curnier in contact mechanics, we simultaneously model contact and exact Coulomb friction as a zero finding problem of a nonsmooth function. A semi-implicit time-stepping scheme is then employed to discretize the dynamics of rods constrained by frictional contact: this leads to a set of linear equations subject to an equality constraint involving a nondifferentiable function. To solve this one-step problem we introduce a simple and practical nonsmooth Newton algorithm which proves to be reasonably efficient and robust for systems that are not overconstrained. We show that our method is able to finely capture the subtle effects that occur when thin elastic rods with various geometries enter into contact, such as stick-slip instabilities in free configurations, entangling curls, resting contacts in braid-like structures, or the formation of tight knots under large constraints. Our method can be viewed as a first step towards the accurate modeling of dynamic fibrous materials.

The numerical integration of switching circuits is known to be a tough issue when the number of switches is large, or when sliding modes exist. Then, classical analog simulators may behave poorly, or even fail. In this paper, it is shown on two examples that the nonsmooth dynamical systems (NSDS) approach, which is made of: 1) a specific modeling of the piecewise-linear electronic devices (ideal diodes, Zener diodes, transistors); 2) the Moreau's time-stepping scheme; and 3) specific iterative one-step solvers, supersedes simulators of the simulation program with integrated circuit emphasis (SPICE) family and hybrid simulators. An academic example constructed in [Maffezzoni,, IEEE Trans. CADICS, vol 25, no. 11, Nov. 2006], so that the Newton–Raphson scheme does not converge, and the buck converter are used to make extensive comparisons between the NSDS method and other methods of the SPICE family and a hybrid-like method. The NSDS method, implemented in thesiconosplatform developed at INRIA, proves to be on these two examples much faster and more robust with respect to the model parameter variations.