Railway rails accumulate large plastic deformations due to cyclic rolling contact loading. The plastic deformations alter the rail geometry, affect material behavior, and cause crack formation and growth. The complex interactions between these phenomena require high fidelity simulations to be understood. 3d finite element simulations are accurate, but their computational cost limits the possible number of simulated cycles. We propose a cyclic finite element simulation in which the wheel and rail remain in contact throughout the simulation. It uses periodic boundary conditions, shadow elements, and model reductions. Compared to previous work, it is 25 times faster. The method is available as an open-source plugin to Abaqus, enabling other researchers to study rolling contact loading coupled with large plastic deformations.