Influence of “J”-Curve Spring Stiffness on Running Speeds of Segmented Legs during High-Speed Locomotion

Both the linear leg spring model and the two-segment leg model with constant spring stiffness have been broadly used as template models to investigate bouncing gaits for legged robots with compliant legs. In addition to these two models, the other stiffness leg spring models developed using inspiration from biological characteristic have the potential to improve high-speed running capacity of spring-legged robots. In this paper, we investigate the effects of “J”-curve spring stiffness inspired by biological materials on running speeds of segmented legs during high-speed locomotion. Mathematical formulation of the relationship between the virtual leg force and the virtual leg compression is established. When the SLIP model and the two-segment leg model with constant spring stiffness and with “J”-curve spring stiffness have the same dimensionless reference stiffness, the two-segment leg model with “J”-curve spring stiffness reveals that (1) both the largest tolerated range of running speeds and the tolerated maximum running speed are found and (2) at fast running speed from 25 to 40/92 m s−1 both the tolerated range of landing angle and the stability region are the largest. It is suggested that the two-segment leg model with “J”-curve spring stiffness is more advantageous for high-speed running compared with the SLIP model and with constant spring stiffness.


Introduction
Owing to the elastic elements (muscles, tendons, ligaments, and other soft tissues) of legged systems, in fast animal locomotion spring-like leg behavior is discovered to represent bouncing gaits like running, hopping, and trotting [1,2]. On the one hand, the leg compliance can reduce the impact of the ground contact phase [3] and recycle kinetic energy by using of elastic strain energy storage and release [4], which can lead to low energy costs. For instance, due to storage and utilization of elastic strain energy, the reduction rate of the required muscle work is up to 40 percent in a trotting or galloping horse [5]. On the other hand, compliant legged systems exhibiting self-stability in response to internal (speed variations) and external perturbations (ground surface irregularities) can simplify the dynamic control of bouncing motion [6]. Therefore, the elastic elements play an important role in fast animal locomotion [7].
According to spring-like leg behavior of bouncy gaits, biomechanists [8] propose a simple spring-mass model to describe animal and human locomotion, which is also called spring-loaded inverted pendulum (SLIP), consisting of a point mass representing the body and a massless linear spring describing the leg. For the purpose of achieving stable running in the SLIP model, an increase in leg stiffness is required for increased running velocities when angle of attack is unchanged [9]. At present, the SLIP model is extensively applied to the study of bouncing gaits for legged robots, including one-legged hopping robots, bipedal robots, quadruped robots, and hexapod robots. Sato and Buehler propose a hopping robot with one leg based on the SLIP model [10].
However, biological limbs are not the telescopic linear leg model, rather, they are made up of multiple joints; their compliance is situated at the joint level [11,12]. According to this concept, Rummel and Seyfarth [6] study the effects of the two-segment leg model with constant spring stiffness on running stability during low-speed running, and results from simulations and experiments show that adjustment of joint stiffness is required to support stable running at different speeds. Furthermore, this segmented leg supports self-stable running at an enlarged range of speeds (lower minimum speed) compared with the SLIP model. At present, the twosegment leg model with constant spring stiffness is widely applied to the design of compliant legged robots such as the intersegmental joint configuration of "JenaHopper" designed by Rummel et al. [13].
It is well known that the SLIP model and the two-segment leg model with constant spring stiffness can be seen as an effective tool to research on bouncing gaits for legged robots. In addition to these two models, Karssen and Wisse [14] also propose the running model with nonlinear leg spring and study the effect of the nonlinear leg springs on disturbance rejection behavior. Some results of this optimization reveal that the push (push forward and backward) disturbances rejection with the optimal nonlinear leg spring is much better than with the optimal linear spring. What is more, it can be seen from the size of the basin of attraction that the range of running speeds of the optimal nonliner spring is largely increased compared with the optimal linear spring. But this running model is analyzed in the low-speed case. At present, a number of prominent high-speed spring-legged robots, such as bipedal robot MABEL [15] and Cheetah-cub quadruped robot [16], are developed. Here, it is interesting to note that a speed record of 3.06 m s −1 is implemented on MABEL by adjusting its effective leg stiffness [17].
Although there are many running robots and running models at present, we focus only on the two-segment leg model in this study. This is because it is the most reduced leg configuration; in spite of the low complexity of this twosegment leg model, it is still suitable to solve the question of how leg segmentation and joint stiffness influence the stability of running at different speeds [6]. On the other hand, so far, the potential effects of joint stiffness on running speeds during high-speed locomotion remain an unresolved issue. Therefore, in the present paper, we propose the twosegment leg model with "J"-curve spring stiffness. Here, "J"curve spring stiffness is inspired by biological materials [18]. Subsequently, we employ this proposed model to investigate the potential role of "J"-curve spring stiffness on running speeds during high-speed locomotion. Here, in this paper running models focus only on the SLIP model, the twosegment leg model with constant spring stiffness and with "J"-curve spring stiffness. This is to develop a deeper understanding of the benefits and drawbacks of the two-segment leg model with "J"-curve spring stiffness compared with the other two models regarding high-speed running capacity.
In this study, we not only hope that the two-segment leg model with "J"-curve spring stiffness will show the largest range of running speed for self-stable high-speed running in all three models but also expect that results of our work will be regarded as a promising concept for the design of bioinspired high-speed robots. Figure 1(a), the SLIP model is modeled as a point mass attached to a massless spring leg with linear stiffness of and rest length of 0 . This model for running can be represented as the flight and stance phases alternatively, which is shown in Figure 1(b). During the flight phase the system dynamics is determined by the point-mass gravity, which results in a ballistic trajectory of the point mass, and then the equation of motion can be expressed as

SLIP Model. As shown in
where F flight is the total force vector during flight and g = [0, − ] T denotes the gravitational acceleration vector. In addition, considering that the total energy total is assumed to be conserved, during flight total is total = where V 0 is the constant horizontal speed of the point mass and apex, denotes the apex height (the point mass has the maximum vertical height at the beginning of the flight phase of running). Afterwards, the transition from flight to stance phase occurs when the spring leg strikes the ground at a given angle of attack 0 . This transition event can be formulated by where ( ), which is the function of time , denotes the current spring leg length. The next moment, the stance phase starts. With the tip of foot regarded as pivot without slipping, the stance phase can be split into the compression and decompression subphases. Here, the transition between the above-described two subphases takes place when the leg reaches its maximum compression. Note that leg stiffness is a fixed value of and the direction of leg force is from the tip of foot to the point mass during these two subphases. Thus, the motion equations of stance phase are written as where F leg is the leg force vector and F stance denotes the total force vector. Finally, when the tip of foot leaves the ground the flight phase starts, and then the point mass reaches the subsequent apex variable apex, +1 . Consequently, a step, called a cycle, can be defined as the movement between apex, and apex, +1 .

Two-Segment Leg Model.
The following section describes the configuration of the two-segment leg model and its dynamics of running with spring-like legs. As illustrated in Figure 2, the two-segment leg model is described by a point mass attached to a rotating segmented leg, and this segmented leg is represented by massless upper and lower leg linked by the intersegmental joint, with the joint angle of 0 at rest and the radius of cable pulley of . What is more, this model can be considered as hip movement actively and knee movement passively. Note that knee joint elasticity originates from spring compliance, and spring property has a significant influence on running stability. Thus, the two-segment leg model can be divided into the two-segment leg model with constant spring stiffness and with nonlinear spring stiffness in terms of the linear or nonlinear characteristics of tensile spring, and its joint stiffness (Δ ) is where denotes instantaneous joint angle, Δ represents the amount of joint flexion, and spring is the tensile force of the spring. In addition, it is evident that both upper leg of length 1 and lower leg of length 2 affect the dynamics of running. Therefore, to make the analysis of the mathematical model easily, the two segment lengths are defined as 1 = 2 = ; in order to facilitate the comparison between the SLIP model and the two-segment leg model, rest length of the virtual leg is also defined as 0 . Figure 3 illustrates the two-segment leg model during a running period. It is worth noting that the two-segment leg model can be conceived as an equivalent SLIP model whose leg stiffness is nonlinear. Thus, similar to the SLIP model, the two-segment leg model for running is also composed of flight and stance phases, and stance phase can also be divided into the compression and decompression subphases. Similarly, the total energy, transition conditions, the total force during flight, and the direction of leg force of the stance phase are identical with those of the SLIP model, respectively.

Two-Segment Leg Model with Constant Spring Stiffness.
As can be seen in Figures 2 and 3, if tensile spring has a constant stiffness, the system can be regarded as the twosegment leg model with constant spring stiffness reported in [6], and the equation of motion for this running system is given by the following equation (based on [6]):

Two-Segment Leg Model with "J"-Curve Spring Stiffness.
In this section we illustrate a nonlinear "J"-curve spring force-elongation relationship of the proposed model and its dynamics of running. A schematic diagram of a joint of large mammals, presented in [19], describes a knee or an ankle; joint configuration shown in this schematic diagram is similar to the two-segment leg model. Here, owing to the similar role of elastic elements, muscles and tendons shown in this schematic diagram can be considered as tensile spring of the two-segment leg model. it is very interesting to note that muscles [20], tendons [21], and muscle-tendon complex [22] force-elongation curves resemble a "J"-shape, and Steven Vogle also notices that biological tissues (ligaments, skin, etc.) force-length curves resemble a "J"-shape [23]. Thereby, we imitate a joint configuration of large mammals and adopt "J"curve spring stiffness inspired by biological materials [18] to establish the two-segment leg model with "J"-curve spring   stiffness. By doing this, we hope that the proposed twosegment leg model with "J"-curve spring stiffness will be capable of realizing high-speed locomotion, resulting in the imitation of biological high-speed running capacity. Here, five points shown in Table 1, presented in [18], can be used to describe "J"-curve spring properties, and then mathematical formulation of the relationship between the "J"-curve spring force ( ) and the "J"-curve spring elongation is obtained by using of the interpolation method, and this equation is represented by where 4 , 3 , 2 , 1 , and 0 denote the undetermined coefficients; subsequently, substituting these five points into (7) can yield where max represents the maximum elongation of "J"-curve spring and max denotes the "J"-curve spring force corresponding to max .
The virtual leg force vector F leg (Δl) relates gravitational force vector g to the total force vector F stance during the stance phase with in [6] the virtual leg force leg (Δ ) is a function of the virtual leg compression Δ with where the joint torque can be determined according to ( ) with = ( ) ; (11) the relationship between the amount of the joint flexion Δ and can be given by Δ can be then represented by Seeing that the maximum compression of the virtual leg rarely exceeds 30% of the rest virtual leg length among running animals [24], both the spring's maximum compression in the SLIP model and the maximum compression of the virtual leg in the two-segment leg model are defined as 0.3 0 in this study, and then the corresponding maximum amount of the joint flexion Δ 30 in the two-segment leg model is formulated by Finally, the radius of cable pulley can be obtained based on the amount of the joint flexion Δ 30 with Applied Bionics and Biomechanics 5

Analysis Methods.
The steps-to-fall analysis and the apex return map, which are reported in [6,9], are adopted to analyze the running dynamics as represented by the abovementioned three models. In the following section, we briefly review these two methods. The steps-to-fall analysis can record the maximum number of steps to fall for a given system parameters ( , 0 , V 0 , etc.), and simulation calculation is not stopped until the number of predefined steps is reached.
Here, the maximum number of successive steps is defined as 50 [6]. Although the steps-to-fall analysis can provide a method for counting the number of successive steps, the system might still fall after a finite threshold. Therefore, in order to solve this question, the apex return map is adopted. The second approach, the apex return map, can identify the fixed point. Furthermore, the stability of the fixed point * = apex, +1 = apex, can be described by here if is smaller than 1, this condition means that the system is stable.

Simulation Parameters Setup.
For the purpose of facilitating the comparison of all three models, simulation parameters can be defined as follows. (1) Seeing that a representative value for humans running is leg compression at 10% of rest leg length [25], the reference stiffness 10% can be defined by the following equation (based on [6]): where 10% denotes a reference leg compression at 10% of rest leg length and 10% is the corresponding leg force; (2) a dimensionless reference stiffness̃can be expressed by the following equation (based on [26]): (3) the dimensionless reference stiffness̃1 0% can be given by the following equation (based on [6]): three models parameters, = 80kg, 0 = 1m, 30 ∘ ≤ 0 ≤ 90 ∘ , and 0 ≤̃1 0% ≤ 50, and the initial apex height of 1 m are defined [6]. Additionally, in the two-segment leg models the small ( 0 = 115 ∘ ) and large nominal joint angles ( 0 = 165 ∘ ) are adopted to analyze the effects of different nominal joint angles on running speeds, respectively; in our model ( 0 = 115 ∘ and 165 ∘ ), the maximum speed of stable running is 92 and 40 m s −1 , respectively. Thus, two speed ranges (5 to 92 m s −1 and 5 to 40 m s −1 ) are utilized to investigate the advantages and disadvantages of our model compared with the other two models with respect to highspeed running capacity, respectively. Here, the speed of 5 m s −1 is the minimum speed of stable running in our model at 0 = 115 ∘ .

Normalized Force-Length Relationships of the Two-Segment Leg Model with "J"-Curve Spring Stiffness.
In this section, we not only analyze normalized force-compression relationships of the two-segment leg model with "J"-curve spring stiffness but also investigate the effects of different nominal joint angles 0 on these relationships. For ease of understanding, we divide the range of Δ / 0 into two subintervals. Here, the first one is 0 ⩽ Δ / 0 < 10, and the second is 10 ⩽ Δ / 0 ⩽ 30. For 0 ⩽ Δ / 0 < 10, results shown in Figure 4 reveal that large variations in Δ / 0 result in small changes in leg force leg / 10% . In contrast, for 10 ⩽ Δ / 0 ⩽ 30, large changes in Δ / 0 lead to large variations in leg force. Furthermore, at the beginning of the second subinterval, a given variations in Δ / 0 lead to increasingly large changes in leg force with increasing virtual leg compression. Next, we consider roles of nominal joint angles. The larger the nominal joint angle, the slightly smaller the leg force for a given Δ / 0 in the first subinterval. In the second subinterval, the smaller the nominal joint angle, the greater the rise in leg stiffness at a given virtual leg compression.

Comparison of Normalized Force-Length Relationships between the Two-Segment Leg Model with "J"-Curve Spring
Stiffness and the Other Two Models. Considering that in [6] comparison of normalized force-length relationships between the SLIP model and the two-segment leg model with constant spring stiffness has already been discussed, this section focuses on comparison of relationships between the two-segment leg model with "J"-curve spring stiffness and the other two models. In the first subinterval, leg force of the proposed model is smaller than those of the other two models for a given Δ / 0 , and leg stiffness of our model becomes increasingly small with increasing 0 . In the second subinterval, for a given Δ / 0 our model has the maximum leg force in all three models, and leg force becomes increasingly large with increasing 0 in our model.

Stability Analysis at Different Running Speeds V 0 and
Nominal Joint Angles 0 . In this section, we concentrate on the effects of a running speed and the dimensionless reference stiffness on each other. Again, to gain better insights into advantages and disadvantages of our model compared with the other two models regarding high-speed running capacity, we also analyze the effects of these speeds on angle of attack and the stability region in all three models, respectively. Here, for ease of understanding, some typical examples, which are the regions of stable running at different running velocities (V 0 = 7 m s −1 , V 0 = 29 m s −1 , and V 0 = 36 m s −1 ) and nominal joint angles ( 0 = 115 ∘ and 0 = 165 ∘ ), are used to analyze self-stabilizing behavior of the above-represented three models according to the normalized force-compression relationships depicted in Figure 4. Interestingly, the speed of 29 m s −1 , which is the maximum running speed of the cheetah, is the highest running speed recorded from land animal [27]. In all three models, the proposed model has the tolerated maximum range of the dimensionless reference stiffness at a running speed from 7 to 40/92 m s −1 ; the smaller the nominal joint angles, the larger the stability regions for given combinations of V 0 and̃1 0% , as shown in Figure 5(a). In other words, we can see from Figure 5(a) that for a given dimensionless reference stiffness the proposed model exhibits not only the largest tolerated range of running speeds but also the tolerated maximum running speed in all three models; for example, for a giveñ1 0% = 41.2698, our model at 0 = 115 ∘ is capable of accomplishing stable running behavior at a speed range from 5 to 76 m s −1 . However, at the same dimensionless reference stiffness, the tolerated speed range is from 5 to 27 m s −1 in the SLIP model; worse still, in the two-segment leg model with constant spring stiffness this range is only from 5 to 19 m s −1 . In addition, a tolerated minimum̃1 0% is required in order to guarantee the stability of running in a region (one of the twelve regions in Figure 6). Interestingly, although the values of thesẽ1 0% become increasingly large with increasing running speeds, this tolerated minimum̃1 0% in our model is smaller than those of the other two models for a given speed from 7 to 40/92 m s −1 ; the smaller the nominal joint angle, the smaller the rate of increment in this tolerated minimum̃1 0% .
At fast running speed form 25 to 40/92 m s −1 , our model has the maximum range in angle of attack ( 1 and 2 ). For instance, at high running speed (29 m s −1 ), 1 and 2 (see Figure 5(b)) are 11.43 and 7.62 ∘ in our model at 0 = 115 ∘ , respectively. In contrast, 1 = 2 = 3.8 ∘ in the SLIP model and 1 = 2 = 0 ∘ in the two-segment leg model with constant spring stiffness are found. Additionally, the range of landing angle becomes smaller with increased 0 ; for example, in our model at V 0 = 40 m s −1 and 0 = 115 ∘ , 1 = 9.52 ∘ and 2 = 7.62 ∘ are found. However, at the same speed, both 1 and 2 at 0 = 165 ∘ are only 0.95 ∘ . As for landing angle sensitivity, all three models are sensitive for angle of attack variations because of 1 = 2 , but our model is sensitive for angle of attack variations at higher speeds. Again, when our model is sensitive for landing angel variations, it is necessary to stable running at higher velocities for decreased 0 .
At fast running speed form 17 to 40/92 m s −1 , regions of stable running of the proposed model are larger than those of the other two models, as illustrated in Figure 5(c). For instance, a region representing the number of successful running steps, shown in Figure 6, is made up of the 64 × 64 equidistant grid; for a given running speed of 29 m s −1 the stability region consists of the 247 equidistant grids in our model ( 0 = 115 ∘ ). However, at the same running speed the region of stable running of the SLIP model is only the 16; worse still, there is no stable region in the two-segment leg model with constant spring stiffness. In contrast, at low running speed from 5 to 10 m s −1 , our model demonstrates the minimum stability region in all three models. Again, in the two-segment leg model, an increase in nominal joint angle leads to a decrease of the stability region during high-speed locomotion. For instance, at high running speed (36 m s −1 ), the stability region of our model is decreased from the 205 at 0 = 115 ∘ to the 10 equidistant grids at 0 = 165 ∘ .

Return Map of the Apex Height.
In the following section we analyze stability of the proposed model by using of this single apex return map. For a given total energy = 34424.8 J and system parameters ( 0 = 50 ∘ and̃1 0% = 18.75), results from the effects of different nominal joint angles on stable fixed points (the intersections of between three curves and the diagonal) are shown in Figure 7 Figure 7(b) shows the effects of another factor-different running velocities-on stable fixed points with system parameters ( 0 = 47 ∘ ,̃1 0% = 23.8, and 0 = 115 ∘ ).
We can obtain that the lower the running speeds, the higher the values of fixed points, and the high running velocities (31 m s −1 ) result in the small basin of attraction containing all apex heights.

Discussion
In this paper, we discuss the effects of the two-segment leg with "J"-curve spring stiffness on running speeds during high-speed running. Two methods, the steps-to-fall analysis and the apex return map [6,9], can be adopted to exploit high-speed running capacity of the model.
Compared with the other two models, during fast running it reveals that (1) system can provide the larger regions of stable running; (2) the tolerated range of̃1 0% for selfstable running is even larger for a given running speed; (3) the proposed model shows the larger tolerated speed range and running speeds. In addition, when the proposed model is stable running, the small nominal joint angle can lead to the low tolerated minimum̃1 0% , with the large tolerated speed range and running speed.  Figure 5: Properties of regions of stable running for given combinations of̃1 0% and 0 at different running speeds in the SLIP model, the two-segment leg model with constant spring stiffness and with "J"-curve spring stiffness, respectively. Here, the speed range of stable running for a giveñ1 0% from 0 to 50 is shown in (a) and (b), 1 denotes the difference between the minimum and maximum landing angles, and 2 is the maximum tolerated range of landing angle on a region of stable running.   These characteristics mentioned above are due to the two-segment leg structural configuration and "J"-curve spring stiffness properties, resulting in the nonlinear forcecompression relationships depicted in Figure 4. In addition, owing to the passive compliance identified in the proposed model, this elastic two-segment leg configuration can take advantage of simple control strategies to guarantee the steadystate running behavior, with little or even no sensory information. We believe that our model can be seen as a template to analyze the high-speed dynamic locomotion for animals and robots. However, this two-segment leg model with a "J"-shape force-elongation curve cannot mimic adjustment of leg stiffness in fast animal locomotion. This is because limb stiffness is adapted to running speed; for example, Hobara et al. [28] reveal through experiments that in humans leg stiffness is adjusted to different hopping frequencies, with the fact that there are two different force-elongation curves during compression and decompression. Hence, in order to adequately mimic limb compliant locomotion, these aspects should be taken into account in the design and implementation of robotic system.
Additionally, in the two-segment leg model with "J"curve spring stiffness, the running is simulated across smooth and level terrain. Yet, the terrain is not the same in the real world, where the ground surface irregularities must be taken into account for running robots and running models. Therefore, in order to achieve stable running of the proposed model on uneven terrain, it is necessary to adjust apex height with adequate ground clearance in response to disturbances in ground height. Currently, for a given total system energy, the return map of the apex height apex, +1 ( apex, ) can be adopted to analyze running stability when the system is perturbed by a change in ground height. Furthermore, as the proposed model is conservative, the total system energy is distributed to the vertical energy and the forward kinetic energy by the adjustment of angle of attack; in other words,  different angles of attack can be used to control a desired apex height or, correspondingly, the forward speed. Here, the return map of the apex height is adopted to investigate the potential effects of different landing angles on the apex height in our model (V 0 = 29 m s −1 ,̃1 0% = 40, and 0 = 115 ∘ ). Results shown in Figure 8 reveal that the apex height apex, +1 is dependent on the preceding apex height apex, and the selected angle of attack 0 . What is more, for a given apex height apex, , both the value of fixed point and the subsequent apex height apex, +1 become increasingly small with increasing 0 , and the apex height is sensitive for angle of attack variations. This means that a higher apex height can be implemented by tuning a smaller angle of attack, with the fact that only one step or a few more steps are required to achieve higher steady-state heights in our model. Again, it is interesting to note that Seyfarth and Geyer [29] introduce a generalized return map and derive an optimal control strategy of the apex height in the SLIP model. This method could be utilized to study the control strategy of the apex height of the proposed model, which will be subject of further investigations.
Running at High Speed. In Figure 5(a), the SLIP model indicates that the tolerated minimum leg stiffness has to increase with running speed to achieve stable running. This is in agreement with a simulation study in which increasing leg stiffness is required to guarantee the stability of a galloping quadruped whose single leg consists of a linear spring attached to a prismatic leg when running speed is increasingly high [30]. In the two-segment leg model with constant spring stiffness, an increase in the tolerated minimum leg stiffness is also necessary to accomplish stable running behavior at higher speeds [6]. This finding is supported by robotic trials. Here, hopping robot with one leg, which is made up of a segmented leg with constant spring stiffness, becomes unstable with increasing speed [31]. Therefore, leg stiffness is sensitive for speed variation in these two models. The situation is not all the same in the two-segment leg model with "J"curve spring stiffness. Here, in our model the amount of increment in the tolerated minimum dimensionless reference stiffness, for the same amount of increment in running speed, is largely decreased compared with the other two models. This characteristic of our model (leg stiffness is insensitive for speed variation) enlarges tolerated range of speeds. On the other hand, for a given running speed, the proposed model has the minimum value of leg stiffness to achieve selfstable running in all three models, as depicted in Figure 5(a). For instance, at high running speed (29 m s −1 ), the tolerated minimum dimensionless reference stiffness in our model at 0 = 115 ∘ is 18.25 compared with 44.44 in linear leg model and no stability solutions in the two-segment leg model with constant spring stiffness. As a result, in all three models, for a given nominal joint angle and leg stiffness, the proposed model at high speed has not only the maximum speed range of stable locomotion but also the tolerated maximum running speed. This means that in all three models the twosegment leg model with "J"-curve spring stiffness is most advantageous for high-speed running.
Practically, several researchers have designed their robotic legs similar to our model configuration and leg stiffness. For example, the robotic leg proposed by Schmiedeler and Waldron [32] used seven parameters to accomplish biomimetic feature. With these parameters, we here in this paper have calculated their virtual leg force-compression relationship, as shown in Figure 9. Obviously, this curve is similar to the curves as represented by our model depicted in Figure 4, leading to the following common characteristics: (1) the virtual leg spring is capable of being soft at the beginning of touch down, which can produce a relatively small impact force [32]; (2) it becomes stiffer with increased compression, resulting in a large force, which can prevent the leg to fold over on itself [32]. With the facts that their prototype leg exhibits the potential for fast locomotion in experiment [32] and then the corresponding KOLT quadruped robot can implement gallop and highspeed turning [33], we conclude that our simulation results have been confirmed in practice.

Conclusion
In this paper, we have presented the two-segment leg model with "J"-curve spring stiffness and have analyzed the effects of "J"-curve spring stiffness on running speeds of the proposed model during high-speed running. According to simulation results of all three models, for a given dimensionless reference stiffness, we have demonstrated that (1) our model has not only the largest tolerated speed range but also the tolerated maximum running speed, with the fact that the smaller the nominal joint angle, the better the running capacity mentioned above; (2) at fast running speed from 25 to 40/92 m s −1 our model has not only the largest range of angle of attack but the largest region of stable running. It has already been successfully applied to quadruped robot, such as KOLT, for high-speed running. Next, we will investigate the effects of a broad ratio of leg length and nominal joint angle range on high-speed running performance, respectively. Again, we will analyze the whole stability zone of segmented leg in order to gain better insights into advantages and disadvantages of our model at fast running.