Stable periodic motion of a controlled segmented leg model of pedal locomotion with inelastic ground-foot collision

Human and robotic legged locomotion can be described with complex multi-degree-of-freedom dynamic models, whose bifurcation or parameter analysis may explain some features of typical patterns during motion. In this paper, we focus on the effect of kinematic parameters and foot placement techniques on the ground-foot impact intensity. The work is based on a multibody dynamic model of a segmented leg, which possesses some fundamental characteristics of locomotion systems: (a) distinct topology in the flight and ground phases; (b) kinetic energy absorption due to partially/fully inelastic ground-foot collision; (c) an active control strategy for maintaining a prescribed mechanical energy level; and (d) different control strategies for the flight and ground phases. We obtain a quantitative measure for the foot collision intensity by analytic calculations. The pre-impact velocity conditions are obtained from a hopping three-segmented planar leg model that imitates pedal locomotion. The single-legged model contains the foot, the shank, the thigh and a reaction wheel attached to the hip, which models the effect of the upper body. The existence of stable periodic motion associated with hopping is shown for a wide range of parameters by means of finding suitable control torques in the ankle, the knee and the hip joint. The parameters of the linear feedback controller are tuned to optimize different cost functions, such as running speed, energy efficiency and impact intensity. We also investigate how the stability of the periodic motion depends on the control gains.

control strategy for maintaining a prescribed mechanical energy level; and (d) different control strategies for the flight and ground phases. We obtain a quantitative measure for the foot collision intensity by analytic calculations. The pre-impact velocity conditions are obtained from a hopping three-segmented planar leg model that imitates pedal locomotion. The singlelegged model contains the foot, the shank, the thigh and a reaction wheel attached to the hip, which models the effect of the upper body. The existence of stable periodic motion associated with hopping is shown for a wide range of parameters by means of finding suitable control torques in the ankle, the knee and the hip joint. The parameters of the linear feedback controller are tuned to optimize different cost functions, such as running speed, energy efficiency and impact intensity. We also investigate how the stability of the periodic motion depends on the control gains.

Introduction
Analyzing the dynamics of legged locomotion is a challenging task for several reasons. The associated mechanical models involve multiple degree-offreedom (DoF) systems, varying topology, alternating over-and under-actuated [1] motions related to the different phases of walking or running, nonlinearities related to geometrical configurations [2], and non-smoothness at foot impacts [3][4][5][6][7][8][9]. Moreover, such models often involve some kind of motor pattern generator [10]. Underlying feedback mechanisms introduce additional components into models, such as reaction time delay, sensory uncertainties and the saturation of the forces exerted by the muscles [11]. Furthermore, the control concept employed by the nervous system to keep the body segments on the desired path is not exactly known.
Still, there are high DoF models in the literature for legged locomotion, which are used for the interpretation of experimental results, such as the calculation of joint forces derived from motion capturing [12,13]. These models take the geometry of the human body into account precisely; however, important aspects such as balance, motion pattern generation and control are generally not included.
Another group of models includes the balancing and control strategy of the brain; however, their geometry is reasonably simplified. Such models help one understand how the nervous system keeps the body in balance, how the energy level is maintained and how the motion patterns are generated during locomotion. The balancing process can be understood better by including a controller that mimics the operation of the nervous system [14]. For instance, the spring-loaded inverted pendulum (SLIP) model in [15] involves a sophisticated control algorithm that takes stability and control issues into consideration, while the other elements of the model are excessively simplified compared to the human body. The model in [16], on the other hand, features a geometry very similar to the human leg, but the model still neglects the mass of leg segments. Consequently, the important issues of ground-foot collision and impact-induced energy absorption are missed. In fact, a considerable amount of energy-demand in locomotion is related to ground-foot impact. Furthermore, foot impact is the main source of injuries when walking and running.
In this paper, we present the stability analysis of the three-segmented controlled leg model developed in [17]. The model describes a stable periodic motion associated with hopping. The goal is to present a dynamic analysis of a moderately complex model, for which the number of parameters is still small enough to be treated with qualitative mathematical methods. The new contribution of the present paper is that the local behavior of the model in [17], namely the stability of periodic solutions, is explored as a function of the most important control parameters. The model includes: (1) human-like geometry; (2) energy absorption induced by ground-foot impact; and (3) control for stabilization. Mathematical tools for piecewise-smooth periodic systems are applied to demonstrate the feasibility of the model and to analyze qualitatively the global behavior of different types of running. More specifically, we present stable forward, zero velocity and backward hopping locomotion. We then perform an analysis of how the forward hopping locomotion and its stability depends on parameters and how it can be optimized in terms of locomotion velocity, apex height, energy efficiency and ground-foot impact intensity. Finally, going well beyond what has been reported in [17], we address biomechanical issues: the effect of tiptoe positioning on locomotion velocity and impact intensity.

Motivation and background
Important elements of our model and analysis are ground-foot interaction and stability control, which are both crucial for efficient running.
Depending on the form of running, the ground-foot interaction may lead to high-intensity impacts or collision. This may result in high kinetic energy absorption, which should be avoided in order to minimize the risk of injuries and also in order to increase energy efficiency during running [5,6,9,18,19]. There is a debate in the literature about the effect of running kinematics and foot placement patterns on the ground-foot collision intensity. For instance, Souza [20] states: "At this time, there is limited evidence that any foot strike pattern is more or less likely to cause a runner to sustain an injury". Another example is Jungers [18] who says "More studies ... are required to provide data instead of opinion, and testable models and scientific explanation instead of anecdotes". Therefore, a reliable mathematical model of legged locomotion could strongly contribute to the understanding of legged locomotion.
Despite considerable research efforts [21,22], the neural processes behind human balancing tasks are still not completely understood. An interesting approach is to compare the required mechanical energy consumption when standing still and during locomotion-when walking, running and hopping. Postural sway observed during standing still clearly demonstrates that certain energy input is required for the stabilization of the otherwise unstable upright standing position. Locomotion, on the other hand, is typically a periodic-like motion, which requires energy from the muscles to both overcome losses and to ensure stability, i.e., prevent falling. A long-term goal is to answer the question whether the energy required for stabilization is smaller in the case of locomotion compared to standing still.

Hopping models and ground-foot collision
Many aspects of periodic motion of legged systems were detailed in [15], where several walking, running and hopping models are reviewed. The classical SLIP model provides a fundamental background for hopping leg models. It represents the body inertia as a point mass which bounces along on a single elastic massless leg. The supporting role of the leg in each stance phase is characterized by the inverted pendulum. However, geometric features of a human leg and the groundfoot collision are not described by the SLIP model. The stable static configurations of a three-segmented leg model were presented in [2] for a case where the inertia and dynamic effects of the body segments are neglected and the body weight is modeled by a point mass located at the hip. Self-stable running-like locomotion for this model was presented in [16], showing that the stable domain in terms of landing angle and horizontal velocity is larger than for the SLIP model. Apart from the elimination of inertial forces, masslessness of the segments leads to smooth dynamics, i.e., there are no impulsive forces present when the foot comes into contact with the ground. This smoothness property allows one to use well-established analytical and numerical approaches for smooth system.
A mechanical model containing the foot and the shank was introduced in [19] for the investigation of the impulsive forces and the effect of the foot strike pattern. The foot-shank system falls vertically and the foot hits a fixed point; the contact point position is characterized by the strike index s = l s /l, where l s indicates the location of the resultant force F p of foot pressure and l is the foot length; see Fig. 1. An extension of the two-segmented model of [19] was introduced in [5] by involving a variable ankle angle and a horizontal velocity. Both [5] and [19] report that the effective mass and, therefore, the impact intensity is lower for forefoot-strike than for rear-foot-strike (heel- Results in [6] showed that, in addition to the strike pattern, the angle of the shank also strongly affects the properties of the ground-foot impact. The model in [6] involves the thigh and the total mass of the human body; the segments perform a rigid-body-like motion in the pre-impact phase with no relative motion of the body segments. However, this is not realistic and, in order to perform further biomechanical analysis of human running performance, we introduced a dynamic model with torque control in [17,23] as an extension of the model from [16]. This model includes the masses of the segments so that the effects of the inertial forces can be determined. Moreover, it provides measures for the ground-foot collision, including kinetic energy absorption and effective mass, and it also provides a realistic pre-impact configuration.
The advantage of the presented numerical study over an experiment is that in the model the control parameters can be tuned to optimize different cost functions, such as running speed, energy efficiency and impact intensity. In this paper, we demonstrate the feasibility of this model by reporting results of a detailed analysis, including of the stability of different types of hopping motion; the mathematical details of the model are given in 3.1.

Dynamic model of the controlled hopping leg
Hopping and running motion can be divided into flight (airborne) and ground (stance) phases, which are connected by ground-foot collisions [24]. For running and hopping robotic or biomechanical systems, some segments of the limbs periodically come into contact with the ground; therefore, the motion of these body segments are constrained during the duration of the stance phase. The intermittent constraints of the limb(s) lead to the absorption of a certain portion of the kinetic energy, referred to as constrained motion space kinetic energy (CMSKE) [3,5,6]. Since the energy loss of passive walkers [25] is caused mostly by the foot impact, CMSKE is often used as an energy efficiency indicator. It can also be used as an indicator of foot impact intensity, since it is proportional to the contact reaction force impulse and also to the peak reaction force [5,26]. During the stance phase, kinetic energy accumulates in the rest of the body while the muscles provide further mechanical power; finally, the ground-foot contact phase terminates.
The model used in the present study was developed to include and analyze the foot-ground impact described above. The body segments move downward before the end of the flight phase. The flight phase terminates when the tiptoe, denoted by point A in Fig. 2, comes into contact with the ground. During the ground phase, the active spring-damper system increases the total kinetic energy of the system. As a result, the height of the next hop is not lower than the height of the previous hop, in spite of the energy lost. The controller (detailed in Sect. 3.2) aims to maintain a nominal mechanical energy level E 0 . Figure 2 presents a diagram of the mechanics of the model. Segments 1, 2 and 3 correspond to the foot, shank and thigh, respectively. Similarly, the points A, B, C and D correspond to the tiptoe, the ankle joint, the knee joint and the hip joint, respectively. The reaction wheel plays the role of the upper body, and the torque M D is acting between the reaction wheel and the thigh. The mass of the reaction wheel and its moment of inertia with respect to the y-axis through point D are m r and J r , respectively. The three homogeneous prismatic bars are characterized by their masses m i and lengths l i , i = 1, 2, 3; they are connected by torsional springs of stiffness k B and k C . The overall center of mass (CoM) is located at point G. Control torques are applied at joints B, C and D, as it is shown in the right panel of Fig. 2. Actuating torques M B and M C implement the motion described by the control law, which is introduced below in Sect. 3.2.

Mechanical structure
The model has a total of 6 DoFs and the general coordinates are chosen to be q = [x A , z A , θ 1 , θ 12 , θ 23 , θ r ] T , where x A and z A are the Cartesian coordinates of the tiptoe, θ 1 is the foot angle measured from the horizontal, θ 12 and θ 23 are the relative angles of the ankle and the knee, respectively, and θ r is the angle of the reaction wheel. The tiptoe is fixed to the ground during the ground phase.
The equation of motion assumes the general form: whereq ∈ R n andq ∈ R n are the general velocity and the acceleration, respectively, with n = 6. The general mass matrix is denoted by H ∈ R n×n , and C ∈ R n represents the inertial force vector. The general force Q ∈ R n contains the gravitational forces, the torques of the torsional springs and the control torques. Nonlinear spring characteristics due to the muscletendon dynamics, as mentioned in [2], is not considered here in order to keep the number of parameters manageable. The torques exerted by the linear torsional springs of stiffness k B and k C are given by and where α 12 and α 23 denote the ankle and knee joint angles corresponding to the unstretched springs, respectively.

Torque control
The control law determines the required active torques in Eqs. (3)-(5) and (7)-(9) based on the motion of the segments and information about the actual phase. Control concepts based on trajectory tracking and those based on feedforward inputs form the central pattern generators are clearly distinguished and described in greater detail in [4] and in [15]. Our control approach rather fits into the latter category in that we consider a type of feedback for the control; the strategy of keeping the leg segments on some prescribed trajectories is not employed as part of the work presented here.
In the flight phase, the control torques M B and M C behave as damping torques and suppress the possible oscillations of the relative segmental angles as (3) and (4) show. The horizontal position x A of the tiptoe relative to horizontal position x G of the CoM has a large effect on stability, on locomotion velocityv and on impact intensity. Approaching and reaching the nominal horizontal position of point A is facilitated by a proportional-derivative (PD) controller defined in (5) below. We consider this as the foot touchdown preparation. The torques M F B , M F C and M F D are given by where D B , D C , P and D are tuned control gains. The nominal tiptoe position is modified by where Π A is the angular momentum about the point A, P Π is an associated control gain and K v is a control parameter responsible for tuning the nominal locomotion speedv. During the ground phase, CMSKE is recovered by means of the control torques (7) and (8) where P E , P r and D r are control gains.

Flight and ground phase transitions
The system is described by the state variable vector x = q,q T , where x ∈ R 2n . The system dynamics is given bẏ where and (12) are smooth vector fields that describe the flight and the ground phase, respectively. Here,q red = [θ 1 ,θ 12 , θ 23 ,θ r ] T is a reduced general velocity vector; similarly, the reduced mass matrix H red is obtained by truncating the first two rows and columns of H, while C red and Q red are obtained by truncating the first two rows of C and Q, respectively. Equation (13) in effect imposes x A = 0 andz A = 0 for the tiptoe coordinates.
The overall system features switching and discontinuities of the solutions due to the transitions between the two phases of motion, which also need to be defined to determine the dynamics; one speaks of a hybrid system [3,27,28]. The transition between (10) and (11) occurs at codimension-one switching (hyper) surfaces Σ G2F and Σ G2F ; these are illustrated in Fig. 3.

Flight to ground
The flight phase is terminated by the foot touchdown. In the presented model, the ground contact is modeled as a single point contact at point A. The contact force at any other body segments is not considered. In certain parameter regions, this leads to stable hopping solutions while point B penetrates into the ground, which is physically unfeasible. This modeling issue arises commonly in single contact point models, see, e.g., [5,19]. Physical feasibility is monitored, and unfeasible solutions are detected in Sect. 4.
The flight to ground (F2G) transition occurs when Σ F2G , defined by is reached. The F2G transition involves an impact, which corresponds to a discontinuity of the solution. We assume that the ground-foot collision is instantaneous and completely inelastic so that there is no rebound [5,6,19]. Furthermore, we assume that the friction coefficient is large enough so that the foot does not slip on the ground; this means that it satisfies where γ q,n and γ q,t are the gradients of the constraints in the normal and tangential directions, respectively. Under these assumptions, the tiptoe remains constrained to a single point on the ground during the ground phase, as in [5,6,17,23,29]. Consequently, the velocitiesẋ A andż A become zero and the angular velocities change abruptly when the solution reaches Σ F2G . The jump function maps the state to the specified location [30], by projecting the velocities into the admissible direction of the constraints γ (q) [6,31,32]. Here is the projection matrix onto the subspace of constrained motion.
The instantaneous ground-foot collision leads to infinitely large instantaneous forces, over infinitesimal time duration, so that the net impulse due to the impact force is finite [5,19]. Since CMSKE is proportional to the impulse of the contact reaction force and also to the peak reaction force [5,26], it is used for characterizing the foot impact intensity in this work. CMSKE, which vanishes at foot touchdown, is calculated with the pre-impact general velocity as in [26] as The papers [5] and [26] showed that foot strike intensity can be characterized by the CMSKE, which depends on the pre-impact configuration and velocity, and the effective mass matrix H e = P T c HP c . The effective mass concept for foot impact is introduced in [33] for a one DoF model of legged locomotion.

Ground to flight
The ground to flight (G2F) transition at the surface Σ G2F , defined by is a switch from (11) to (10) without a discontinuity of the solution. Here, the Lagrange multiplier λ z is the vertical ground contact force acting on the tiptoe, and a negative and positive sign of λ z refers to a pushing and a pulling force, respectively. The contact force is calculated by solving the linear algebraic system adopted from [31]. The second time derivativeγ = γ qq +γ qq of the geometric constraint vector γ (q) appears in (20), where is the gradient of the constraint. The ground-foot contact is represented by γ (q) = 0 with Control torques included in Q(q,q) in (20) are calculated with the control law (7)-(9) in the ground phase.

Finding stable relative periodic orbits
Our focus is on hopping motion, which corresponds to a relative periodic orbit with a translational symmetry in the lateral ground direction, given by the variable x. Such an orbit is shown in the phase space sketch in Fig. 3.
A relative periodic orbit is found with a shooting approach as a root x * 0 of where u 0 (x 0 ) is a subvector of a point x 0 ∈ Σ F2G and u e (x 0 ) is the corresponding subvector at the end point of the solution trajectory of the overall system that starts at x 0 and ends again in the section Σ F2G by reaching Σ G2F . We remark that u 0 has lower dimension than x 0 , because some of the state variables (both coordinates and velocities) are known when the tiptoe is in contact with the ground. This reduction in the number of unknowns considerably reduces the computation time.

Monodromy matrix and Floquet multipliers
To determine the stability of the relative periodic orbit given by x * 0 (see Fig. 3), we proceed as described in [3,30,34]. During the flight and ground phases we integrate from t = 0 for the two vector fields f F and f G , respectively, the first variational equatioṅ where I is 2n × 2n identity matrix and f x is Jacobian of the vector field. This results in the fundamental solution matrices Φ F (t) and Φ G (t) [34]. The Jacobian or monodromy matrixΦ(t) over the entire relative periodic orbit is then obtained as the compositioñ Here, S G2F and S F2G describe the contribution of the mappings at the switchings, which are known as the saltation matrices and are related to the Jacobian of the transition condition [35,36]. A saltation matrix can be determined [3,30] as where g − x is the gradient of the jump function, and h − x is the gradient of the indicator function [see (14) and (19)] evaluated at the state just before the transition. For the G2F transition, S G2F is obtained from (27) by applying the substitutions For the F2G transition, the discontinuity mapping g F2G (16) is taken into account, and substitutions are, respectively, the state just before and just after the impacting transition (see Fig. 3).
The stability of the relative periodic orbit is determined from the eigenvalues of the monodromy matrix Φ(t G2F ), which are known as the Floquet multipliers [34]. There is always a Floquet multiplier equal to 1, because of periodicity; moreover, due to the translation in the x-direction there is a second Floquet multiplier 1. The relative periodic orbit is stable when all other Floquet multipliers have modulus less than 1.

Tuning control parameters
Tuning of the physical and control parameters is crucial for the study of the model. Our goal here is to find the control parameter set associated with a stable periodic motion for which some cost function, such as energy use, is minimized. We collect 10 parameters of (3)- (9) in the vector p p = [α 12 , α 23 , D B , D C , P, D, An initial parameter set p, which provides stable hopping, was found using a trial-and-error method as reported in [17]. To tune the parameters p, the associated solution x * 0 of (23) was found by continuation for a parameter sweep during which minima of the given cost function were detected. Note that if the cost function is chosen as the norm of the largest eigenvalue ofΦ, then it characterizes stability. Another possible choice is to use CMSKE [denoted T c in (18)] as the cost function, which characterizes impact intensity and energy efficiency.

Results on stable locomotion
We consider now the test-case example of a 24-yearold male subject with bodyweight of 73 kg and body height of 173.1 cm, adopted from [37] with further anatomical data as given in Table 1. The parameters of the torsional springs at the joints are independent of human  Table 2 Parameters of the torsional springs placed at the ankle and the knee joints anthropometry and are shown in Table 2. The control parameters used in the study presented here are given in Table 3. For these choices of parameters, relative periodic motion (i.e., hopping) was found. In what follows, the parameters E 0 and K v are considered as the main parameters, because they are responsible for the tuning of the apex height z A max and the average horizontal locomotion velocityv of the hopping motion, respectively. The rest of the parameters remain fixed. The model may display stable forward hopping, stable hopping motion on the spot and stable backward hopping, with a locomotion velocity that depends on the choice of K v . Setting K v to 0.4, 0 and − 0.4, respectively, gives the hopping motions denoted A, B and C and are illustrated in Fig. 4 over one complete period of the relative periodic orbit. Panels a1 and a2 for case A show forward locomotion, panel b for case B shows hopping in one place, and panels c1 and c2 for case C show backward hopping locomotion. Here the time history of the tiptoe elevation is plotted together with the stroboscopic view of the model in panels a1, b and c1. The tilt angle of the line between the tiptoe and the CoM at the beginning of the period is in correlation with the direction of locomotion: tilting to the right yields forward locomotion; no tilting yields hopping on the spot; and tilting to the left yields backward locomotion. The paths of the tiptoe and the CoM are shown in panels a2 and c2 for the forward and backward locomotion, respectively. We remark that tiptoe and the CoM paths are qualitatively similar to that of the SLIP model in [15].
Since forward locomotion is the most relevant from a practical point of view, we focus on case example A in the rest of the section. It is illustrated in more detail in Fig. 5. Panel a shows the relative periodic solution in projection onto the subspace (z A , θ r ,ż A ). It starts at the point G2F, which is located at the intersection of the switching surfaces Σ G2F and Σ F2G . After the flight phase, the solution has a jump within Σ F2G . Finally, the ground phase stays on Σ F2G until it reaches Σ G2F again. Different planar projections of the phase space are shown in Fig. 5b. As can be seen, all of the general coordinates are periodic in time, except for x A , which illustrates that the corresponding motion is indeed represented by a relative periodic orbit. Also note that the F2G transition generates a jump in each general coordinate, except in θ r . The reason for the continuous solution for θ r is the lack of its inertial coupling with the other general coordinates. Figure 6 shows, in panels a and b, time evolutions of the kinetic energy T , the potential function U and the total mechanical energy T +U . The constrained motion space kinetic energy (CMSKE) T c is indicated in both cases; this portion of the kinetic energy is absorbed during the F2G transition and then compensated for during the ground phase. tion until the boundary of stability is reached. Figure 7 shows the induced locomotion velocityv as a function of K v , which is approximately a linear function of K v for this particular value of E 0 . Consequently, the locomotion velocity can be predicted by simple linear extrapolation.
The parameter E 0 is responsible for the tuning of the apex height of hopping, and the K v -range of stability can be computed for different values of E 0 . In this way, the stable region in the (K v , E 0 )-plane in Fig. 8 is determined, where parameter points for case examples A, B and C are also indicated. Moreover, the continuation shown in Fig. 7 is indicated in Fig. 8 by the vertical dashed line at E 0 = 1000 J. We have already shown that the locomotion velocity is zero for case example B. The locus of zero locomotion velocity can be computed as well, and it is shown in Fig. 8, where it divides the region of stability in subregions of forward and back- ward hopping. Notice also the narrowing of the stable region at about a nominal energy level of E 0 = 500 J. This happens near the dash-dotted line in Fig. 8 indicating the potential energy level U s = 490.4 J of standing still (ankle and knee angles are α 12 and α 23 , respectively). This energy level is important, because a lower energy level cannot be reached permanently. On the other hand, the controller is able to ensure stable operation for a target energy level of E 0 < U s as well.
In the region highlighted with red color in Fig. 8, the ankle point B of the mechanical model shown in Fig. 2 does not touch the ground and the minimum of the vertical position z B is larger than zero. Hence, the corresponding solutions are physically feasible. As is the case for the entire diagram in Fig. 8, size and shape of this feasibility region depend on the parameters of the model. Note that case examples A and B are physically feasible for the chosen parameter values. Figure 9 shows the properties of the locomotion as a function of E 0 , which were derived from continuation in E 0 from case A of forward locomotion and from case B of hopping on the spot. Panel a shows that the apex height z A max of the hopping motion is not only affected by the nominal energy level E 0 , but also by K v . For E 0 > U s , the apex height z A max can be approximated a b c d Fig. 9 The dependence of some key indicators on the target energy level E 0 . a Apex height of the trajectory with different parameters; b ratio of CMSKE to the total kinetic energy T ; c time period and flight phase time to time period ratio; and d angle range θ r max − θ r min of the rotation wheel as a linear function of E 0 (light blue lines). Moreover, the hopping height for cases A and B, i.e., the elevation of the CoM, is realistic according to [38,39], showing that the generated motion is indeed a good replication of human hopping.
The ratio of CMSKE to the total kinetic energy T is plotted in Fig. 9b as the function of E 0 . It is possible to find the optimal parameter set where the ratio of CMSKE is the smallest. The ratio T c /T can therefore be used as a cost function during the controller design of artificial running and walking systems. We can also deduce that the model can be applied for studying the energetic costs of human locomotion. Figure 9c shows the time period t per and the ratio t F /t per of the flight phase time duration and the full period. This ratio t F /t per is about 70% and 75% for cases A and B, which are close to typical values for human running [24]. Finally, panel d of Fig. 9 shows the range of angular oscillation of the reaction wheel, which plays the role of the torso of the runner. For both cases A and B, the angle θ r max − θ r min is smaller than 12 • , which again agrees with realistic human running.

Biomechanical effect of foot positioning
To demonstrate that model and analysis methods presented in this paper can be used to analyze the effect of biomechanical parameters on the locomotion performance, we focus on the effect of pre-impact tiptoe positioning on locomotion velocity and impact intensity. To this end, the range of investigated control parameters is extended to higher energy levels in order to have a more global view of the applicable parameter set.
The dependence of locomotion velocityv on the tiptoe positioning, characterized by the distance x A − x G , is depicted in panel a of Fig. 10 for five different levels of the energy E 0 . It can be seen that the relation betweenv and the tiptoe positioning is almost linear for each value of E 0 . This observation indicates that footlanding in front of the CoM is beneficial during running, because it results in a higher locomotion velocity. This case, when x A − x G > 0, is referred to as overstriding in the literature [20]. Figure 10b shows the effect of tiptoe positioning on the CMSKE, which is an indicator of ground-foot impact intensity and impacting forces [5,26]. For negative values of x A −x G the kinetic energy ratio T c /T is small (around 1%), while it tends to be higher (up to 4%) for positive values (i.e., when overstriding). The general consequence for athletes is that foot-landing in front of the CoM induces high kinetic energy loss and, therefore, high intensity ground-foot collisions and low energy efficiency. In other words, overstriding is disadvantageous in this respect, which is the opposite outcome from the velocity consideration represented in Fig. 10a. Hence, a runner has to find the balance between energy efficiency and running speed: short-distance sprinters may touch down their tiptoe in front of the CoM in order to maximize speed, while a b Fig. 10 Average locomotion velocityv and CMSKE (as given by T c ) as functions of the horizontal tiptoe position x A − x G relative to the body CoM long distance runners may avoid overstriding in order to reduce energetic costs.

Conclusions and further steps
A planar three-segmented leg model was analyzed, which takes into account the inelastic ground-foot impact-induced energy absorption and an active controller, which maintains a uniform energy level. We showed that there exist relative periodic orbits that corresponds to hopping, and that they are stable over a large and realistic parameter range. This was achieved by means of computing the monodromy matrix with a well-established method for piecewise-smooth dynamical systems. We found that both forward and backward motion is possible and that the hopping height can be tuned. Locomotion velocities in the range (−7, 3) m/s and tiptoe elevations in range (0.01, 4) m were discovered by tuning the main control parameters E 0 and K v only, while the remaining parameters we have chosen to guarantee stable operation in a wide region of E 0 and K v . We further showed that it is possible to tune the parameters to achieve optimal motion, not only in the sense of lowest foot impact intensity (i.e., energy loss), but also in the sense of apex height, stability and robustness. Overall, this demonstrates the usefulness of the model and our analysis approach for obtaining insights into the biomechanics of running.
Future work will concentrate on tuning control parameters to achieve better performance and a larger region of stability of forward hopping. In this way, we intend to demonstrate that the presented methodology and results can be applied also to more complex models of legged locomotion. The longer term goal is to develop the model further toward more humanlike motion; in particular, two legs and a more accurate model for the upper body can be considered. Furthermore, ground contact of multiple points of the foot will be included in the model. Such a further developed model can also be used to investigate the effect of the terrain on the optimal motion, for instance, by considering inclination of the ground.
Although some conclusion regarding human motion can be drawn from mathematically generated trajectories of a model, a comparison with laboratory experiments with human subjects remains essential. After refinement of the proposed mechanical model and the controller, the determination of other performance measures of human running may become feasibleenabling fine-tuning against measured data.

Compliance with ethical standards
Conflict of interest The authors declare that they have no conflict of interest.
Funding This study was funded by MTA-BME Lendület Human Balancing Research Group, Hungary.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.