This technical brief revisits the method outlined in Tasora and Anitescu 2011 [“A Matrix-Free Cone Complementarity Approach for Solving Large-Scale, Nonsmooth, Rigid Body Dynamics,” Comput. Methods Appl. Mech. Eng., **200**(5–8), pp. 439–453], which was introduced to solve the rigid multibody dynamics problem in the presence of friction and contact. The discretized equations of motion obtained here are identical to the ones in Tasora and Anitescu 2011 [“A Matrix-Free Cone Complementarity Approach for Solving Large-Scale, Nonsmooth, Rigid Body Dynamics,” Comput. Methods Appl. Mech. Eng., **200**(5–8), pp. 439–453]; what is different is the process of obtaining these equations. Instead of using maximum dissipation conditions as the basis for the Coulomb friction model, the approach detailed uses complementarity conditions that combine with contact unilateral constraints to augment the classical index-3 differential algebraic equations of multibody dynamics. The resulting set of differential, algebraic, and complementarity equations is relaxed after time discretization to a cone complementarity problem (CCP) whose solution represents the first-order optimality condition of a quadratic program with conic constraints. The method discussed herein has proven reliable in handling large frictional contact problems. Recently, it has been used with promising results in fluid–solid interaction applications. Alas, this solution is not perfect, and it is hoped that the detailed account provided herein will serve as a starting point for future improvements.