There are mainly two kinds of methods to model the unilateral constraints. The first kind is based on smooth contact dynamics, including methods using Hertz's models, penalty methods, and some regularization force models, while the second kind is based on the non-smooth contact dynamics, which models the system with unilateral contacts as variational inequalities.