Vibro-impact crawling robots driven by the emerging dielectric elastomer actuators (DEAs) feature the advantages of reduced system complexity and bidirectional locomotion capability. However, due to the lack of systematic dynamic models and in-depth investigations, the fundamental locomotion principles of the robots are still unclear and therefore their applications in real-world scenarios remain hindered. In this paper, a comprehensive dynamic model of this robot is developed by considering the complex interactions between the electro-mechanical coupling, impact mechanism, and multiple nonlinear friction characteristics. By incorporating extensive modeling and experimental studies, we explain the fundamental principles that lead to the bidirectional locomotion of the robot. The actuation strategies for both forward and backward locomotion are characterized in-depth, which include the actuation frequencies, relative phases, and waveforms. Three typical contact surface cases (i.e. rigid & dry, rigid & viscous, compliant & viscous) are considered in this paper, where we report the changes in the optimal actuation strategies for bidirectional locomotion with the contact surfaces and reveal the influences of key factors of friction coefficient, viscosity, and material compliance. The key findings reported in this work can build foundations for developing a highly robust and efficient DEA-driven vibro-impact crawling robot for broad applications.