The linear spring models, by virtue of their linearity, were unable to correctly model the dependence of CR on the impact speed. Nevertheless, their integrability allows fast numerical solution, which might make them preferable in some cases where speed takes precedence over accuracy. The logarithmic spring model is new, derived by analysis of the dynamics at very small Weber numbers and extension to the parameter regime of interest. It has three free parameters, which can be determined by measuring TC and CR for two different impact speeds. Only the logarithmic spring model was found to be consistent with the experimental data for Weber numbers smaller than 1. Even the bouncing threshold, which in general requires understanding of the intervening layer dynamics, could be almost entirely rationalized by considering mode transitions and contact time. For applications requiring a highly accurate representation of the temporal dependence of the reaction force, as will be the case for the model of walking drops treated in our subsequent paper (MBII), the logarithmic spring model will be invaluable.