Humanoid Robot: Balance, Joint Torque and Runtime

Estimate humanoid robot battery life, ZMP balance margins and hip/knee stance torque. Animate a walking step and compare battery, gait and duty-cycle assumptions.

#battery runtime #cost of transport #humanoid robot #joint torque #LIPM #ZMP

Estimate humanoid robot battery life, inspect a ZMP walking step and compare simplified hip and knee loads in one calculator. It is intended for engineers choosing early design targets, students learning walking dynamics and anyone checking whether a proposed battery and duty cycle are plausible. All robot specifications are unofficial user inputs; the defaults describe a generic scenario.

Start with robot mass including its battery, then add payload and enter the combined center-of-mass height. Step length is the distance between consecutive opposite-foot contacts. Nominal speed sets the step duration. Enter the dimensions of the planted sole, the hip height and the two leg lengths. Unreachable poses and swing-knee floor penetration are flagged across the sampled step. Press ▶ beside step progress to watch the body move over the stance foot while the other leg swings forward. The footprint panel distinguishes the required ZMP from the moving CoM projection.

The balance model follows the constant-height linear inverted pendulum equations described by Kajita and colleagues. A symmetric nominal step is combined with optional initial acceleration biases. These biases illustrate how an open-loop motion can drift even while its required ZMP stays fixed. The signed sole margin is the distance to the nearest rectangular contact edge: a negative value means the requested contact is infeasible. A separate friction check compares horizontal reaction force with available sole friction. Neither check establishes that a complete gait is stable or recoverable.

Use the joint torque calculator results as preliminary stance-load estimates. Ground reaction moments are evaluated at the hip and knee of a massless two-link leg. The cards show sampled peak magnitudes multiplied by your sizing allowance. Lowering the hip generally increases the knee lever arm. These values exclude swing inertia, link weight, impacts, gearbox friction, motor speed limits and continuous thermal duty; they cannot select a complete actuator by themselves.

For runtime, enter pack energy, usable fraction, incremental electrical cost of transport and always-on power. Cost of transport converts weight and distance into battery energy; it already includes the motion drive losses. Always-on power covers computing, sensing, cooling and stationary holding. If a published CoT already includes those loads, set always-on power to zero. The walking fraction weights motion power across the mission. Manipulation and charging time are excluded. The speed sweep holds CoT fixed, so it is a sensitivity study rather than a measured efficiency map. Changing leg posture does not automatically recalibrate CoT.

TrendForce reported typical two-to-four-hour operation and packs generally below 2 kWh on 28 January 2026. The 1X NEO page, checked on 27 October 2026, lists 842 Wh and four hours without defining an equivalent continuous walking test. These figures provide context, not calibration. Replace assumptions with measured power and contact data before using this tool for a particular robot.

How it is calculated

The model combines one single-support balance step, sagittal stance load moments and an independent electrical energy budget.

Kajita et al. (2003), attached source s1, page 2, equations (6)–(10), gives px=x−hx¨/gp_x=x-h\ddot{x}/g and py=y−hy¨/gp_y=y-h\ddot{y}/g for constant CoM height and negligible centroidal angular-momentum rate. Here ω=g/h\omega=\sqrt{g/h} and T=S/vT=S/v. The symmetric nominal trajectory is

xn(t)=S2sinh⁡[ω(t−T/2)]sinh⁡(ωT/2).x_n(t)=\frac{S}{2}\frac{\sinh[\omega(t-T/2)]}{\sinh(\omega T/2)}.

With initial position offsets and acceleration biases bx,byb_x,b_y,

x=x0+xn+bxω2(cosh⁡ωt−1),y=y0+byω2(cosh⁡ωt−1).x=x_0+x_n+\frac{b_x}{\omega^2}(\cosh\omega t-1),\qquad y=y_0+\frac{b_y}{\omega^2}(\cosh\omega t-1).

Thus the required ZMP is constant: px=x0−bx/ω2p_x=x_0-b_x/\omega^2, py=y0−by/ω2p_y=y_0-b_y/\omega^2. The nominal initial forward velocity is Sω/[2tanh⁡(ωT/2)]S\omega/[2\tanh(\omega T/2)]; lateral initial velocity is zero. Only the unbiased step has mean forward speed equal to the nominal mission speed.

Accelerations are evaluated as analytic second derivatives: x¨=x¨n+bxcosh⁡(ωt)\ddot{x}=\ddot{x}_n+b_x\cosh(\omega t) and y¨=bycosh⁡(ωt)\ddot{y}=b_y\cosh(\omega t). An independent four-state ode45 solution integrates x˙=vx\dot{x}=v_x, v˙x=g(x−px)/h\dot{v}_x=g(x-p_x)/h and the corresponding lateral equations from the prescribed initial state. Property checks compare the full trajectories, accelerations, and initial/final velocities, sweeping nonzero forward and lateral biases and nominal speed. The agreement tolerances are 2 μm for position, 5×10⁻⁵ m/s² for acceleration and 2 μm/s for endpoint velocity.

The sole margin is min⁡(L/2−∣px∣,W/2−∣py∣)\min(L/2-|p_x|,W/2-|p_y|). The model requires Fz=mgF_z=mg and Fx=mx¨F_x=m\ddot{x}. Available friction is compared with max⁡x¨2+y¨2/g\max\sqrt{\ddot{x}^2+\ddot{y}^2}/g.

A forward-bending knee is obtained by two-circle intersection using the thigh and shank lengths. At each joint jj, the sagittal ground-reaction load moment is

Mj=(px−xj)Fz+zjFx.M_j=(p_x-x_j)F_z+z_j F_x.

The balancing actuator torque has the opposite sign. Peak cards use magnitudes across 201 evenly spaced stance samples, multiplied by the user sizing allowance. Swing-knee floor clearance is checked over all 201 step samples and at the current animated pose; a negative clearance blocks results. No full multibody inverse dynamics is claimed.

Collins et al. (2005), s2 page 2, defines dimensionless electrical CoT as energy divided by weight and distance, and distinguishes it from positive mechanical work. Here em=Cwmge_m=C_w mg, Pm=emvP_m=e_m v, Pwalk=Pm+PauxP_{walk}=P_m+P_{aux}, and Pmission=dPm+PauxP_{mission}=dP_m+P_{aux}. Available energy is Eu=EpackfuE_u=E_{pack}f_u; runtime is Eu/PE_u/P and mission walking distance is vdtvdt. Use explicit kWh/Wh and hour/second conversions. The CoT input is incremental: a total published electrical CoT requires Paux=0P_{aux}=0.

Assumptions

  • Generic, unofficial robot inputs; no manufacturer presets or inferred specifications.
  • Flat, rigid rectangular stance sole, ankle centered in it, no slip and no vertical CoM acceleration. Single support carries the full combined weight.
  • Robot plus payload is represented by one mass at constant height; the hip is directly below it. Both leg links are massless. Lateral moments and link inertias are omitted.
  • The illustrated swing ankle follows a smooth clearance path proportional to leg length. Its inverse-kinematic pose is schematic and contributes no torque or energy. Footprint lateral spacing is schematic.
  • Replaying the animation resets a single step. Impact, double support, push-off, support exchange and whole-gait feedback are not simulated. The lateral initial condition is not a periodic gait controller.
  • Electrical CoT stays constant for the mission and the speed sweep. Joint geometry and electrical CoT are separate estimates; posture changes require a measured or otherwise justified CoT update.
  • Always-on power persists during stationary time and walking. Reserve, pack aging and other energy derating are combined in the usable fraction. No extra efficiency is applied to electrical CoT.

Where it applies

Use for preliminary comparisons and teaching, on level ground with moderate lateral motion. ZMP inside the sole is a necessary contact-feasibility condition under these assumptions; it is not proof of stability, capturability or a feasible whole-body gait. Tedrake's Underactuated Robotics, chapter 5, explains these limits.

Correct all reach and knee-clearance errors before interpreting joint loads. Swing-knee floor clearance is checked at all 201 step samples and at the current animated pose. A longer support interval amplifies bias exponentially; the tool warns when ωT>3\omega T>3. The sagittal load estimate is flagged for lateral excursions above 20% of CoM height. Friction is a separate necessary contact check.

The sampled stance moments exclude link gravity/inertia, swing loads, impact peaks, compliance, ankle sizing, torso angular momentum, joint speed and motor/gearbox limits. A sizing multiplier cannot replace full inverse dynamics, structural analysis or actuator thermal evaluation. The results do not establish a safe robot design.

The energy budget excludes manipulation, slopes, starts/stops and charging. Battery sag, discharge-rate limits, cold-temperature performance and changing compute loads require measured corrections. At zero walking fraction and zero always-on power, a finite runtime is undefined and the tool asks for a positive load.

The Cornell and ASIMO checks compare published electrical estimates only; their generic leg geometry is not a reconstruction of those robots. The Cornell paper's rounded CoT predicts 10.20 W versus its measured 11 W. For ASIMO, its approximate CoT 3.2 gives 0.529 h versus about 0.5 h; the reported raw battery/speed numbers imply CoT about 3.39. These differences are documented in the reference cases rather than removed by changing physics.

Industry context is dated: TrendForce, 28 January 2026; 1X NEO manufacturer page checked 27 October 2026. A vendor's quoted runtime lacks the duty-cycle definition needed to validate a continuous walking prediction.

References

  • Kajita et al. (2003). Biped Walking Pattern Generation by using Preview Control of Zero-Moment Point. ICRA, pp. 1620–1626. DOI: https://doi.org/10.1109/ROBOT.2003.1241826 . Open-access author paper: https://people.csail.mit.edu/katiebyl/kb/DW2008/papers_of_tangential_interest/kajita03.pdf
  • Collins, Ruina, Tedrake & Wisse (2005). Efficient bipedal robots based on passive-dynamic walkers. Science 307, 1082–1085. DOI: https://doi.org/10.1126/science.1107799 . Author-hosted paper: https://groups.csail.mit.edu/robotics-center/public_papers/Collins05.pdf
  • Tedrake. Underactuated Robotics, chapter 5: Highly-articulated Legged Robots; especially ZMP Stability Metric. Free course text: https://underactuated.mit.edu/humanoids.html
  • TrendForce (28 January 2026). Humanoid Robots Move Toward Commercialization; Solid-State Batteries Set to Break Through Power Bottlenecks. https://www.trendforce.com/presscenter/news/20260128-12902.html
  • 1X. NEO Home Robot — battery and hardware specifications. Manufacturer claims checked 27 October 2026; runtime duty unspecified. https://www.1x.tech/neo

Inputs

InputSymbolUnitDefault
Robot mass, including battery
Unofficial user specification. Include the battery, exclude the added payload.
m_robotkg60
Added payload
The combined center of mass must be updated separately when payload position changes.
m_loadkg5
Combined center-of-mass height
Constant height above the floor for robot plus payload.
h_comm0.85
Step length
Distance between consecutive opposite-foot contacts; a stride is two steps.
S_stepm0.3
Nominal walking speed
Sets step time as length / speed. It is the actual mean CoM speed only when the acceleration biases are zero.
v_walkm/s0.6
Stance sole length
Flat rectangular contact; ankle is assumed at its center.
L_footm0.26
Stance sole width
Only the planted foot contributes to the support polygon.
W_footm0.14
Fore-aft CoM offset
Shift of the nominal mid-stance CoM relative to the stance ankle.
x_midm0
Lateral CoM offset
CoM lateral position at the beginning of this single-support step, relative to the stance ankle.
y_midm0.01
Fore-aft acceleration bias at step start
Added to the nominal LIPM acceleration at step start. Held-ZMP dynamics amplify this bias during the step.
a_bias_xm/s^20
Lateral acceleration at step start
Starts with zero lateral velocity. This single-step model has no lateral foot-placement controller.
a_bias_ym/s^20
Available sole/floor friction coefficient
Assumed dry Coulomb friction; compare with the peak horizontal/vertical ground-force ratio.
mu_floor0.6
Hip height
Hip stays directly below the CoM in the sagittal view. Use an appropriate crouched height.
h_hipm0.74
Thigh length
Hip joint center to knee joint center.
L_thighm0.44
Shank length
Knee joint center to ankle joint center. Ankle height is idealized as zero.
L_shankm0.44
Torque sizing multiplier
Applied to the sampled stance peak. It is a user allowance, not a substitute for inverse dynamics or thermal rating.
K_size1.5
Nominal battery energy
Generic scenario, not a specification of a named robot.
E_packkWh1
Usable battery fraction
Fraction remaining after reserve, state of health and operating derating. Enter a fraction, e.g. 0.85.
f_use0.85
Incremental electrical cost of transport
Battery-side motion energy / (combined weight × distance), excluding the always-on load below. If using a published total electrical CoT, set always-on power to zero.
C_walk0.7
Always-on battery power
Compute, sensing, cooling and stationary holding, averaged over the mission. Include any idle joint losses here; avoid counting them again in CoT.
P_auxW100
Fraction of mission spent walking
Walking occurs at the specified speed. Remaining time uses only the always-on power; manipulation energy is excluded.
f_walk0.7
Step progress · ▶ play
One step from contact to contact; replay resets the initial state. Swing-leg motion is a schematic clearance path.
phase0.65

Results

With the default inputs above, the calc gives:

ResultValue
Mixed-mission runtime2.96 h
Walking distance per mission4.47 km
Signed ZMP sole margin60 mm
Sole-contact checkZMP inside sole
Hip sizing torque · stance only18.6 N·m
Knee sizing torque · stance only230.4 N·m
Average mission battery power287.4 W
Continuous walking runtime2.31 h
Walking energy per distance170.2 Wh/km
Required peak friction coefficient0.176

Checked against sources

Each case below was calculated by this calc and agrees with the value from its source (within the tolerance where one is shown).

Cornell walker: published electrical power

Source: Collins et al., Science 307 (2005), attached page 2 (printed p. 1083), DOI 10.1126/science.1107799. The paper reports about 0.2 electrical CoT, 13 kg, 0.4 m/s and measured 11 W. Rounded CoT gives 10.20 W; the 1 W absolute tolerance admits that published rounding. P_aux=0 because this published CoT includes electronics. Checks only energy; generic leg geometry is not a reconstruction of the Cornell robot.

Inputs: Robot mass, including battery = 13 kg; Added payload = 0 kg; Nominal walking speed = 0.4 m/s; Incremental electrical cost of transport = 0.2; Always-on battery power = 0 W; Fraction of mission spent walking = 1; the others as above.

ResultSourceThis calc
Total power while walking11 ± 1 W10.199 W

ASIMO: published battery/runtime energy estimate

Source: Collins et al. (2005), attached page 4 (printed p. 1085), note 15, DOI 10.1126/science.1107799: 510 N weight, 1.6 km/h, 38.4 V × 10 Ah, about 30 min, estimated electrical CoT ≈3.2. Those battery/speed numbers imply CoT ≈3.39, not exactly 3.2. Keeping the physical equations and reported CoT gives 0.5294 h; 0.04 h tolerance reflects approximate, non-simultaneous published figures. P_aux=0 avoids double counting. This historical case validates the energy budget only, not ASIMO kinematics.

Inputs: Robot mass, including battery = 52.006 kg; Added payload = 0 kg; Nominal walking speed = 0.4444 m/s; Incremental electrical cost of transport = 3.2; Always-on battery power = 0 W; Nominal battery energy = 0.384 kWh; Usable battery fraction = 1; Fraction of mission spent walking = 1; the others as above.

ResultSourceThis calc
Vertical single-support ground force510 N510 N
Continuous walking runtime0.5 ± 0.04 h0.5294 h
Available battery energy384 Wh384 Wh

Related Resources

Responses

Your email address will not be published. Required fields are marked *

L ading...