Legged locomotion of robots has advantages in reducing payload in contexts such as travel over deserts or in planet surfaces. A recent study (Li et al. 2013) partially addresses this issue by examining legged locomotion over granular media (GM). However, they miss one extremely significant fact. When the robot’s wheels (legs) run over GM, the granules are set into motion. Hence, unlike the study of Li et al. (2013), the viscosity of the GM must be included to simulate the kinematic energy loss in striking and passing through the GM. Here the locomotion in their experiments is re-examined using an advanced Navier-Stokes framework with a parameterized granular viscosity. It is found that the performance efficiency of a robot, measured by the maximum speed attainable, follows a six-parameter sigmoid curve when plotted against rotating frequency. A correct scaling for the turning point of the sigmoid curve involves the footprint size, rotation frequency and weight of the robot. Our proposed granular response to a load, or the ‘influencing domain’ concept points out that there is no hydrostatic balance within granular material. The balance is a synergic action of multi-body solids. A solid (of whatever density) may stay in equilibrium at an arbitrary depth inside the GM. It is shown that there exists only a minimum set-in depth and there is no maximum or optimal depth. The set-in depth of a moving robot is a combination of its weight, footprint, thrusting/stroking frequency, surface property of the legs against GM with which it has direct contact, and internal mechanical properties of the GM. If the vehicle’s working environment is known, the wheel-granular interaction and the granular mechanical properties can be grouped together. The unitless combination of the other three can form invariants to scale the performance of various designs of wheels/legs. Wider wheel/leg widths increase the maximum achievable speed if all other parameters are unchanged.