This paper focuses on the development of a model-based feedback controller to realize high versatility of fully actuated planar bipedal robotic walking. To conveniently define both symmetric and asymmetric walking patterns, we propose to use the left and the right legs for gait characterization. In addition to walking pattern tracking error, a biped's position tracking error in Cartesian space is included in the output function in order to enable high-level task planning and control such as multi-agent coordination. A feedback controller based on input–output linearization and proportional–derivative control is then synthesized to realize exponential tracking of the desired walking pattern as well as the desired global position trajectory. Sufficient stability conditions of the hybrid time-varying closed-loop system are developed based on the construction of multiple Lyapunov functions. In motion planning, a new method of walking pattern design is introduced, which decouples the planning of global motion and walking pattern. Finally, simulation results on a fully actuated planar biped show the effectiveness of the proposed walking strategy.

## Introduction

There are a variety of approaches to realize bipedal robotic walking. One of the most frequently used methods is based on the zero moment point (ZMP) balance criterion [1,2]. Enforcing the ZMP criterion turns a biped into a fixed-based robot manipulator as it requires flat-footed walking, which results in limited walking speed.

The concept of viability has been introduced to walking control [3,4], and capturability has been proposed as a computationally feasible approximation of viability [5–8]. Because capturability is defined based on the viability of walking instead of the state of balance, the associated walking strategies may be less conservative than the ZMP-based approach.

Another approach to realizing bipedal robotic walking is based on nonlinear control theories and the concept of hybrid zero dynamics (HZD) [9–13]. With the HZD framework, the walking stabilization problem becomes a stabilization problem of the closed-loop control system. Orbitally, asymptotic stabilization of an underactuated walking biped with hybrid dynamics was first achieved based on feedback control under the assumption of finite time convergence during the continuous walking phase [9]. Later on, finite time stabilization was relaxed to sufficiently fast exponential stabilization, and the concepts of virtual constraints and HZD were introduced [12–14]. Because the HZD framework only specifies a general structure of walking dynamics, it can be applied to a variety of biped models. In addition to provable closed-loop stability of the control system, another advantage of the HZD-based approach is the achievable high walking speed. Furthermore, the walking strategy based on HZD has been extended to rough terrain walking [15], three-dimensional walking [16], and neutrally stable walking [17]. Besides the HZD framework, orbital stabilization of underactuated walking has also been realized based on transverse linearization [18–20].

For fully actuated walking, the HZD framework has been extended to achieve exponential tracking of the desired walking pattern as well as the desired walking speed [11,21,22]. However, because orbital stabilization cannot control a biped's global position trajectory in Cartesian space, only velocity tracking has been addressed. Therefore, satisfactory position tracking in Cartesian space is not guaranteed, and the walking versatility is limited. Controlling a biped's global position trajectory is desirable because it enables high-level task planning and control such as multi-agent coordination and obstacle avoidance.

Also, previous studies on orbital stabilization typically utilize the traditional support-swing gait characterization that describes a bipedal gait using the support and the swing legs. This gait characterization can be used to conveniently define symmetric gaits but not asymmetric gaits. Therefore, planning and control of asymmetric bipedal robotic gaits have not been fully investigated in previous studies. However, asymmetric gait tracking is potentially meaningful for related research areas such as motion planning and control of prosthetic devices and exoskeletons because asymmetric gaits are common in injured human walking.

In this study, exponential stabilization of fully actuated planar bipedal walking is realized with enhanced versatility as compared with the previous studies on orbital stabilization. The left–right gait characterization is utilized to describe a bipedal gait so as to conveniently define both symmetric and asymmetric gaits. Although this type of gait characterization has been previously adopted [23], planning and control of asymmetric walking have not been fully studied. Another benefit of the left–right gait characterization is that the states representing joint positions become continuous and well defined upon a swing-foot landing. In control design, input–output linearization is utilized to synthesize a controller to achieve exponential tracking of the desired global motion and walking pattern, and the closed-loop stability conditions are derived based on the construction of multiple Lyapunov functions [24]. In motion planning, a new method of walking pattern design is presented, which not only guarantees that the planned motion respects the impact events but also enables decoupled planning of the desired global position trajectory and the desired walking pattern. The proposed walking strategy is validated through simulated walking of a planar biped model with three revolute joints. Simulation results show that the proposed control design can achieve exponential tracking of the desired global position trajectory in Cartesian space as well as the desired walking pattern in the configuration space. A comparison with previous studies on orbital stabilization is also presented to validate the improved versatility.

In Sec. 2, the left–right gait characterization and the full-order hybrid walking dynamics are presented along with the tracking error of the desired position trajectory in Cartesian space and the desired walking pattern. A model-based feedback control law based on input–output linearization is presented in Sec. 3. In Sec. 4, the closed-loop stability is analyzed based on the construction of multiple Lyapunov functions. A new method of walking pattern design is proposed in Sec. 5, which enables decoupled planning of global motion and walking pattern. In Sec. 6, improved versatility of the proposed walking strategy as compared with previous studies is illustrated through simulations, and the effects of control gains on the closed-loop stability are analyzed.

Preliminary results of control design and stability analysis presented in this paper were initially reported in Ref. [25]. The present paper includes important results of motion planning and fully developed theorems on stability analysis, which were omitted from Ref. [25]. Also, a comparison with previous work on versatility is presented through simulations.

## Problem Formulation

The objective of this study is to achieve exponentially stable walking with high versatility so that a fully actuated planar bipedal robot can exponentially track the desired walking pattern, both symmetric and asymmetric, as well as the desired global position trajectory in Cartesian space. In order to reach this goal, a bipedal gait is characterized by the left and the right legs instead of the traditional characterization by the support and the swing legs, the full-order hybrid walking dynamics are modeled under the left-right gait characterization, and the expressions of the walking pattern tracking error and the global position tracking error are derived.

### Gait Characterization.

The traditional gait characterization describes a bipedal gait using the support and the swing legs [9,20,21]. Under the support-swing gait characterization, the states that represent positions and velocities of the support and the swing legs are always discontinuous at the end of a step because of the role switching of the support and the swing legs. Here, we utilize the left–right gait characterization to describe a bipedal gait. There are two main advantages of using the left and the right legs for gait characterization. First, the states that represent joint positions will be continuous and well defined upon an impact [23,26,27]. Although the states that represent joint velocities may still experience sudden jumps due to landing impacts, the definition of these states is at least consistent. Second, the left–right gait characterization enables planning and control of asymmetric gaits, which will be detailed in Secs. 3 and 5.

### Hybrid Walking Dynamics With Impulse Effects.

In this study, a fully actuated planar biped with three revolute joints, identical legs, and massless thin feet is considered (see Fig. 1). There are two actuators at the hip and one at each ankle, and three of them are active at any moment except for the one at the swing ankle. Assume that the swing foot always lands flat and remains in static, full contact with the walking surface until the next landing occurs. Under this assumption, the biped is considered fully actuated during continuous phases because it has three degrees-of-freedom during continuous phases and three independent active actuators. Also, assume that the landing impact is a rigid-body contact with an impulse effect and that the double-support phase is instantaneous [9]. Without loss of generality, suppose that the swing leg length is adjustable. Otherwise, when the swing leg passes the support leg, the planar biped with a compass gait can only avoid scuffing the walking surface exactly at the vertical configuration of both legs [9].

**q**denote the joint position vector of the bipedal robot

where *q*_{1}, *q*_{2}, and *q*_{3} represent the joint angles of the left leg, the right leg, and the trunk with respect to the world coordinate frame *O _{w}X_{w}Z_{w}*, respectively.

where *u*_{1}, *u*_{2}, *u*_{3}, and *u*_{4} are defined as in Fig. 1.

Based on the left–right gait characterization introduced in Sec. 2.1, a complete walking cycle (or, a stride) consists of two successive steps, which include two single-support phases (SSP) and two swing-foot landings. The two SSPs are the left-in-support SSP and the right-in-support SSP. The two swing-foot landings are the left-to-right-support landing and the right-to-left-support landing. System dynamics during the two SSPs are continuous time, but the swing-foot landings are modeled as impulse effects.

*i*∈ {

*L*,

*R*}) as

where the subscript *i* ∈ {*L*, *R*} indicates whether the left (*L*) or the right (*R*) leg is in support, $hL(q):=l\u2009cos\u2009(q1)\u2212l\u2009cos\u2009(q2)$ is the swing foot height when the left leg is in support, and $hR(q):=l\u2009cos\u2009(q2)\u2212l\u2009cos\u2009(q1)$ is the swing foot height when the right leg is in support (*l* is the leg length).

where the derivation of $Ii:Q\u2192\mathbb{R}3\xd73$ is given in Ref. [28]. Equation (5) indicates that the states representing joint velocities may experience a sudden jump at the switching event, but the states representing joint positions are always continuous because of the left–right gait characterization.

### Tracking Error of the Desired Global Position Trajectory.

To accomplish complex tasks such as multi-agent coordination and obstacle avoidance, it is necessary for a biped to follow the desired travel path with the desired motion in Cartesian space, which can be formulated as a contouring control problem in general. By constructing an orthogonal global task coordinate frame along the desired travel path, the contour error and the motion along the desired travel path can be separately represented in two sets of coordinates, based on which contour error minimization and desired motion following along the contour can be decoupled into a stabilization problem and a trajectory tracking problem [29]. In this study, we want to solve this contouring control problem for the fully actuated planar biped model in Fig. 1, which reduces to the problem of position trajectory tracking along the *X _{w}*-axis because the

*X*-axis is the only feasible travel path for the planar biped. The complete problem of contouring control will be considered in our future work on three-dimensional walking of fully actuated bipedal robots.

_{w}*X*-axis of the world coordinate frame (see Fig. 1). Let (

_{w}*x*,

_{l}*z*) and (

_{l}*x*,

_{r}*z*) be the left foot and the right foot positions in the world coordinate frame, respectively. Let the biped's hip position

_{r}*s*along the

*X*-axis represent its global position in Cartesian space. With reference to Fig. 1,

_{w}*s*can be expressed as

Let *s _{d}*(

*t*) denote the desired trajectory of

*s*along the

*X*-axis. Then, the tracking error of the desired global position trajectory

_{w}*s*(

_{d}*t*) can be expressed as

*s*−

*s*(

_{d}*t*).

*s*(

_{d}*t*) as

Then, the tracking error of *q*_{st} − *q*_{std}(*t*) can be used to indicate the global position tracking error *s* − *s _{d}*(

*t*). If a control law is synthesized such that

*q*

_{st}−

*q*

_{std}(

*t*) is driven to zero exponentially, then the exponential tracking of the desired global position trajectory

*s*(

_{d}*t*) is realized.

### Tracking Error of the Desired Walking Pattern.

A walking pattern is defined as the relative evolution of a biped's joint positions with respect to a reference (or, encoding) variable in a complete walking cycle [9]. Tracking a preplanned walking pattern is desirable partly because the joint motion can be synchronized with respect to the reference variable. Furthermore, a proper walking pattern design can be utilized to construct the hybrid invariance of the desired motion, which will be detailed in Sec. 5.

where the functions $\varphi ji(s\xaf)\u2009(i\u2208{L,R},\u2009j\u2208{1,2})$ will be determined with a new method of walking pattern design in Sec. 5.

Hence, the tracking error of the desired walking pattern is simply $gi(s\xaf,qsw,q3)$. If $gi(s\xaf,qsw,q3)$ is exponentially driven to zero with a controller design, then exponential tracking of the desired walking pattern is realized.

*q*

_{st}is introduced to represent $\varphi ji(s\xaf)$ as

*i*∈ {

*L*,

*R*} and

*j*∈ {1, 2}. Thus, under the left–right gait characterization, the desired walking pattern of the swing-leg angle

*q*

_{sw}with respect to the support-leg angle

*q*

_{st}can be expressed as

Such a walking pattern is illustrated in Fig. 2. From the solid lines in Fig. 2, one can see that an asymmetric gait can be conveniently defined by differentiating the left and the right legs because $\varphi \u03031L(qst)$ and $\varphi \u03031R(qst)$ can be chosen independently. However, as shown by the dashed line in Fig. 2, the traditional walking characterization based on the support and the swing legs can at most represent a symmetric walking pattern where $\varphi \u03031L(qst)=\varphi \u03031R(qst)=\varphi \u03031(qst)$.

## Feedback Control Based on Input–Output Linearization

The control objective of this study is to realize exponential tracking of the desired walking pattern in the configuration space as well as the desired position trajectory in Cartesian space. As the first step of our ongoing research, it is assumed that there are no modeling errors or disturbances, and input–output linearization is utilized to synthesize the needed controller to achieve the control objectives. Specifically, the swing leg and the trunk are driven to exponentially converge to the desired walking pattern encoded by the support leg angle, and the support leg is driven to realize exponential tracking of the desired motion in Cartesian space.

Note that the output function **y*** _{i}* (

*i*∈ {

*L*,

*R*}) can be designed as three-dimensional because the biped has three independent active actuators. By exponentially driving the output function

**y**

*(*

_{i}*i*∈ {

*L*,

*R*}) to zero, exponential tracking of the desired global motion

*s*(

_{d}*t*) and the desired walking pattern $gi(s\xaf,qsw,q3)=0$ (

*i*∈ {

*L*,

*R*}) can be realized simultaneously.

From Eq. (15), it can be known that **P*** _{i}*(

*q*

_{st}) is always invertible.

where $Ni=PiMi\u22121Bi$ is proved to be invertible and $Li=PiMi\u22121hi+zi$.

where *χ* is the full state space, $03\xd73\u2208\mathbb{R}3\xd73$ is a zero matrix, and $I3\xd73\u2208\mathbb{R}3\xd73$ is an identity matrix.

**K**

*and*

_{Pi}**K**

*are chosen such that*

_{Di}**A**

*is Hurwitz, then there exists a real positive-definite-symmetric matrix*

_{i}**W**

*such that $Vi(x)=xTWix$ is a Lyapunov function candidate for the continuous-phase dynamics in Eq. (20) and there exist positive constants*

_{i}*c*

_{1}

*,*

_{i}*c*

_{2}

*, and*

_{i}*c*

_{3}

*(*

_{i}*i*∈ {

*L*,

*R*}) such that

*V*(

_{i}**x**) satisfies

**x**during continuous phases [30].

## Stability Analysis

where the expressions of $\Delta L\u2192R$ and $\Delta R\u2192L$ can be derived from $\Delta i$ and **y*** _{i}* and the expressions of the switching surfaces

*S*

_{L}_{→}

*(*

_{R}*t*,

**x**) and

*S*

_{R}_{→}

*(*

_{L}*t*,

**x**) can be obtained from

*S*and

_{i}**y**

*.*

_{i}Without loss of generality, suppose that the walking process begins with the left leg in support and the right leg in the air at *t* = *t*_{0}. Without loss of generality, suppose *t*_{0} = 0. Let $TLk$ and $TRk,\u2009k\u2208{1,2,\u2026}$ represent the actual moments of the *k*th left-to-right-support and the *k*th right-to-left-support impacts, respectively. Without loss of generality, assume $TR0=t0$. Let $\tau Lk$ and $\tau Rk,\u2009k\u2208{1,2,\u2026}$ denote the desired moments of the *k*th left-to-right-support impact assuming $x(t)=0\u2200t>TRk\u22121$ and the *k*th right-to-left-support impact assuming $x(t)=0\u2200t>TLk$, respectively.

Properties of $TLk$ and $TRk$ are summarized in Theorem 1, which is introduced based on Lemma 2 in Ref. [31].

Theorem 1. *Consider the fully actuated walking system in Eq. (24). Let the following conditions hold:*

(A1) There is no beating effect at impacts.

(A2) The desired global position trajectory

*s*(_{d}*t*) is continuously differentiable and monotonically increasing for*t*> 0.(A3) The function $\varphi \u0303ji(qst)\u2009(i\u2208{L,R},\u2009j\u2208{1,2})$ in Eq. (11) is continuously differentiable in

*q*_{st}.

Then, there exists a small neighborhood *U* of the point $(\tau ik,0),\u2009i\u2208{L,R},\u2009k\u2208{1,2,\u2026}$ such that $Tik(\tau ik,pik(\tau ik))$ is a unique continuously differentiable function in *U*, where $pLk(t)$ is the solution of $x\u02d9=ALx$ with initial condition $pLk(TRk\u22121+)=x(TRk\u22121+)$ for $t\u2208(TRk\u22121,+\u221e)$ and $pRk(t)$ is the solution of $x\u02d9=ARx$ with initial condition $pRk(TLk+)=x(TLk+)$ for $t\u2208(TLk,+\u221e)$. Also, $Tik$ has the following properties:

(P1) $TLk(\tau Lk,0)=\tau Lk,\u2009TRk(\tau Rk,0)=\tau Rk$; and

(P2) there exists a positive number

*L*such that $|Tik(\tau ,w)\u2212Tik(\tau ,u)|\u2264L\tau \Vert w\u2212u\Vert ,\u2009\u2200(\tau ,w),\u2009(\tau ,u)\u2208U$._{τ}

*Proof.* By the definitions of $Tik$ and $\tau ik\u2009(i\u2208{L,R},\u2009k\u2208{1,2,\u2026})$, the property (P1) holds. From the conditions (A1)–(A3) and Eqs. (3), (13), and (24), it can be known that the functions that define the continuous dynamics, the reset maps, and the switching surfaces are all continuously differentiable in *t* and **x**. Then, by the implicit function theorem, the property (P2) holds. *▪*

We are now ready to present the main theorem.

Theorem 2. Let the conditions (A1)–(A3) hold. Assume that $\Delta L\u2192R(\tau Lk\u2212,0)=0$ and $\Delta R\u2192L(\tau Rk\u2212,0)=0$ hold for any *k* ∈ {1, 2,…}. Then, there exist positive-definite diagonal matrices **K*** _{Pi}* and

**K**

*(*

_{Di}*i*∈ {

*L*,

*R*}) and a positive number

*δ*such that the hybrid time-varying closed-loop system in Eq. (24) is locally exponentially stable for all $x(0)\u2208B\delta (0):={x\u2208\chi :\Vert x\Vert <\delta}$.

*Proof.* Without loss of generality, suppose that the walking process begins with the left-in-support continuous phase.

Let *V _{L}*(

**x**) and

*V*(

_{R}**x**) be the Lyapunov functions associated with the left-in-support and the right-in-support phases, respectively. When the walking process begins with the left-in-support continuous phase, the

*K*th step (

*K*∈{1, 3, 5…}) is a left-in-support phase and the

*K*th switching is a left-to-right-support switching. Let $VR|K+$ and $VL|K+1+$ (

*K*∈ {1, 3, 5,…}) denote the values of Lyapunov functions right after the

*K*th and the (

*K*+ 1)th switchings, respectively. By stability analysis via multiple Lyapunov functions [24], the overall system is exponentially stable if

*V*(

_{L}**x**) and

*V*(

_{R}**x**) are exponentially decreasing in the left-in-support and the right-in-support phases, respectively, and if ${VR|1+,VR|3+,VR|5+\u2026}$ and ${VL|2+,VL|4+,VL|6+\u2026}$ are both strictly decreasing sequences.

**K**

*and*

_{Pi}**K**

*(*

_{Di}*i*∈ {

*L*,

*R*}) are chosen such that

**A**

*is Hurwitz, then the continuous-phase subsystems are exponentially stabilized. Therefore, the remaining task is to derive stability conditions to guarantee that the sequences ${VR|1+,VR|3+,VR|5+\u2026}$ and ${VL|2+,VL|4+,VL|6+\u2026}$ are both strictly decreasing. This requirement can be rewritten as*

_{i}First, we prove that there exist positive-definite diagonal matrices **K*** _{Pi}* and

**K**

*(*

_{Di}*i*∈ {

*L*,

*R*}) and a positive number

*δ*such that $VR|K+2+<VR|K+$ for any $x(0)\u2208B\delta R(0)$.

_{R}*K*th impact and

*K*+ 1)th impact.

where *T _{K}*

_{+1}is the moment of the actual (

*K*+ 1)th impact, $x|K+1\u2212$ and $x|K+1+$ represent the values of

**x**right before and after the (

*K*+ 1)th impact, respectively, and

*τ*

_{K}_{+1}is the moment of the desired (

*K*+ 1)th impact assuming

**x**(

*t*) = ∀

*t*>

*T*.

_{K}*t*and

**x**, it is locally Lipschitz continuous in

*t*and

**x**. Hence, there exists $rR*>0$ such that for any $x(0)\u2208BrR*(0)$, one has

where $L\Delta R:=max(2L\Delta Rx2,2L\Delta Rt2L\tau 2)$.

hold.

where $\Delta \tau K:=\tau K+1\u2212TK$ and $\Delta \tau K+1:=\tau K+2\u2212TK+1$. Note that Δ*τ _{K}* is the desired duration of the (

*K*+ 1)th step, which is known right after the

*K*th actual impact occurs.

*c*

_{3}

*(*

_{i}*i*∈ {

*L*,

*R*}) is determined by

**K**

*and*

_{Pi}**K**

*. Hence, if the PD gains are sufficiently large such that*

_{Di}**A**

*is Hurwitz and that there exists a positive number $\delta =min(\delta L,\delta R)$ such that*

_{i}*K*∈ {1, 3, 5,…}, then $VR|K+2+<VR|K+$ and $VL|K+3+<VL|K+1+$ hold for any $x(0)\u2208B\delta (0)$ and any

*K*∈ {1, 3, 5,…}; i.e., the closed-loop system in Eq. (24) is locally exponentially stable.

*▪*

The stability conditions in Eqs. (46) and (47) indicate that the output function should converge to zero sufficiently fast so as to diminish the possible divergence caused by reset maps. In previous studies, rapidly exponential convergence of output function has been proposed to deal with the possible expansiveness of a landing impact [21], which can also be applied here to further increase the convergence rate.

In Sec. 5, it will be shown that the assumption of $\Delta L\u2192R(\tau K\u2212,0)=0$ and $\Delta R\u2192L(\tau K+1\u2212,0)=0$ in Theorem 2 will always hold for any *K* ∈ {1, 3, 5,…} if the desired walking pattern is designed properly.

## Desired Walking Pattern Design

The desired walking motion **q*** _{id}*(

*t*,

*q*) (

_{st}*i*∈ {

*L*,

*R*}) in Eq. (13) is completely defined by the desired global position trajectory and the desired walking pattern. Suppose that the desired global position trajectory

*s*(

_{d}*t*) is determined by the high-level task planner, which is continuously differentiable and monotonically increasing for

*t*> 0. Then, the remaining task of motion planning is walking pattern design, which should guarantee that the desired motion

**q**

*(*

_{id}*t*,

*q*

_{st}) will satisfy the following conditions:

(C1) $\Delta L\u2192R(\tau K\u2212,0)=0$ and $\Delta R\u2192L(\tau K+1\u2212,0)=0,\u2009\u2200\u2009K\u2208{1,3,5,\u2026}$

(C2) forward walking direction;

(C3) ground-contact constraints including the friction cone and the unilateral constraint;

(C4) joint position and velocity limits; and

(C5) joint torque limits.

Since the last four conditions (C2)–(C5), can be easily met through trajectory optimization, they are not further discussed in this paper. The first condition (C1) essentially states that the desired gait should respect the reset map; i.e., if $x(\tau K\u2212)=0$, then $x(\tau K+)=0$ should always hold. As presented in Sec. 4, the first condition (C1) is important because it can greatly simplify the stability analysis of the hybrid time-varying closed-loop control system. However, the satisfaction of (C1) is not straightforward as it involves both the desired global position trajectory and the desired walking pattern.

In this section, a new method of walking pattern design is proposed, which guarantees that (C1) is always satisfied for any feasible *s _{d}*(

*t*) that is continuously differentiable and monotonically increasing for

*t*> 0. This is advantageous because the high-level planning of the desired global position trajectory

*s*(

_{d}*t*) and the low-level planning of the desired walking pattern represented by $\varphi \u0303ji(qst)$ are decoupled for the satisfaction of (C1).

### Hybrid Invariance of Desired Motion.

The condition (C1) can be decomposed into two parts. One part requires that the desired walking pattern should respect the reset map, which can be satisfied based on the same method of walking pattern design for constructing HZD [13]. The other part is tricky to meet, which requires that the desired position trajectory of the support leg should respect the reset map. As indicated in Eq. (7), the desired support-leg trajectory is updated at the beginning of each actual step, and thus it depends on when and where the last actual swing-foot touchdown occurs.

When the walking process begins with the left-in-support continuous phase, the *K*th step (*K* ∈ {1, 3, 5…}) is a left-in-support step starting at $t=TK\u22121+$, and the *K*th desired switching is a left-to-right-support switching at $t=\tau K\u2212$. Suppose that the desired walking pattern $gL(l\u2009sin\u2009(q1),q2,q3)=0$ has at least one intersection with the switching surface $SL(q,q\u02d9)$. Also, suppose that one of these intersections has left leg position at $q1*$. In the following, we will develop the conditions on $q1*$ and the desired walking pattern $gi(l\u2009sin\u2009(q1),q2,q3)=0$ (*i* ∈ {*L*, *R*}) such that $\Delta L\u2192R(\tau K\u2212,0)=0$ in condition (C1) is satisfied for any *K* ∈ {1, 3, 5,…}.

From Eq. (7), the desired support-leg trajectory *q*_{1}* _{d}*(

*t*) will increase monotonically within the

*K*th step if

*s*(

_{d}*t*) is planned as feasible and monotonically increasing. Then, the

*K*th desired landing moment $t=\tau K\u2212$, if exists, can be uniquely determined from the equation $q1d(\tau K\u2212)=q1*$.

where $I\u0303L(q1*):=IL(HqL(q1*))=IL(q(\tau K\u2212))$.

*K*∈ {1, 3, 5,…}. Similarly, suppose that the desired walking pattern $gR(l\u2009sin\u2009(q2),q1,q3)=0$ has at least one intersection with the switching surface $SR(q,q\u02d9)$ and that one of these intersections has right leg position at $q2*$. If $\varphi \u0303ji\u2009(i\u2208{L,R},\u2009j\u2208{1,2})$ is designed to satisfy

### Walking Pattern Parameterization With Beziér Curves.

*M*th-order Beziér curves as

where *q*_{10} and *q*_{20} are the initial support-leg angles of the left-in-support and the right-in-support phases, respectively, determined by the desired walking pattern and the switching surfaces.

The functions $\varphi \xafji(\lambda i)\u2009(i\u2208{L,R},\u2009j\u2208{1,2})$ have the following properties [12]:

- (1)
$\varphi \xaf1i(0)=\alpha i0,\u2009\varphi \xaf2i(0)=\beta i0$;

- (2)
$\varphi \xaf1i(1)=\alpha iM,\u2009\varphi \xaf2i(1)=\beta iM$;

- (3)
$d\varphi \xaf1i/d\lambda i(0)=M\alpha i1\u2212\alpha i0,\u2009d\varphi \xaf2i/d\lambda i(0)=M\beta i1\u2212\beta i0$; and

- (4)
$d\varphi \xaf1i/d\lambda i(1)=M\alpha iM\u2212\alpha i(M\u22121),\u2009d\varphi \xaf2i/d\lambda i(1)=M\beta iM\u2212\beta i(M\u22121)$.

These properties can be utilized to construct hybrid invariance of the desired motion.

### An Example of Walking Pattern Design With Third-Order Beziér Curves.

*α*and

_{ip}*β*(

_{ip}*i*∈ {

*L*,

*R*},

*p*∈ {0, 1,…,

*M*}) of $\varphi \u0303ji\u2009(i\u2208{L,R},\u2009j\u2208{1,2})$ in Eq. (60) is explained. Suppose

*M*= 3. From Eqs. (56)–(61) and the properties of $\varphi \u0303ji\u2009(i\u2208{L,R},\u2009j\u2208{1,2})$, one has

There are 16 unknown parameters of the four third-order Beziér curves, and there are ten equations and two inequality constraints in Eqs. (64)–(67). Therefore, the function $\varphi \u0303ji\u2009(i\u2208{L,R},\u2009j\u2208{1,2})$ can be determined through numerical search in order to satisfy the first condition (C1).

## Simulation Results

In this section, a fully actuated planar biped with three revolute joints is simulated to show the validity of the proposed walking strategy. We first compare our proposed walking strategy with previous studies to show that we can achieve exponential position tracking in Cartesian space but the previous study can only achieve exponential velocity tracking. Then, a symmetric walking pattern is simulated with two different desired global position trajectories *s _{d}*(

*t*)–one with a constant walking speed and the other with a time-varying walking speed—to show the versatility of the proposed walking strategy. The simulation results also show that the desired gait respects the reset map regardless of the choice of

*s*(

_{d}*t*) when the desired walking pattern is designed as introduced in Sec. 5. Furthermore, exponential tracking of an asymmetric walking pattern is illustrated through simulations. Finally, effects of the continuous-phase convergence rate on the closed-loop stability are analyzed through simulations, which validates the stability conditions in Sec. 4.

### Comparison With Previous Work on Orbital Stabilization.

In the previous work on orbital stabilization [11], the bipedal gait is characterized by the support and the swing legs. Thus, only a symmetric walking pattern can be exponentially tracked. Besides walking pattern tracking, another control objective of the previous work is velocity tracking in Cartesian space.

In order to compare our proposed walking strategy with the previous orbital stabilization [11], the desired walking pattern is chosen as the symmetric walking pattern in Fig. 4(a). The desired global position trajectory *s _{d}*(

*t*) is defined as monotonically increasing with a constant speed. Because the previous walking strategy focuses on velocity tracking in Cartesian space, its desired global motion is defined as $s\u02d9d(t)$.

The simulation results corresponding to the previous work are shown in Fig. 5. From Fig. 5, it is clear to see that exponential tracking of the desired symmetric walking pattern is achieved. However, there is always a nonzero steady-state tracking error of the desired global position trajectory *s _{d}*(

*t*), although the desired global velocity trajectory $s\u02d9d(t)$ can be exponentially tracked.

In contrast, with our proposed walking strategy, we can realize exponential walking pattern tracking and exponential global position tracking, as shown in Fig. 6.

This comparison clearly illustrates that our proposed walking strategy can greatly improve walking versatility as compared with the previous work on orbital stabilization.

### Exponential Tracking of Different Global Position Trajectories.

*s*(

_{d}*t*). Without loss of generality, the control gains are chosen the same for both cases:

**K**

*=*

_{PL}**K**

*= diag[28, 28, 28] and*

_{PR}**K**

*=*

_{DL}**K**

*= diag[11, 11, 11]. These control gains are chosen such that the matrices*

_{DR}**A**

*and*

_{L}**A**

*in Eq. (20) are both Hurwitz and that the conditions in Eqs. (46) and (47) are met under relatively large initial tracking errors. The two sets of bipedal walking share the same initial conditions*

_{R}where **W** is the solution of the Lyapunov equation **AW** + **WA**^{T} + **Q** = **0** with **Q** = **I**_{6 × 6} [30].

Figure 7 shows the results with *s _{d}*(

*t*) = 0.6

*t*− 0.1(m), and Fig. 8 corresponds to $sd(t)=2.3e\u22120.3(t+0.5)+0.6t\u22122.1(m)$. From the plots, we can see that exponential tracking of the desired hip trajectory

*s*(

_{d}*t*) with either a constant or a time-varying walking speed is achieved under the same desired walking pattern. It shows that a walking pattern generated through the proposed gait design method in Sec. 5 can be automatically incorporated with an arbitrary hip trajectory

*s*(

_{d}*t*) that is differentiable and monotonically increasing. Note that the Lyapunov function plot in Fig. 7 shows a relatively large jump at the first impact at

*t*= 0.5 s while the one in Fig. 8 shows no significant jump at the first landing impact at

*t*= 2 s. Because the desired global position trajectory in Fig. 8 has a much lower velocity in the first few seconds than that in Fig. 7, the duration of the first step is much longer in Fig. 8, which results in the much smaller trajectory tracking error right before the first impact and the much smaller jump of the Lyapunov function right after the first impact in Fig. 8.

### Stable Asymmetric Walking.

As stated earlier, the left-right gait characterization enables planning and tracking of an asymmetric walking pattern, which is illustrated with simulation results in this subsection. Figure 9 shows simulated bipedal walking with the initial conditions in Eq. (68), the desired asymmetric walking pattern in Fig. 4(b), and the desired global position trajectory *s _{d}*(

*t*) = 0.6

*t*− 0.1(m). The control parameters are chosen as

**K**

*=*

_{PL}**K**

*= diag[12, 12, 12] and*

_{PR}**K**

*=*

_{DL}**K**

*= diag[7, 7, 7], and the Lyapunov functions are defined as Eq. (69).*

_{DR}As shown in Fig. 9, the desired asymmetric walking pattern as well as the desired global position trajectory is exponentially tracked, which validates the high versatility of our proposed walking strategy.

### Effects of Continuous-Phase Convergence Rate on Closed-Loop Stability.

Theorem 2 introduced in Sec. 4 indicates that the continuous-phase convergence rate determined by the PD gains should be sufficiently fast to guarantee the stability of the closed-loop hybrid dynamical system in Eq. (24). In this subsection, two sets of simulated bipedal walking (see Figs. 10 and 11) are presented under different PD gains. They share the same initial conditions, desired walking pattern, desired hip trajectory, and Lyapunov function definitions as in Fig. 9. However, the PD gains corresponding to Fig. 10 are larger than those in Fig. 9, and the PD gains corresponding to Fig. 11 are less than those in Fig. 9. Accordingly, it can be seen from Figs. 9–11 that the actual motion converges to the desired motion faster in Fig. 10 but slower in Fig. 11 as compared with Fig. 9. Therefore, this comparison validates Theorem 2 and shows that a higher convergence rate during the continuous phases results in faster closed-loop convergence.

## Conclusion

In this paper, provably stable, fully actuated, planar bipedal robotic walking has been achieved with improved versatility as compared with previous studies. In order to define both symmetric and asymmetric walking patterns, the left and the right legs were used to characterize a bipedal gait. Under the left–right gait characterization, the full-order hybrid walking dynamics were modeled. With the output function designed as the global position tracking error and the walking pattern tracking error, an input–output linearizing controller was then synthesized to exponentially drive the output function to zero. Closed-loop stability conditions were analyzed based on the construction of multiple Lyapunov functions, which essentially requires that the continuous-phase convergence rate of the output function should be sufficiently fast in order to overcome the possible divergence caused by landing impacts. A new method of walking pattern design was proposed, which guarantees that the low-level planning of the desired walking pattern can be decoupled from the high-level planning of the desired global motion. Provable exponential stabilization and high versatility of the proposed walking strategy were validated with simulated bipedal walking.

## Acknowledgment

Any opinion, findings, and conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the National Science Foundation.

## Funding Data

Division of Information and Intelligent Systems, National Science Foundation (Grant No. IIS-0916807).