The nonlinear dynamic behavior of superparamagnetic beads exposed to a periodic array of micromagnets and an external rotating field is simulated as a function of the relative size of the bead with respect to the micromagnet size and the strength of the external field relative to the pole density of the substrate. For large bead sizes, it is confirmed that the motion of the beads corresponds to the dynamics of an overdamped nonlinear harmonic oscillator. For lower bead sizes, additional subharmonic locking effects are observed along with the emergence of bounded orbits. These results qualitatively support previous experimental investigations of traveling-wave magnetophoresis and provide guidelines for achieving nearly infinite separation resolution between differently sized beads.