In the process of modelling of the Celtic stone rotating and rolling on a plane surface, different versions of simplified models of contact forces between two contacting bodies are used. The considered contact models take into account coupled dry friction force and torque, as well as rolling resistance. They were developed using Padé approximations, their modifications, and polynomial functions. Before the use of these models, some kinds of regularisations have been made, allowing to avoid singularities in the differential equations. The models were tested both numerically and experimentally, giving some practical guesses of the most essential elements of contact modelling in the Celtic stone numerical simulations. Since the tested contact models do not require the space discretisation, they can find application in relatively fast numerical simulations of rigid bodies with friction contacts.