---
title: Angular Momentum Linear Inverted Pendulum Model
url: https://www.emergentmind.com/topics/angular-momentum-linear-inverted-pendulum-model-amlip
type: topic
---

# Angular Momentum Linear Inverted Pendulum Model

Angular Momentum Linear Inverted Pendulum Model (AMLIP) is a reduced-order model for bipedal locomotion in which the horizontal center-of-mass (CoM) displacement and the angular momentum about the contact point are taken as the principal dynamic states. Under constant CoM height, pinned or point-foot contact, and negligible ankle-torque or centroidal-angular-momentum terms, AMLIP yields a linear unstable pendulum model whose continuous dynamics are formally identical to those of the classical Linear Inverted Pendulum (LIP), but with contact-point angular momentum replacing linear velocity as the state used for prediction and feedback. In the recent locomotion literature, closely related formulations also appear under the abbreviation ALIP, especially when the same state substitution is embedded in hybrid walking, constrained footholds, moving rigid surfaces, stair climbing, and full-body momentum regulation [2008.10763, 2105.08170].

## 1. Historical emergence and terminology

The formulation entered the recent legged-locomotion literature through the argument that angular momentum about the contact point is a better indicator of robot state than linear velocity. “Angular Momentum about the Contact Point for Control of Bipedal Locomotion: Validation in a LIP-based Controller” reformulated a standard LIP controller in terms of angular momentum and validated the resulting feedback law on Cassie Blue [2008.10763]. “Zero Dynamics, Pendulum Models, and Angular Momentum in Feedback Control of Bipedal Locomotion” then placed the same idea in the broader context of pendulum models and virtual-constraint-induced zero dynamics, arguing that an approximate model parameterized by angular momentum provides better predictions for foot placement on a physical robot than a related approximate model parameterized in terms of linear velocity [2105.08170].

The cited literature uses both “AMLIP” and “ALIP.” In both cases, the defining feature is that angular momentum about the stance or contact point is elevated to a first-class reduced-order state. In the planar sagittal setting, the state is commonly written as
\[
\xi = \begin{bmatrix} x_n \\ L \end{bmatrix},
\]
where \(x_n\) is the horizontal CoM displacement measured from the stance foot and \(L\) is the sagittal-plane angular momentum about the stance foot [2105.08170]. In 3D ALIP formulations, the same idea is applied separately in sagittal and frontal planes, with states \(\mathbf{x}=[p_x,\;L_{c,y}]^\top\) and \(\mathbf{y}=[p_y,\;L_{c,x}]^\top\) [2408.05308].

The original reduced-order assumptions are narrow and explicit: the CoM is constrained to move at a constant height \(H\) or \(h\), the stance contact is pinned or point-foot and non-slipping, ankle torques are neglected or set to \(u_a \simeq 0\), and other degrees of freedom are assumed to be regulated by high-bandwidth virtual-constraint or whole-body controllers so that their effect enters only through small neglected terms [2105.08170, 2104.10367]. Much of the subsequent literature can be read as a sequence of controlled relaxations of these assumptions.

## 2. Canonical continuous-time dynamics

The starting point is the exact contact-point angular momentum relation
\[
L = L_c + m(x_n \dot z_n - z_n \dot x_n),
\]
where \(L_c\) is the angular momentum about the CoM, \(m\) is total mass, and \((x_n,z_n)\) are CoM coordinates in the stance-foot frame [2105.08170]. Under constant CoM height \(z_n=H\) and \(\dot z_n=0\), this reduces to
\[
L = L_c + mH v_n.
\]
Neglecting the small CoM-about-CoM angular momentum term \(L_c\) gives the AMLIP kinematic equation
\[
\dot x_n = \frac{L}{mH}.
\]
Taking moments about the contact point and setting the ankle torque to zero yields
\[
\dot L = mg x_n.
\]
The resulting linear time-invariant reduced model is therefore
\[
\dot x_n = \frac{1}{mH}L,
\qquad
\dot L = mg x_n,
\]
or, in second-order form,
\[
\ddot x_n = \frac{g}{H}x_n = \omega^2 x_n,
\qquad \omega^2 = \frac{g}{H}.
\]
This is the canonical AMLIP system [2105.08170].

The formal relation to the traditional LIP is immediate. In the standard LIP, the state is \([x_n;\,v_n]\) and the dynamics are \(\ddot x_n=(g/H)x_n\). In AMLIP, the second state variable is replaced by \(L/(mH)\), so the continuous dynamics are formally identical after the substitution \(v_n \mapsto L/(mH)\) [2105.08170, 2008.10763]. In planar underactuated formulations, a mass-normalized momentum coordinate is also used:
\[
L_y = h\,\dot x,
\qquad
\begin{pmatrix}\dot x\\\dot L_y\end{pmatrix}
=
\begin{pmatrix}
0 & \tfrac{1}{h}\\
g & 0
\end{pmatrix}
\begin{pmatrix}x\\L_y\end{pmatrix},
\]
which is algebraically equivalent after scaling by mass [2104.10367].

For 3D ALIP, the same constant-height approximation leads to decoupled planar systems,
\[
\dot{\mathbf{x}} = A\mathbf{x},
\qquad
\dot{\mathbf{y}} = -A\mathbf{y},
\qquad
A=
\begin{pmatrix}
0 & \tfrac{1}{mH}\\
mg & 0
\end{pmatrix},
\]
with closed-form step solutions expressed through hyperbolic functions [2408.05308]. This decoupled structure underlies much of the subsequent foot-placement and template-tracking design.

## 3. Hybrid step-to-step regulation and stability

AMLIP becomes a walking controller when the unstable continuous-time template is closed through discrete foot placement. In the fixed-step-time formulation with \(T\) the step duration and \(\ell=\sqrt{g/H}\), the one-step-ahead predictor for contact-point angular momentum is
\[
\hat L(T_k^-,t)
=
mH\ell \sinh\!\bigl(\ell(T_k^- - t)\bigr)x_n(t)
+
\cosh\!\bigl(\ell(T_k^- - t)\bigr)L(t).
\]
Using impact continuity of \(L\) on flat ground and the reset of the horizontal coordinate by the landing location, the desired swing-foot landing position that achieves a target one-step-ahead momentum is
\[
p^{\,\mathrm{des}}
=
\frac{
L^{\mathrm{des}}
-
\cosh(\ell T)\,\hat L(T_k^-)
}{
mH\ell \sinh(\ell T)
}.
\]
This is the core AMLIP foot-placement law [2105.08170, 2008.10763].

A more general asymptotic law introduces a gain \(\alpha\in[0,1)\), producing a step-to-step error recursion
\[
e_{k+1}=\alpha e_k.
\]
The corresponding discrete Poincaré map has eigenvalues \(\{\alpha,0\}\), and the fixed point
\[
x_n^* = \frac{(1-\cosh(\ell T))L^{\mathrm{des}}}{mH\ell \sinh(\ell T)},
\qquad
L^* = L^{\mathrm{des}}
\]
is exponentially stable for any \(\alpha\in[0,1)\) [2105.08170]. In simulation on the planar Rabbit robot, the dominant eigenvalue of the full-order Poincaré map closely matches \(\alpha^2\) predicted by the ALIP model across \(\alpha\in[0,0.9]\), supporting the interpretation of AMLIP as an approximation of the slow zero dynamics rather than only a heuristic template [2105.08170].

The physical interpretation is that AMLIP does not stabilize walking by canceling the unstable continuous-time pendulum; it stabilizes walking by making the impact map and step placement reshape the state from one support phase to the next. In the canonical formulation, no ankle torque is used to actively change \(L\) during a single support phase; all corrections to \(L\) occur via choosing the next foot placement [2105.08170].

## 4. Extensions to non-stationary, constrained, and time-varying walking

Several later variants retain the angular-momentum state but relax the assumptions that made the original model time-invariant and point-contact based.

**Constrained footholds and impact momentum injection.** “Bipedal Walking on Constrained Footholds: Momentum Regulation via Vertical COM Control” derives an underactuated AMLIP in which the discrete impact map couples pre-impact vertical CoM velocity to post-impact angular momentum:
\[
x^+ = x^- - x_{\rm sw}^-,
\qquad
L_y^+ = L_y^- + x_{\rm sw}^- \dot z^- - z_{\rm sw}^- \dot x^-.
\]
When the swing foot strikes at ground level, \(z_{\rm sw}^-=0\), this reduces to
\[
\Delta L_y = x_{\rm sw}^- \dot z^-.
\]
The control strategy chooses a desired post-impact momentum from a viable LIP orbital energy, computes the required pre-impact vertical velocity, and realizes it through a small CoM-\(z\) QP combined with a task-space whole-body QP [2104.10367]. Under the reported settings, AMBER achieved stable periodic walking speeds from \(0.4\,\mathrm{m/s}\) up to \(1.6\,\mathrm{m/s}\), traversing both fixed-pattern and randomly varying stepping stones for thousands of steps with no failures, while Cassie in simulation robustly ascended and descended staircases and crossed random stones at \(\sim 0.8\,\mathrm{m/s}\) [2104.10367].

**Moving rigid surfaces and linear time-varying ALIP.** “Time-Varying ALIP Model and Robust Foot-Placement Control for Underactuated Bipedal Robot Walking on a Swaying Rigid Surface” extends ALIP from stationary to horizontally moving surfaces with known trajectory \(x_S(t)\). The continuous reduced-order dynamics become
\[
\dot x_{SC} = \frac{L_S}{mH} - \dot x_S(t),
\qquad
\dot L_S = mg x_{SC},
\]
which is a linear time-varying, non-homogeneous ODE. Foot placement enters through the impact reset
\[
\Delta \mathbf{x} = \begin{pmatrix}-u\\0\end{pmatrix},
\]
and the step-to-step controller takes the form
\[
u = u^* + \mathbf{K}\bigl(\mathbf{x}^- - \mathbf{x}^*\bigr).
\]
With the monodromy operator \(M=[I+B]\Phi(T,0)\), the cited theorem states that if all eigenvalues of \(M\) lie strictly inside the unit circle, then the zero solution of the homogeneous hybrid system is exponentially stable and the unique \(T\)-periodic solution of the non-homogeneous hybrid system is exponentially stable [2210.13371].

**Variable pendulum length and stair climbing.** “Stair Climbing using the Angular Momentum Linear Inverted Pendulum Model and Model Predictive Control” replaces constant \(z\) by a known nominal \(z(t)\) or \(r(t)\), preserving linearity in the state while making the system linear time-varying. The key equations are
\[
\ell = m(z\dot x - x\dot z),
\qquad
\dot \ell = mg(x-p) + \tau,
\qquad
\ddot x = \frac{g}{z(t)}(x-p) + \frac{1}{m z(t)}\tau,
\]
or, in angular form,
\[
\dot\theta = \frac{\ell}{m r(t)^2},
\qquad
\dot\ell = mg\,r(t)\,\theta + \tau.
\]
The paper combines virtual-constraint-based posture control with an MPC that acts on the stance ankle torque, using \(\Delta t \approx 1\,\mathrm{ms}\), prediction horizon \(N=T/\Delta t \approx 400\), and ankle-motor limits \(u_{\min/\max}=\pm 23\,\mathrm{Nm}\) [2307.02448]. In simulation on a 20 degree-of-freedom model of Cassie, turning off ankle torque causes falls within two steps on flat ground with fixed-step length, whereas with MPC ankle torque the robot remains stable despite large \(20\%\) torque-drop perturbations and climbs stairs for \(10+\) steps while \(\ell\) and \(\theta\) track nominal periodic profiles [2307.02448].

## 5. Embedding AMLIP in full-body humanoid control

A central question is whether AMLIP survives once the robot has flat feet, non-passive ankles, nontrivial torso and limb inertia, and non-negligible centroidal angular momentum. “Moving past point-contacts: Extending the ALIP model to humanoids with non-trivial feet using hierarchical, full-body momentum control” addresses exactly this regime [2408.05308]. The paper starts from the floating-base rigid-body dynamics
\[
M(q)\ddot q + C(q,\dot q) + G(q) - J_c(q)^\top \lambda = S^\top \tau,
\]
with no-slip contact constraints
\[
J_c \ddot q + \dot J_c \dot q = 0,
\]
and the centroidal momentum mapping
\[
h_{\mathrm{CoM}}=
\begin{pmatrix}
L_{\mathrm{CoM}}\\
m\dot p_{\mathrm{CoM}}
\end{pmatrix}
=
A_{\mathrm{CoM}}(q)\dot q.
\]
The control architecture uses a hierarchical task-space whole-body controller to regulate centroidal momentum so that \(L_{\mathrm{CoM}}\to 0\), \(p_z\to H\), and the horizontal CoM follows ALIP trajectories, while the desired foot placement is tracked as a lower-priority task [2408.05308].

The significance of this construction is not that the original ALIP assumptions become true; rather, the full-order robot is coerced to match the reduced ALIP dynamics through contact wrench modulation and momentum regulation. In MuJoCo simulations on the Sarcos Guardian XO robot, which has mass \(150\,\mathrm{kg}\), flat offset feet, and substantial limb and torso inertia, the measured CoM closely follows the ALIP-prescribed trajectory with RMS error \(<0.02\,\mathrm{m}\), centroidal angular momentum is driven to near zero of order \(10^{-3}\,\mathrm{kg\cdot m^2/s}\), predicted end-of-step contact angular momentum matches measured values within \(<5\%\), and the \(1\,\mathrm{s}\) moving-average of \(\dot p_x,\dot p_y\) stays within \(\pm 0.02\,\mathrm{m/s}\) of the commanded speeds for stepwise changes up to \(\pm 0.45\,\mathrm{m/s}\) forward and \(\pm 0.225\,\mathrm{m/s}\) lateral [2408.05308].

A related, but distinct, development expands the Divergent Component of Motion framework to include angular coordinates by combining classical 3D linear DCM with a 1D angular DCM derived from a flywheel-approximation SRBM. The resulting \(3\)D linear and \(1\)D angular DCM framework has open-loop unstable dynamics and decoupled exponentially stable closed-loop subsystems under VRP and VRO feedback. MATLAB simulations and TORO humanoid experiments report angular error \(\|\theta-\xi_{\rm a}\|_\infty<0.02\,\mathrm{rad}\), CoP excursion \(<4\,\mathrm{cm}\), ground reaction torque matching the predicted \(\tau_{\rm ext}=I\omega_a^2(\theta-\phi_{\rm vro})\) within \(5\%\), and simulation-to-hardware agreement within \(10\%\) over a \(5\,\mathrm{s}\) trial [2409.12796]. This suggests one route by which AMLIP-style angular-state reasoning can be connected to spatial DCM planning.

## 6. Advantages, common misconceptions, and limitations

A recurring misconception is that AMLIP is only the LIP written in different coordinates. Formally, the continuous equations are identical after the substitution \(v\mapsto L/(mH)\), but the cited works emphasize several structural consequences of the momentum coordinate. In the planar constant-height setting, \(\dot L = mgx\) depends only on CoM position; \(L^y\) is invariant across flat impacts when \(v_z=0\); and the angular-momentum output has relative degree \(3\), yielding weak coupling to joint-torque peaks and swing-leg disturbances [2008.10763]. On Cassie, one-step-ahead angular-momentum prediction error remained within \(\pm 2\%\) over each step at \(2\,\mathrm{m/s}\), whereas the standard LIP-based one-step velocity prediction would have exhibited prediction errors in excess of \(20\%\) under the same conditions in simulation comparisons [2105.08170].

A second misconception is that AMLIP necessarily implies direct ankle-torque regulation of the reduced state. In the canonical step-to-step AMLIP of [2105.08170], no ankle torque is used to actively change \(L\) during a single support phase; all regulation is accomplished through foot placement. Other papers deliberately break that restriction when the task requires it: stair climbing uses stance-ankle torque as the control input in an LTV ALIP-MPC controller [2307.02448]; constrained-foothold walking uses pre-impact vertical CoM velocity to inject the desired angular momentum at impact [2104.10367]; and flat-foot humanoid locomotion uses full contact wrenches and centroidal-momentum regulation to force a full-order robot to behave like an ALIP template [2408.05308].

The model’s limitations are equally explicit in the literature. The original formulations assume constant CoM height, zero or negligible angular momentum about the CoM, pinned or point-foot contact, and zero vertical CoM velocity at impact; they neglect ankle torques in prediction and require reliable IMU or EKF estimation of \(L^y\) [2008.10763, 2105.08170]. Stairs violate both the constant-height assumption and unrestricted sagittal foot placement, motivating the variable-length and torque-actuated LTV model [2307.02448]. Moving support surfaces introduce explicit time variation and non-homogeneous forcing into the reduced dynamics [2210.13371]. Flat feet and non-centralized limb inertia invalidate the negligible-centroidal-momentum and zero-contact-moment assumptions unless a whole-body controller actively regulates them [2408.05308]. Impact-driven oscillations and spring dynamics can also generate brief torque spikes not predicted by the simplified ALIP model [2307.02448].

Within those bounds, AMLIP has served as a compact analytical template that links continuous pendulum dynamics, hybrid impact maps, foot-placement control, whole-body momentum regulation, and zero-dynamics interpretations. The common thread across its variants is not a single fixed equation set, but the use of angular momentum about the contact or stance point as the reduced coordinate through which walking balance, velocity regulation, and disturbance rejection are organized [2008.10763, 2105.08170].

Source: https://www.emergentmind.com/topics/angular-momentum-linear-inverted-pendulum-model-amlip