Pendular models for walking over rough terrains
Abstract
The NewtonEuler equations that govern multibody systems are not integrable in general. However, they become so in the pendular mode, a specific way of moving where conservation of the angular momentum is enforced. This property was successfully showcased for locomotion over horizontal floors (2D locomotion) by walking pattern generators based on the LIPM and CARTtable models. In this talk, we will see how to generalize these two models to 3D locomotion while taking into account both friction and tilted contacts, resulting into the FIP (3D version of LIPM) and COMaccel (3D version of CARTtable) models. We will demonstrate both approaches in live simulations with the HRP4 humanoid model.
Content
Slides of the UniRoma version (October 2017)  
Slides of the JNRH version (June 2017) 
References
Paper on the COMaccel model (Humanoids 2016)  
Paper on the FIP model (IROS 2017) 
Discussion
Thanks to all those who have contributed to the conversation so far. Feel free to leave a reply using the form below, or subscribe to the Discussion's atom feed to stay tuned.

Attendee #1
Posted on
In this presentation, you linearized friction cones using the innerpyramid approximation. What is the cost of this approximation? Why not using the full secondorder friction cones?

Stéphane
Posted on
To answer the first part of your question, in my opinion the cost of this approximation is small: the forces we miss in the corner volumes (between the secondorder cone and its inner pyramids) are extreme values anyway, and we want to avoid them in practice. This is to be balanced with the benefits of working with polytopes:
 The class of optimization problems becomes LP or QP, and these problems are usually solved faster than the SOCPs you get with secondorder cones. (It is not true in general that "SOCP is slower than LP/QP", but in my experience it is the case in force optimization.)
 Polyhedral manipulations allow you to compute more complex objects, such as torquelimited friction cones.
By the way, secondorder friction cones are themselves an outer approximation of a more complex phenomenon, see e.g. CoulombOrowan pseudocodes.
The debate is not closed at any rate. As a matter of fact it is still ongoing in my team as well, where Hervé Audren uses SOCPs in his computation of robust staticequilibrium polytopes.


Attendee #2
Posted on
You added friction to your model but for walking on flat floors it is usually neglected. What is the motivation for that?

Stéphane
Posted on
At the end of the day the motivation is always experimental: with common friction coefficients (0.7 and more), the most stringent contactstability constraint in practice is the ZMP one. When you don't need to model friction to get your experiments to work, it makes sense not to model it. Now, this observation becomes less and less true when your contacts are tilted (especially when they are tilted outward since angularmomentum conservation tends to get forces pointing to the CoM), so for general walking over rough terrains you want to consider modelling friction.
The second good reason for doing that is that modelling friction is easy in terms of polyhedra (look e.g. at the condition in the FIP paper). The hard part is always the ZMP condition, which you cannot do without at any rate.


Attendee #3
Posted on
How do you implement swing foot trajectories in your solution?

Stéphane
Posted on
We interpolate Hermit curves with parameters selected to avoid sharp foot accelerations (HOUBA curves, detailed in Section IV of this paper. There is nothing special here, actually I believe a better way to do these polynomial interpolations is the one given in Section III of the paper Footguided Agile Control of a Biped Robot through ZMP Manipulation by Sugihara and Yamamoto (IROS 2017).


Attendee #4
Posted on
In the pendular model you keep a constant angular momentum, while other motion generation techniques [1, 2] also optimize over kinematics to take it into account. What's missing in your solution to make the angular momentum a proper control?

Stéphane
Posted on
Angular momentum is definitely one of the next big questions. It is an open question even for 2D walking, while here we are doing 3D walking over rough terrains! There are a number of properties that explain why it is not as straightforward as COM control: its nonholonomy (see Wieber (2009)) and dependence on the full kinematic state (joint angles and velocities). Hence the idea of a kinematics optimization coupled with centroidal dynamics.
To the best of my knowledge today, these solutions have been applied so far to motion generation (i.e. with more computation time, not in a closed control loop). When developing the two walking pattern generators presented here, the first problems were to be (1) general enough for 3D walking and (2) fast enough for the control loop. Now that we've made a first step here, I'll be definitely looking forward to the angular momentum in the future.
