Breaking speed scaling in quadrupedal robots via Huygens' coupled-pendulum dynamics
Authors: Yucheng Tao, Yongbin Jin, Shaowen Cheng, Xianwei Liu, Yanyan Yuan, Yanhong Liang, Chengkai Su, Chaojie Fu, +3 more
Organizations: Center for X-Mechanics, Zhejiang University, Hangzhou, China. · ZJU-Hangzhou Global Scientific and Technological Innovation Center, Hangzhou, China. · Mirrorme Technology Co., Ltd., Shanghai, China.
Achieving biological-level running speeds has largely been pursued through advances in control algorithms, which improve the utilization of existing hardware. However, the ultimate speed limits remain governed by the underlying force and torque requirements of rapid locomotion, which are typically addressed through increased actuator capacity. Inspired by Huygens' coupled pendulums, we demonstrate that superior locomotion can emerge from principled exploitation of intrinsic dynamics rather than brute-force hardware scaling. Inter-limb inertial coupling redistributes energy across the gait cycle and reduces peak joint torque required for rapid periodic motion, thereby expanding the achievable speed without proportional increases in actuator capability. Incorporating hardware parameters as additional design variables further extends this analysis into a co-optimization framework, enabling the systematic utilization of inertial coupling in robot design. Guided by this framework, a quadruped robot achieves a running speed of 10.74 m/s (Froude number 21.4) and completes a 100-meter sprint in 12.2 seconds, representing the first legged robot to surpass 10 m/s. These results establish inertial coupling as an underlying mechanism governing high-speed legged locomotion and highlight its role in reducing force requirements, offering new insights into the design of agile robotic systems.
Figures & tables
Figure 1: Statistics of the maximum speed and body mass of mammals and quadrupedal robots in logarithmic scales. There are total 171 mammalian species of mammals with size ranging from an adult mouse to an elephant. Data for mammalian species are obtained from Hirt et al. [ 17 ] , in which body mass and running speed adhere to a general scaling law: v=25.5m0.26(1−e−22m−0.6) . The lower-left inset summarizes peak locomotion performance and actuator capability for representative quadrupedal robots. Detailed data are provided in Table S1 .
Figure 2: Mechanical structure diagram of BP2. BP2 has a standing height of 55 cm, a body length of 80 cm, a width of 44 cm, and a total mass of 36 kg with battery. The shank structure is enlarged to clarify the assembly of the spring components and foot.
Figure 3: Analysis of the coupled dynamic model. ( A ) Experimental setup for Huygens’ coupled pendulums. ( B ) The simplified coupled dynamic model. ( C ) Dimensionless per-pendulum energies E~i=Ei/(mgl) vs. normalized time t/T0 (with T0=2πl/g ). Light traces are instantaneous, bold traces are trailing time-averages ⟨E~1,2⟩T . ( D ) Poincaré section analysis of the coupled dynamic system. ( E ) Active torque requirements under specific kinematic conditions ( R=1,A=1rad,ω=5Hz,ϕ=41π ).
Figure 4: Force decomposition and contribution analysis of the coupled dynamic model. ( A-C ) Stacked force profiles at 10 m/s for key degrees of freedom. ( D ) Optimized robot trajectory at 10 m/s, with the red arrow indicating the direction of forward motion. ( E-H ) Mean magnitude of each force/torque component in different running speeds. ( I-L ) Relative contribution of each dynamic term in different running speeds.
Figure 5: Experimental results of high-speed locomotion. ( A ) Representative keyframes of three consecutive gait cycles. Intermediate frames are faded and spatially shifted to illustrate the temporal evolution of leg configurations. ( B ) Bird’s-eye view of the experimental environment. ( C ) GPS-recorded global velocity profile. ( D ) Poincaré section of body angular-velocity dynamics during steady high-speed running. ( E ) Motor torque-speed operating points over three cycles during steady high-speed running.
Figure 6: Hundred-Meter Sprint Benchmark. ( A ) Field snapshot of the race between BP2 and Lyles [ 39 ] . ( B ) Speed–distance profiles for BP2 and world-class athletes. The kinematic data for Noah Lyles and Usain Bolt are derived from the official 10 m split times reported for the 2024 Paris Olympics [ 42 ] and the biomechanical report of the 2009 Berlin World Championships [ 14 ] , respectively, and the corresponding continuous curves are obtained through curve fitting. ( C ) Speed–distance profiles for BP2 and trained athletes. ( D ) Speed–distance profiles for BP2 and untrained adult participants.
Figure 7: Peak Speed Endurance Test. ( A ) Experimental setup, snapshot of the robot at peak velocity. ( B ) Thermal images captured at the initial pose and during peak speed operation. ( C ) Torque-speed profile of the key joint motor during the endurance test. ( D ) Temperature curve of the key actuator during the endurance test.
Figure S1 : Hybrid dynamic model of the spring-legged robot with a complete motor inertia model. Schematic of the sagittal-plane spring-legged model used for deriving the hybrid dynamics. The model is used for the forelimbs and adapted for the hindlimbs by setting L0 and l0 to negative values. In the derived equations of motion, the motor- and transmission-related parameters enter only through the rotor inertia Im and the gear ratio r .
Figure S2 : Co-optimization results for the sagittal-plane model at 10 m/s. Throughout, blue and red curves represent the left and right legs, respectively. ( A ) Robot trajectory, with the red arrow indicating forward motion. ( B-E ) Velocity profiles for horizontal center-of-mass (CoM) motion (B), vertical CoM motion (C), hip rotation (D), and knee rotation (E). ( F, G ) Ground reaction forces. ( H, I ) Hip (H) and knee (I) joint torques. Gray dashed lines indicate the upper and lower torque limits of the hind leg motor at the corresponding joint velocities.
Robot Model
Body Mass
Motor Peak Torque
Motor Max Speed
Peak Speed
m (kg)
τmax (Nm)
ωmax (rad/s)
v (m/s)
MIT Mini Cheetah [ 26 ]
9
18
40
3.9
BlackPanther [ 25 ]
10
18
42
5
BlackPanther2
36
120
21.4
10.74
StarlETH [ 20 ]
23
24
10.5
0.75
ANYmal [ 22 ]
30
40
12
1.6
Table S1 : Performance specifications of representative legged robots. The table summarizes the body mass, motor peak torque, maximum motor speed, and achievable peak locomotion speed of different robotic platforms.
Category
Level
Event
Qualifying Time
Units
Male
Level-1
100 m
10.93
s
Male
Level-2
100 m
11.74
s
Male
Level-3
100 m
12.64
s
Female
Level-1
100 m
12.33
s
Female
Level-2
100 m
13.04
s
Female
Level-3
100 m
14.04
s
Table S2 : Performance standards for China National Level athletes in the 100 m sprint. The table summarizes the qualifying times for male and female athletes under the official athletics grading system.
Symbols
Physical Meaning
Typical Values
Units
m0
Mass of torso
6.15
kg
m1
Mass of thigh
2.2
kg
m2
Mass of shank
0.1
kg
m3
Mass of toe
0.5
kg
I0
Inertia of torso
0.10
kg ⋅ m 2
I1
Inertia of thigh
0.014
kg ⋅ m 2
Table S3 : Typical physical parameters of the sagittal-plane model. Values correspond to the nominal configuration.
Table S4 : Reward terms used for reinforcement learning controller training. The imitation stage emphasizes reference-motion guidance, whereas the relaxed stage continues training from the imitation policy with reweighted objectives for high-speed locomotion. Base coefficients are given as imitation/relaxed. Terms below the separator are added for the higher-speed BP2 training.
Symbol
Description
qref , q˙ref
Reference joint position and velocity generated by the kinematic gait generator.
q , q˙
Measured joint position and velocity of the robot.
vcmd , ωcmd
Commanded linear and angular velocities.
v , ω
Measured body linear and angular velocities.
τn
Normalized joint torque.
τn,last
Normalized joint torque at the previous control step.
Achieving high-speed locomotion in quadrupedal robots remains highly challenging, as actuators operate near their physical limits and exhibit pronounced nonlinearities. However, many existing methods neglect actuator nonlinearities and physical constraints during training, leading to a significant sim-to-real gap under highly dynamic motions and limiting achievable performance. To address this issue, we propose a high-speed locomotion framework that reduces sim-to-real discrepancies and stabilizes learning over a wide command distribution. A refined actuator model explicitly captures high-speed voltage coupling and magnetic saturation, enabling a more accurate representation of the torque-speed envelope. In addition, a reinforcement learning framework incorporating a two-stage curriculum and adaptive command scheduling (ACS) ensures stable training. Experiments on the 36.5 kg quadruped BlackPanther2 (BP2) demonstrate speeds of up to 13.2 m/s on a treadmill and 11.65 m/s outdoors, establishing a new state-of-the-art and, to the best of our knowledge, a world record for quadrupedal robot locomotion. The results further highlight the importance of accurate actuator modeling in preventing non-physical policy exploitation, and show that ACS improves robustness without sacrificing performance.
Yucheng Tao, Shaowen Cheng, Guorong Lan +3
Center for X-Mechanics and Institute of Applied Mechanics, Zhejiang University, Hangzhou 310012, China · ZJU-Hangzhou Global Scientific and Technological Innovation Center, Hangzhou 311200, China · MirrorMe Robotics Co., Ltd., Hangzhou 311215, China
Legged robots operate across a wide range of physical scales, but how their designs should be adapted as size changes remains unclear. Here, we tackle this question in two ways. First, we survey existing legged robots to provide a broad context for the key scaling variables, robot mass m and leg length L. We find the surprising result that bipedal robot mass generally scales with the length squared, L^2, rather than the isometric prediction L^3. Then, to reduce the variance in design choices, we focus on a pair of previously developed bipeds that share the same quasi-passive morphology but differ by a factor of six in leg length, use different feet and controllers, and achieve different relative speeds. We reconstruct both robots in a common 3-D simulation environment and scale each design over leg lengths from 0.02 to 1.2 meters under both mass models (mass is proportional to L^2 and is proportional to L^). The controlled comparison shows that velocity follows dynamic similarity, velocity is proportional to L^{1/2}, across designs and mass models, while the torque needed to sustain walking follows that tau is proportional to mL. Consequently, torque scales approximately with L^3 when m is proportional to L^2 and L^4 when mass is proportional to L^3. A 3-D foot-shape sweep further shows that foot dimensions that permit walking scale approximately linearly with leg length, but the speed-maximizing shape and the mechanism by which each robot moves do not transfer by scaling alone. Overall, the results provide practical insights for rescaling legged systems that leverage natural body dynamics.
Naomi Oke, Aja M. Carter, Ben Gu +4
Department of Mechanical Engineering, Carnegie Mellon University, Pittsburgh, PA 15213, USA · Robotics Institute, Carnegie Mellon University, Pittsburgh PA 15213, USA
The spine plays a crucial role in the dynamic locomotion of quadrupedal animals, improving the stability, speed, and efficiency of their gait, especially for fast-paced and highly agile movements. Therefore, the spine is also a promising and natural way to extend the capabilities of quadruped robots. This paper empirically investigates the benefits of an actuated spine for learning agile quadruped locomotion. We evaluate whether the use of the spine brings benefits in terms of high-speed running, climbing stairs, climbing high-angle slopes, hurdling, and crawling scenarios. We conducted an empirical study in MuJoCo simulation using the Silver Badger robot from MAB Robotics with an actuated 1-DOF spine in the sagittal plane. The obtained results show that the use of the spine provides the robot with increased agility and allows it to overcome higher stairs, steeper slopes, higher obstacles, and smaller passages.
Nico Bohlinger, Piotr Kicki, Davide Tateo +2
Department of Computer Science, Technical University of Darmstadt, 64289 Darmstadt, Germany · IDEAS Research Institute, 00-060 Warsaw, Poland · Institute of Robotics and Machine Intelligence, Poznan University of Technology, 60-965 Poznan, Poland