Modeling and simulation of a planar rigid multibody system with multiple revolute clearance joints based on variational inequality