The derailment of trains is a complex phenomenon that requires an elaborate contact model in simulation to better understand its mechanism. The CONTACT program is a well-known reference for wheel-rail contact modeling due to its high accuracy. However, its low computational efficiency restricts its applications especially in the context of a multi-body simulation. Therefore, a high computational efficient, simplified and approximate non-Hertzian contact is preferred in derailment simulation. The aim of this research is to verify the efficiency of a recently developed non-Hertzian wheel-rail contact model in derailment simulation, which is a combination of the Kik-Piotrowski model and the KBTNH that is a fast creep force solver for non-Hertzian contacts. To assess the performance of the non-Hertzian model in derailment simulation, the derailment coefficient for steady-state and quasi-steady conditions, the wheel/rail contact forces during flange contact, and the dynamics behaviors of the wheelset prior to the derailment are compared with the state of the art contact methods representing different levels of modeling complexity, accuracy and efficiency, namely the classical approach (Hertz theory+FASTSIM algorithm) and the ‘exact’ solver CONTACT.