The tank-treading rotation of red blood cells (RBCs) in shear flows has been studied extensively with experimental, analytical, and numerical methods. Even for this relatively simple system, complicated motion and deformation behaviors have been observed, and some of the underlying mechanisms are still not well understood. In this study, we attempt to advance our knowledge of the relationship among cell motion, deformation, and flow situations with a numerical model. Our simulation results agree well with experimental data, and confirm the experimental finding of the decrease in frequency/shear-rate ratio with shear rate and the increase of frequency with suspending viscosity. Moreover, based on the detailed information from our simulations, we are able to interpret the frequency dependency on shear rate and suspending viscosity using a simple two-fluid shear model. The information obtained in this study thus is useful for understanding experimental observations of RBCs in shear and other flow situations; the good agreement to experimental measurements also shows the potential usefulness of our model for providing reliable results for microscopic blood flows.