The lateral resistance of ballast bed is an important parameter to prevent track expansion and maintain track stability. The invasion of sand particles can cause the change of lateral resistance of ballast bed and affect the stability of track structure, but little attention has been paid to the change characteristics of nonlinear lateral resistance of sandy ballast bed. In this paper, the field tests on the lateral resistance of windblown sand ballast bed were carried out to establish a multiscale three-dimensional discrete element model of sleeper-ballast bed. A systematic analysis on the evolution of lateral resistance, resistance to lateral deformation, micro-contact characteristics and lateral stability of ballast bed is performed. The results show that sand intrusion can increase the lateral resistance of ballast bed, which is approximately 40 % higher than that of clean ballast bed. In nonlinear strengthening stage and yield stage, the enhancement effect of sand particles on the lateral resistance of ballast bed is relatively weaker in comparison with the linear growth stage. With the increase in sand intrusion depth, the lateral resistance and resistance work of ballast bed both gradually go up, and the contribution of ballast shoulder to lateral resistance tends to play a leading role. Sand intrusion can increase the lateral stiffness of ballast bed and reduce the elasticity of track structure. Therefore, the maintenance operation should be carried out in time for the section with severe sandstorm.