Influence of the swing ankle angle on walking stability for a passive dynamic walking robot with flat feet

Transcription

1 Special Issue Article Influence of the swing ankle angle on walking stability for a passive dynamic walking robot with flat feet Advances in Mechanical Engineering 016, Vol. 8(3) 1 13 Ó The Author(s) 016 DOI: / aime.sagepub.com Xizhe Zang, Xinyu Liu, Yixiang Liu, Sajid Iqbal and Jie Zhao Abstract To achieve high walking stability for a passive dynamic walking robot is not easy. In this article, we aim to investigate whether the walking performance for a passive dynamic walking robot can be improved by just simply changing the swing ankle angle before impact. To validate this idea, a passive bipedal walking model with two straight legs, two flat feet, a hip joint, and two ankle joints was built in this study. The walking dynamics that contains double stance phase was derived. By numerical simulation of the walking in MATLAB, we found that the walking performance can be adjusted effectively by only simply changing the swing ankle angle before impact. A bigger swing ankle angle in a reasonable range will lead to a higher walking stability and a lower initial walking speed of the next step. A bigger swing ankle angle before impact leads to a bigger amount of energy lost during impact for the quasi-passive dynamic walking robot which will influence the walking stability of the next step. Keywords Walking dynamics, passive dynamic walking robot, biped robot, numerical simulation, walking stability, step length, step duration Date received: 11 November 015; accepted: 31 December 015 Academic Editor: Mario L Ferrari Introduction Bipedal walking robot research is a popular area in robotics research due to its high resemblance to humans. During the past decades, several bipedal robots have been built. The famous bipedal walkers, ASIMO 1 and NAO, exhibit versatile gaits when walking and walk with a high stability and robustness. However, such robot features high energy consumption 3 and unnatural gait during walking because of their zero moment point (ZMP) 4 -based static walking control strategy. The concept of passive dynamic walking 5,6 provides a possible way to develop more human-like robot. However, the fully passive dynamic walking robot is very sensitive to its initial condition and can only walk successfully within a range of slopes. In order to achieve level ground walking, several quasi-passive dynamic walking robot 7 9 have been built, in which some joints are driven by motors or pneumatic muscles. The passive dynamic walker is more sensitive to its mechanical construction and mechanical parameters than the active robot because not all the joints are controlled. To improve this situation, several studies have been reported to explore more reasonable mechanical State Key Laboratory of Robotics Institute, Harbin Institute of Technology, Harbin, China Corresponding author: Xinyu Liu, State Key Laboratory of Robotics Institute, Harbin Institute of Technology, Science Park Robotics Institute 03, Yikuang Street no., Nangang District, Harbin , Heilongjiang Province, China. wo-shiliuxinyu@163.com Creative Commons CC-BY: This article is distributed under the terms of the Creative Commons Attribution 3.0 License ( which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages ( open-access-at-sage).

2 Advances in Mechanical Engineering construction and optimized mechanical parameters to obtain better walking stability and efficiency. It is found that passive robot with flat feet with springs can obtain better walking stability than that with curved feet or with the feet without ankle compliance. 10,11 It is also found that adding an upper body to the passive robot can enhance the walking stability to the robot. 1 Different types of passive models were optimized to obtain a better walking performance based on the local stability analysis or global stability analysis with cell mapping method. 16 The optimal mechanical structure of the passive robot in Collins et al. 17 is obtained based on several times trail, error and correction process in experiments. In addition, because only few joints are actuated, the common effective control idea, such as the hybrid zero dynamics 18 that is well used on the under-actuated dynamic biped robot such as RABBIT, 19 is hard to apply to the passive dynamic robot. So, usually the control for passive dynamic walkers is applied to prevent falling forward or backward in a catch-up way to keep dynamic equilibrium during walking but not to define the exact joint trajectories. It is found that the limit cycle walker will never fall forward if the swing leg swings forward fast enough and will not fall backward if the stance leg angle is not too big, 0 which has been proved effective to guide control method design in some other works. 1, MW Spong et al. 3 used virtual gravity principle to control a compass gait robot walking on level ground in simulation. The energy consumption is very low, but it is hard to apply it on a real robot. In Schuitema et al., 9 a reinforcement learningbased controller is applied to a five-link passive robot which can deal with 1 cm step-down during walking. Hobbelen and Wisse 1 applied local ankle feedback control and ankle push-off to control a limit cycle walker which can walk over a 3 cm step-down successfully. It is also found that the posture of the robot has a big effect on the energy lost at impact, 4 which could affect the walking performance definitely. The swing leg retraction method to improve the waking stability for a passive dynamic walker proposed in Hobbelen and Wisse 5 is based on this idea. So, this study aims to investigate whether we can improve the walking performance for a passive dynamic walking robot by just simply changing the swing ankle angle so that to adjust the robot posture before impact. In this study, a fully passive dynamic walking model that consisted of two straight legs, two flat feet, two ankle joints, and a hip joint was built. A detailed walking dynamics that contains double stance phase was derived based on the virtual power principle. 14,6 Numerical simulations were performed in MATLAB. The influence of the swing ankle angle on walking stability, step length and step duration was analyzed. The effect of different values of the swing ankle angle on walking stability of the quasi-passive dynamic walker was also analyzed by observing the amount of lost energy during impact. The article is organized as follows. In section Walking dynamics, we derive the walking dynamics. In section Numerical simulation and results, we introduce the numerical simulation process and the simulation results. In section Effect of swing ankle angle s value on walking performance, detailed analysis based on the simulation results is discussed. We conclude in section Conclusion. Walking dynamics Walking phases To study the influence of the swing ankle angle on walking performance, a passive walking model that consisted of two straight legs, two flat feet, two ankle joints, and a hip joint was built, as shown in Figure 1, where u 1, u, u 3, and u 4 denote the trailing leg, leading leg, trailing foot, and leading foot angle, respectively. C 1, C, C f1, and C f denote the position of the center of mass (COM) for the trailing leg, leading leg, trailing foot, and leading foot, respectively. K h and k a denote the joint spring stiffness at the hip and ankle joints, respectively. We only focus on the influence of the swing ankle just before impact, so straight leg without knee joints is enough for our analysis. Double stance phase (i.e. both feet contact with the ground at the same time) exists in human and animal walking. 7 So there is double stance phase during walking in this study. One typical walking phase contains first impact phase, first double stance phase, second impact phase, second double stance phase, and swing phase, as depicted in Figure. The red cycles in the feet represent the constraints exerted by the ground contact forces in Figure. Phase a: First impact phase. The swing foot s heel impacts with the ground and the old swing leg becomes the new stance (i.e. leading) leg (dotted line in Figure ), and the old stance leg becomes the new trailing leg (solid line in Figure 1. Passive dynamic walking model. θ4 y x ka C θ Cf β kh θ1 H C1 ka Cf1 θ3

3 Zang et al. 3 Table 1. Mechanical parameters and their default values. Parameter Description Dimensionless value Figure. Walking phases in one walking step. Figure ). Typically, the instant just after impact is chosen as the initial condition of the walking. However in this article, to decrease the number of independent initial state variables, the instant just before heel impact is chosen as the start of the walking. Phase b: First double stance phase. The trailing leg s toe and the leading leg s heel keep contact with the ground at the same time. Whether the trailing leg s heel keeps contact with the ground depends on where the ground reaction force under the trailing leg happens. In this article, we can find that the trailing leg s heel loses contact immediately after the first impact, which will be discussed in detail in the following section. So the robot rotates forward around the trailing foot s toe and the leading foot s heel. Phase c: Second impact phase. When the leading foot rotates to a proper angle, the second impact happens. The leading toe impacts with the ground and then the leading foot keeps full contact with the ground. Phase d: Second double stance phase. The trailing foot s toe and the whole leading foot keep contact with the ground at the same time. The robot rotates forward around the trailing foot s toe and the leading foot s ankle joint. Phase e: Swing phase. When the trailing foot s toe loses contact with the ground, the walking comes to the swing phase. The leading leg rotates forward around the ankle joint and the swing leg rotates forward around the hip joint. The leading foot is kept at a constant angle with respect to the leading leg in this phase by the constraint in equation (0). This phase ends before the instant that the swing heel impacts with the ground again. Then, the robot will repeat these phases to keep walking. Walking dynamics The parameters and their default values are listed in Table 1. The values are determined based on the l Leg length 1 l f Foot length 0.5 l f1 Frontal foot length 0. (from toe to ankle) l f Hinder-foot length 0.05 (from ankle to heel) b l Leg COM 0.5 (from COM position to hip) b f Foot COM 0.1 (from COM position to heel) m h Hip mass 0.4 m f Foot mass 0.05 m l Leg mass 0.5 J h Hip moment of inertia 0 J f Foot moment of inertia J l Leg moment of inertia b Slope angle 3 ( ) K h Hip joint spring stiffness 0.4 k a Ankle joint spring stiffness 0.6 COM: center of mass. analogy to human body parameters. The values for k h and k a are determined by several times tuning in simulation to find the fixed point. 8 To make the simulation results more universal, all the parameters are transferred to their dimensionless form. The parameters related to length is obtained by dividing l, the parameters related to mass is obtained by dividing M which denotes the total mass of the robot, the parameters related to moment of inertial is obtained by dividing Ml, the time is obtained by dividing (l/g) 1/ and the spring stiffness is obtained by dividing Mg/l. Usually, we use the Poincare mapping method 8 to describe the periodic walking for a passive robot. We can choose a Poincare section through which the walking trajectory passes periodically at any instant of the Poincare mapping process, and the state on this section is chosen as the initial condition for the walking. If the robot can converge to one fixed initial condition after several steps, we deem that the walking is locally stable, and the converged initial condition is the so-called fixed point. 8 Usually, the instant just after foot ground impact is chosen as the Poincare section. However, in this study, the instant just before foot ground impact is chosen as the Poincare section (i.e. initial condition). If we choose the instant just before impact as the initial condition, the number of initial independent state variables can be reduced to three which are u 1, _u 1, and _u, and all other state variables can be described by the three initial independent state variables, as shown in equations (1) and (). Fewer initial independent state variables enhance the chance of obtaining fixed point greatly. As

4 4 Advances in Mechanical Engineering mentioned above, the trailing foot is kept at a constant angle with respect to the trailing leg before impact, which means that _u 3 = _u 1. Like human-beings, we assume that the stance foot keeps full contact with the ground before the swing heel touches the ground as in equation (1). In equation (), u w is defined as the swing ankle angle that denotes the angle between the swing leg and the foot and the superscript denotes the instant just before impact. In this article, the y-axis is set as the reference axis for u 1 and x-axis is that for u 3, as shown in Figure 1. So the value of u w is 0 when the actual physical angle between the swing leg and the foot is p=. By changing the value of u w, we can obtain different angle values for u 3 as in equation (). The value of u is obtained by the physical constraint of the robot just before impact as in equation () _u 3 = _u 1 _u 4 = u 4 = 0 ð1þ 8 u 3 = u 1 + u w p >< 180 u = arccos l f 1 sin u 3 + l cos u ðþ 1 >: l First impact phase. We assume that the foot ground impact is instantaneous and fully inelastic, and the angles of each part cannot change suddenly at the instant just after impact. According to the angular momentum theorem, the angular momentum of the entire robot is conserved around the leading heel F; the angular momentum of the leading leg, trailing leg, trailing foot, and hip is conserved around the leading ankle A1; the angular momentum of the trailing leg and trailing foot is conserved around the hip joint H; and the angular momentum of the trailing foot is conserved around the trailing leg ankle A, as in equation (3) 8 _x + htt = _x+ hlh _y + htt = _y+ hlh >< ~L + F = ~L F +(x F3 r x F ) Iy (y F3 y F ) Ix ~L + A1 = ~L A1 +(x F3 r x A1 ) Iy (y F3 y A1 ) Ix ~L + H = ~L H +(x F3 r x H ) Iy (y F3 y H ) Ix >: ~L + A = ~L A +(x F3 r x A ) Iy (y F3 y A ) Ix ð3þ In equation (3), ~L + denotes the angular moment about a point just after impact and ~L denotes that just before impact, respectively. The left side of the first two equations is calculated based on the trailing leg s toe ground contact position, and the right side is calculated based on the leading leg s heel ground contact position, respectively. These two equations can guarantee the constraint that the trailing leg s toe and the leading leg s heel must keep contact with the ground without sliding and lifting during impact. During impact, there are impulses both under the leading and trailing foot as the trailing leg will not lift up immediately. Ix and Iy denote the component of impulse along x-axis and y-axis under trailing foot, respectively in equation (3). Because the position where the impulse generates is unknown, a new parameter r that denotes the distance from the trailing foot s toe to the position that the impulse generates is introduced here. To solve r, we assume that the trailing foot s heel loses contact with the ground immediately after impact (i.e. the impulse generates at the toe point) first, and so r = 0. There are six unknown parameters that are _u + 1, _u +, _u + 3, _u + 4, I x, and I y, respectively, and can be solved by the six equations in equation (3). After we solve this equation, we examine the value of Iy. If the value of Iy is greater than 0, the trailing toe keeps contact with the ground indeed and thus we deem that the assumption about the value of r is reasonable. If Iy 0, we deem that r 6¼ 0, and thus the trailing foot keeps full contact with the ground during impact (i.e. _u + 4 = 0). We can also solve all the unknown parameters in equation (3) with the additional function _u + 4 = 0. First double stance phase. During this phase, the degree of freedom (DOF) of the system is, but the number of the generalized coordinates is 4. So the system is a holonomic system that contains redundant coordinates, and thus the Euler Lagrange equation 9 cannot be applied here. The virtual power principle 14,6 is suitable in this phase, which can handle holonomic system that contains redundant coordinates or non-holonomic system. The walking dynamics can be described by the rectangular coordinates x or the generalized coordinates q as follows x = ½x h, y h, x c1, y c1, x c, y c, x cf 1, y cf 1, x cf, y cf Š T q = ½x h, y h, u 1, u, u 3, u 4 Š T ð4þ In equation (4), the parameters in x denote the position of the COM of the hip, trailing leg, leading leg, trailing foot, and leading foot along x-axis and y-axis, respectively. For q, please see Figure 1. We add two extra redundant generalized coordinate x h and y h to the system, which will be discussed below. The rectangular coordinates can be converted to the generalized coordinates as follows x = x(q) ð5þ Convert equation (5) into its variation form and second differential form as well, and we can get equation (6)

5 Zang et al. 5 d_x = J d_q x = J q + J9 _q ð6þ where J = x= q, J9 = (J _q)= q. The virtual power approach follows the following form 14 as follows d_x T ½f M xš = 0 ð7þ Equation (7) means that the sum of the virtual power of the active forces and the inertial forces is zero at any instant for an ideal mass points system, so we convert it to the following form as DP + DP = 0 ð8þ Equation (8) is a general form for the virtual power principle, where DP denotes the active force s virtual power and DP denotes the inertial force s virtual power, respectively. The detailed forms of them are as follows 8 < DP = P f i d_x i + P t i d_q i i i DP = P m i x i d_x i P ð9þ : I i q i d_q i i i In equation (9), x i and q i denote each element in equation (4), respectively. Based on equation (6), we can convert equation (9) to its matrix form in generalized coordinates as follows DP = d_q T J T f + d_q T t DP = d_q T f½j T MJ + IŠ q + ½J T ð10þ MJ9Š_qg In equation (10), M denotes the mass matrix for each mass point (i.e. each COM of the robot), I denotes the moment of inertia matrix for each part of the robot, f denotes the force matrix which can be induced by gravity or other external forces, and t denotes the torque matrix to each generalized coordinate which can be induced by passive springs or active driving sources. In this article, all the torques exerted to each joints are induced by passive springs. The details of them are listed as follows M = Diag(m h, m h, m l, m l, m l, m l, m f, m f, m f, m f ) I = Diag(J h, J h, J l, J l, J f, J f ) f = Mg½sin (b), cos (b), sin (b), cos (b), sin (b), cos (b), sin (b), cos (b), sin (b), cos (b)š T t = ½t hx, t hy, t sw, t st, t f 1, t f Š T ð11þ The descriptions and values for all elements in equation (11) are listed in Table 1. Based on equation (10), we can transfer equation (8) to the following form d_q T fj T ½f MJ9 _qš + t ½J T MJ + IŠ qg = 0 ð1þ Because the generalized coordinates are not all independent variables (i.e. the number of generalized coordinates is bigger than that of the DOF), we introduce l multiplier and the relative constraint function into the system, based on the Lagrange equation of first kind. 15,6 The instant when the trailing toe loses contact with the ground depends on the ground reaction force exerted on the toe and the direction of the acceleration of the toe ground contact point. 30 So in order to solve the ground reaction force, at least two extra generalized coordinates x h and y h should be added to the system. During this phase, the introduced constraint function mentioned above can be described as follows 3 x h + l sin (u 1 )+l f 1 : cos (u 3 ) l f y e(q)= h l cos (u 1 )+l f 1 : sin (u 3 ) x h + l sin (u ) l f cos (u 4 ) L 5 = 0 ð13þ y h l cos (u ) l f sin (u 4 ) In equation (13), L denotes the step length and can be handled as a constant value here. The first two constraint functions guarantee that the trailing leg s toe must keep contact with the ground (i.e. the first two functions are the calculated position of the trailing toe along x-axis and y-axis, respectively). The second two constraint functions guarantee that the leading leg s heel must keep contact with the ground (i.e. the second two functions are the calculated position of the leading heel along x-axis and y-axis, respectively). Transfer equation (13) into its variation form and multiply by the l multiplier as follows d_q T T T l = 0 ð14þ In equation (14), T = e(q)= q and l is the multiplier that denotes the constraint force induced by each constraint functions in equation (13). 15,6 In this phase, the first two elements of l denote the ground contact force under trailing toe along x-axis and y-axis, respectively, and the last two elements of l denote the ground contact force under leading heel along x-axis and y-axis, respectively. On rearranging equations (1) and (14), we can get ½J T MJ + IŠ q = J T ½f MJ9 _qš + t + T T l ð15þ Equation (13) can be converted into its second differential form T q + T _q _q = 0 ð16þ q On rearranging equations (15) and (16), we can obtain the walking dynamics in this phase as follows

6 6 Advances in Mechanical Engineering J T MJ + I T T q = J T ½f MJ9 _qš + t T 0 l ( (T _q)= q) _q ð17þ Second impact phase. When the value of the leading toe along y-axis is equal to 0, the walking comes to the second impact phase. The assumption for this phase is the same as it is in the first impact phase. The angular momentum of the entire robot is conserved around the tailing toe F3; the angular momentum of the leading leg, leading foot, trailing leg, and hip is conserved around the trailing ankle joint A; the angular momentum of the leading leg and leading foot is conserved around the hip joint H; the angular momentum of the leading foot is conserved around the leading leg ankle joint A1, as in equation (18) 8 x + htt = x+ hla y + htt = y+ hla _u >< + 4 = 0 ~L + F3 = ~L F3 +(x F1 r x F3 ) Iy (y F1 y F3 ) Ix ~L + A = ~L A +(x F1 r x A ) Iy (y F1 y A ) Ix ~L + H = ~L H +(x F1 r x H ) Iy (y F1 y H ) Ix >: ~L + A1 = ~L A1 +(x F1 r x A1 ) Iy (y F1 y A1 ) Ix ð18þ In equation (18), the first two equations guarantee the trailing toe and leading ankle keep contact with the ground in this phase as mentioned in the first impact phase. r* denotes the distance from the leading toe to the position where the impulse happens. We assume that the leading foot keeps full contact with the ground as mentioned above, which means that _u + 4 = 0. There are seven unknown parameters _u + 1, _u +, _u + 3, _u + 4, Ix, I y, and r, which can be solved by equation (18) directly. Second double stance phase. The walking dynamics for this phase shares the same form as in the first double stance phase, as shown in equation (17). However, the constraint function for this phase is different from that in the first double stance phase. During this phase, the trailing toe still keeps contact with the ground, and the leading foot keeps full contact with the ground. So the constraint function in this phase is as follows 3 x h + l sin (u 1 )+l f 1 : cos (u 3 ) l f y h l cos (u 1 )+l f 1 : sin (u 3 ) e(q)= x h + l sin (u ) l f cos (u 4 ) L y h l cos (u ) l f sin (u 4 ) 5 = 0 ð19þ y h l cos (u )+l f 1 sin (u 4 ) The first two constraint functions guarantee that the trailing leg s toe must keep contact with the ground. The second two constraint functions guarantee that the leading leg s heel must keep contact with the ground. The last function guarantees that the leading leg s toe must keep contact with the ground (i.e. the last function is the calculated position of the leading toe along y- axis). The leading toe shares the same constraint function with the leading heel along x-axis as they are part of the same rigid leading foot so that three constraint functions are enough for the leading foot. When the value of l multiplier under the second constraint function equals to 0, the walking robot comes to the swing phase. In other words, when the ground contact force under the trailing toe reduces to 0, the trailing foot loses contact with ground. Swing phase. The walking dynamics for this phase shares the same form as in the first double stance phase, as shown in equation (17). However, the constraint function for this phase is different from that in the first double stance phase. During this phase, the leading foot keeps full contact with the ground. So the constraint function in this phase is as follows 3 x h + l sin (u ) l f cos (u 4 ) L y e(q)= h l cos (u ) l f sin (u 4 ) y h l cos (u )+l f 1 sin (u 4 ) 5 = 0 ð0þ u 3 u 1 u w p=180 In equation (0), the last constraint function controls the angle between the swing leg and swing foot as mentioned in equation (1). When u w = 0, this constraint function means that the swing foot keeps the same angle with the swing leg all over the swing phase. When the value of the trailing heel along y-axis is equal to 0, the walking comes to the first impact phase again and repeat the walking periodically. Numerical simulation and results Numerical simulation process We perform the numerical simulation process in MATLAB based on the derived walking dynamics in section Walking dynamics. We can test whether the derived walking dynamics is correct by the simulation process and make further analysis based on the simulation results. To find out the fixed point, the Poincare mapping 8 function should be built first. The Poincare mapping function can be described in the following form as x n + 1 = P(x n ) ð1þ In equation (1), P is the Poincare function which describes the mapping from (nth) system initial state to the ((n + 1)th) system initial state. After several times of mapping, if the system state can converge to itself,

7 Zang et al. 7 (17). The function in the swing phase is solved by the ODE45 algorithm in MATLAB as well. If the result of this function satisfies the condition 3, the robot comes to the first impact phase again. By now, one loop for the Poincare mapping function is ended. If not, we return to the start of the walking and give another initial condition for the walking and perform the same Figure 3. Details of Poincare mapping function. loop. Condition 3 can be described as follows 8 < u \u 1 vya 1 = l f 1 cos (u 4 ) _u 4 l sin (u ) _u + l sin (u 1 ) _u 1 \0 ð5þ : yf = l f 1 sin (u 4 )+l cos (u ) l cos (u 1 ) l f sin (u 3 )=0 we can say that the fixed point 8 is obtained, as in equation () x = P(x ) ðþ In other words, the Poincare mapping method can describe the walking process for one step. The input of this function is the initial condition for this step, and the output of this function is the ending of this step. The ending of this step can be used as the initial condition for the next step which forms a loop for the walking. The details of Poincare mapping function are described in Figure 3. As shown in Figure 3, we first use the first impact function as in equation (3). The result of this function is used as the initial condition for the first double stance function as in equation (17). This function is solved by the ODE45 algorithm in MATLAB. If the result of this function satisfies the condition 1, the robot comes to the second impact phase. If not, we return to the start of the walking and give another initial condition for the walking and perform the same loop. Condition 1 can be described as follows yf 1 = l f sin (u 4 )=0 ð3þ This function means that when the position of the leading toe along y-axis is equal to 0, the robot comes to the second impact phase. After this phase, the robot comes to the second double stance phase. The function in the second double stance phase is solved by the ODE45 algorithm in MATLAB as well. If the result of this function satisfies the condition, the robot comes to the swing phase. If not, we return to the start of the walking and give another initial condition for the walking and perform the same loop. Condition can be described as follows l F3 = 0 ð4þ This function means that when the ground reaction force under the trailing leg s toe is equal to 0, the robot comes to the swing phase. It can be solved by equation The first equation guarantees that the swing leg must be at the front of the stance leg before impact. The second equation guarantees that the leading toe is in the state of falling-down, but not lifting-up. Taking the leading ankle to calculate the velocity gives a more precise prediction about whether the swing leg is swinging upward or downward than that of the leading heel. The last equation guarantees that the leading heel is going to impact with the ground immediately. We use the Newton Raphson iteration algorithm 14 to obtain the fixed point based on the Poincare mapping function described above. After the fixed point has been obtained, we use it as the initial condition for the walking and perform the simulations based on the Poincare mapping function. Simulation results First, we set u w = 0 in equation (1) to perform the simulation process and find out the fixed point. After trying several times, we obtained the fixed point at x = [ ; ; ]. We took this as the initial condition of the walking and performed five steps simulation. As mentioned above, by taking the instant just before first impact as the start of the walking, the number of the independent initial state variables can be reduced to three, so there are only three elements for the fixed point. The walking gait analysis based on the simulation results is shown in Figure 4. Figure 4 shows the values of the legs and feet angle for two steps. In Figure 4(a), the leg angle seems to be nearly the same as the passive walking without a double stance phase. But, in Figure 4(b), we can find obvious evidence for the existence of the double stance phase. The second step starts from the point A as shown in Figure 4(b). After first impact phase, it comes to point B. Then, the first double stance phase starts and the leading toe impacts with the ground at point C. This phase spends only about 0.15 dimensionless time. Then, the second impact and double stance phases start. At point D, the trailing toe loses contact with the ground and the swing

8 8 Advances in Mechanical Engineering Figure 5. (a) Angular velocity of the swing and stance legs and (b) Angular velocity of the swing and stance feet. Figure 4. (a) Swing and stance leg angle values and (b) Swing and stance feet angle values. phase starts. After the instant of losing contact, the swing foot is set to have the same angle with the swing leg immediately. So there is a sudden change in the value of trailing foot s angle at point D. The leading foot s angle is equal to 0 after the second impact phase until the next first impact phase starts again. If there was no double stance phase, the swing foot will have totally the same curve with the swing leg. So just from this figure, we can conclude that there must be double stance phase during the walking. Figure 5 shows the values of the angular velocity for the legs and feet for two steps. In Figure 5(a), we found that the second impact phase has only tiny influence on the angular velocity of the legs. In Figure 5(b), we found that the leading foot s angular velocity changed greatly after the first impact, because the ground reaction impulse is very big but the mass of the foot is very small. The angular velocity of the trailing foot changed three times in one step. The first two times are caused by the impact. The last change is because that the trailing foot is constrained to swing together with the swing leg immediately after it has lost contact with the ground. Figure 6 shows the limit cycles for the legs and the feet for five steps. The two legs of the walking model show two limit cycles during walking which implies that the model can converge to its initial condition for every step. From points A to E and points A# to E#, the model walks stably from the first impact phase to the swing phase. The two feet of the walking model show two limit cycles as well, as shown in Figure 6(b). From points A to E and points A# to E#, the model walks stably from the first impact phase to the swing phase. For the stance (i.e. leading) foot, its angle and angular velocity value is kept at 0 after the second impact phase, so A#, B#, and E# share the same point. The stick diagram of the walking process is depicted in Figure 7, in which the robot shows a stable walking with periodical gait. The main features for a walking model with flat feet that contains double stance phase in this study are (1) the double stance phase accounts for 5.5% of one step time and () there are two times of impact during walking and the second impact has tiny influence on the walking dynamics. Effect of swing ankle angle s value on walking performance Fully passive dynamic biped walking robot In this section, we introduce the influence of swing ankle angle on walking performance for the fully

9 Zang et al. 9 Figure 6. (a) Limit cycles of the swing and stance legs and (b) Limit cycles of the swing and stance feet. Figure 7. Stick diagram of the passive walking model. passive walking model in this study. The walking performance refers to that of walking stability, step length, and step duration. We try to change the swing ankle angle before first impact based on the following equation u 3 = u 1 + u w p ð6þ 180 The value of the ankle joint angle u w is set immediately to a value and kept constant during the swing phase by the last constraint function listed in equation (0) at the start of the swing phase in simulation, which means that the control of this joint is finished immediately. However, in practice, the change in the ankle joint angle should be realized by an active controlled spring that is connected to a motor. When u w = 0, the ankle is in its equilibrium position. So we change the value for u w in a reasonable range and try to obtain the fixed point for different ankle angle first. We also consider that the results can be varied greatly with the change of the hinder-foot length, so we change the hinder-foot length from 0.05 to 0.1 to investigate this effect. Figure 8 shows the difference of the two limit cycles for the swing ankle angle at 15 and 35, respectively. From points A to E and points A# to E#, the model Figure 8. Limit cycles under different swing ankle angles. Figure 9. Walking stability under different swing ankle angles.

10 10 Advances in Mechanical Engineering Figure 10. Step length under different swing ankle angles. Figure 11. Step duration under different swing ankle angles. walks successfully from the first impact phase to the swing phase. Figure 9 shows the influence of different ankle angle values on walking stability under different hinder-foot lengths. In Figure 8, we found that the limit cycles under different swing ankle angle values are quite different which implies that the change in ankle angle value has an obvious influence on the walking performance. When u w = 15, we can hardly find the first double stance phase which implies that the swing foot is nearly parallel to the walking surface at the instant just before first impact. We cannot find the fixed point when the swing ankle angle is smaller than 17.7 or bigger than 38 when the hinder-foot length is equal to 0.05, and we cannot find the fixed point when the swing ankle angle is smaller than 17 or bigger than 10 when the hinderfoot length is equal to 0.1. Figure 9 shows the biggest eigenvalues of the Jacobian matrix for the Poincare mapping function under different ankle angle values which can be used as a criterion to determine whether the passive robot is locally stable. When the biggest eigenvalue is smaller than 1, we deem that the robot is locally stable. Analytical solution for the Jacobian matrix is extremely hard to obtain, so we obtain the Jacobian matrix based on the numerical method as described in Kim et al. 31 We take the module of the biggest eigenvalue if it is a plural. From the figure, we can see that the robot is getting more stable when the ankle angle value is changed from its equilibrium position in most area within the changing range. When the hinder-foot length is equal to 0.1, this influence is more obvious. In general, the different ankle angle values can affect the walking stability greatly which means that we can obtain a better walking stability by just changing the swing ankle angle before impact. Figure 10 shows the effect of different swing ankle angle values on step length under different hinder-foot lengths. Figure 11 shows the effect of different swing ankle angle values on step duration under different hinder-foot lengths. The step length is calculated as follows L = l f + l sin (u + 1 ) l f cos (u + 3 ) l sin (u + ) l f 1 cos (u + 4 ) ð7þ In equation (7), u + means the instant just after first impact. From Figure 10, we can see that the step length changes slowly when the swing ankle angle is smaller than 0 and changes rapidly when it is bigger than 0. The step length can be changed from to 0.486, which is about 5% of its initial value. So we can change the step length by only changing the swing ankle angle before impact. When the hinder-foot length is equal to 0.1, this change happens rapidly. In Figure 11, we can see that the step duration changes slowly when the swing ankle angle is smaller than 5 and changes rapidly when it is bigger than 5. The step duration can be changed from.1864 to , which is about 66.8% of its initial value. When the hinder-foot length is equal to 0.1, this change is more obviously, which is about 94.4% of its initial value. So the changing of swing ankle angle has a great influence on the step duration. In other words, the changing of swing ankle angle value has a great influence on the walking speed. In general, we can change the walking stability, step length, and step duration by just changing the swing ankle angle before first impact, especially when we want a smaller step length and a slower walking speed. It should be mention that the changing of the walking performance is mainly caused by the change in of the fixed point which is caused by the change in of the swing ankle angle (i.e. changing of the swing ankle angle is the direct factor but changing of the fixed point is the primary factor).

11 Zang et al. 11 Quasi-passive dynamic walking robot For a quasi-passive dynamic walking robot prototype, it is extremely hard to obtain the fixed point in the experiments as the passive robot is usually initialed manually. So it is hard to compute the biggest eigenvalues of the Jacobian matrix in experiments. In practice, the step length or step duration (i.e. walking velocity) is usually used as the criterion to determine whether the robot is walking stably, like the gait sensitivity norm 3 criterion. If the robot has nearly the same step length or step duration for every step, it is considered that the walking is stable. Humans try to keep the walking balance by changing the step length or step duration, 33 and usually a faster walking speed needs a bigger step length and vice versa. 33 So in this section, we try to investigate whether we can adjust the step duration of the next step by only changing the swing ankle angle before impact. We investigate this by observing the amount of total kinetic energy of the robot before impact and the amount of lost kinetic energy during impact. More remained kinetic energy after impact implies a higher initial walking speed (i.e. less walking duration) of the next step and vice versa. The kinetic energy can be calculated based on the following equation E = m l(v c1 + v c ) + J l(u 1 + u ) + m h v h + J f (u 3 + u 4 ) + m f (v cf 1 + v cf ) ð8þ We also use the passive walking model here to make the analysis, but we give the same initial condition for one walking step with different swing ankle angle values. This can guarantee that the change in the walking performance is only caused by the change in the ankle angle. Figure 1(a) shows the total kinetic energy change under different swing ankle angle values when the slope angle value is equal to 3. In Figure 1(a), the top block denotes the kinetic energy before first impact and the bottom block denotes the remained kinetic energy after impact. The vertical line denotes the amount of the lost kinetic energy during impact, which means that the longer the line the more kinetic energy lost in the impact. From the figure, we found that when the swing ankle angle is getting bigger before impact, the amount of total kinetic energy before impact is getting smaller. When the swing ankle angle is getting bigger before impact, the amount of lost kinetic energy during impact is getting bigger. So the amount of remained kinetic energy for the next step is getting smaller with the increase in the swing ankle angle before impact. The amount of remained kinetic energy after impact can be reduced from to 0.051, about 5.8% of its Figure 1. (a) Changing of the kinetic energy of the robot when the slope angle is equal to 3 (degree) and (b) Changing of the kinetic energy of the robot when the slope angle is equal to 5 (degree). initial value when the swing ankle angle changes from 0 to 35, and can be enhanced from to , about 9.9% of its initial value when the swing ankle angle changes from 0 to 15. There are two main reasons for this effect: (1) the increase in the swing ankle angle value makes the swing heel impact with the ground at an earlier time, which decrease the angular velocity of each part of the robot and thus the amount of kinetic energy and () the increase in the swing ankle angle value changes the impact posture of the robot, which influence the amount of lost kinetic energy during impact. So we can obtain a more proper initial walking speed of the next step by just simply changing the swing ankle angle, especially when the walking speed of this step is too fast. This can thus enhance the walking stability of the next step. We also want to investigate how the changing of the swing ankle angle can influence the initial walking speed of the next step when the walking is faster. So we change the slope angle from 3 to 5 which can give a higher speed of the walking. The bigger slope angle value has the same effect with a bigger driving torque.

12 1 Advances in Mechanical Engineering As shown in Figure 1(b), the remained kinetic energy changes more than that under the slope angle value of 3. The remained kinetic energy can be reduced from to 0.064, which means that when the robot is walking faster the influence is getting bigger. In Hobbelen and Wisse, 1 the authors adjust the walking speed by adjusting the amount of ankle pushoff torque. This is more effective when the walking is too slow. In Hobbelen and Wisse, 5 the authors adjust the initial walking speed of the next step by swing leg retraction, which is an effective way to get better walking stability of the next step. However, changing the state of the swing leg costs more energy than changing the state of the swing ankle angle as the mass of the foot is much lower than the leg. Our method also provides a way to adjust the initial walking speed of the next step but with low energy consumption. And the method is more effective when the walking is too fast, which is a complementary method to that in Hobbelen and Wisse, 1 so combining the method in Hobbelen and Wisse 1 and in this article may obtain a better walking stability to resist disturbances during walking. Conclusion In this article, a passive walking model that consisted of two straight legs, two flat feet, two ankle joints, and a hip joint was built. A detailed walking dynamics was derived which contains double stance phase that makes the walking gait more natural. The instant just before impact is chosen as the initial condition for the walking, which reduces the number of the independent initial state variables from 5 to 3. The main features for a walking model with flat feet that contains double stance phase in this study is that (1) the double stance phase accounts for 5.5% of one step time and () there are two times of impact during walking and the second impact has tiny influence on the walking dynamics. A detailed analysis based on the numerical simulation for the influence of the swing ankle angle on the walking performance was introduced. For a fully passive walking robot, we can obtain a better walking stability in a reasonable range when changing the swing ankle angle from its equilibrium position. The step length could be changed about 5% and the step duration could be changed about 66.8% by just changing the swing ankle angle before first impact. A bigger swing ankle angle will lead to a lower initial walking speed of the next step. When the swing hinder-foot length is getting bigger, this influence is getting bigger. When it comes to the quasi-passive dynamic walking robot, we can reduce the amount of remained kinetic energy after impact about 5.8% when the swing ankle angle changes from 0 to 35 and enhance that about 9.9% when the swing ankle angle changes from 0 to 15. So we can obtain a more proper initial walking speed of the next step by just simply changing the swing ankle angle, especially when the walking speed of this step is too high. Thus, this might enhance the walking stability of next step when the robot has encountered disturbances that cause speed change of this step. When the robot is walking faster, this influence is getting bigger. Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article. Funding The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research is funded and supported by National Magnetic Confinement Fusion Science Program Multi- Purpose Remote Handling System with Large-Scale Heavy Load Arm (01GB10004). References 1. Sakagami Y, Watanabe R, Aoyama C, et al. The intelligent ASIMO: system overview and integration. In: IEEE/RSJ international conference on intelligent robots and systems, Lausanne, 30 September 4 October 00, vol. 3, pp New York: IEEE.. Strom J, Slavov G and Chown E. Omnidirectional walking using ZMP and preview control for the NAO humanoid robot. In: Baltes J, Lagoudakis MG, Naruse T, et al. (eds) RoboCup 009: robot soccer world cup XIII. Berlin, Heidelberg: Springer, 010, pp Collins SH and Ruina A. A bipedal walking robot with efficient and human-like gait. In: Proceedings of the 005 IEEE international conference on robotics and automation, ICRA 05, Barcelona, 18 April 005, pp New York: IEEE. 4. Kajita S, Kanehiro F, Kaneko K, et al. Biped walking pattern generation by using preview control of zeromoment point. In: Proceedings of the ICRA 03 IEEE international conference on robotics and automation, Taipei, Taiwan, September 003, vol., pp New York: IEEE. 5. McGeer T. Passive dynamic walking. Int J Robot Res 1990; 9: McGeer T. Passive walking with knees. In: Proceedings of the 1990 IEEE international conference on robotics and automation, Cincinnati, OH, May 1990, pp New York: IEEE. 7. Collins S, Ruina A, Tedrake R, et al. Efficient bipedal robots based on passive-dynamic walkers. Science 005; 307: Hobbelen D, De Boer T and Wisse M. System overview of bipedal robots flame and tulip: tailor-made for limit cycle walking. In: IEEE/RSJ international conference on intelligent robots and systems, IROS 08, Nice, 6 September 008, pp New York: IEEE.

13 Zang et al Schuitema E, Hobbelen DGE, Jonker PP, et al. Using a controller based on reinforcement learning for a passive dynamic walking robot. In: 5th IEEE-RAS international conference on humanoid robots, Tsukuba, Japan, 5 December 005, pp New York: IEEE. 10. Wisse M, Hobbelen DGE, Rotteveel RJJ, et al. Ankle springs instead of arc-shaped feet for passive dynamic walkers. In: 6th IEEE-RAS international conference on humanoid robots, Genova, 4 6 December 006, pp New York: IEEE. 11. Narukawa T, Yokoyama K, Takahashi M, et al. A simple 3D straight-legged passive walker with flat feet and ankle springs. In: IEEE/RSJ international conference on intelligent robots and systems, IROS 08, Nice, 6 September 008, pp New York: IEEE. 1. Wisse M, Hobbelen DGE and Schwab AL. Adding an upper body to passive dynamic walking robots by means of a bisecting hip mechanism. IEEE T Robot 007; 3: Hass J, Herrmann JM and Geisel T. Optimal mass distribution for passivity-based bipedal robots. Int J Robot Res 006; 5: Schwab AL and Wisse M. Basin of attraction of the simplest walking model. In: Proceedings of the ASME design engineering technical conferences and computers and information in engineering conference, Pittsburgh, Pennsylvania, 9 1 September, vol. 6, pp New York: ASME. 15. Huang Y, Wang Q, Xie G, et al. Optimal mass distribution for a passive dynamic biped with upper body considering speed, efficiency and stability. In: 8th IEEE- RAS international conference on humanoid robots, Daejeon, South Korea, 1 3 December 008, pp New York: IEEE. 16. Hsu CS. Global analysis by cell mapping. Int J Bifurcat Chaos 199; : Collins SH, Wisse M and Ruina A. A three-dimensional passive-dynamic walking robot with two legs and knees. Int J Robot Res 001; 0: Westervelt ER, Grizzle JW and Koditschek DE. Hybrid zero dynamics of planar biped walkers. IEEE T Automat Contr 003; 48: Chevallereau C, Gabriel A, Aoustin Y, et al. Rabbit: a test bed for advanced control theory. IEEE Contr Syst Mag 003; 3: Wisse M, Schwab AL, Van der Linde RQ, et al. How to keep from falling forward: elementary swing leg action for passive dynamic walkers. IEEE T Robot 005; 1: Hobbelen DGE and Wisse M. Ankle actuation for limit cycle walkers. Int J Robot Res 008; 7: Hobbelen DGE and Wisse M. Controlling the walking speed in limit cycle walking. Int J Robot Res 008; 7: Spong MW, Holm JK and Lee D. Passivity-based control of bipedal locomotion. IEEE Robot Autom Mag 007; 14: Hu rmüzlu Y and Moskowitz GD. The role of impact in the stability of bipedal locomotion. Dynam Stabil Syst 1986; 1: Hobbelen DGE and Wisse M. Swing-leg retraction for limit cycle walkers improves disturbance rejection. IEEE T Robot 008; 4: Huang Y, Wang QN, Gao Y, et al. Modeling and analysis of passive dynamic bipedal walking with segmented feet and compliant joints. Acta Mech Sinica 01; 8: Ruina A, Bertram JEA and Srinivasan M. A collisional model of the energetic cost of support work qualitatively explains leg sequencing in walking and galloping, pseudoelastic leg behavior in running and the walk-to-run transition. J Theor Biol 005; 37: Goswami A, Thuilot B and Espiau B. Compass-like biped robot part I: stability and bifurcation of passive gaits. INRIA research report no. 996, October Spong MW, Hutchinson S and Vidyasagar M. Robot modeling and control. New York: Wiley, Duindam V. Port-based modeling and control for efficient bipedal walking robots, 006, free.fr/bibliographie_net/duindam_phdthesis.pdf 31. Kim J, Choi CH and Spong MW. Passive dynamic walking with knee and fixed flat feet. In: IEEE international conference on systems, man, and cybernetics (SMC), Seoul, South Korea, October 01, pp New York: IEEE. 3. Hobbelen DGE and Wisse M. A disturbance rejection measure for limit cycle walkers: the gait sensitivity norm. IEEE T Robot 007; 3: Kuo AD. A simple model of bipedal walking predicts the preferred speed step length relationship. J Biomech Eng 001; 13: