Predicting wheel‐terrain interaction with semiempirical models is of substantial importance for developing planetary wheeled mobile robots (rovers). Primarily geared toward the design of manned terrestrial vehicles, conventional terramechanics models do not provide the sufficient fidelity required for application on autonomous planetary rovers. To develop a high‐fidelity interaction mechanics model, in this study the physical effects of wheel lug, slip sinkage, wheel dimension, and load are analyzed based on experimental results, including wheel sinkage, drawbar pull, normal force, and moment, which are measured on a single‐wheel test bed. The mechanism of lug‐terrain interaction is investigated systematically to clarify the principle of increasing shear stress, conditions of forming successive shearing among adjacent lugs, and the influence on shear displacement of soil. A mathematical model for predicting the concentrated forces and torque of rigid wheels with lugs for planetary rovers moving on sandy terrain is derived by integrating the improved models of normal and shearing stress distributions. In addition to the wheel parameters, terrain parameters, and motion state variables, wheel‐terrain interaction parameters, such as the linear varying sinkage exponent, the soil displacement radius, and load effect parameters, were proposed and explicitly included in the model. In the single‐wheel experiments, the slip ratio was increased approximately from 0.05 to 0.6, and the relative errors of the predicted results using the proposed model are less than 10% for all the wheels when compared with the experimental data. The proposed model has been used in the simulation of a four‐wheeled rover, and its effectiveness is evaluated by comparing the simulation results with experimental results.