In this article, the existing nonsmooth generalized- ? method for simulation of mechanical systems with frictionless contacts, modeled as unilateral constraints, is extended to frictional contacts. On that account, we complement constraints set-valued Coulomb-type friction laws. Moreover, devise a set benchmark systems, which can be used validate numerical schemes Finally, benchmarks numericall...