Omnidirectional Walking Pattern Generator Combining Virtual Constraints and Preview Control for Humanoid Robots

Abstract

This paper presents a novel omnidirectional walking pattern generator for bipedal locomotion combining two structurally different approaches based on the virtual constraints and the preview control theories to generate a flexible gait that can be modified on-line. The proposed strategy synchronizes the displacement of the robot along the two planes of walking: the zero moment point based preview control is responsible for the lateral component of the gait, while the sagittal motion is generated by a more dynamical approach based on virtual constraints. The resulting algorithm is characterized by a low computational complexity and high flexibility, requisite for a successful deployment to humanoid robots operating in real world scenarios. This solution is motivated by observations in biomechanics showing how during a nominal gait the dynamic motion of the human walk is mainly generated along the sagittal plane. We describe the implementation of the algorithm and we detail the strategy chosen to enable omnidirectionality and on-line gait tuning. Finally, we validate our strategy through simulation experiments using the COMAN + platform, an adult size humanoid robot developed at Istituto Italiano di Tecnologia. Finally, the hybrid walking pattern generator is implemented on real hardware, demonstrating promising results: the WPG trajectories results in open-loop stable walking in the absence of external disturbances.

1 Introduction

The algorithm results in a light-weight, omnidirectional WPG¹ based on a simple template model that can be easily deployed on any humanoid robot (as shown in Figure 1). This method allows for a more dynamic motion without dropping the advantages of the PC and the ZMP criterion. Furthermore, redesigning the strategy for synchronization not only allows to drop a heuristic approach, but provides for the implementation of omnidirectionality in the WPG. The template model is mapped on the robot using the whole-body inverse kinematics (IK) solver by specifying a set of desired Cartesian tasks sorted by priority, so as to impose a desired body posture besides directly controlling the model variables as high priority tasks, i.e., the feet trajectories and the composite motion of the sagittal and the lateral components of the CoM. Finally, we discuss the advantages and the limitations of this approach and we propose future improvements.

FIGURE 1

2 Related Works

The first paradigm formulated for humanoid walking is known as static stability, which restricts the projection on the ground of the center of mass (CoM) to lay inside the support polygon at all times. In this way, a static walk can be simply synthesized by slowly moving the CoM from one foot to another (Kato et al., 1974). The seminal work of Kajita (Kajita et al., 2003) demonstrated how, using the Preview Control theory and leveraging on the notion of the ZMP, a more dynamic, stable gait can be achieved. The effectiveness of this approach made it widespread among the humanoid locomotion community and various improvement were proposed: the LIPM model was extended to reduce the modeling error and generate a robust pattern locomotion (Shimmyo et al., 2013). Application that enhance the preview control reacting against external perturbations (Nishiwaki and Kagami, 2006; Nishiwaki and Kagami, 2009) or adapting to uneven terrains (Kajita et al., 2006). The simple yet successful strategy is usually framed into an Model Predictive Control (MPC) formulation, which is still the most widely adopted: a method to guarantee its intrinsic stability was presented in Scianca et al. (2016), and it was extended for disturbances rejection (Wieber, 2006), push recovery (Stephens and Atkeson, 2010) and automatic footstep placement (Herdt et al., 2010). The new strategy which released the constraints of the ZMP by introducing a different notion of stability based on limit cycles was first proposed in Hürmüzlü and Moskowitz (1986). The theory was further developed in Westervelt et al. (2007) using Hybrid Zero Dynamics (HZD) to actively control the humanoid and enlarge the basin of attraction of a limit cycle gait. This approach yielded successful results, but it is usually limited to robots with particular body structures (Rezazadeh et al., 2015; Gong et al., 2019). Some strategies were explored to generalize this method to the 3D case, such as the functional Routhian reduction introduced in Chevallereau et al. (2014), or the self-synchronization exploiting the symmetry of the system, formalized in Razavi et al. (2015). In fact, while producing dynamical and energy-efficient gaits, this method is not trivially generalized to the 3D case, where hybrid invariant manifolds are more challenging to describe. Nonetheless, some techniques were developed to implement 3D walking on different robots (Da et al., 2016; Reher et al., 2016), most of them decoupling the walking gait into forward and lateral motion. Furthermore, most of the strategies solely consider straight walking gaits, as steering cannot be trivially integrated. While some elegant solution can be found in literature, they are mostly theoretical results, usually verified only in simulation (Motahar et al., 2017; Chevallereau et al., 2010).

The proposed hybrid WPG relies on the synchronization of the sagittal and lateral planes of the robot, as in the previous works. However, it merges results from the VC and the PC theory: it exploits a more dynamical approach that only constrains the ZMP trajectory to the support polygon along the lateral direction. Furthermore, as opposed to the classical implementation of the VC, our algorithm does not heavily depend on the kinematics of the robot and does not require burdensome parameter tuning: template trajectories can be designed and realized on-line on the real robot, making it a computationally lightweight and flexible tool independent of the structure of the walker.

3 Sagittal Plane: Virtual Constraints

Virtual constraints are employed to combine the kinematics of many joints into a single template motion. This simplifies their coordination and reduces the dimension of the problem: the desired motion for a n-dimensional system can be encoded by controlling a smaller set of state variables n−n_m, where n_m is the number of constrained variables. On the other hand, its inherent limitation resides in the cyclic nature of this template movement, which evolves along a periodic orbit. While VCs are useful for the generation of repetitive gaits, they are unsuited for any action that alters their evolution, such as turning steps or aperiodic manoeuvres.

To illustrate this, let be the state vector of a n-link bipedal mechanism and a collection of m internal variables bound together by the VCs , to encode a specific locomotion. One VC is imposed on the mechanical system by zeroing the output of the function:so that evolves along the periodic orbit h_d as function of , known as phase variable. The phase variable, in turn, is defined as a function of a set of r variables , which, oppositely to the controlled variables, remains free. While monotonically increases from to , the controlled variables are forced along a desired trajectory that translates into a specific stepping motion, as shown in Figure 2. Usually the output in (1) is zeroed by a feedback controller, which asymptotically brings the state variables towards the desired orbit h_d. In the work (Westervelt et al., 2007) is shown how suitable VCs can produce cyclic motions for a walker: the robot is treated as a point-foot template model, and the evolution of the internal joints is constrained as a function of the angle between the ground and the virtual leg connecting the stance foot and the CoM. Our strategy is based on the same formulation, which we detail in the following paragraph.

FIGURE 2

3.1 Template-Space Constraints

The template model of the humanoid robot is the 2-link compass walker, with its CoM located at the joint connecting the two links of equal length, as shown in Figure 1. The configuration vector is , where and are the angles describing respectively the stance and swing leg. The selected VC exploits the symmetry of the model by forcing the swing leg to behave as a mirror image of the stance leg. Its formulation for the compass geometry is the following:

According to (2), is selected as the phase variable , while is the controlled variable constrained to the evolution of . Finally, z is the output to be driven to zero.

3.2 Task-Space Constraints

The control of our humanoid robot relies on a whole-body hierarchical inverse kinematics framework: multiple Cartesian tasks can be defined and organized in a Stack of Tasks (SoT) fashion. Therefore, instead of defining the VC in the joint-space, we describe it as a function of Cartesian variables. Similarly to the formulation in the template-space, we are interested in binding relevant quantities of the robot (i.e., the CoM and the feet) with the phase variable α². In the real robot, we choose α as the angle between the normal of the ground and the virtual leg connecting the CoM to the stance foot, as shown in Figure 1, which corresponds to the tilt angle. As increases, the robot tilts on his stance foot in the sagittal direction, resulting in a CoM displacement c:where h is the height of the CoM w.r.t. the ground. By formulating the step displacement s as function of the configuration vector we obtain:

Finally, substituting in (2), yields s as a function of the tilt angle only:

This VC restricts the sagittal displacement of the swing foot to keep the CoM at a fixed distance from each leg: when the robot tilts forwards, the phase variable α increases, constraining the swing foot to perform a stepping motion until it impacts the ground, resetting the cycle. The VC does not constrain the trajectory of the leg, but only its overall displacement: thus, it can be treated as a task for the whole-body IK solver. This allow to tune the parameters of the stepping: We design the motion as a polynomial primitive with a peak of a desired height, which corresponds to the step clearance.

Following the theory in (Westervelt et al., 2007), the template model assumes fully inelastic collision, i.e. the landing foot does not experience slips or bounces at the contact point. For this reason, we impose zero-velocity of the swing leg at the beginning and at the end of the step, to reduce high impacts resulting in undesired rebounds of the foot on the ground that could disrupt the walking cycle. Finally, the impact is instantaneous: the stance leg is lifted at the same time the swing foot touches the ground, without any double stance phase during the gait.

3.3 Ankle Actuation

Once a suitable VC is applied, the template model can be treated as a simplified system with only one DoF, the tilt angle. Its evolution, if not controlled, injects the free-fall dynamics of the stance leg, modeled as a linear inverted pendulum (LIP), into the system. This free evolution exploits the natural dynamics to reduce energy consumption during the stepping motion, but limits the flexibility of the gait, which is bound to a non-controlled variable. The feet of our humanoid are composed of flat soles and fully actuated ankles. Differently from the point feet of the template model, this allows to directly control the tilt angle of the robot and, as a consequence, the phase variable α. The CoM can be fully controlled as long as the ZMP remains inside the support polygon (the implications of this choice are discussed in Appendix). When on the edge, the robot becomes underactuated and starts tilting with its natural dynamics. Explicitly controlling the sagittal tilt precludes the possibility to exploit the natural dynamics of the robot in the forward motion to reduce the power consumption of the gait. On the other hand, dealing with a fully controllable system is advantageous, since it can be easily regulated to impose a desired behaviour. For instance, when steering, changing the feet distance or increasing the step stride. Finally, by fully controlling the CoM, it is possible to define a consistent system along the two planes to guarantee that the lateral plane can rely on the LIPM: without violating the Cartesian constraints between the sagittal displacement of the CoM and the swing leg, we impose a constant height to the CoM trajectory, realized by the HIK solver, so that the corresponding system in the lateral plane can be modeled as a LIPM.

4 Lateral Plane: Preview Control

The strategy exploited for the sagittal plane is simple to define in the planar case and allows a more dynamical motion, but it is not trivial to extend for 3D robots. Furthermore, the work in (Sardain and Bessonnet, 2004) demonstrates how the dynamic component of the walking is found along the sagittal plane of the foot: during a nominal gait, the ZMP travels mainly in the sagittal direction until it reaches the edge of the foot, without moving considerably along the lateral axis. For these reasons, we choose a conservative approach, as we rely on the well-grounded theory of preview control (Kajita et al., 2003): by treating the robot as a LIPM, we can impose a desired ZMP and track it by controlling the trajectory of the CoM.

In particular, the cart-table model describes a dependency between the projection along the lateral plane of the CoM and the ZMP :where is the CoM height w. r. t the ground and g is the magnitude of the gravity vector.

Eqn. 6 is used to formulate the desired output-tracking problem. A common choice is to select the state vector as:while setting the control input as the lateral CoM jerk . By discretizing via zero-order-hold, We obtain a discrete-time LTI system, with time step T:whereand the output, according to (6) is:

As shown in (Kajita et al., 2003), the solution to this output tracking problem requires knowledge of the desired ZMP over a finite future horizon. We realize the anticipative action by framing the problem as a finite-horizon LQR. In particular, the optimal controller is formulated as follow:where is the desired ZMP at the kth future time step, and r > 0 is a weight that penalizes the control action. The QP problem (11) consists in minimizing the state and the input trajectory over the horizon of length N > 0, constrained to the dynamics of the system (8).

A solution is found via the KKT equations, which take the form of a highly sparse set of n = 4N + 3N equations in as many unknowns, solved through a sparse LU decomposition from the Eigen3 library. Finally, at each iteration the first control input u_t is applied to the system (8); the resulting control law is linearly dependent on the current state x_t and the future reference ZMP trajectory :where is a state feedback matrix, whereas is a reference feed-forward matrix. The state vector x is not updated with the real state of the robot, hence its future value x_t+1 is obtained solely by integrating the current over a chosen period T, subject to the constant control input u_t. The length N of the prediction horizon affects the tracking of the ZMP: a high N slows down the computation without greatly improving the performance, while a low N generates overshoots in the ZMP tracking. The system is kept stable thanks to the feedback term, hence computing the closed-loop eigenvalues magnitude is useful to assess the stability of the system:if the stability criterion (13) is not met, more of the future ZMP reference should be fed to the controller, increasing the length of the preview horizon, while ensuring that the resulting computational burden introduced by the length of the horizon can be handled within one control cycle of the system.

5 Synchronization

During a nominal walk, the motions along the two planes are strictly coupled together, since a sagittal displacement corresponds to a lateral movement of the CoM. In our previous work (Ruscelli et al., 2019), a strategy for synchronization between the lateral plane and the sagittal plane was proposed, showing how two structurally different control strategies can generate feasible references for walking. However, while being a viable method, it implied jumps or temporal freezes of the receding preview window: while the sagittal motion depended on the phase variable α, the ZMP window was dependent on time, and was reset at each step to synchronize with the stepping motion. Furthermore, it presented significant limitations for the generation of the future ZMP, as the receding window had a fixed ZMP pattern with predetermined values for the left and right foot. In this work, we propose a novel method for the synchronization of the planes. In the remaining of this Section, we outline the behaviour of the two independent component and we illustrate how they are coupled.

Sagittal plane. The sagittal displacement is dependent on the phase variable according to (3) and (5). A full cycle, from the take off of the swing foot to its landing on the ground entails increasing monotonically from to , as shown in Figure 2.

Lateral plane. One complete cycle in the lateral plane corresponds to the advancement of the preview window: the ZMP reference resides inside the sole of the stance foot until the impact with the ground, when it switches to the opposite leg. We remove the dependency of the preview window on time and bind its advancement to the phase variable . In particular, in one cycle the receding window slides from to , as shown in Figure 3. The WPG automatically generates the values of the ZMP reference: the first segment of the preview window is the y-component of the position vector of the current stance foot, while the second segment is the y-component of the goal position of the swing foot. The remaining tail of the preview window is filled with a repeating pattern based on the ZMP values of the current and next steps.

FIGURE 3

Synchronization. The phase variable is measured at each iteration, and its evolution coupled with . Specifically, is directly driven by , as its progression coincides with the monotonic increment of the sagittal phase variable, as depicted in Figure 4: the two planes are automatically synchronized, and the advancement of the CoM and the swing foot in the sagittal direction drives the lateral swing of the CoM along the full stepping cycle.

FIGURE 4

6 Implementation

The proposed algorithm cyclically generates reference trajectories for the CoM and feet during a single step: their evolution on each plane, from taking off of the swing foot to its landing, is described by the phase variables and . The structure of the WPG is shown in Figure 5.

FIGURE 5

6.1 Online Walking Pattern Generator

The CoM is directly controlled through the phase variable : corresponds to the maximum inclination of the robot, which translates into the length of the step according to the VC in (5): by changing the interval a desired stride length can be imposed on the gait. If a new is issued, the CoM velocity during the step changes accordingly, as shown in Figure 6.

FIGURE 6

Notice that, due to the imposed VC, the CoM lies always in the middle of the two feet in the sagittal direction, and when , the CoM is directly above the stance foot, as shown in Figure 2.

These enhancements entail two other necessary improvements: first, footsteps are no more planned beforehand, but automatically generated on-line during the previous step given the maximum inclination and the angle of curvature ϕ. Secondly, the preview window is continuously updated with a suitable ZMP reference to accommodate the gait, as shown in Figure 4. The proposed WPG is enclosed in a state machine to manage the walking phases, such as starting and stopping. Finally, the state machine exposes a simple interface for the communication between user and robot.

6.2 Whole-Body Inverse-Kinematics

The trajectory are generated on-line by a whole-body IK solver in a SoT fashion (Rocchi et al., 2015; Laurenzi et al., 2019). This allows to run low priority task in the null space of higher priority ones, so as to guarantee the execution of critical tasks above the others. The following SoT was designed for the WPG:

The / symbol is used to impose strict hierarchy among sets of tasks. Among the same level, tasks are listed via the + symbol, which impose soft hierarchy. The description define a task of the frame B expressed w.r.t. the frame A. If A is not specified, task is expressed in joint space. In the chosen stack, refers to a generic task B w.r.t. the world frame. The Torso task reduces the swing of the robot at each step by keeping the torso perpendicular to the ground. The contact tasks and are at the highest priority level. All the other tasks, including the CoM , act in the null-space of the contacts to guarantee a consistent solution within the base under-actuation. The pose of the joints not involved in the walking motion is maintained thanks to the Posture task, located in the null space of all the other defined tasks. Finally, the symbol is used to impose constraints such as joint limits and joint velocity limits.

7 Omnidirectional Walking

The previous WPG proposed in (Ruscelli et al., 2019) generates a forward walk without the ability to steer to another direction or walk backwards: due to the inherent decoupling between lateral and sagittal planes, the hybrid WPG is not trivially generalized for omnidirectional walking. In this section, we describe the strategy devised to extend the WPG for steering in any given direction, as shown in Figure 7.

FIGURE 7

For a three-dimensional robot, we define the tilt angle as:where and are the x and z components of , i.e., the stance foot ankle w.r.t. a reference frame centered at the CoM and oriented with the x-axis aligned to the direction of motion. When a steering angle ϕ is issued, a steering step is performed to change the direction of motion from the heading angle ψ to the new one , and the reference frame at the CoM is rotated by ϕ along the z-axis. The position of the stance ankle is recomputed according to the rotation:where is the rotation matrix representing a rotation by ϕ on the horizontal plane (Figure 8). This cause the tilt angle to change during the steering motion: once the rotation (16) is applied, the new is computed at the beginning of the steering step by applying (15).

FIGURE 8

the CoM moves forward in the sagittal plane until α reaches and the foot impacts the ground. At the impact, the new angle of heading ψ is updated with the angle of steering ϕ, and the walk resumes in the new direction.

The distance d between the right and the left sole in the lateral plane and the tilt in the sagittal plane (which corresponds to the step stride) are kept constant before and after the curve. During the steering, due to the rotation of the reference frame, both d and are adjusted to generate the desired step: the ZMP pattern fed to the preview control and the phase variable are updated as shown in Figure 9.

FIGURE 9

8 Experiments and Discussion

The walking pattern generator was tested in simulation using Gazebo and validated on COMAN+, a humanoid robot developed at Istituto Italiano di Tecnologia. COMAN + has 28 DoFs, weights 70 kg and is 1.7 m tall. A particular four-bar mechanism controls the pitch and the roll at the ankle level of the robot (Ruscelli et al., 2018). The WPG prescribes position references that are tracked on the robot by a whole-body IK framework named CartesI/O (Laurenzi et al., 2019) based on OpenSoT Rocchi et al. (2015), which allows to specify the desired SoT solved by quadratic programming (QP) optimization. COMAN+ is powered by the XBotCore real-time software architecture (Muratore et al., 2017). A ROS pipeline allows communication between user WPG and lower level. The WPG runs at 100 Hz in a ROS node and exposes a simple ROS interface for issuing commands to the robot.

8.1 Experiments

A set of experiments was carried out both in simulation and on real hardware for the robot COMAN+, in order to validate the WPG. We detail here one meaningful example in simulation. COMAN + walks 32 steps, with a 0.05 m step clearance, a 0.4 s step duration and a maximum tilt of 0.08 rad, which corresponds to a 0.28 m step stride. Figures 10,11 shows the projection of the walking pattern on the x-y plane and the evolution of the soles in space, respectively: the CoM oscillates laterally between the legs, from the stance towards the swinging foot. The faster the gait velocity, the smaller the lateral swings of the CoM. During step 5 and 13 a new is issued, while at step 19 and 27 the robot steers back of to the straight direction. The WPG can steer up to 30 deg in a single step. The maximum steering angle is limited by the hybrid nature of the proposed architecture: during the steering step the preview control only works in the previous lateral direction, and a component of the lateral CoM trajectory will not be guaranteed to belong to the support polygon. The maximum speed of the robot is 1.25 m/s in simulation when moving forward, while the maximum step length change is 0.4 m (from 0.1 to 0.5 m). Figure 12 depicts the time evolution of the ZMP reference in the lateral plane and the phase variable along the sagittal direction. In particular, the preview controller tracks the reference ZMP, which is continuously updated during the stepping motion, and outputs the lateral CoM trajectory. The phase variable constrains the evolution of the step, spanning from to at each cycle after an impact occurs, and it is synchronized with the evolution of the phase variable in the lateral plane. A sequence of frames of the experiment is depicted in Figure 13. The WPG was deployed on real hardware to assess its effectiveness (Figure 14). The experiment consisted of a 10 step walk with 0.23 m stride length. The selected step clearance was 0.05 m and the step duration of 0.4 s. A steering angle of was issued at step 3. While the algorithm generates feasible trajectories for the desired gait, the WPG does not guarantee a stable walking: the reference trajectories are open-loop, since no stabilization controller is added to react against an unwanted loss of balance or reject external disturbances. Videos of the experiments both in simulation and on the real robot can be found at https://www.youtube.com/playlist?list=PL7c1ZKncPan7yphDxvZDtaqzmgCRtpT5-.

FIGURE 10

FIGURE 11

FIGURE 12

FIGURE 13

FIGURE 14

8.2 Discussion

The proposed WPG is a lightweight algorithm that can be easily deployed on real hardware. The classical implementation of VC entails joint-trajectory optimization to find energy-efficient gaits. Often the optimization is run offline individually for each desired gait that corresponds to one optimal orbit (or a family of similar orbits). Providing a walker with different gaits usually requires a library of optimal parameters collected beforehand. This leaves less margin to tune a walking pattern or change it on the fly, hindering the versatility of the method. Furthermore, the VC usually binds joint variables, making it heavily dependent on the mechanism kinematics. Instead, we treat the VC as Cartesian tasks in a SoT and we design template trajectories that can be tuned and realized on line by the whole-body IK solver. This allows a light and flexible tool for stepping that only depends on a simple template model and a few step parameters, such as length, duration and clearance of a step. The new WPG also allows to modify the heading angle and the distance between feet. According to the belief that a more dynamic approach can be used along the sagittal plane while a more conservative one can be exploited laterally, in our strategy the CoP is guaranteed to lay inside the support polygon only in the lateral direction, while in the sagittal one it can move freely. On the other hand, our strategy has a higher energy consumption with respect to the classical approach, because it doesn’t exploit the natural dynamic of the system during a step. In fact, the trajectory of the CoM is fully controlled. Right now the commanded trajectory is linear, but a possible solution is finding an optimized trajectory of the CoM to minimize a desired cost function. Notice that the CoM trajectory during a step is directly related to the tilt of the robot, which corresponds to the phase variable. Hence, modifying the CoM trajectory amount to changing the behaviour of the stepping. Another limitation of the strategy is the amount of modification that the gait can sustain in one single step: being a cyclic motion, abrupt changes in the direction of movement or in the step length will result in the robot losing balance. Since the work was only focused in the development of a WPG, a feedback term of the measured state of the robot was not included: hence no information from the actual robot is exploited in the algorithm.

The major downside of this choice is reducing the dynamic component of the walking motion in the sagittal direction. While part of the walking can be controlled directly imposing a CoM trajectory, the portion of the gait which entails a “controlled” fall can be exploited only given a proper state estimation: sensing the evolution of the tilt angle injects into the system the free-fall dynamics of the stance leg, used by a continuous feedback controller to restrict the stepping motion on a periodic orbit. Finally, while simulation experiments show successful results, their hardware counterpart is not as satisfying: this discrepancy is mainly due the lack of a stabilization effect that becomes necessary for real robot experiments. While the open-loop WPG plans feasible trajectories, as discussed in Appendix and demonstrated in simulation where the factors acting as disturbances can be limited, a control layer is essential when deploying the WPG in a real-word scenario. Nonetheless, including hardware results highlights how this strategy is viable, displaying promising results towards a WPG that integrates a planning and a control layer.

9 Conclusions and Future Work

The proposed walking pattern generator, due to its hybrid nature, is not trivially generalizable for omnidirectional gaits: in order to do so, we improved the synchronization of the lateral and sagittal planes, enabling on-line modifications to the walking gait. The resulting algorithm is a lightweight, omnidirectional WPG which generates feasible whole-body trajectories according to a few user-parameters. Besides off-line specifications such as duration, length and stride of the desired step, the proposed WPG allows on-line commands to change heading angle, step length and feet distance. The sagittal motion of the CoM is directly commanded as a simple linear trajectory. In future works, we will use optimization techniques to select a trajectory that minimizes a desired function such as power consumption. Simulation experiments demonstrate the effectiveness of our strategy: the robot walks and changes the nominal gait according to user-inputs without falling. Furthermore, preliminary experiments on the humanoid COMAN + shows how the proposed hybrid WPG is promising, successfully allowing the robot to change its heading direction. However, without a stabilization controller, any external disturbance (e.g., uneven or sloped terrain) is enough to make the robot lose balance after a few steps.

Similarly, the maximum steering angle that can be issued on the real robot is limited, as the weaknesses highlighted in the straight walk are sharpened when changing direction.

In future works, we will close the loop with the external environment by implementing stabilization techniques to add robustness to the gait.

References

• 1

  ChevallereauC.GrizzleJ. W.Ching-Long ShihC. (2010). Steering of a 3d Bipedal Robot with an Underactuated Ankle. IEEE/RSJ International Conference on Intelligent Robots and Systems, Taipei, Taiwan, October 18–22, 2010. 1242–1247. 10.1109/IROS.2010.5648801

  • CrossRef
  • Google Scholar

• 2

  ChevallereauC.SinnetR. W.AmesA. D.UniversitéL. (2014). Models, Feedback Control, and Open Problems of 3D Bipedal Robotic Walking. Automatica.50(8), 1955–1988. 10.1016/j.automatica.2014.04.021

  • CrossRef
  • Google Scholar

• 3

  DaX.HaribO.HartleyR.GriffinB.GrizzleJ. W. (2016). From 2d Design of Underactuated Bipedal Gaits to 3d Implementation: Walking with Speed Tracking. IEEE Access4, 3469–3478. 10.1109/ACCESS.2016.2582731

  • CrossRef
  • Google Scholar

• 4

  GongY.HartleyR.DaX.HereidA.HaribO.HuangJ.et al (2019). Feedback Control of a Cassie Bipedal Robot: Walking, Standing, and Riding a Segway. In American Control Conference, Macau, China, November 3–8, 2019. 4559–4566.

  • Google Scholar

• 5

  HerdtA.DiedamH.WieberP.-B.DimitrovD.MombaurK.DiehlM. (2010). Online Walking Motion Generation with Automatic Footstep Placement. Adv. Robotics24, 719–737. 10.1163/016918610x493552

  • CrossRef
  • Google Scholar

• 6

  HürmüzlüY.MoskowitzG. D. (1986). The Role of Impact in the Stability of Bipedal Locomotion. Dyn. Stab. Syst.1, 217–234. 10.1080/02681118608806015

  • CrossRef
  • Google Scholar

• 7

  KajitaS.KanehiroF.KanekoK.FujiwaraK.HaradaK.YokoiK.et al (2003). Biped Walking Pattern Generation by Using Preview Control of Zero-Moment point. IEEE Int. Conf. Robotics Automation2, 1620–1626.

  • Google Scholar

• 8

  KajitaS.MorisawaM.HaradaK.KanekoK.KanehiroF.FujiwaraK.et al (2006). Biped Walking Pattern Generator Allowing Auxiliary Zmp Control. IEEE/RSJ International Conference on Intelligent Robots and Systems, Beijing, Taiwan, October 9–15, 2006. 2993–2999. 10.1109/IROS.2006.282233

  • CrossRef
  • Google Scholar

• 9

  KatoI.OhteruS.KobayashiH.ShiraiK.UchiyamaA. (1974). Information-power Machine with Senses and Limbs. On Theory and Practice of Robots and Manipulators. Tokyo, Japan. Springer, 11–24. 10.1007/978-3-7091-2993-7_2

  • CrossRef
  • Google Scholar

• 10

  LaurenziA.HoffmanE. M.MuratoreL.TsagarakisN. G. (2019). Cartesi/o: A Ros Based Real-Time Capable Cartesian Control Framework. 2019 International Conference on Robotics and Automation (ICRA), Montreal, QC, Canada, May 20–24, 2019. 591–596. 10.1109/ICRA.2019.8794464

  • CrossRef
  • Google Scholar

• 11

  MotaharM. S.VeerS.PoulakakisI. (2017). Steering a 3d Limit-Cycle walker for Collaboration with a Leader. IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Vancouver, BC, Canada, September 24–28, 2017. 10.1109/IROS.2017.82064165251–5256.

  • CrossRef
  • Google Scholar

• 12

  MuratoreL.LaurenziA.HoffmanE. M.RocchiA.CaldwellD. G.TsagarakisN. G. (2017). XBotCore: A Real-Time Cross-Robot Software Platform. First IEEE International Conference on Robotic Computing, Taichung, Taiwan, April 10–12, 2017. 77–80. 10.1109/IRC.2017.45

  • CrossRef
  • Google Scholar

• 13

  NishiwakiK.KagamiS. (2006). High Frequency Walking Pattern Generation Based on Preview Control of Zmp. Proceedings 2006 IEEE International Conference on Robotics and Automation, Orlando, FL, United States, May 15–19, 2006. 2667–2672. 10.1109/ROBOT.2006.1642104

  • CrossRef
  • Google Scholar

• 14

  NishiwakiK.KagamiS. (2009). Online Walking Control System for Humanoids with Short Cycle Pattern Generation. Int. J. Robotics Res.28, 729–742. 10.1177/0278364908097883

  • CrossRef
  • Google Scholar

• 15

  PrattJ.CarffJ.DrakunovS.GoswamiA. (2006). Capture point: A Step toward Humanoid Push Recovery. In 2006 6th IEEE-RAS International Conference on Humanoid Robots, Genova, Italy, December 4–6, 2006. 200–207. 10.1109/ICHR.2006.321385

  • CrossRef
  • Google Scholar

• 16

  RazaviH.BlochA. M.ChevallereauC.GrizzleJ. W. (2015). Restricted Discrete Invariance and Self-Synchronization for Stable Walking of Bipedal Robots. Proceedings of the American Control Conference, Chicago, IL, United States, July 1–3, 2015. 4818–4824. 10.1109/ACC.2015.7172088

  • CrossRef
  • Google Scholar

• 17

  ReherJ.CousineauE. A.HereidA.HubickiC. M.AmesA. D. (2016). Realizing Dynamic and Efficient Bipedal Locomotion on the Humanoid Robot DURUS. In IEEE International Conference on Robotics and Automation, Stockholm, Sweden, May 16–21, 2016. 1794–1801.

  • Google Scholar

• 18

  RezazadehS.HubickiC.JonesM.PeekemaA.Van WhyJ.AbateA.et al (2015). Spring-mass Walking with Atrias in 3d: Robust Gait Control Spanning Zero to 4.3 Kph on a Heavily Underactuated Bipedal Robot. ASME 2015 dynamic systems and control conference, Columbus, OH, United States, October 28–30, 2015. American Society of Mechanical Engineers Digital Collection). 10.1115/dscc2015-9899

  • CrossRef
  • Google Scholar

• 19

  RocchiA.HoffmanE. M.CaldwellD. G.TsagarakisN. G. (2015). Opensot: A Whole-Body Control Library for the Compliant Humanoid Robot Coman. In IEEE International Conference on Robotics and Automation (ICRA), Seattle, WA, United States, May 26–30, 2015. 6248–6253. 10.1109/ICRA.2015.7140076

  • CrossRef
  • Google Scholar

• 20

  RuscelliF.LaurenziA.HoffmanE. M.TsagarakisN. G. (2019). Synchronizing Virtual Constraints and Preview Controller: a Walking Pattern Generator for the Humanoid Robot Coman+. In IEEE/RSJ International Conference on Intelligent Robots and Systems, Macau, China, November 3–8, 2019. 3876–3881.

  • Google Scholar

• 21

  RuscelliF.LaurenziA.Mingo HoffmanE.TsagarakisN. G. (2018). A Fail-Safe Semi-centralized Impedance Controller: Validation on a Parallel Kinematics Ankle. In IEEE/RSJ International Conference on Intelligent Robots and Systems, Madrid, Spain, October 1–5, 2018, 1–9. 10.1109/IROS.2018.8594112

  • CrossRef
  • Google Scholar

• 22

  SardainP.BessonnetG. (2004). Zero Moment Point-Measurements from a Human Walker Wearing Robot Feet as Shoes. IEEE Trans. Syst. Man. Cybern. A.34, 638–648. 10.1109/TSMCA.2004.832833

  • CrossRef
  • Google Scholar

• 23

  SciancaN.CognettiM.De SimoneD.LanariL.OrioloG. (2016). Intrinsically Stable MPC for Humanoid Gait Generation. IEEE-RAS International Conference on Humanoid Robots, Cancun, Mexico, November 15–17, 2016. 601–606. 10.1109/HUMANOIDS.2016.7803336

  • CrossRef
  • Google Scholar

• 24

  ShimmyoS.SatoT.OhnishiK. (2013). Biped Walking Pattern Generation by Using Preview Control Based on Three-Mass Model. IEEE Trans. Ind. Electron.60, 5137–5147. 10.1109/TIE.2012.2221111

  • CrossRef
  • Google Scholar

• 25

  StephensB. J.AtkesonC. G. (2010). Push Recovery by Stepping for Humanoid Robots with Force Controlled Joints. IEEE-RAS International Conference on Humanoid Robots, Nashville, TN, United States, December 6–8, 2010. 52–59.

  • Google Scholar

• 26

  VukobratovićM.BorovacB. (2004). Zero-moment point-thirty Five Years of its Life. Int. J. humanoid robotics1, 157–173.

  • Google Scholar

• 27

  WesterveltE. R.GrizzleJ. W.ChevallereauC.ChoiJ. H.MorrisB. (2007). Feedback Control of Dynamic Bipedal Robot Locomotion, CRC Press.Vol. 26. 10.1201/9781420053739

  • CrossRef
  • Google Scholar

• 28

  WieberP. (2006). Trajectory Free Linear Model Predictive Control for Stable Walking in the Presence of strong Perturbations. IEEE-RAS International Conference on Humanoid Robots, Geneva, Italy, December 4–6, 2006. 137–142.

  • Google Scholar