Ignoring the non linearties the contact force can be easily and accurately can be calculated using Lagrangian Method. In this method you have to calculate the position vector of contact point at wheel/rail from fixed frame of reference.
In this you can account linear Kalker theory for formulation. It is linear analysis, valid for small creepages, therefore may not be exact. However it is a convenient way of calculation.