Assessing the dynamic behavior of railways vehicle by numerical simulations requires a consistent wheel rail contact model that gives accurate normal and tangential forces. The current contact theories implemented on most multi-body systems (MBS) consider the adhesion coefficient as constant in time during the simulation. This means that the wheel and the rail cleaning due to axles sliding and leading to regenerate adhesion is generally neglected. However, such an assumption is invalid for many applications especially for traction and braking where longitudinal creepage could be relatively high leading to important effect of cleaning and high adhesion recovery. In this paper, the frictional contact model of the multibody software VOCO is extended to take into account the adhesion recovery due to axles sliding in a braking context. The adhesion time-dependent model is based on the calculation of the power dissipated by friction at the wheel-rail contact level and was validated by comparisons with measurements provided from on-track braking tests performed with the Regio 2N ONO vehicle. The model has improved predictions of braking distance, axle speeds and train deceleration when performing braking simulations. A virtual Wheel Slide Protection (WSP) system that co-simulates with VOCO the dynamic response of the vehicle inside an exchange procedure was used in this work. The exchange procedure was designed to run in real-time in order to be further implemented inside a Hardware In the Loop (HIL) architecture with a real WSP device.