Static friction, differential algebraic systems and numerical stability

We show how Differential Algebraic Systems (Ordinary Differential Equations with algebraic constraints) in mechanics are affected by stability issues and we implement Lubich’s projection method to reduce the error to practically zero. Then, we explain how the “numerically exact” implementation for static friction by Differential Algebraic Systems can be stabilized. We conclude by comparing the corresponding steps in the “Contact mechanics” introduced by Moreau.