**Applied Mathematics**

Vol.06 No.03(2015), Article ID:54411,9 pages

10.4236/am.2015.63043

Influence of Slope Angle on the Walking of Passive Dynamic Biped Robot

Nita H. Shah^{1}, Mahesh A. Yeolekar^{2}

^{1}Department of Mathematics, Gujarat University, Ahmedabad, India

^{2}Amiraj College of Engineering & Technology, Gujarat Technological University, Ahmedabad, India

Email: nitahshah@gmail.com, mahesh_mca07@yahoo.com

Copyright © 2015 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

Received 10 February 2015; accepted 28 February 2015; published 5 March 2015

ABSTRACT

In this paper, we modeled a simple planer passive dynamic biped robot without knee with point feet. This model has a stable, efficient and natural periodic gait which depends on the values of parameters like slope angle of inclined ramp, mass ratio and length ratio. The described model actually is an impulse differential equation. Its corresponding poincare map is discrete case. With the analysis of the bifurcation properties of poincare map, we can effectively understand some feature of impulse model. The ideas and methods to cope with this impulse model are common. But, the process of analysis is rigorous. Numerical simulations are reliable.

**Keywords:**

Biped Robot, Limit Cycle Walking, Passive Dynamic Walking, Poincare Map, Orbital Stability, Bifurcation

1. Introduction

Human walking can be considered as a complex mechanical process controlled by Newton’s law of motion which inspired to model the biped robot. Furthermore, the biped robot is more capable to work in the hazardous atmosphere than the wheeled vehicles. The fundamental problem of biped robot is to achieve the stable gait but it is not a trivial task. So, most of the researchers have focused on locomotion while working on bipedal robots till today. We will mainly focus on the walking pattern of robots. This paper will discuss about the effect of slope on the behavior of robot while walking on the inclined ramp.

This paper organized in several steps. Section 2 presents mathematical model of passive dynamic knee-less biped robot, powered only by gravity while walking on the inclined ramp. The knee-less biped robot is built of two rods where rods are considered as the legs of biped while joint of rods is considered as a hip. The gait of biped is dividing into two phase: swing phase and impact phase. The motion of the swing phase is equivalent to the motion of double inverted pendulum [1] . The model has two basic parts: one is a differential equation for the swing motion and second is the impact function for the heel-strike when foot of swing leg touching the ground. The model build here is a representative of many biped models developed in the literature [2] - [7] .

Section 3 explains the basic definitions for the solution of discrete system which developed in Section 2 [8] - [10] . It also describes the method of poincare section to the system with consideration of impact effect [11] - [14] . The discrete-time map defines from the plane of Poincare section to itself, that is, physically define from one heel-strike to next successive heel-strike which is called poincare return map [15] - [17] . The fixed point analysis of the Poincare return map will generate the periodic motion of the system. One of difficult task to generate the periodic motion is the selection of initial condition which described in this section.

Section 4 describes the bifurcation route to chaos as an influence of slope angle [18] . The steps convergent sequences and phase space diagrams simply explain bifurcation theory. Section 5 gives the conclusion and scope of future work.

2. Modeling

This section presents the dynamic model of simplest passive biped robots, consisting two equal length rigid legs pivoted at the hip with no ankles and no knees. It is a structure of three point-masses: one at the hip and other two at the centre of mass of legs. These point masses are not independent due to the distance constraints imposed by the stance and the swing legs. So, the system has four degrees of freedom. The walk of robot is considered only in the sagittal plane and on a level surface. Furthermore, the a periodic motion consists of successive swing phases, means one leg (swing leg) is free for moving and other leg (stance leg) is touching the ground, and transition phase, means the transition from one leg to another taking place in an infinitesimal length of time [3] [8] . With this consideration, the model of the biped robot made up of two parts: the differential equations which describe the dynamics of the biped robot during the swing phase and the impact model which describe the contact event and transition phase. Such models are very common in the field of biped locomotion. The linearized form of this model will be used for further analysis. It is assumed that the biped has no foot scuffing and no bounce during the walking.

2.1. Motion Equation

The angular coordinates, masses of the legs and hip, length parameters are indi-

cated in Figure 1. Note that, the positive angles are calculated anticlockwise with respect to the indicated vertical lines in Figure 1. Moreover, it walks only under the gravity and all the external forces are absents. Using the Lagrange method [3] , the passive dynamic model of the robot between successive impacts is derived in a standard second order system

(1)

Figure 1. Schematic representing the generalized coordinates of a biped robot with two point masses of mass m at the centre of mass of legs and one at hip of mass M. Two legs are of equal length .

where: and parameterize the angles of the stance leg and the swing leg respectively, parameterize slope of a walking ramp. The matrix represents the inertia matrix, the matrix

includes terms of centrifugal and coriolis forces, contains the gravity term are given as follows:

;

;

.

The effect of centrifugal and coriolis forces on the walk is very low so it can be ignored to find the walking pattern of a periodic motion of a biped robot. Neglecting the above forces, the non-linear Equation (1) formed into liearized equation,

(2)

where is the equilibrium point. The state space form of the above linerized equations model can be written as

(3)

where the state space variable,. The matrix is given as follows:

where

.

2.2. Impact Model

For the biped with equal length and without a knee, the identification of the point of contact of the swing leg with the walking surface would seem to be physically ambiguous because the swing leg must scuff along the ramp when it passes through the stance leg. It is assumed that the walker is prevented from the scuffing and the impact of the swing leg with the walking surface has no rebound and no slipping and stance leg naturally lifting from the ground without integration [2] [11] [12] . The impact between the swing leg and the ground is considered as a fully elastic collision between two rigid bodies [2] . At the time of impact, the vertical distance between foot of swing leg and the walking surface will become a zero, so it met the constraint. Considering the elastic collision between the swing leg and the ground, the angular momentums are conserved at the foot of a swing leg and a hip indicated in Figure 2.

The conservation laws of the angular momentum gives to the following compressed equation between the pre- and post-impact angular velocities after the impact [3] which denoted by and respectively. Moreover, “-” and “+” represent the left and right limits respectively.

Figure 2. Orientation of the angular momentums of the swing leg’s foot before and after the impact are displayed by the blue lines and for the hip, it displayed by the green lines.

where,; and .

After the impact the swing leg will be the new stance and the stance leg will be the new swing leg for the next subsequent step, it can be expressed by the following equation:

.

The state space form of the impact model can be written as follows:

where

(4)

2.3. The Complete Hybrid Model

The complete model of a biped passive dynamic robot is a hybrid of the motion equation and impact model which described as:

(5)

where. Figure 3 shows the walking steps of biped robot on the inclined ramp.

Regarding this model, the evolution process of periodic motion of a walking robot is described in the following steps: first, a trajectory of the robot is evaluated by the solution of motion equation until an impact occurs. The impact occurs when attains the impact condition and generates a contact point on the walking ramp. As the result of impact, there is a very fast change in the velocity components of the state vector of the contact point which is done by the impact model. This instantaneous change in the velocity components is resulting in a discontinuity of the velocities which make the biped system discontinuous. The last result of the impact model is a new initial state of a robot for the evolution of next step until the next impact. Due to the instantaneous change at the time of impact, the robot have two different state space positions at the same time which not to be obliged, so the impact event is explained with two notation: first, the state vector “” at impact time “” for the state just before impact and second, the state vector “” at impact time “” for the state just after impact.

3. Poincare Map and Analytic Solution

In this section, we define the poincare map for discontinuous system and describe the stability in sense of a fixed point of poincare map and describe the analytic solution of the impulse model (5).

Figure 3. Schematics show the steps of a periodic motion of a walking robot: Step I indicates the beginning of the walk, Step II indicates that swing leg swings until the next impact occur and Step III shows the impact of swing leg with the walking ramp.

3.1. Poincare Map

If we assume that the function is continuous where S is given the subset topology from Z, and then is an open subset of S. So The Poincare map define by where

is the solution of a continuous linear system (3) and represents impact time with respect to initial state, is well define and continuous. With this notion of Poincare map the following statements hold for the linear sys-

tem (5). (i) If is a periodic orbit of (5) that is transversal to S, then there exists a point that generates. (ii) is as table fixed point of if, and only if, the orbit of is stable in the sense of Lyapunov. (iii) is an asymptotically stable fixed point of if, and only if, the orbit of is asymptotically stable in the sense of Lyapunov [8] .

3.2. Analytic Solution of Linearized Equation of Motion and Proper Initial Condition

The analytic solution of linearized Equation (5) of motion is where _{ }is the initial condition of

walking. To get the stable walking, the initial condition should be proper otherwise random selection of initial conditions leads to walking falling forward or backward. To find the proper initial condition, considered the Poincare map of linearized model which is given by:

.

The eigenvector y corresponding to unity eigenvalue of the matrix of Poincare map for some values of such that and the corresponding slope angle can be found using the relation.

This y could be the proper choice of initial conditions for some range of slope angle [3] .

4. Influence of Slope Angle on the Periodic Orbit

The state variable y, the half inter leg angle, the step period are focused to examine the influence of slope angle on the periodic orbit. The bifurcation diagram is the graphical evolutions of such variables as the function of slope angle. We analyzed that the one periodic gaits turns periodic when slope angle is increased and further it directed to chaos. This research found that deg is the bifurcation point where the bifurcation occurs. The continuous spreading of points in the bifurcation diagrams represent the chaotic gaits, shown in Figures 4-6. We observed that the behavior of variables is monotonic during the period doubling flow. The following table details the behavior of the biped robot during a period doubling flow as the influenced of slope angle α.

Analyzing Table 1, the intervals of slope angle are gradually decreasing for periodic orbit. 8-periodic orbits are very sensitive so it difficult to observed, shown in Figure 10 and after that number of bifurcations are occurred for small change in the value of slope angle but those bifurcations are difficult to separate individually.

(a) (b)

Figure 4. Bifurcation graph: (a) angle of stance leg as a function of slope angle; (b) angular velocity of stance leg as a function of slopeangle.

(a) (b)

Figure 5. Bifurcation graph: (a) angle of swing leg as a function of slope angle; (b) angular velocity of swing leg as a function of slopeangle.

(a) (b)

Figure 6. Bifurcation graph: (a) impact time as a function of slope angle; (b) half inter leg angle as a function of slope angle.

The phase plane diagram can be used to show the period doubling flow. The 1 periodic gait occurs when the sequence of impact time converses to one point after some finite number of steps with respect to the initial condition, shown in Figure 7(a). In other words, after some finite steps from the initial conditions, the state va-

riables repeat themselves after every heel-stride i.e., so the phase plane diagram has only one

closed limit cycle which is identical to both legs, shown in Figure 7(b) [3] . For the 2-periodic gaits, the sequence of impact time has two convergent points after some finite number of steps with respect to the initial condition, shown in Figure 8(a), it means, after some finite steps, the state variables repeat themselves after two

successive heel-strike:. As a result, the phase plane diagram generates two closed limit cycles

each one is associated to one leg, shown in Figure 8(b) [3] . In case of 4-periodic gaits, the sequence of impact time has four convergent points, shown in Figure 9(a) and the state variables repeat themselves after 4 steps. Consequently, the phase plane has four closed limit cycles where each leg is associated two closed limit cycles, shown in Figure 9(b) [3] . Also in case of 8-periodic gaits, the sequence of impact time has eight convergent points Figure 10(a) and the phase space plane has eight closed limit cycles where four for each leg, shown in Figure 10(b). In general, for -periodic gaits, each leg is associated to closed limit cycles in the phase plane [3] . In case of chaotic gait, difficult to find number of convergent points for the sequence of impact time, shown in Figure 11(a), so hard to separates the closed limit cycle associate to legs in the phase plane, shown in Figure 11(b) [3] .

5. Conclusion and Future Scope

In this work, we analyzed the effect of slope angle on the stability of walking of passive dynamic biped robot on linearized model. From the results of bifurcation diagram, the stable walking range is obtained with respect to slope angle. Also we observed that the stability of walking of passive biped robot is deceased as the slope is increased. This linear model can be helpful for guessing the initial states of the robot for the stable walking of its nonlinear model. These results can be helpful for finding effect of slope angle on non-linearized model and also for designing a biped robot.

Table 1. Influence of slope angle on periodic orbit.

(a) (b)

Figure 7. 1-periodic limit cycles: (a) a diagram of impact time vs step numbers; (b) phase space diagram for deg. (---leg 1 and ---leg 2)

(a) (b)

Figure 8. 2-periodic limit cycles: (a) a diagram of impact time vs step numbers; (b) phase space diagram for deg. (---leg 1 and ---leg 2)

(a) (b)

Figure 9. 4-periodic limit cycles: (a) a diagram of impact time vs step numbers; (b) phase space diagram for deg. (---leg 1 and ---leg 2)

(a) (b)

Figure 10. 8-periodic limit cycles: (a) a diagram of impact time vs step numbers; (b) phase space diagram for deg. (---leg 1 and ---leg 2)

(a) (b)

Figure 11. Chaotic limit cycles: (a) a diagram of impact time vs step numbers; (b) phase space diagram for deg. (---leg 1 and ---leg 2)

Acknowledgements

This research is partially funded by DST-FIST2014 file number: MS1-097.

References

- Shah, N.H. and Yeolekar, M.A. (2013) Pole Placement for Controlling Double Inverted Pendulum. Mechanical and Mechanics Engineering, Global Journal of Researches in Engineering (USA), 13, 17-23.
- Wisse, M., Schwab, A.L. and van der Helm, F.C.T. (2004) Passive Dynamic Walking Model with Upper Body. Robotica, 22, 681-688. http://dx.doi.org/10.1017/S0263574704000475
- Goswami, A., Thuilotz, B. and Espiauy, B. (1998) A Study of the Passive Gait of a Compass-Like Biped Robot: Symmetry and Chaos. International Journal of Robotics Research, 17, 1282-1301. http://dx.doi.org/10.1177/027836499801701202
- Hurmuzlu, Y., Genot, F. and Brogliato, B. (2004) Modeling, Stability and Control of Biped Robots―A General Frame- work. Automatica, 40, 1647-1664. http://dx.doi.org/10.1016/j.automatica.2004.01.031
- Shah, N.H. and Yeolekar, M.A. (2014) Graphical Analysis of Passive Dynamic Biped Robot. International Journal of Robotics Research and Development, 4, 1-8.
- Shah, N.H. and Yeolekar, M.A. (2014) Local Stability Analysis of Passive Dynamic Biped Robot. International Journal of Robotics Research and Development, 4, 9-16.
- Collins, S.H., Wisse, M. and Ruina, A. (2001) A Three-Dimensional Passive-Dynamic Walking Robot with Two Legs and Knees. International Journal of Robotics Research, 20, 607-615. http://dx.doi.org/10.1177/02783640122067561
- Grizzle, J., Abba, G. and Plestan, F. (2001) Asymptotically Stable Walking for Biped Robots: Analysis via Systems with Impulse Effects. IEEE Transactions on Automatic Control, 46, 51-64. http://dx.doi.org/10.1109/9.898695
- Garcia, M., Chatterjee, A., Ruina, A. and Coleman, M.J. (1998) The Simplest Walking Model: Stability, Complexity, and Scaling. Journal of Biomechanical Engineering, 120, 281-288. http://dx.doi.org/10.1115/1.2798313
- Hsu, C.S. (1980) A Theory of Cell-to-Cell Mapping Dynamical Systems. Journal of Applied Mechanics, 47, 931-939. http://dx.doi.org/10.1115/1.3153816
- Mochon, S. and McMahon, T.A. (1980) Ballistic Walking: An Improved Model. Mathematical Biosciences, 52, 241- 260. http://dx.doi.org/10.1016/0025-5564(80)90070-X
- McGeer, T. (1990) Passive Dynamic Walking. International Journal of Robotics Research, 9, 62-82. http://dx.doi.org/10.1177/027836499000900206
- McGeer, T. (1990) Passive Walking with Knees. Proceeding of IEEE International Conference on Robotics and Automation, Cincinnati, 13-18 May 1990, 1640-1645. http://dx.doi.org/10.1109/ROBOT.1990.126245
- Khraief, N., M’sirdi, N.K. and Spong, M.W. (2003) Nearly Passive Dynamic Walking of a Biped Robot. Proceedings of European Control Conference, University of Cambridge, UK, 1-4 September 2003.
- Aroudi, A., Debbat, M. and Martinez-Salamero, L. (2007) Poincaré Maps Modeling and Local Orbital Stability Analysis of Discontinuous Piecewise Affine Periodically Driven Systems. Nonlinear Dynamics, 50, 431-445. http://dx.doi.org/10.1007/s11071-006-9190-1
- Iribe, M. and Osuka, K. (2006) A Designing Method of the Passive Dynamic Walking Robot via Analogy with the Phase Locked Loop Circuits. Proceedings of the IEEE International Conference on Robotics and Biomimetics, Kunming, 17-20 December 2006, 636-641.
- Collins, S.H., Ruina, A., Tedrake, R. and Wisse, M. (2005) Efficient Bipedal Robots Based on Passive-Dynamic Walk- ers. Science, 307, 1082-1085. http://dx.doi.org/10.1126/science.1107799
- Lin, X., Du, H. and Li, S. (2013) Parameter Influence on Passive Dynamic Walking of a Robot with Flat Feet. Kybernetika, 49, 792-808.