Numerical simulation of the dynamics of boats by a variational inequality approach