For improving the energy efficiency of hopping robot, an asymmetric spring loaded inverted pendulum hopping model with leg mass is considered. The period orbit problem of twolegged hopping robot is investigated. Firstly, the hybrid dynamic model is constructed. Then the passive hopping gaits are found using quasinewton optimization method. Secondly, a PD controller is implemented to track the desired pitch trajectory of the body. Through applying control during stance phase, period orbits of the robot with offset body mass is obtained. Finally, the effect of the location of the leg mass and the body mass on hopping performances is investigated.
Keywords: Hopping robot, asymmetric, passive dynamics, periodic motion, energy efficiency
Hopping robots use isolated footholds, so it can negotiate uneven terrain. Due to these advantages, it has been studied by many researchers in recent years. Raibert and his coworkers had developed onelegged and twolegged hopping robots [1] based on a springloaded inverted pendulum (SLIP) model. These robots had a massless telescopic leg and a point foot. They use a threepart algorithm to control the robot. Since then, many jumping robots such as robots with spring leg [25] and robots with multijoint leg [68] are constructed.
The SLIP model is the simplest mathematical model for hopping robot and it captured the basic dynamicbehaviors of the running of many animals [9]. The control of the SLIP model is very simple because it doesn’t consider the dynamic coupling between the body and the legs. The control strategy based on SLIP model is necessary to be modified in order to successful apply in the control of the real robots [10], so it is necessary to construct a more realistic model similar to the real hopping robots.
In the hopping control, the energy consumption is especially crucial because it directly affects the size of the actuator and in return influences the weight of the robots. So, it is necessary to achieve the hopping gait with minimal power consumption. In this field, there are some remarkable researches on passive hopping. It has been found that using the proper initial conditions, the hopping robots can excuates periodic motion without input torque [1112]. Papadopoulos and Cherouvim proved that the passive gait is the least energyconsuming gait [2]. Ahmadi and Buehler used hip and leg compliance to realize passive motion [34]. Schammass et al exploited the passive dynamics of the system to reduce the energy consumption [13]. The simulation results shown that the robot including hip compliance reduced approximately 67% of the total energy compared to the robot without hip compliance. Hyon et al used a new gaitsearching algorithm to find the passive orbit and proposed two kinds of controllers to realize orbital stabilization [1416].
All of these hopping robots have symmetric configuration with the center of mass (COM) along the geometrical axis of a telescopic leg. In reality, the COM of most animals offsets the hip joint and the weight of the leg can’t be ignored. The offset bodymass offers greater flexibility in the design of hopping robot because the actuator can be placed within the body freely. The asymmetrical body mass configurations have been considered by several scholars [1723]. Kuswadi et al. proposed an asymmetric model of a singlelegged robot and realized a sustained hopping motion. The center of gravity of the model isn’t on the line of action of the actuator [19]. Sayyad et al analyzed a springylegged offsetmass (SLOM) configuration and proposed an energy compensated model to explore the possibility of a periodic motion [20]. But all these paper doesn’t consider the rotation of the body. Poulakakis et al used a SLIPembedding control law to control an asymmetric spring loaded inverted pendulum (ASLIP) hopping model, but the leg is assumed to be massless [2223].
In this paper, an asymmetric spring loaded inverted pendulum hopping model (ASLIP) is proposed, where the leg mass is considered and the center of mass of the body is not located in the hip joint. The goal is to study of the effect of the offset mass of the body and the mass of the leg on the robot’s hopping performance.
This paper is organized as follows. In Section 2, the dynamic model is presented along with its mechanical properties and motion equations. In Section 3, the Poincare mapping is used to analyze the passive periodic hopping motion of the hopping robot and the fix point is obtained using optimization, and the effects of the offset mass and the leg mass on the hopping motion are presented in Section 4.
Figure 1 is the simplified model of the ASLIP hopping robot. The hopping robot consists of two parts, a body and a springy telescope leg. The body and the leg is connected by a torsional spring on the hip joint . In the model, the hip joint doesn’t coincide with the center of mass of the body. The mass of the leg is also taken into account. A coupling between the leg swing motion and the body rotation is introduced, which made the analysis of the dynamics of the robot more complicated. Supposing the position of the center of mass of the leg and the position of the center of mass of the body are given by and , where . The variable denotes the length of the leg, represents the leg angle with respect to the vertical axis indicates the body angle with respect to the horizontal axis. Other physical parameters and its values are listed in Table 1.
Figure 1. Model of robot 
The hopping motion of the robot is composed of two phases: stance phase and flight phase. The switch from the stance phase to the flight phase is called liftoff, triggered by . The switch from the flight phase to the stance phase is called touchdown, triggered by . We use the eventbased subscript to describe the motion during different phases. The subscript represents the stance phase, represents the flight phase, is the moment at touchdown and is the moment at lift off.
Notation  Units  Value 

body mass  kg  9 
leg mass  kg  3 
body length  m  0.22 
body inertia  kg.m^{2}  0.5 
leg inertia  kg.m^{2}  0.11 
leg spring stiffness  N/m  3000 
hip spring stiffness  Nm/rad  10 
natural leg length  m  0.5 
gravity acceleration  m/s^{2}  9.81 
During the stance phase the leg is in contact with ground at a toe point (end of the leg). Supposing that the robot does not bounce back nor slip relative to the ground, the leg rotates around the toe. Then the toe can be modeled as a pivot. We choose polar coordinates , as the generalized coordinates. Using the Lagrange method [22], the equations of motion is obtained as:

(1) 
where is the symmetric positive definite inertia matrix during the stance phase, which deponds on the generalized coordinates ; is the matrix containing centrifugal force and Coriolis force during the stance phase, which is a function of and ; is the gravity vector during the stance phase, which is a function of ; is the equivalent driving torque vector of torsion spring during the stance phase, which is related to ; ; is the driving force to the leg, while denotes the hip torque during the stance phase.
During the flight phase the leg loses contact with the ground. The COM of the robot moves along the ballistic trajectory and the angular momentum of the robot about the COM is conserved. The length of the leg is assumed to keep constant. We select Cartesian coordinates as the generalized coordinates, is the position of the hip joint. The dynamics is given by:

(2) 
where is the symmetric positive definite inertia matrix during the flight phase, which deponds on the generalized coordinates ; is the matrix containing centrifugal force and Coriolis force during the flight phase, which is a function of and ; is the gravity vector during the flight phase, which is a function of ; is the equivalent driving torque vector of torsion spring during the flight phase, which is related to ; ; is the hip torque during the flight phase.
At touchdown, the leg touches down the ground and inelastic impulse occurs. The generalized displacements keep constant, but the generalized velocities change instantaneously

(3) 
where the subscript and represent the moment before and after touchdown respectively.
At liftoff, there are no discontinuous changes of states except for:

(4) 
It can be seen that during the landing phase, the leg compresses and releases the spring to realize the energy conversion from the above equations. Since the loss caused by friction is not considered, the mechanical energy is conserved. In the flight phase, the mechanical energy of the system is also conserved. At touchdown, the robot legs collide with the ground, which will lead to energy loss. In order to simplify the analysis, the energy loss at touchdown is not considered in this paper. The robot is a complex switching system with multi variables, strong coupling and nonlinearity, so it is impossible to obtain the analytical solution of the trajectory of the system, and can only be solved by numerical calculation.
To improve the energy efficiency, firstly we consider the completely passive hopping which means that there is no driving force applied on the robot. A very useful and common method to study the existence and stability of periodic motions is the Poincare map. According to the trajectory of the state variables in the Poincare section, we can find the periodic motion. Through Poincare map, we can transform the existence of the periodic orbits of the hopping robot to the fixed point in the Poincare section.
The robot has 4 DOF during flight and 3 DOF during stance. For reducing the number of the independent variables, we choose the moment just after touchdown as the Poincare section and state vector is defined as . For at touchdown, the number of the independent variables reduces to 5. In the process of iteration, let the state vector in the step is , we can get the state vector in the step after one mapping. This process can be described as a Poincare map .
For a certain state vector , the point is called a fixed point if it maps itself through one map,

(5) 
Then we build the Poincare map of the hopping robot. The dynamic Eqs. (1) and (2) are rewritten into a standard statespace form.

(6) 

(7) 
where is the state vector during the stance phase and during the flight phase. Because the number of the state variables between the stance phase and the flight phase is different, we need to construct two mapping between the transitions of these two phases. The mapping is represented by the following equations

(8) 

(9) 
Combining these four submaps, then the whole Poincare mapping is .
Since the Poincare map of the hopping robot have been constructed, now we consider the method of finding the fixed point. The NewtonRaphson’s method [15] is commonly used to find the fixed point

(10) 
where, is the identity matrix, is the jacobian matrix of the Poincare map at .
Because the NewtonRaphson’s method demands the initial value of the state vector to be very close to the accurate value of the fixed point, it’s difficult to give the initial value. If the initial value is selected improperly, the algorithm will diverge soon.
In this paper we used the optimization algorithm to solve the nonlinear equations. Firstly we transformed the problem into a minimum value problem. Eq. (5) can be rewritten as

(11) 
Then the problem becomes an unconstrained optimization problem. We can use the quasiNewton method to find the minimum value of Eq.(11). If the minimum value of is zero, the point is considered as the fixed point. In reality, due to calculation error, the point is considered as the fixed point if .
In this paper, quasi Newton method is used to solve the minimum problem

(12) 
where , , , .
Using the proposed algorithm, passive hopping gaits are found. According to the Hyon [14], all passive running gaits of the ARL monopod robot were unstable except for a vertical hopping. The obtained gaits of our hopping robot are also quasiperiod ones, not stable limit cycle. The distance between the COM of the leg and the toe is , and the distance between the COM of the body and the hip joint is . The location of the COM of the leg and the body is depend on and , respectively, we examine the influence of and on the jumping performance of the robot.
Firstly the effect of α is examined when the COM of the body is coincide with hip joint ( is zero). The base parameters given in Table 1. The results are shown in Figure 2. The location of the COM of the leg is displaced along the leg from the toe to the hip joint, which means that the value of is between 0 and 1. When , the center of the leg is located in the toe and the model is similar to a twomass model; when , the center lies in the hip joint and the model is the ARL monopod. In Figure 2 it can be seen that when α increases, the apex height of the COM of the robot and the compression of the leg spring decreases, while the step length increases. The step length reached at the maximum value when . The following conclusions can be drawn. If the robot is wanted to have a long jump, the COM of the leg should be close to the hip joint. If the robot is wanted to have a high jump, it should be far away from the hip joint and at the same time the robot should provide greater power by increasing the compression of the spring.
Figure 2. The step length , the apex height and the minimum value of the leg length during stance phase 
In Figure 3 with the increase of α, the angle of the center of mass at the takeoff time (takeoff angle) decreases, the horizontal velocity of the center of mass increases, and the vertical velocity decreases. With these velocities are approximately equal. The decrease of vertical velocity means the decrease of flight time, but the step size of the robot increases with the increase of horizontal velocity. Due to the dual influence of flight time and horizontal velocity, the robot . The maximum step size is reached at 0.7.
Figure 3. The horizontal and vertial velocity of the COM at takeoff 
During the change of leg centroid, the stance time changed very little and remained around 0.2s (Figure 4). However, as the COM of the leg gets closer to the toe, the total jump time increases and the duty cycle decreases, which means that the time spent by the robot in the takeoff phase increases and the gait frequency decreases. Before the next jump, the swing of the leg is controlled more easily.
Figure 4. Bounce period, stance time v.s. 
In Figures 5 and 6, it can be seen that with the increase of α, the center of mass of the leg is closer to the hip joint, and the body angle amplitude decreases, but the leg angle amplitude increases. The closer the center of mass of the leg is to the hip joint, the closer the leg swing curve is to the sine curve. When the center of mass of the body is located at the hip joint, the leg swing curve is a sine curve.
Figure 5. Trajectory of the leg 
Figure 6. Trajectory of the body 
Then we investigate the effect of the COM of the body on the hopping performance, it can be expressed by . The periodic gait of the hopping robot without active control is in Figure 7 when the COM of the body doesn’t locate in the hip joint. It shows that there was little change in the angle of the leg and the body, the motion of the robot is a vertical hopping and the robot doesn’t leave the ground, so it is necessary to control the robot in order to get a forward hopping. The control scheme is using the hip actuator to control the pitch angle to track the desired trajectory during the stance phase. A PD (ProportionalDerivative) controller is used to track the desired pitch trajectory [2426].
Figure 7. Periodic passive hopping gait without active control 
Using the pitch angle of passive hopping of ARL robot as the desired trajectory and computing the limit cycle with the optimization method, the passive gait of the asymmetrical hopping robot is obtained. Figures 8 and 9 are the posture and the limit cycle during one hoping cycle. It is can be seen that with active control, a satisfactory periodic motion is obtained when the robot has an asymmetric body mass. It can be seen that the toes of the robot collided with the ground in the flight phase. In order to prevent the toes of the robot from contacting the ground in the flight stage, the jumping height of the robot should be increased.
Figure 8. Active periodic gait of the hopping robot when , 
Figure 9. Limit cycle of the hopping robot 
Although in theory can be taken in [0,1], is less than 1 for the following two reasons. The first reason is that it is unable to find fixed point. The second reason is that the minimum length of the leg length at the landing stage is 0, which means that the compression of the spring reaches 100%. In fact, due to the structure of the spring itself, the compression cannot be pressed to the bottom. The above two cases can be increased by improving the control method and increasing the spring stiffness to increase the range of . As can be seen from Figure 10, with the increase of , the maximum allowed value of is small. It shows that the closer the robot leg is to the hip joint, the smaller the allowable eccentricity of the body.
Figure 10. The maximum allowable value of v.s. 
Figure 11 shows that with the increase of , the angle of the center of mass at the takeoff decreases, the vertical velocity remains almost unchanged, while the horizontal velocity increases rapidly, and the step size of the robot becomes larger. The farther the body centroid is from the hip joint, the better. Because the change of bounce cycle time is very small, most of the energy provided by the leg spring is transformed into forward movement. However, due to the low centroid height of the robot at the takeoff, the maximum height in the flight phase is also low, and the robot toe is easy to collide with the ground. So cannot be too large.
Figure 11. The horizontal and vertical velocity of the center of mass at takeoff 
As can be seen from Figure 12, with the increase of , the minimum length of the leg in the landing stage decreases rapidly, and the compression amount of the spring changes from 38% to 69%, indicating that the farther the center of mass of the body is from the hip joint, the greater the force required to be provided by the leg spring. It should be considered to increase the stiffness of the spring.
Figure 12. Minimum length of legs in landing phase with different of 
Figure 13 shows that with the increase of , the center of mass of the body deviates more and more from the hip joint, and the rotation amplitude of the body and legs increases. The closer the body centroid is to the hip joint, the closer the leg swing curve is to the sinusoidal curve. The angle of the body and leg increases at the takeoff, correspondingly the takeoff angle of the body centroid decreases. Body eccentricity provides a larger takeoff range for robot motion.
 
Figure 13. Leg angle and the body angle vs. time 
Figure 14 shows that the greater the body eccentricity, the greater the driving torque required by the robot to drive the body pitching motion in the landing phase, because the swing range of the body in the landing phase also becomes larger.
Figure 14. Influence of body eccentricity ratio on driving torque 
The period orbit problem of an asymmetric spring loaded inverted pendulum is addressed. The dynamic equations of the two link robot are established. Then the Poincare mapping method is used to turn the solution of the periodic motion of the hopping robot into the solution of the fixed point, and then the optimization method is used to obtain the fixed point of the hip robot. With a PD controller, a period orbit can be obtained. Through investigating the effect of the location of the leg mass and the body mass on the hopping performance, if the robot is wanted to jump longer, the center of mass of the leg and the body should be far away from the hip joint.
[1] Raibert M.H. Legged robots that balance. The MIT Press, Cambridge, pp. 233, 1986.
[2] Zeglin G. The bow leg hopping robot. Carnegie Mellon University, pp. 161, 1999.
[3] Ahmadi M., Buehler M. A control strategy for stable passive running. Proc. IEEE Conf. Intelligent Robots and Systems, 3:152157, Pittsburgh, PA, 1995.
[4] Ahmadi M., Buehler M. Stable control of a simulated onelegged running robot with hip and leg compliance. IEEE Transactions on Robots and Automation, 13(1):96104, 1997.
[5] Zeglin G.J. Unirooa one legged dynamic hopping robot. Massachusetts Institute of Technology, Thesis, 1991.
[6] Recovering, storing and releasing energy based on a natural model. 2014. http://www.festo.com/cms/en_corp/13704.htm
[7] Hutter M., Remy C.D., Hoepflinger M.A., Siegwart R. ScarlETH: Design and control of a planar running robot. IEEE/RSJ International Conference on Intelligent Robots and Systems, San Francisco, CA, USA, pp. 562567, 2011.
[8] Blickhan B. The springmass model for running and hopping. J.Biomechanics, 22(1112):12171227, 1989.
[9] Cherouvim N., Papadopoulos E. Single actuator control analysis of a planar 3DOF hopping robot. Robotics: Science and Systems I, Massachusetts Institute of Technology, Cambridge, Massachusetts, June 811, 2005.
[10] Papadopoulos E., Cherouvim N. On increasing energy autonomy for a onelegged hopping robot. Proceedings of the 2004 IEEE International Conference on Robotics and Automation, pp. 46454650, 2004.
[11] François C., Samson C. A new approach to the control of the planar onelegged hopper. International Journal of Robotics Research, 17(11):11501166, 1998.
[12] Thompson C.M., Raibert M.H. Passive dynamic running. The First International Symposium on Experimental Robotics I, pp. 7483, 1989.
[13] Schammass A., Caurin G.A.P., Valente C.M.O. Control of a onelegged robot with energy savings. Journal of the Brazilian Society of Mechanical Sciences & Engineering, 23(1):4148, 2001.
[14] Hyon S., Emura T. Quasiperiodic gaits of passive onelegged hopper. IEEE/RSJ International Conference on Intelligent Robots and Systems, 3:26252630, Lausanne, Switzerland, 2002.
[15] Hyon S., Emura T. Energypreserving control of a passive onelegged running robot. Advanced Robotics, 18(4):357381, 2004.
[16] Hyon S., Emura T., Ueta T. Delayed feedback control of onelegged passive running robot. Sice 2004 Annual Conference, pp. 949954, Sapporo, Japan, 2004.
[17] Lei B., Wang P., Zha F., Li M., Guo W. Stability analysis of quasiperiodic hopping of a passive onelegged robot with compliant hip joint. IEEE International Conference on Mechatronics and Automation, pp. 623628, Tianjin, China, 2014.
[18] Wei T.E., Nelson G.M., Quinn R.D., Verma H., Garverick S.L. Design of a 5cm Monopod Hopping Robot. Proceedings 2000 ICRA, Millennium Conference, IEEE International Conference on Robotics and Automation, 3:28282833, San Francisco, CA, USA, 2000.
[19] Kuswadi S., Ohnishi A., Takahashi A., Sampei M., Nakaura S. A one linear actuator hopping robot：modeling and control. Advanced Robotics, 17(8):709737, 2003.
[20] Sayyad A., Seth B., Issac K.K. Dynamics and control of a onelegged 2D SLOM popping robot. 12th IFToMM World Congress, Besançon, 2007.
[21] Seth B., Seshu S., Shanmuganathan P.V., Vichare V.V., Raj P. Search for initial conditions for sustained hopping of passive springyleg offsetmass hopping robot. Journal of Dynamic Systems, Measurement, and Control, 129(4):522526, 2007.
[22] Poulakakis I., Grizzle J.W. The spring loaded inverted pendulum as the hybrid zero dynamics of an asymmetric hopper. IEEE Transactions on Automatic Control, 54(8):17791793, 2009.
[23] Poulakakis I. Spring loaded inverted pendulum embedding: Extensions toward the control of compliant running robots. 2010 IEEE International Conference on Robotics and Automation (ICRA), pp. 10504729, Anchorage, Alaska, USA, 2010.
[24] Terry P., Piovan G., Byl K. Towards precise control of hoppers: Using high order partial feedback linearization to control the hopping robot FRANK. IEEE 55th Conference on Decision and Control (CDC), pp. 66696675, Las Vegas, USA, 2016.
[25] Zou H., Schmiedeler J.P. The effect of asymmetrical bodymass distribution on the stability and dynamics of quadruped bounding. IEEE Transactions on Robotics, 22(4):711723, 2006.
[26] Nie H., Sun R., Hu L., Su Z., Hu W. Control of a cheetah robot in passive bounding gait. Journal of Bionic Engineering, 13(2):283291, 2016.
Published on 13/01/23
Accepted on 01/01/23
Submitted on 15/10/22
Volume 39, Issue 1, 2023
DOI: 10.23967/j.rimni.2023.01.002
Licence: CC BYNCSA license
Are you one of the authors of this document?