- 1Institute for Robotics and Intelligent Machines, Georgia Institute of Technology, Atlanta, GA, United States
- 2George W. Woodruff School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, GA, United States
- 3School of Physics, Georgia Institute of Technology, Atlanta, GA, United States
- 4The Wallace H. Coulter Department of Biomedical Engineering, Georgia Institute of Technology, Atlanta, GA, United States
Much of the Earth and many surfaces of extraterrestrial bodies are composed of non-cohesive particulate matter. Locomoting on such granular terrain is challenging for common robotic devices, either wheeled or legged. In this work, we discover a robust alternative locomotion mechanism on granular media-generating movement via self-vibration. To demonstrate the effectiveness of this locomotion mechanism, we develop a cube-shaped robot with an embedded vibratory motor and conduct systematic experiments on granular terrains of various particle properties and slopes. We investigate how locomotion changes as a function of vibration frequency/intensity on such granular terrains. Compared to hard surfaces, we find such a vibratory locomotion mechanism enables the robot to move faster, and more stably on granular surfaces, facilitated by the interaction between the body and surrounding grains. We develop a numerical simulation of a vibrating single cube on granular media, enabling us to justify our hypothesis that the cube achieves locomotion through the oscillations excited at a distance from the cube’s center of mass. The simplicity in structural design and controls of this robotic system indicates that vibratory locomotion can be a valuable alternative way to produce robust locomotion on granular terrains. We further demonstrate that such cube-shaped robots can be used as modular units for vibratory robots with capabilities of maneuverable forward and turning motions, showing potential practical scenarios for robotic systems.
1 Introduction
Many terrestrial and extraterrestrial landmasses are composed of soft flowable material, like sand and snow. Among them, granular media consists of discrete solid particles that dissipate energy upon interaction. Such terrain poses locomotor challenges for conventional wheeled and tracked devices. To surmount this challenge, diverse robotic systems have been developed to produce controllable movement in granular materials. These systems encompass legged configurations (Liang et al., 2012; Qian et al., 2015; Zhong et al., 2018; Lee et al., 2022), limbless structures (Maladen et al., 2011; Marvi et al., 2014), and rover designs (Knuth et al., 2012; Shrivastava et al., 2020). A common issue faced by such robotic systems attempting locomotion within granular media lies in the potential transition of the granular terrain from a solid to a flowing state, which occurs when the force per unit area exerted by the robot exceeds its yield stress. The interaction between the robot body and the terrain often engenders a coupled relationship encompassing the robot’s movement, the resistive force the robot experiences, and the dynamic changing of the terrain state before, during, and after the robot’s movement. Particularly in cases involving shear-based locomotion, the granular material can aggregate into formations that act as obstacles, hindering the robot’s progress (Shrivastava et al., 2020).
A vibration-driven robot leverages periodic oscillations or rotations of internal components to generate motion throughout the entire body (Becker et al., 2011; Reis et al., 2013; Calisti et al., 2017; Golitsyna, 2018). In contrast to conventional wheeled, legged, or elongated limbless robots, vibration-driven robots offer a simpler and more compact body design, without the need for body deformation and exposed actuators or joints. This intrinsic feature effectively prevents potential damage to the robot arising from particle infiltration or obstruction of the actuation mechanisms. Further, such design allows for a greater contact area between the robot’s body and the terrain. Therefore, the yield strain and stress for supporting body weight and producing propulsion can be reduced; consequently, the risk of terrain collapse is mitigated.
Previous research has designed and studied vibration-driven robots that can locomote on various terrains, such as hard surfaces (Rubenstein et al., 2012; Notomista et al., 2019; Li et al., 2021), fluid surfaces (Cocuzza et al., 2021; Tarr et al., 2024), and amphibious environments (Tang et al., 2018; Wang et al., 2022). These vibratory robots typically utilize a consistent actuation principle, wherein vibratory actuators induce multi-directional movement through alterations in force orientation, and the accumulated movement direction is dictated by anisotropic friction or resistive forces. These forces stem from specific directional components of periodic internal forces or inherent body features (Calisti et al., 2017; Golitsyna, 2018). However, the performance of such a vibration-driven mechanism on granular media remains unexplored.
Due to its design simplicity, ease of control, and enlarged body-terrain contact area, a vibratory robot presents a promising avenue for achieving robust and effective locomotion on granular media. This work introduces the development of a novel vibrating robotic system (Figure 1) to investigate the locomotion capabilities of such vibratory mechanisms within different granular media and identify potential practical applications.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1298676/frobt-11-1298676-HTML-r1/image_m/frobt-11-1298676-g001.jpg)
Figure 1. Vibratory single cube and bi-cube fabrication and motion. Panels (A,C,E) show the appearance, perspective drawing, and the turning maneuver of a single cube. By combining two single cubes together, the bi-cube is fabricated. Panels (B,D,F) show the appearance, perspective drawing, and the forward-moving maneuver of a bi-cube. The locomotion of single cube and bi-cube are recorded in Supplementary Video S1.
2 Materials and methods
2.1 Robot design and maneuver manipulation
2.1.1 Single cube
For design simplicity, we used Solidworks to create a cubic box with an open side for the motor and a press-fit lid to close the cube (as shown in Figure 1). The 4-cm-side-length cube was then 3D printed using a LulzBot TAZ Workhorse 3D Printer and polylactic acid (PLA) as the printing material. To generate the vibration, we selected an eccentric rotating mass vibration motor (ERM), specifically the VJQ24 from Vybronics Inc. (motor weight 31 g, 2,550 rpm at 5 V DC). The rotary axis of the vibration motor lies horizontally in the x direction, parallel to the cube’s bottom surface. When power is applied, the uneven mass starts rotating, which leads to rotary oscillation about the x-axis. By switching the voltage from positive to negative, the vibration motor rotation direction converts from counter-clockwise to clockwise (view in the positive x direction), which allows the single cube and bi-cube to generate maneuvers. The position of the vibration motor is adjusted to the cube’s center, guaranteeing alignment between both the center of the cube’s mass and cube’s geometric center. Inside the cube, we implemented a structure to securely press fit the motor in place to ensure that it remains stationary while vibrating.
Influenced by the simultaneous lateral oscillatory inertia generated by the single cube system, the cube’s forward locomotion is accompanied by turning (shown in Figure 1). When the input voltage is positive, the counter-clockwise-rotating vibration motor induces a leftward turning. When the input voltage is negative, the clockwise-rotating vibration motor induces a rightward turning.
2.1.2 Bi-cube and maneuver control
The bi-cube robot is comprised of two firmly bonded identical cubes by combining the cubes in the orientation in which the ERM’s rotary axes are parallel to each other. Such a bi-cube robot cannot only move forward but can also turn.
Based on the maneuver test shown in Figure 2, when the left cube (marked as #1 in Figure 2) vibration motor rotates counter-clockwise and the right cube (marked as #2) rotates clockwise, the lateral oscillation of each motor can be counterbalanced, thus the bi-cube robot performs a forward motion. When #1 cube turns off and #2 cube rotates counter-clockwise, the bi-cube robot performs a left turn. Similarly, when #2 cube turns off and #1 cube rotates clockwise, the bi-cube robot performs a right turn.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1298676/frobt-11-1298676-HTML-r1/image_m/frobt-11-1298676-g002.jpg)
Figure 2. Bi-cube steering mechanism. The bi-cube can execute three maneuvers by controlling the rotational state of vibration motors #1 and #2 inside the left and right cubes: forward motion, left turning, and right turning, from the left column to the right column, with each row illustrating the motor state from the perspective view, back view, and the corresponding maneuver.
2.1.3 Vibration frequency
To correlate input voltage with the cube’s vibration frequency, we recorded acceleration data from a single cube hanging from a thread at various voltage inputs using an IMU sensor, sampled at 200 Hz. The vibration frequency data was derived from the acceleration data through Fast Fourier Transform (FFT) analysis in Figure 3. As the voltage increases, the vibration frequency increases.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1298676/frobt-11-1298676-HTML-r1/image_m/frobt-11-1298676-g003.jpg)
Figure 3. Vibration frequency of a cube hanging in air as a function of input voltage. The frequency increases from 0 Hz to 65 Hz when the voltage increases from 2 V to 8 V.
2.2 Experiment environment
2.2.1 Air-fluidized testbed
To ensure consistent initial conditions for our experiments on granular media, we implemented a terrain creation and locomotion testing system (as shown in Figure 4), following the approach described in previous work (Qian et al., 2013). The system uses an air-fluidized bed, a container (60 cm long, 30 cm wide) filled with granular material (5 cm deep) such that the surface can be flattened at the beginning of each experiment by blowing air through a porous rigid plastic layer. A vacuum (Vacmaster) is used to force air through the bottom of the bed. The flow distribution layer, which evenly distributes the air across the bed, is approximately 0.5 cm thick and has randomly distributed pores with a diameter of 50 μm. This allows for precise control over the bed’s surface and ensures that the initial conditions of each experiment are consistent; in these experiments all states are of loosely packed granular media (Qian et al., 2013).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1298676/frobt-11-1298676-HTML-r1/image_m/frobt-11-1298676-g004.jpg)
Figure 4. The automated terrain creation and locomotion testing system. (A) Diagram of the air-fluidized bed for robot locomotion testing. (B) The process of preparing the granular media before each experiment and creating a loosely packed state (left: before, mid: blowing air to the fluidized bed, right: loosely packed state achieved after cessation of air flow).
2.2.2 Granular materials
We expect that the locomotor capabilities of our robotic cube can vary depending on the type of granular material it is tested on. To investigate this, we tested robot speed and energy use on three different types of granular material: glass particles with an average diameter of approximately 200 μm (Figure 6A), sand with particle sizes ranging from 500 to 700 μm (Figure 6B, fine sand), and 1,000–1,200 μm (Figure 6C, coarse sand). The corresponding angles of repose for these materials are 18.6°, 28.8°, and 36.7°, which indicates different values of inter-grain friction. As a comparison, we also tested cube locomotion on a hard wooden board (Figure 6D).
2.3 Experimental protocol
For each of the four terrains we considered, we conducted a series of experiments with the bi-cube robot. We varied the input voltage from 0 to 8 V to test forward speed and cost of transport (CoT). Specifically, we limited the robot’s rotational movement using two parallel walls that were 10 cm away from each other (slightly larger than the robot’s width 8 cm). For each voltage input value, we carried out three repeated experiments. In each experiment, we placed the robot at one end of the tank (
2.4 Discrete element method (DEM) simulation of a single cube
To enhance our understanding of cube movement on granular materials, we utilized Discrete Element Method (DEM) simulations, a methodology used in prior investigations into the physical interactions between robots and granular terrains Zhang et al. (2013); Pullinz et al. (2013). We used the open-source software Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) Thompson et al. (2022) to represent the motion of a vibrating single cube. To capture normal forces in granular media interactions during simulation, we used spherical particles in conjunction with the Hertzian contact model, where the normal force between two overlapping particles is proportional to the overlapping area. To model tangential forces between particles we used Mindlin No-slip solution. We incorporated rolling friction between particles, which is computed using a spring-dashpot-slider model. The parameters defining the granular particles in the simulations are detailed in Table 1, where we align the particle density and size distribution with the coarse sand particles. The simulation box of particles was created by randomly generating and pouring them into the bed. To construct the cube in simulation, we utilized 1 mm diameter spherical particles and connected them rigidly to form a hollow cube with the same dimensions as the actual cube Figure 11A). Inside the cube, a single particle is positioned at a distance
Here,
3 Result and discussion
3.1 Locomotor performance
3.1.1 Forward velocity
We measured the average speed under various applied voltages to quantify the bi-cube robot’s locomotion performance.
Under specific voltage input conditions, as illustrated in Figure 5 with the 8 V input on fine sand, the bi-cube robot may unpredictably become stuck in a collapsed pit (bubble remnants from rapid solidification), particularly during the initial phase of each test trial. Once freed from the pit, the robot exhibits steady-state motion. Utilizing tracked data, we computed the average forward speed (measured in cm/s) achieved by the bi-cube robot during each experiment’s steady-state motion. The values were averaged over three separate trials, as depicted in Figure 6, where error bars represent the standard deviation.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1298676/frobt-11-1298676-HTML-r1/image_m/frobt-11-1298676-g005.jpg)
Figure 5. Steady-state motion and a sticking case in bi-cube locomotion on find sand. In various test conditions and trials, the bi-cube robot demonstrates smooth movement, as depicted by the blue curve, illustrating a motion tracking trial under a 6 V input. In challenging scenarios, such as navigating loosely compacted terrain pits (illustrated by the red curve under an 8 V input), the robot may encounter temporary impediments before successfully escaping and achieving a stable, steady-state motion. The averaged velocity in Figure 6 is computed by the steady-state motion.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1298676/frobt-11-1298676-HTML-r1/image_m/frobt-11-1298676-g006.jpg)
Figure 6. Bi-cube speed test on (A) glass particles, (B) fine granular media (fine sand), (C) coarse granular media (coarse sand), and (D) hard wooden surface. Each figure in the left column shows the averaged velocity of three experiment trials as a function of input voltage from 0 to 8 V, with the error bar showing the standard derivation. The local maximum velocity on each terrain is marked by the blue dot. The particle pile image shows the angle of repose of glass particles (18.6°), fine sand (28.8°), and coarse sand (36.7°). The right column shows an example locomotion of the robot operated with the optimal input voltage in each terrain, in which the frames are recorded at t = 0 and t = 4s. The black scale bars indicate 5 cm. The performance is shown as videos in Supplementary Video S2.
In glass particles, the bi-cube robot cannot generate effective locomotion when the input voltage is less than 3 V. With increasing input voltage, the forward locomotion speed increases, reaching the local maximum of
Through these experimental findings, we demonstrate that the bi-cube robot is capable of achieving effective locomotion across various granular media types (at rates greater than 1 body length per second). Furthermore, we observed that the optimal operating voltage increases with the grain size and coefficient of friction. This optimal voltage may be influenced by factors such as robot dimensions or hardware limitations. These limitations encompass actuator output power and vibration intensity, as well as the unsteady desynchronization between the two vibration motors in the high-voltage input region.
3.1.2 Locomotion efficiency
In addition to evaluating the forward locomotion speed, we conduct real-time power consumption measurements for each experiment. Subsequently, we calculate the cost of transport (CoT) following
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1298676/frobt-11-1298676-HTML-r1/image_m/frobt-11-1298676-g007.jpg)
Figure 7. Bi-cube locomotion cost of transport (CoT) test on (A) glass particles, (B) fine granular media (fine sand), (C) coarse granular media (coarse sand), and (D) a hard wooden surface. Each figure shows the averaged cost of transport (CoT) of three trials as a function of input voltage with the error bars showing the standard derivation. Data points with extremely large CoT values
Experimental results reveal that the local minima of CoT vary corresponding to grain sizes and coefficient of friction. Notably, these local minimum cost of transport values remain below 10 across all granular terrains:
In comparison, we surveyed available literature values of the cost of transport (CoT, including all work input to the system) or the mechanical cost of transport (mCoT, ignoring internal wasted work) for some other living creatures and robots locomoting on dry granular terrain, such as sandfish (CoT: 100), shovel-nosed snake (CoT: 50) (Hosoi and Goldman, 2015), limbless robot (mCoT: 5-15) (Wang et al., 2023), bipedal rover (mCoT: 20-200), quadrupedal rover (mCoT: 10-100) (Bagheri et al., 2023), and Bipedal walking (mCoT: 0.67) (Hosoi and Goldman, 2015). The vibrating cube’s CoT falls within a reasonable range compared to other systems at the same scale. Furthermore, as we measure the overall cost of transport rather than the mechanical cost of transport, the vibrating cube’s mCoT value is anticipated to be smaller.
3.1.3 Voltage difference impact on forward maneuver velocity
The forward maneuver necessitates maintaining an equal input voltage magnitude for both Cube #1 and Cube #2. Any divergence in voltage between the two sides adversely affects maneuver performance and, with increasing divergence, eventually leads to turning. In Figure 8, we explore the impact of input voltage divergence on the forward motion velocity component. Specifically, we applied a fixed voltage input (6.5 V) to one side’s vibration motor while varying the voltage on the other side. As the voltage difference increases, the forward motion velocity component consistently decreases.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1298676/frobt-11-1298676-HTML-r1/image_m/frobt-11-1298676-g008.jpg)
Figure 8. The influence of a voltage difference to the two motors of the bi-cube on the forward velocity. During each test, we fixed motor 2 input at 6.5 V and varied motor 1 input from 6.5 V to 2 V. The corresponding voltage difference increases from 0 V to 4.5 V. Each velocity is tested 3 times repeatedly, with bars representing standard deviation.
3.2 Motion recording via high speed imaging
To gain insight into the mechanisms governing the bi-cube robot’s translation through vibration and subsequently develop either a kinematic or dynamic model to elucidate its motion, we utilized a high-speed 240 fps camera to capture the rigid body orientations of the robot. Subsequently, we extracted specific frames to calculate the instantaneous displacements in the forward
As discussed in Section 3.1.1, the inherent unsteady desynchronization issues identified in the bi-cube system can randomly disrupt the forward direction trajectory. Consequently, we opt to constrain the lateral motion of the single cube, focusing solely on tracking the forward direction component.
In Figure 9A, we depict a sequence of the robot’s postures during a motion period in the quasi-2D setup under a 5 V input on fine sand. Figure 9B illustrates the
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1298676/frobt-11-1298676-HTML-r1/image_m/frobt-11-1298676-g009.jpg)
Figure 9. Single cube motion tracking. (A) A sequence of the robot posture in the abstracted 2D workspace over one complete period, the red dot represents the center of geometry. The corresponding video is provided in Supplementary Video S3. (B) Tracked trajectories of
Our observations revealed that the single cube exhibited a higher velocity than the bi-cube system. Specifically, in the context of the coarse sand forward speed component, the single cube (Figure 11E) demonstrated a velocity ranging from 8–14 cm/s as the voltage increased from 4 V to 8 V (where the excitation frequency simultaneously varies from 35 to 55 Hz). In contrast, the bi-cube system (Figure 6C) exhibited a lower speed, ranging from 4–9 cm/s. This discrepancy can be attributed to the lateral forces induced by the single cube, resulting in lateral shaking and, consequently, additional forward propulsion as its back corner anchors in the sand, complementing the turning effect. Conversely, the bi-cube system counterbalances such lateral forces, emphasizing its capability for maneuvering.
Moreover, we observed enhanced stability in the single cube, particularly in the high-voltage region (6.5–8 V). This increase in stability can be attributed to the avoidance of unsteady desynchronization issues inherent in the bi-cube system, as discussed in Section 3.1.1.
3.3 Locomotion mechanism hypothesis
The
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1298676/frobt-11-1298676-HTML-r1/image_m/frobt-11-1298676-g010.jpg)
Figure 10. A simplified model for understanding single cube locomotion. (A) The 3D cube model. The vibration motor pendulum mass
In the single cube 3D model (Figure 10A), the pendulum mass
As force
During the swing-up phase in each oscillation cycle, the cube elevates its front part and undergoes a counter-clockwise rotation (as shown in Figure 10B, left half part body diagram). The cube tends to push the sand beneath it in the positive x-direction and generate the negative x-direction force
In contrast, during the touch-down phase, the front part of the cube drops down and rotates clockwise (illustrated in Figure 10B right half part body diagram). The cube tends to push the sand around its back bottom corner to the negative x-direction which leads to a positive x-direction force
During the swing-up phase, the cube is lifted by the upward direction force
3.4 Simulation analysis
We used DEM simulation to test the picture above and to determine if we could quantitatively reproduce the behavior of the single cube on coarse sand (see Figure 11A). By applying a rotating interaction force
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1298676/frobt-11-1298676-HTML-r1/image_m/frobt-11-1298676-g011.jpg)
Figure 11. Single cube system DEM simulation. (A) DEM model of a single vibrating cube. External forces
To further test our model’s capabilities we compared the forward locomotion speed of the experimental cube with that of simulation across different excitation amplitude and frequencies on coarse sand as illustrated in Figure 11E. For comparison purposes, we used excitation frequency for both simulation and experiment where the voltage applied to the experimental cube is converted to excitation frequency as explained in Section 2.1.3. The simulations and experiments exhibit monotonically increasing speeds as the excitation frequency increases, suggesting that increased energy input to the cube yields faster locomotion. This result for a single cube is in contrast with the bi-cube system speed performance, depicted in Figure 6. We note that the performance mismatch between single and bi-cube systems might be partly due to the phase difference and slight difference in mechanics between the vibration motors in the bi-cube system. We leave investigating this disagreement in a future study. Nevertheless, the simulation results not only yielded displacement dynamics similar to those of the experimental cube but also exhibited similar speeds across various excitation profiles.
3.5 Payload test
As payload capability is crucial for a robotic locomotion system, we assessed the bi-cube robot’s capacity by introducing additional weights to its body, as illustrated in Figure 12. To maintain the center of mass position, we evenly attached magnetic blocks to both sides of the robot. The original mass of the bi-cube is 146 g. As the payload increased, the forward velocity steadily decreased, eventually ceasing forward motion when the payload reached 100 g. The 146 g bi-cube demonstrated the ability to carry up to 75 g (approximately 50% of its original weight) payload while still moving, albeit at a slower speed. The payload capacity could potentially be enhanced with a more robust vibration motor.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1298676/frobt-11-1298676-HTML-r1/image_m/frobt-11-1298676-g012.jpg)
Figure 12. Payload test for the bicube during forward locomotion. The payload is set on two sides of the bi-cube evenly, ranging from 25, 50, 75, to 100 g. Compared to the 0-payload condition, normalized velocity is calculated with a standard deviation based on 3 repeated trials. The experiment is conducted on coarse sand at 6.5 V.
3.6 Surface friction
According to the locomotion hypothesis in Section 3.2, the forward propulsion of the bi-cube motion primarily arises from the interaction with granular terrain and the particle jamming at the bi-cube’s bottom rear part and the inertial swing forward; the cube’s stoppage during each cycle should be dominated by the energy dissipated during swing down (due to multiple inelastic and frictional grain collisions). We thus expected that for a range of cube belly surface friction values, the locomotor performance should be insensitive to surface-grain friction. To test this, we conducted experiments measuring forward motion velocity with various surface friction coefficients. Different sandpapers were attached to the robot’s bottom, as illustrated in Figure 13. Since it was challenging to measure the friction between the cube and the environment, we characterized the belly resistance by measuring the friction coefficients between the bi-cube’s surface and the wooden board. The results (Figure 13) indicate that for the coarse sand, varying belly surface friction does not significantly impact velocity performance. We expect that for sufficiently low friction surfaces, performance could improve although this could be masked by the large dissipation due to the collision.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1298676/frobt-11-1298676-HTML-r1/image_m/frobt-11-1298676-g013.jpg)
Figure 13. Robot-substrate contact surface friction test. The influence of cube surface friction on forward velocity is tested by attaching different sandpaper sheets to the bottom. The friction coefficients between the bi-cube surface and the wooden board are measured, which are 0.17, 0.36, and 0.58, respectively. The normalized velocity is calculated with a standard deviation based on 3 repeated trials. The experiment is conducted on coarse sand at 6.5 V.
3.7 Slope climbing
Slope climbing remains challenging for many robots designed to navigate granular environments (Marvi et al., 2014; Shrivastava et al., 2020), as steep granular slopes tend to be sensitive to stress and shear forces disturbances (Tegzes et al., 2003; Gravish and Goldman, 2014), which can cause avalanche. The intricate nature of these terrains becomes particularly evident during climbing maneuvers, as even minor disturbances introduced by robot motion can trigger avalanches and instigate the transition of the terrain from a solid to a fluid-like state, resulting in failure.
To investigate the slope-climbing capabilities of the bi-cube robot, we conducted tests on the testbed with coarse sand, with one side of the testbed elevated, thereby forming granular slopes with 4, 8, and 12 degrees of inclined angle (Figure 14B). Figure 14A shows the average speed and standard deviation of the cube’s motion, powered by 6 V input voltage, averaged over three trials. With increasing slope angles, there was a continuous decrement in the cube’s velocity. At a 12-degree slope angle, we observed that the granular slope would sometimes collapse during robot locomotion. Consequently, the robot would become stuck shortly after climbing a few body lengths.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1298676/frobt-11-1298676-HTML-r1/image_m/frobt-11-1298676-g014.jpg)
Figure 14. Bi-cube slope climbing. (A) The averaged velocity of three trials versus slope angle with standard derivation bars. The cube is tested on coarse granular media (coarse sand) at 6 V. At 12-degree slope, we recorded the average speed before the bi-cube got stuck. (B) The bi-cube slope climbing motion (4, 8, 12°), recorded at frames of
3.8 Escaping from entrapment
Given that the testbed was fluidized before each experiment, the loosely compacted granular terrain could unpredictably collapse during robot locomotion (due to the inherent bubbling in fluidized beds of this particle size Anderson et al., 1995). This led to the formation of small pits in the terrain which could cause the robot to get stuck. Thus we conducted experiments to test the bi-cube robot’s ability to escape from pits. Through experiments, we verified that with appropriate input voltages (e.g., 5 V in fine granular media), the bi-cube robot could extricate itself from a pit: The robot first engaged in a process of crawling, gradually moving the sand pile from its front to rear; and eventually, this enabled the robot to successfully escape entrapment. We provide a demonstration in Supplementary Video S6, in which the robot first became stuck in a pit in fine granular media (fine sand) when actuated by 3 V, and escaped with a 5 V input.
3.9 Maneuverability test
The bi-cube robot exhibits the capacity to execute forward, left, and right turning maneuvers, enabling effective navigation across 2D granular terrains. We illustrate this capability through a maneuver demonstration conducted on a coarse granular media surface. Notably, two cylindrical obstacles have been rigidly placed within the terrain, as depicted in Figure 15. In this demonstration, the robot’s maneuvering actions are manually switched among forward moving, left turning and right turning as shown in Figure 2. The cube follows an “
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1298676/frobt-11-1298676-HTML-r1/image_m/frobt-11-1298676-g015.jpg)
Figure 15. Bi-cube maneuver demonstration. The cube moves in an “
4 Conclusion
In this paper, we systematically tested the capability of vibratory locomotion on granular media in experiments and simulations. The vibration cube exhibits the capability to navigate across granular terrains of various particle sizes as well as solid ground. Notably, this vibratory locomotion method performs better on granular media surfaces than on hard ground, excelling in both velocity and efficiency. We posit that the flowable nature of the granular terrain plays a pivotal role in stabilizing motion and attenuating extraneous vibration energy, thereby amplifying locomotive efficacy. The inherent vibratory locomotion mechanism showcases an inherent affinity for granular terrains, suggesting a harmonious alignment between the mechanism and such environments. We employed DEM simulations to replicate the locomotion patterns of a vibrating single cube on granular particles. The simulation outcomes corroborated the experimental observations, demonstrating that the cube attains forward motion through interplay of centrifugal forces, torques and rotations generated by the vibration motor and its interaction with the granular media.
A single cube only demonstrates the capacity to execute left and right turns. However, through the fusion of two individual cubes, a bi-cube configuration is achieved, facilitating forward, left, and right turning. This amalgamation imparts a notable enhancement to maneuverability. Moreover, the inherent simplicity of the vibration cube underscores its potential to exhibit swarming capabilities on granular terrain.
Future work includes revealing the granular vibratory locomotion mechanism via both comprehensive theoretical modeling and experimental validation. Further, we intend to conduct investigations into the influence of particle size and density on locomotion efficiency based on simulation. Additionally, we will upgrade the vibration cube into a swarm robotic system, integrating self-feedback loops and inter-unit communication, for potential future application in exploration.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
BL: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Software, Validation, Visualization, Writing–original draft, Writing–review and editing. TW: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Software, Validation, Visualization, Writing–original draft, Writing–review and editing. DK: Data curation, Formal Analysis, Investigation, Methodology, Software, Validation, Writing–review and editing. VK: Data curation, Investigation, Methodology, Validation, Writing–original draft, Writing–review and editing. FH: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Validation, Writing–original draft, Writing–review and editing. DG: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Validation, Writing–original draft, Writing–review and editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. The authors are grateful for funding from Army Research Office Grant (W911NF-11-1-0514), SitS NSF-UKRI Award (1935548), and Google Cloud Research Credits program with the award code 0H64-QMAN-2BBH-0FMG. Google was not involved in the study design, collection, analysis, interpretation of data, the writing of this article, or the decision to submit it for publication.
Acknowledgments
DG would like to thank the Dunn Family Professorship.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/frobt.2024.1298676/full#supplementary-material
References
Anderson, K., Sundaresan, S., and Jackson, R. (1995). Instabilities and the formation of bubbles in fluidized beds. J. Fluid Mech. 303, 327–366. doi:10.1017/s0022112095004290
Bagheri, H., Jayanetti, V., Burch, H. R., Brenner, C. E., Bethke, B. R., and Marvi, H. (2023). Mechanics of bipedal and quadrupedal locomotion on dry and wet granular media. J. Field Robotics 40, 161–172. doi:10.1002/rob.22121
Becker, F., Minchenya, V., Zeidis, I., and Zimmermann, K. (2011). “Modeling and dynamical simulation of vibration-driven robots,” in 56th internat. Scientific colloquium innovation in mechanical engineering: shaping the future USA, (Universitätsbibliothek Ilmenau).
Calisti, M., Picardi, G., and Laschi, C. (2017). Fundamentals of soft robot locomotion. J. R. Soc. Interface 14, 20170101. doi:10.1098/rsif.2017.0101
Cocuzza, S., Doria, A., and Reis, M. (2021). Vibration-based locomotion of an amphibious robot. Appl. Sci. 11, 2212. doi:10.3390/app11052212
Golitsyna, M. (2018). Periodic regime of motion of a vibratory robot under a control constraint. Mech. Solids 53, 49–59. doi:10.3103/s002565441803007x
Gravish, N., and Goldman, D. I. (2014). Effect of volume fraction on granular avalanche dynamics. Phys. Rev. E 90, 032202. doi:10.1103/physreve.90.032202
Hosoi, A., and Goldman, D. I. (2015). Beneath our feet: strategies for locomotion in granular media. Annu. Rev. fluid Mech. 47, 431–453. doi:10.1146/annurev-fluid-010313-141324
Knuth, M. A., Johnson, J., Hopkins, M., Sullivan, R., and Moore, J. (2012). Discrete element modeling of a mars exploration rover wheel in granular material. J. Terramechanics 49, 27–36. doi:10.1016/j.jterra.2011.09.003
Lee, K., Ryu, S., Kim, C., and Seo, T. (2022). A compact and agile angled-spoke wheel-based mobile robot for uneven and granular terrains. IEEE Robotics Automation Lett. 7, 1620–1626. doi:10.1109/lra.2022.3141204
Li, S., Dutta, B., Cannon, S., Daymude, J. J., Avinery, R., Aydin, E., et al. (2021). Programming active cohesive granular matter with mechanically induced phase changes. Sci. Adv. 7, eabe8494. doi:10.1126/sciadv.abe8494
Liang, X., Xu, M., Xu, L., Liu, P., Ren, X., Kong, Z., et al. (2012). “The amphihex: a novel amphibious robot with transformable leg-flipper composite propulsion mechanism,” in 2012 IEEE/RSJ International Conference on Intelligent Robots and Systems, China, 7-12 Oct. 2012 (IEEE), 3667–3672.
Maladen, R. D., Ding, Y., Umbanhowar, P. B., and Goldman, D. I. (2011). Undulatory swimming in sand: experimental and simulation studies of a robotic sandfish. Int. J. Robotics Res. 30, 793–805. doi:10.1177/0278364911402406
Marvi, H., Gong, C., Gravish, N., Astley, H., Travers, M., Hatton, R. L., et al. (2014). Sidewinding with minimal slip: snake and robot ascent of sandy slopes. Science 346, 224–229. doi:10.1126/science.1255718
Notomista, G., Mayya, S., Mazumdar, A., Hutchinson, S., and Egerstedt, M. (2019). “A study of a class of vibration-driven robots: modeling, analysis, control and design of the brushbot,” in IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), USA, 06 December, 2022 (IEEE), 5101–5106.
Pullinz, A., Fearing, R. S., Goldman, D. I., Masarati, P., Hoover, A. M., Birkmeyer, P., et al. (2013). Walking and running on yielding and fluidizing ground. Robotics Sci. Syst. VIII, 345–352. doi:10.7551/mitpress/9816.003.0049
Qian, F., Daffon, K., Zhang, T., and Goldman, D. I. (2013). “An automated system for systematic testing of locomotion on heterogeneous granular media,” in Nature-inspired mobile robotics (China: World Scientific), 547–554.
Qian, F., Zhang, T., Korff, W., Umbanhowar, P. B., Full, R. J., and Goldman, D. I. (2015). Principles of appendage design in robots and animals determining terradynamic performance on flowable ground. Bioinspiration biomimetics 10, 056014. doi:10.1088/1748-3190/10/5/056014
Reis, M., Yu, X., Maheshwari, N., and Iida, F. (2013). Morphological computation of multi-gaited robot locomotion based on free vibration. Artif. life 19, 97–114. doi:10.1162/artl_a_00084
Rubenstein, M., Ahler, C., and Nagpal, R. (2012). “Kilobot: a low cost scalable robot system for collective behaviors,” in 2012 IEEE international conference on robotics and automation (IEEE), China, 14-18 May 2012, 3293–3298.
Shrivastava, S., Karsai, A., Aydin, Y. O., Pettinger, R., Bluethmann, W., Ambrose, R. O., et al. (2020). Material remodeling and unconventional gaits facilitate locomotion of a robophysical rover over granular terrain. Sci. Robotics 5, eaba3499. doi:10.1126/scirobotics.aba3499
Tang, C., Li, B., Fang, H., Li, Z., and Chen, H. (2018). A speedy, amphibian, robotic cube: resonance actuation by a dielectric elastomer. Sensors Actuators A Phys. 270, 1–7. doi:10.1016/j.sna.2017.12.003
Tarr, S. W., Brunner, J. S., Soto, D., and Goldman, D. I. (2024). Probing hydrodynamic fluctuation-induced forces with an oscillating robot. Phys. Rev. Lett. 132, 084001. doi:10.1103/physrevlett.132.084001
Tegzes, P., Vicsek, T., and Schiffer, P. (2003). Development of correlations in the dynamics of wet granular avalanches. Phys. Rev. E 67, 051303. doi:10.1103/physreve.67.051303
Thompson, A. P., Aktulga, H. M., Berger, R., Bolintineanu, D. S., Brown, W. M., Crozier, P. S., et al. (2022). LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm. 271, 108171. doi:10.1016/j.cpc.2021.108171
Wang, D., Liu, Y., Deng, J., Zhang, S., Li, J., Wang, W., et al. (2022). Miniature amphibious robot actuated by rigid-flexible hybrid vibration modules. Adv. Sci. 9, 2203054. doi:10.1002/advs.202203054
Wang, T., Pierce, C., Kojouharov, V., Chong, B., Diaz, K., Lu, H., et al. (2023). Mechanical intelligence simplifies control in terrestrial limbless locomotion. Sci. Robotics 8, eadi2243. doi:10.1126/scirobotics.adi2243
Zhang, T., Qian, F., Li, C., Masarati, P., Hoover, A. M., Birkmeyer, P., et al. (2013). Ground fluidization promotes rapid running of a lightweight robot. Int. J. Robotics Res. 32, 859–869. doi:10.1177/0278364913481690
Keywords: vibration, granular media, robot, locomotion, modular robot, robophysics, DEM simulation
Citation: Liu B, Wang T, Kerimoglu D, Kojouharov V, Hammond FL III and Goldman DI (2024) Robust self-propulsion in sand using simply controlled vibrating cubes. Front. Robot. AI 11:1298676. doi: 10.3389/frobt.2024.1298676
Received: 22 September 2023; Accepted: 11 July 2024;
Published: 30 August 2024.
Edited by:
Murat Reis, Bursa Uludağ University, TürkiyeReviewed by:
Vincent Lebastard, IMT Atlantique Bretagne-Pays de la Loire, FranceGuillermo Javier Amador, Wageningen University and Research, Netherlands
Copyright © 2024 Liu, Wang, Kerimoglu, Kojouharov, Hammond and Goldman. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Frank L. Hammond III, frank.hammond@me.gatech.edu; Daniel I. Goldman, daniel.goldman@physics.gatech.edu
†These authors have contributed equally to this work