a matrix - free cone complementarity approach for solving large - scale , nonsmooth , rigid body dynamics preprint anl / mcs - p 1692 - 1109 .
large - scale parallel multibody dynamics with frictional contact on the graphical processing unit .