# Stability Analysis of Bipedal Robots Using the Concept of Lyapunov Exponents

Authors: Liu Yunping; Wang Lipeng; Mei Ping; Hu Kai
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

The dynamics and stability of passive bipedal robot have an important impact on the mass distribution, leg length, and the angle of inclination. Lyapunov’s second method is difficult to be used in highly nonlinear multibody systems, due to the lack of constructive methods for deriving Lyapunov fuction. The dynamics equation is established by Kane method, the relationship between the mass, length of leg, angle of inclination, and stability of passive bipedal robot by the largest Lyapunov exponent. And the Lyapunov exponents of continuous dynamical systems are estimated by numerical methods, which are simple and easy to be applied to the system stability simulation analysis, provide the design basis for passive bipedal robot prototype, and improve design efficiency.


---

## Body

## 1. Introduction

The walking gait of active bipedal robot is achieved by control and drive system tracking the joint angle trajectories. The huge energy constrains the development of bipedal robot [1]. The passive bipedal robot can walk naturally under the drive of gravity without the outside force. It is important to study the bipedal robot walking in low energy from the bionics [2, 3]; some research by the universities of Cornell, Mit, and Delft was published in the Science [4].However, the stability of the passive bipedal robot has high sensitivity to its structure parameters; any change of structure parameter will lead to its gait characteristics being fluctuated evidently [5]. Therefore, it is very important to optimize the configuration for the stability of the system [6, 7]. The effect of structure parameters to the stable fixed point and speed of periodic motion was analyzed under the given initial conditions in literature [5]. Asymptotical stability of passive bipedal robot in three dimensional space was researched by the method that combined extended virtual constraints and hybrid zero dynamics [8].The dynamics equation of passive bipedal robot is modeled using the Kane method in this paper, and the method of Lyapunov exponent is applied to analyze the relationship between the stability of the system and mass, length of leg, and slope angle. This method is simple and reliable for the optimization design of passive bipedal robot prototype and is easy to program.


## 2. Lyapunov Exponents

Consider the two following equations(1)xn+1=f(xn),yn+1=f(yn).Suppose that the error of initial value is|x0-y0|, through the first iteration by
(2)|x1-y1|=|f(x0)-f(y0)|=|f(x0)-f(y0)||x0-y0||x0-y0|≈|dfdx|x0||x0-y0|,
where
(3)|dfdx|x0|=limx0→y0|f(x0)-f(y0)||x0-y0|.After the second iteration, we have(4)|x2-y2|≈|dfdx|x1||x1-y1|≈|dfdx|x1dfdx|x0||x0-y0|.And through then time iteration
(5)|xn-yn|≈|∏n=0n-1df(xn,μ)dx|xn||x0-y0|.From the above, the sensitivity of two systems to initial disturbance is effected by the value of derivative|df/dx| at x0.The sensitivity of overall system mapping to initial value is achieved depending on the average of all initial conditions that needs thenth time iteration; the value of every deviation is
(6)(|∏n=0n-1dfdx|xn|)1/n.If there are small deviations from the initial value in the two systems, the result will be divergent along with the time (or times of iteration). The deviation is measured by Lyapunov exponent, the logarithmic of geometric average in following form:(7)λ=1nln(|∏n=0n-1df(xn,μ)dx|xn|),
where xn is the value of n times iteration. The computational formula of Lyapunov exponent (8) is received with n tending to infinity
(8)λ=limn→∞1n∑n=0n-1ln|df(xn,μ)dx|.The stability of system state is related to the divergence or convergence of two adjacent trajectories through the time evolution, which can be measured by Lyapunov exponent.In phase space, the initial conditions of the system are defined as an infinitesimalm dimensional ball. The ball will naturally deform as a super circle ellipsoid due to the dynamics effects. All the main axes of ellipsoid are arranged according to their length, and the Lyapunov exponent λi can be achieved by the following form:
(9)λi=limn→∞1nεi(n)ε0(n),(i=1,2,…).The Lyapunov exponent is related to the divergence or convergence of system. The trajectory is convergent in the direction of the value of Lyapunov exponent of less than 0, and the system is stable with no sensitivity to the initial conditions [9]. The trajectory is divergent in the direction of the value of positive Lyapunov exponent, and the system is unstable with sensitivity to the initial conditions. Usually, Lyapunov exponent is arranged according to λ1≥λ2≥λ3≥⋯≥λm.Where the convergence and divergence of two adjacent trajectories in phase space are quantitative description by maximum Lyapunov exponent; the motion is stable when the maximum Lyapunov exponent of system is less than zero.


## 3. Lyapunov Exponent of Continuous System

Ann-dimensional differential equation of continuous smooth dynamics system is defined as the following form:
(10)x˙=F(x),
where x˙=dx/dt, x(t)∈Rn is the state vector. Usually, suppose two close points x0, x0+u0, and the initial point x0 located in the basin of attraction; u0 is the disturbance of initial x0 (Figure 1). After a period of time t, the disturbance ut is as the following form:
(11)ut≡ft(x0+u0)-ft(x0)=Dx0ft(x0)·u0,
where ut is tangent vector in (11), which satisfies vary equation in following form [10]:
(12)Φ˙t(x0)=DxF(ft(x0))·Φt(x0),Φ0(x0)=I,
where Φt(x0) is the derivative at x0, that is, Φt(x0)=Dx0ft(x0). In order to calculate the trajectory, (13) needs to be integrated. Consider the following:
(13){x˙Φ˙}={F(x)DxF(x)·Φ},{x(t0)Φ(t0)}={x0I}.Figure 1
Distance evolution of the adjacent trajectory.Then, the average exponent of two trajectories’ divergence or convergence is defined as(14)λ(x0,u0)=limt→∞1t∥ut∥∥u0∥=limt→∞1tln∥Dx0ft(x0)·u0∥,
where ∥ut∥ is the vector of length. If λ>0, the exponent will be divergent nearby track. For x0∈M, (14) can calculate the maximum Lyapunov exponent in a very weak smoothness conditions.


## 4. Dynamics Model of Passive Bipedal Robot

The simplified biped model is showed in Figure2. The mass, length, moment of inertia are, respectively, represented with m, L, and I. The two legs are connect with a passive joint. The distance between the center of mass and hip joint is L0=L/2. The body of leg is rigid. The collision between foot and ground is completely inelastic contact with no friction and slippage. The ground is also rigid. The robot will walk along the slope face automatically under its own gravity and inertia with the leg given an initial speed.Figure 2
The simplified biped model.The process of robot’s motion is divided into parts. (1) The leg 1 will swing around the hip joint after it left off the ground. The total procedure is only in gravity acting, so the total mechanical energy is conserved. (2) When the swing leg contacts with the ground, it will exchange role with the supporting leg. The collision between the foot and the ground is instantaneous, and there is no sliding during the process of collision. The angular momentum of system is conservation.The dynamics equations are modeled by Kane method in following form:(15)[M11M12M21M22]{u˙}=f(u),f(u)=[(mL02+mL0(L0-L))gLmcos[θ]sin[q1(t)]0+mg(L-L0)cos[θ]sin[q1(t)]0+mLL0sin[q1(t)-q2(t)]u220-mgL0cos[θ]sin[q2(t)]0+(mL02+mL0(L0-L))sin[q1(t)-q2(t)]u12],M11=-I0+mL(L-L0)-mL0L+m(L-L0)2,M12=mL0Lcos[q1(t)-q2(t)],M21=(mL0(L-L0)+mL02)cos[q1(t)-q2(t)],M22=-I0-mL02.


## 5. Simulation

SupposingL=0.62 m, θ=31°, and m=14 kg, the Lyapunov exponent spectrum will be received by (14) (showed in Figure 3). The whole calculation is used by software Mathematica. From the Lyapunov exponent spectrum of Figure 3, we can find one of the values of Lyapunov exponent above 0, so the passive biped robot system is unstable. It is greatly important for optimization of the structure parameter and stability of control to analyze the dynamic characteristic of the passive biped robot.Figure 3
Lyapunov exponent spectrum.The time series of all variables is showed in Figure4. The dynamics characteristic of system can be observed by changing one of the structure parameters. From lots of results of simulation, it is concluded that small changes in mass cannot affect the dynamic characteristics of the whole system, but the length of leg and slope angle’s effect are obvious.Figure 4
Time series of the variable.


## 6. Conclusion

It is significant to optimizet mass distribution, leg length, slope angle, and other parameters of passive bipedal robot for improving the stability of the system. The dynamics equation of passive bipedal robot is modeled by Kane method in this paper, and the method of Lyapunov exponent is applied to analyze the relationship between the stability of the system and mass, length of leg, and slope angle. This method is simple and reliable for the optimization design of passive bipedal robot prototype.



---

*Source: 546520-2013-11-28.xml*