In this paper, we examine the question of whether parallel elastic actuation or series elastic actuation is better suited for hopping robots. To this end, we compare and contrast the two actuation concepts in energy optimal hopping motions. To enable a fair comparison, we employ optimal control to identify motion trajectories, actuator inputs, and system parameters that are optimally suited for each actuator concept. In other words, we compare the best possible hopper with parallel elastic actuation to the best possible hopper with series elastic actuation. The optimizations are conducted for three different cost functions: positive mechanical motor work, thermal electrical losses, and positive electrical work. Furthermore, we look at three representative cases for converting rotary motor motion to linear leg motion in a legged robot. Our model featured an electric DC-motor model, a gearbox with friction, damping in the leg spring, and contact collisions. We find that the optimal actuator choice depends both on the cost function and conversion of motor motion to leg motion. When considering only thermal electrical losses, parallel elastic actuation always performs better. In terms of positive mechanical motor work and positive electrical work, series elastic actuation is better when there is little friction in the gear-train. For higher gear-train friction parallel elastic actuation is more economical for these cost functions as well.

## Introduction

Springs play a fundamental role in legged locomotion. In nature, elastic elements are used for energy storage, as return springs, and to cushion impacts [1]. Model-based analyses have shown that compliant legs can explain the dynamics of human walking and running [2], as well as a wide variety of quadrupedal gaits, including walking, trotting, tölting, and galloping [3,4]. In all these cases, elastic energy storage enables the recycling of energy and improves energetic economy. Cavagna et al. [5] and Blickhan [6] observed that during running gaits, animals conserve energy by having the body undergo an elastic bouncing motion. In human in-place hopping, for example, the majority of the energetic fluctuations are generated passively through elastic energy storage in muscles and tendons [7,8]. Motivated by these biological benefits, elastic elements have been incorporated successfully in robotic simulations, e.g., see Refs. [9] and [10], and in hardware prototypes, e.g., see Refs. [11–15].

In robotic hardware, the two primary ways of incorporating elasticity into legs are to place springs in parallel or in series with the actuation source (Fig. 1). In a parallel elastic actuator (PEA), the motor contracts and extends the entire leg and spring. The spring force and motor force, therefore, act additively. The motor inertia moves with the joint and is not isolated from impacts. In a series elastic actuator (SEA), the motor moves the proximal end of the spring. The motor force must overcome the spring force. The motor inertia does not directly add to the joint motion and is isolated from impacts. Both types of actuator have been implemented in prototypes, e.g., see Refs. [16] and [17]. Yet, to date, there is still no detailed understanding of the effect of each actuator type on the energetics and motion characteristics of legged systems. In particular, there is still disagreement about which actuator solution is more energetically economical.

There is a broad range of similar questions that essentially try to determine the benefits of discrete design choices. Can a compliant ankle joint in a robot's leg improve efficiency? Is a quadrupedal robot faster with a spine that is rigid or deformable? Is a set of swinging arms useful for robotic locomotion? These questions are fundamentally difficult to answer. One cannot simply swap a particular feature in a robotic implementation or robotic model and then compare the performances of the two different variations. This is because an optimal robot with one of the two variations will likely have to look and behave very differently than an optimal robot with the other. For example, a robot with PEA might require a different spring stiffness and a different transmission ratio than a robot with SEA, and the two actuator types might need very different ways to power their motion. If we keep the system parameters of the two robots identical or apply the same type of motion in both cases, we will implicitly bias the comparison to favor one variation over the other. After all, we are not interested in the question of whether a *certain* robot with a certain motion strategy and a certain choice of parameters performs better with one design variation or the other. We would like to know if the *best possible* robot, with the best possible set of parameters, and the best possible motion strategy can benefit from the design variation. To answer this question, we have to investigate optimal motion, optimal parameters, and optimal morphology at the same time. Due to the complexity of the problem, it is impossible to do this comparison by conducting an exhaustive search or an analytical evaluation. Instead, we propose to use trajectory optimization to find the best possible actuator inputs and motion trajectories while simultaneously optimizing the system parameters.

This approach extends upon existing morphological optimizations in legged robotics [18–21]. These past studies optimize parameters to optimally adapt a model to a particular task but do not use the technique to compare discrete designs. It also extends on existing discrete structural comparisons [22] by using trajectory optimization on a realistic motor model and including morphological parameters directly in the optimization. In particular, this technique has not been applied in previous studies comparing PEA and SEA. An initial comparison was conducted, for example, by James et al. [23], and found that PEA was more economical. Grimmer et al. [24], on the other hand, found SEA to be more economical for a model of a human ankle joint. The difference in the two studies likely stems from their comparisons of *particular* implementations of each actuator type, rather than the *best* implementation. Furthermore, the studies significantly simplified the problem. James et al. [23] only looks at a single cost function and a single set of parameters for actuator comparisons. Grimmer et al. [24] has similar simplifications, including ignoring the motor inertia and damping in the springs.

In this work, we performed a thorough comparison of PEA and SEA for the exemplary case study of in-place hopping. In a trajectory optimization framework, we considered inputs, motion trajectories, and system parameters simultaneously to understand their mutual effects and their full implications on the choice of PEA or SEA. As a basis for our study, we used a model that, while simple enough to lead to general conclusions, has enough detail to model more realistic motions during hopping, such as leg dynamics during flight. In particular, our models incorporated feet with mass, detailed electric DC motor models, damping in the springs, and gearbox friction. In addition, we enforced realistic constraints on the possible parameter values, motion trajectories, and motor inputs. We considered two different cases of rotary-to-linear transmissions that represent common values in modern legged robotics. We contrasted these two types with the more theoretical scenario of a completely frictionless transmission. Within each of these three cases, we examined a variety of cost functions quantifying work- and force-based efforts. We found that the optimal actuator type and motion was highly dependent on both the type of rotary to linear transmission as well as on the chosen cost function.

As hopping is highly relevant in legged locomotion, which often uses the template of a spring loaded inverted pendulum for bipedal [6] and multilegged [25] locomotion, our results are relevant for both a fundamental understanding of the effect of elastic actuators on legged motion as well as for applied hardware design. In addition to answering the question of which actuator type is better in different situations, our study provides a detailed understanding of what optimal motions look like and can, therefore, be used to “calibrate” trajectory optimization results in real hardware. To this end, the paper represents a culmination of our previous efforts analyzing optimal actuation [26–28].

## Theoretical Discussion

In the past, discussions about which actuation concept is better suited for legged robotics often argued that PEA is better when minimizing motor torque and SEA is better when minimizing mechanical work, e.g., see Ref. [29]. To illustrate this argument and to use it as a departure point for our own work, we can formalize this distinction in terms of two different cost functions: the positive mechanical motor work *C*_{mech} (where the work done by the motor is equal to *T _{m}ω_{m}*) and the thermal electrical losses

*C*

_{therm}(which is proportional to $Tm2$). In these expressions,

*T*is the torque created by the motor, and

_{m}*ω*is the motor velocity. Since

_{m}*C*

_{therm}is only dependent on

*T*, one could think that this cost function would indeed prefer PEA, where motor force and spring force act additively and the motor torque is thus reduced. SEA would be at a disadvantage, as its motor must actively provide the entire force acting on the leg. For

_{m}*C*

_{mech}, the situation is reversed. For SEA, the motor and main body components are well insulated from impacts, and the inertia of the motor does not add directly to the joint motion, potentially leading to reduced mechanical losses. PEA seems to be at a disadvantage, as the rotor inertia now factors into collisions.

*j*that is connected to the joint via a frictionless gearbox with transmission ratio

_{m}*n*. For both the PEA and SEA version of this hopper, we can, in fact, find entirely lossless solutions for

*both*cost functions. PEA can achieve

*C*

_{mech}= 0 by having the motor input no torque (

*T*= 0) and, therefore, having the spring perform the hopping task passively. With this motion strategy, the only possible source of mechanical losses is collision losses in the rotor inertia

_{m}where *m*_{1} is the main body mass and *v*_{foot} represents the velocity of the foot at touchdown [26]. The reflected rotor inertia of the motor *j _{o}* scales with the square of the transmission ratio

*n*: where we assume that the transmission ratio converts rotary to linear motion, giving

*j*units of mass. By setting

_{o}*n*= 0 and thus

*j*= 0, the collision losses can be completely removed, which leads to

_{o}*C*

_{mech}= 0. There are no disadvantages to this choice, since no torque is required from the motor. This strategy also achieves

*C*

_{therm}= 0, as

*T*= 0. In contrast, SEA can yield

_{m}*C*

_{mech}= 0 by setting

*ω*= 0. This fixes the location of the proximal end of the spring, letting the spring perform the hopping task entirely passively. Since in SEA, the rotor inertia is decoupled from the motion, no losses occur and the motion is periodic. This solution, however, requires a nonzero force to keep the motor in place as the leg spring compresses. In order to achieve

_{m}*C*

_{therm}= 0, SEA must, thus, set the transmission

*n*to infinity. Since

*T*=

_{o}*nT*, where

_{m}*T*is the output torque after the gearbox,

_{o}*C*

_{therm}is proportional to $To2/n2$. Using

*n*=

*∞*yields

*C*

_{therm}= 0. Again, there are no disadvantages to this choice, since no motion is required from the motor.

Even though the model in this analysis is a contrived example, it shows two important issues. First, with the right choice of motion and parameters, PEA and SEA can achieve the same performance for *both* work-based and torque-based cost functions. No actuator type is better per se. Second, the two actuator types require almost diametrically different strategies and parameters. PEA applies no force as the motor moves with the spring and utilizes a transmission ratio of *n* = 0. Choosing the opposite extreme of *n* = *∞* drives *C*_{therm} to infinity, as the collision losses in the motor become infinite. SEA keeps the motor at rest and requires a transmission ratio of *n* = *∞*. If SEA were to choose the opposite extreme of *n* = 0 then *C*_{therm} would be driven to infinity, as it would take an infinite *T _{m}* to keep the motor still. It is clear that when trying to answer the question of which actuator type is better, we need to take into account that each actuator requires substantially different motion strategies and parameters.

The extreme values for *n* in this example are a consequence of the assumption that the hopping is otherwise lossless. As soon as we introduce damping and collision losses, these transmission values become extremely nonoptimal. When the motors must move and apply torques to do positive mechanical work to replace the inevitable energetic losses, finite nonzero transmission values *n* are necessary. In particular, SEA can no longer have an infinite *n*, as the resulting infinite reflected rotor inertia would require an infinite *T _{m}* to move the motor to replenish losses, which will drive both cost functions to infinity. PEA can also no longer have

*n*= 0, as it will take an infinite

*T*to transmit force through the gearbox and add energy to the system. That infinite

_{m}*T*will again drive both cost functions to infinity. The infinite costs at both extremes of

_{m}*n*suggest that there exist optimal choices for

*n*in between. For a simplified SEA hopper, for example, such an optimal value has been derived analytically in Ref. [28]. Similar types of dependencies likely exist for other parameters and must be resolved by including the parameters in an optimization formulation.

*T*≈ 0), which minimizes

_{m}*T*and thus

_{m}*C*

_{therm}. Moreover, this oscillation can be timed such that the foot has zero velocity at touchdown, leading to no collision losses even in the presence of a foot mass. This timing is achieved by tuning the natural frequency

*f*of the leg spring

*k*, transmission ratio

*n*, and foot mass

*m*

_{2}. That is, through this natural dynamic oscillation, the choice of motion and parameters is coupled. Such an exploitation of natural dynamics can happen in a variety of ways and must, thus, be accounted for when trying to answer the question of which actuator type is more economical.

The main takeaway of these considerations is that motion and optimal morphology are strongly coupled and will differ greatly between hoppers with PEA and SEA. In order to figure out which actuator type is better, we have to investigate optimal motion and optimal morphology at the same time, such that we can compare the best possible SEA robot to the best possible PEA robot. In addition to the examples mentioned earlier, this question is even further complicated by the presence of friction in the gearbox and limits on the motor force and joint motion. With torque inputs, motion trajectories, system parameters, and limits mutually affecting each other in such complex ways, numerical optimization is the only suitable tool to understand the full implications of the choice of PEA or SEA.

## Methods

In this paper, we used trajectory optimization to find the most energetically economical motions and parameters for each actuator configuration and used these optimal results as a basis for comparison. In particular, we examined optimal, periodic, one-dimensional, in-place hopping, and studied a variety of cost functions in a model-based approach. To this end, we established parametrized models of electrically driven PEA and SEA hoppers that included a dynamical model, a motor model, a detailed transmission model, and realistic limitations on all states and parameters. In this section, we outline these models, discuss our parameter choices, describe the cost functions, and detail our optimization approach. To minimize the number of free parameters in our analysis, all states and parameters were normalized with respect to total mass *m _{o}*, uncompressed leg length

*l*, and gravity

_{o}*g*.

### Dynamical Model.

Our study was based on the simple model of a hopper (Fig. 1). The hopper consisted of a main body with mass *m*_{1} and a point foot segment with mass *m*_{2}. The motion of the hopper was restricted to a pure vertical movement; i.e., hopping in place. Its state was defined by the position and velocity of the main body (given by *y* and $y\u02d9$) and the length and contraction velocity of the leg (*l* and $l\u02d9$). The motion of the foot was coupled to the main body by a spring with stiffness *k* and damping ratio *ζ*. Having damping in the springs and a foot with mass meant that the system was energetically nonconservative. As a result, positive network had to be performed over the course of a stride. This work was created by an electric DC motor with an attached transmission that produced a torque of *T _{o}* and had a reflected inertia of

*j*. The mass of the motor is included as part of

_{o}*m*

_{1}. This actuator was either connected in parallel to the spring (PEA) or in series with the spring (SEA). During flight, the motion of the hopper was governed by the equations of motion $Mq\xa8=h+\tau $, with the mass matrix

**M**, the gravitational terms

**h**, and the generalized forces $\tau $.

*F*, was the sum of the force in the spring and the motor force

*b*was computed from the damping ratio

*ζ*

**q**

_{PEA}= (

*y*,

*ℓ*)

^{T}was only two-dimensional. The reflected rotor inertia

*j*acted along the coordinate of

_{o}*ℓ*, which led to the following mass matrix for PEA:

The system was driven by the generalized forces $\tau PEA=(0,\u2009F)T$. The differentiable force vector was given by **h**_{PEA} = (–(*m*_{1} + *m*_{2})*g*, *m*_{2}*g*)^{T}.

*u*and $u\u02d9$) described the position of the proximal end of the spring (coordinate origins were defined such that for

*u*= 0, and an uncompressed spring, the leg length

*ℓ*was equal to the resting length

*ℓ*). The generalized coordinate vector, thus, had three components:

_{o}**q**

_{SEA}= (

*y*,

*ℓ*,

*u*)

^{T}. Since the motor was in series with the spring, the leg force,

*F*, was equal to the force produced by the spring

The system was driven by the generalized forces $\tau SEA=(0,\u2009F,\u2009To\u2212F)T$. The differentiable force vector was given by $hSEA=(\u2212(m1+m2)g,m2g,0)T$.

For both models, the main body velocity and the leg retraction rate were equal during ground contact $(y\u02d9=\u2113\u02d9)$. That led to the modified equations of motion of $y\xa8=\u2113\xa8=F/(m1+jo)\u2212m1g/(m1+jo)$ for PEA and $y\xa8=\u2113\xa8=(F/m1)\u2212g$ for SEA. This constraint implies that at the moment of touch-down, a collision brought the foot velocity $vfoot=y\u02d9+\u2113\u02d9$ to zero.

### Motor and Transmission Model.

*i*via the motor torque constant

*k*(

_{T}*T*=

_{m}*ik*). The thermal losses in the motor are given by

_{T}where *R* is the armature resistance, *k _{b}* is the motor speed constant, and

*K*is the speed torque gradient of the motor. Here, we made use of the fact that in SI units,

*k*=

_{T}*k*.

_{b}*j*and was modeled to be connected to a transmission system. This transmission had two roles. First, it amplified the motor torque (

_{m}*T*=

_{o}*nT*), while reducing the output speed $(\omega o=(1/n)\omega m)$. Second, it converted the rotational motion of the motor to a linear motion that allowed the hopper's leg to extend and contract. We defined the overall transmission as the combination of a rotary gearbox with transmission ratio

_{m}*n*(a unitless quantity), and a rotary to linear transmission with transmission ratio

_{r}*n*(with normalized units of $rad/\u2113o$). The overall transmission ratio

_{ℓ}*n*was given by the product of

*n*. For simplicity, we neglected the inertia associated with the transmission, assuming that the inertia of the actuator is dominated by the reflected inertia of the rotor

_{r}n_{ℓ}We considered two representative values used for converting rotary to linear motion in modern legged robotics: a large value of $n\u2113=200\u2009rad/\u2113o$ and a small value of $n\u2113=2\u2009rad/\u2113o$.

*n*were chosen because they are roughly representative of the transmission given by a ball screw [31] and a knee joint bent at approximately 45 deg, respectively. For both cases, for simplicity, we considered

_{ℓ}*n*to be a frictionless transmission. This is evident for the knee joint and good approximation for the ball screw, which can have an efficiency of torque conversion of 90% and higher [31]. The only source of transmission friction was, thus, in the rotary gearbox

_{ℓ}*n*. This was an important consideration. A particular choice of

_{r}*n*may be able to better avoid

_{ℓ}*n*frictional losses. We modeled these

_{r}*n*losses as being created by a planetary gearbox with dry friction. The torque output after the transmission is, therefore, given by

_{r}*T*is the torque output by the motor,

_{m}*ω*is the motor speed, and

_{m}*T*is the coulomb frictional torque of the gearbox. We calculated

_{f}*T*in Eq. (10) from the maximum gearbox torque conversion efficiency

_{f}*ε*

_{max}

*T*

_{o,max}is the maximum continuous output torque. Substituting into Eq. (10) and rearranging yields

For comparison, we included a third case in our analysis which assumed a completely frictionless transmission system (with *ε*_{max} = 1). This was done to understand which transmission choices would be optimal in the absence of any penalty associated with larger gearbox friction. Since in this frictionless case, there is no difference between *n _{ℓ}* and

*n*; we simply set

_{r}*n*= 1.

_{ℓ}### Model Parameters.

where *m _{m}* is the mass of the motor in kg,

*T*

_{m}_{,max}is the maximum motor torque,

^{1}and

*ω*

_{m}_{,max}is the maximum permissible speed of the motor. Note that the values in these equations are not yet normalized.

*n*by the exponent

_{r}*γ*. This decision was motivated by the observation that gear ratio scales approximately exponentially with stage number, which, in turns, scales approximately linearly with efficiency [32]. Our efficiency estimate was based on a logarithmic regression over 809 gearboxes from a major motor and gearbox manufacturer [32]. The results of the regression are shown in Fig. 2. We found that

where the final approximation ensured that the efficiency was 100% when there was no gearbox (*n _{r}* = 1).

*T*

_{o,max}, by

All parameters and their limits are reported in Table 1. They are based off of our previous hardware [33], as done in Ref. [28]. In particular, we chose a value of *m _{m}* = 0.6 kg,

*m*= 5 kg, and

_{o}*ℓ*= 0.4 m.

_{o}m_{1} = 0.95 m_{o} | u_{min} = –0.15 ℓ_{0} | n_{r,min} = 1 |

m_{2} = 0.05 m_{o} | u_{max} = 0.15 ℓ_{0} | n_{r,max} = 6285 |

ζ = 0.2 | k_{min} = 0.0001 m/_{o}gℓ_{o} | ℓ_{min} = 0.5 ℓ_{0} |

ℓ_{max} = 1.15 ℓ_{0} | k_{max} = 1000 m/_{o}gℓ_{o} | $K=2.3\xd710\u221204\u20091/mog\u2113o$ |

$jrot=1.48\xd710\u221205\u2009mo\u2113o2$ | T_{m,max} = 0.0157 m_{o}gℓ_{o} | $\omega m,max=748\u2009radg/\u2113o$ |

m_{1} = 0.95 m_{o} | u_{min} = –0.15 ℓ_{0} | n_{r,min} = 1 |

m_{2} = 0.05 m_{o} | u_{max} = 0.15 ℓ_{0} | n_{r,max} = 6285 |

ζ = 0.2 | k_{min} = 0.0001 m/_{o}gℓ_{o} | ℓ_{min} = 0.5 ℓ_{0} |

ℓ_{max} = 1.15 ℓ_{0} | k_{max} = 1000 m/_{o}gℓ_{o} | $K=2.3\xd710\u221204\u20091/mog\u2113o$ |

$jrot=1.48\xd710\u221205\u2009mo\u2113o2$ | T_{m,max} = 0.0157 m_{o}gℓ_{o} | $\omega m,max=748\u2009radg/\u2113o$ |

### Cost Functions.

To measure the energetic efficiency of each actuation variation, we used three cost functions: positive mechanical motor work, thermal electrical losses, and positive electrical work. All cost functions are expressed as an integral over a single hop, from time *t* = 0 until *t* = *T*.

#### Positive Mechanical Motor Work.

#### Thermal Electrical Losses.

*i*

^{2}

*R*losses in the motor windings. It is given by

#### Positive Electrical Work.

As in Ref. [28] we assumed that the robot was unable to recover negative electrical work and store it. Negative motor work could be used, however, to compensate for thermal losses. This cost function provides a close approximation of the total energy discharge from a battery that would be required to power the motors.

### Trajectory Optimization.

A hopping height of *y*(*T*) = 1.3*ℓ*_{0} was chosen to ensure that the hoppers had a flight phase. All bounds are given in Table 1. The motor torque *T _{m}*(

*t*) was parametrized as a piecewise linear function. The aforementioned optimization problem was implemented in the optimization package muscod, which utilizes a fourth/fifth-order Runge–Ketta–Fehlberg numerical integration algorithm [34–36]. We tried multiple initial conditions, all leading to the same results.

Note that the gearbox gear ratio *n _{r}* and the spring stiffness

*k*were free parameters in the optimization. System properties that were unrelated to the actuator configuration were left fixed. These parameters included the masses (

*m*

_{1}and

*m*

_{2}) and the damping ratio

*ζ*.

*x*) function using the following equation:

*σ*= 0.001. We also smoothed the function sign(

*ω*) using the following equation:

_{m}*α*= 100.

## Optimal Configurations, Motions, and Parameters

As a result of the optimization-based motion generation and parameter identification, we found that the ideal actuator type was both dependent on the choice of rotary to linear transmission and on the selected cost function. For positive electrical work, which combines mechanical motor work with thermal losses and reflects a realistic trade-off between these two contributions, we found that SEA was the optimal actuator type for an ideal, frictionless transmission as well as for the $n\u2113=200\u2009rad/\u2113o$ case. For a hopper with $n\u2113=2\u2009rad/\u2113o$, PEA was the optimal actuator type (Fig. 3). In the following, we elaborate on the optimal actuator type for each configuration and detail the resulting motion profiles and parameter choices for each case.

### Frictionless Transmission.

For a completely frictionless transmission, the energetically optimal actuator type depended on the cost function (Fig. 4). For positive mechanical motor work *C*_{mech}, the SEA hopper was 67% more energetically economical than PEA (0.070 *m _{o}gℓ_{o}* versus 0.14

*m*). For both actuator types, energy was primarily lost to damping, accounting for 85% of

_{o}gℓ_{o}*C*

_{mech}for SEA and 96% for PEA. In terms of thermal electrical losses

*C*

_{therm}, the PEA hopper was 72% more economical than SEA (0.017

*m*versus 0.036

_{o}gℓ_{o}*m*). For the electrical work

_{o}gℓ_{o}*C*

_{el}the SEA hopper was 75% more energetically economical than PEA (0.12

*m*versus 0.25

_{o}gℓ_{o}*m*). For both actuator types, the

_{o}gℓ_{o}*C*

_{el}cost was primarily caused by damping and thermal losses in the spring and motor. Negative mechanical motor work was negligible, indicating a preference for following the natural dynamics of the system and exploiting passive storage of excess energy.

When optimizing for *C*_{mech} and *C*_{therm} individually, the results can largely be explained by the inherent advantages of each actuator type. SEA's advantage when optimizing for *C*_{mech} stemmed from the fact that SEA can directly influence the relative rate of spring motion and therefore decrease damping losses. For example, during the second half of stance, the actuator *u* is pushed downwards to inject energy but also to reduce the extension rate of the spring (as discussed in Ref. [38]). Such a strategy is not possible for PEA during stance where the spring motion is inherently coupled to the main body motion. PEA can reduce damping losses during flight, however, by holding the leg at its maximum length. This strategy reduces the effective hopping height and avoids oscillations with their associated damping losses in the air. Additionally, the foot is released slightly before touchdown to have its relative velocity with the ground be approximately zero (i.e., the leg velocity $\u2113\u02d9$ approximately matched the main body velocity $y\u02d9$), minimizing collision losses (we will refer to this as the *clamping strategy*). As for minimizing *C*_{therm}, the main advantage of PEA stemmed from the fact that PEA did not need to support the weight of the robot, which reduced the required motor forces. Furthermore, the parameters for spring stiffness and gear ratio were tuned precisely such that when the leg oscillated nearly freely during flight (without any large peaks in the motor torque *T _{m}*) the leg velocity, $\u2113\u02d9$ approximately matched the main body velocity at touchdown (we will refer to this as the

*oscillation strategy*), which again minimized collision losses (Fig. 5).

What is surprising is that the rather small margins for each of these two cost functions turned into a large benefit for SEA in the combined cost function of positive electrical work *C*_{el}. This clear advantage of SEA can be explained by the fact that the optimal motion profiles for *C*_{mech} and *C*_{therm} were very similar (Fig. 5). The motor trajectory held the leg length nearly constantly at its maximum extension during flight and injected energy into the system during stance. It was, therefore, easy for SEA to find a strategy that could keep both *C*_{mech} and *C*_{therm} small at the same time when optimizing for *C*_{el}. The total *C*_{el} effort is merely a sum of the individual cost functions (Fig. 4). For PEA, in contrast, the motion trajectory differed significantly depending on the chosen cost function. For *C*_{mech}, the motion of the leg was similar to SEA, utilizing the clamping strategy; whereas for *C*_{therm}, the optimal motion is the oscillation strategy. For the combined cost function *C*_{el}, PEA, therefore, had to resolve a severe trade-off between two very different actuation strategies. As a result, for *C*_{el}, PEA had both higher mechanical work (primarily driven by an increase in damping losses) and higher thermal losses than for either *C*_{mech} or *C*_{therm}.

The aforementioned flight motion strategies indicate an important difference between our modeling and previous modeling efforts. The presence of foot mass and actuator inertia in our model allows us to present these optimal actuation strategies during the flight phase. Those strategies differ, for example, from those presented in Refs. [9], [39], and [40] which, due to their massless foot, only look at forces during the stance phase. Furthermore, having massless feet, and no damping, as the aforementioned papers do, allows for entirely lossless hopping strategies. For example, Alexander [9] states that for a SEA “it is possible to find a running gait that requires no work from the telescopic actuators, for any combination of speed and stride length.” For PEA, the same can be said. Without collision losses or damping, the spring can passively store and replace all of the energy during hopping. Having both cases be entirely lossless would make a comparison of the two types of actuators impossible, and therefore justifies the added complexity of our model.

For all cost functions, the optimal PEA hopper had a significantly lower transmission ratio *n* = *n _{r}n_{l}* than SEA (Table 2). Smaller values for

*n*are better for PEA, because the reflected rotor inertia directly adds to the collision losses (Eq. (1)). Since the spring force and motor torque act additively for PEA, the output torques from the motor are much smaller than for SEA, and there is, thus, no large penalty for using a smaller gearbox. Along the same lines of reasoning, SEA required a larger

*n*to reduce the effective motor torque

*T*and the associated thermal losses. Since leg motion and actuator motion were decoupled, there was no penalty for the resulting larger reflected rotor inertia.

_{m}Cost function | Actuator type | k (m/_{o}gℓ)_{o} | n (rad/ℓ)_{o} |
---|---|---|---|

Positive motor work C_{mech} | PEA | 27.5 | 41.0 |

Positive motor work, C_{mech} | SEA | 192 | 789 |

Thermal losses, C_{therm} | PEA | 21.5 | 193 |

Thermal losses, C_{therm} | SEA | 37.8 | 703 |

Positive electrical work, C_{el} | PEA | 7.80 | 103 |

Positive electrical work, C_{el} | SEA | 186 | 797 |

Cost function | Actuator type | k (m/_{o}gℓ)_{o} | n (rad/ℓ)_{o} |
---|---|---|---|

Positive motor work C_{mech} | PEA | 27.5 | 41.0 |

Positive motor work, C_{mech} | SEA | 192 | 789 |

Thermal losses, C_{therm} | PEA | 21.5 | 193 |

Thermal losses, C_{therm} | SEA | 37.8 | 703 |

Positive electrical work, C_{el} | PEA | 7.80 | 103 |

Positive electrical work, C_{el} | SEA | 186 | 797 |

For both actuation types, *n* was much smaller for the positive motor work cost function *C*_{mech} than for the other two cost functions. This result reflects again the fundamental differences that arise when considering *C*_{therm}. For *C*_{therm}, minimizing motor torques plays a more important role. Increasing *n* leads to smaller required *T _{m}*. This consideration is not important for

*C*

_{mech}, which is independent of

*n*.

### Rotary Gearbox With Friction, nℓ=200 rad/ℓo.

Introducing friction in the transmission led to a new trade-off: in addition to creating a larger reflected rotor inertia, larger transmission ratios were now additionally penalized by larger friction values. For a rotary to linear transmission value of $n\u2113=200\u2009rad/\u2113o$, this trade-off is not particularly grave.

As a result, the energetically optimal actuator type was cost function dependent in the same way as for the frictionless transmission (Fig. 6). For *C*_{mech}, the SEA hopper was 65% more energetically economical than PEA (0.096 *m _{o}gℓ_{o}* versus 0.19

*m*). For both actuators, energy was still primarily lost to damping, accounting for 84% of

_{o}gℓ_{o}*C*

_{mech}for SEA and 97% for PEA. Frictional losses in the gearbox were relatively low, as both actuators chose small optimal

*n*values, which, therefore, had high efficiencies (Table 3). For SEA, the optimizer chose an overall transmission ratio of $n=256\u2009rad/\u21130$, while PEA had $n=200\u2009rad/\u21130$ (i.e., no rotary gearbox at all). As a result, for SEA, the frictional losses accounted only for 7% of the SEA losses and for none of the PEA losses. For

_{r}*C*

_{therm}, the PEA hopper was 131% more energetically economical than SEA (0.017

*m*versus 0.080

_{o}gℓ_{o}*m*). For

_{o}gℓ_{o}*C*

_{el}, the SEA hopper was 20% more energetically economical than PEA (0.28

*m*versus 0.35

_{o}gℓ_{o}*m*).

_{o}gℓ_{o}Cost function | Actuator type | k (m/_{o}gℓ)_{o} | n (rad/ℓ)_{o} |
---|---|---|---|

Positive motor work, C_{mech} | PEA | 16.4 | 200 |

Positive motor work, C_{mech} | SEA | 23.3 | 256 |

Thermal losses, C_{therm} | PEA | 22.9 | 200 |

Thermal losses, C_{therm} | SEA | 25.6 | 492 |

Positive electrical work, C_{el} | PEA | 5.93 | 200 |

Positive electrical work, C_{el} | SEA | 23.4 | 500 |

Cost function | Actuator type | k (m/_{o}gℓ)_{o} | n (rad/ℓ)_{o} |
---|---|---|---|

Positive motor work, C_{mech} | PEA | 16.4 | 200 |

Positive motor work, C_{mech} | SEA | 23.3 | 256 |

Thermal losses, C_{therm} | PEA | 22.9 | 200 |

Thermal losses, C_{therm} | SEA | 25.6 | 492 |

Positive electrical work, C_{el} | PEA | 5.93 | 200 |

Positive electrical work, C_{el} | SEA | 23.4 | 500 |

The results for the three cost functions followed the same general trend as in the completely frictionless transmission. The results held largely because of the advantage SEA obtains from utilizing a highly efficient rotary to linear transmission with a large transmission ratio. With the rotary to linear transmission doing most of the reduction $(nl=200\u2009rad/\u21130)$, only small values for *n _{r}* are necessary and there is a very small penalty for SEA to choose a large transmission ratio

*n*=

*n*(Eq. (17)). Therefore, the SEA costs showed only a slight increase due to gear friction. PEA, however, favored smaller transmission ratios for

_{l}n_{r}*C*

_{mech}and

*C*

_{el}. Here, PEA attempted to have as low a transmission ratio as possible, choosing to have

*n*= 1 for all cost functions. Still, with $nl=200\u2009rad/\u21130$, it was forced to have a minimum transmission ratio that was larger than its optimal choice for the frictionless transmission. This nonoptimal choice of

_{r}*n*drives up PEA costs for

*C*

_{mech}and

*C*

_{el}much higher than the slight gear friction for SEA, maintaining the trends from the frictionless transmission while increasing costs. In particular, for PEA, negative mechanical motor work became a significant portion of the

*C*

_{el}losses, likely to avoid excess collisional losses from the larger transmission (Eq. (1)). For

*C*

_{therm}, the PEA cost and transmission ratio are nearly identical to the frictionless transmission.

For both PEA and SEA, the optimal motion was nearly identical to the frictionless transmission for *C*_{mech} and *C*_{therm} (Fig. 7). For *C*_{el}, there was again a trade-off between the oscillation strategy and the clamping strategy. The optimal *C*_{el} motion, however, changed. Whereas for the frictionless case, PEA chose the oscillation strategy, here it chose the clamping strategy. The consequences of the choice to suppress the natural oscillatory motion of the spring can be seen in the largely increased negative mechanical motor work for *C*_{el}.

### Rotary Gearbox With Friction, nℓ=2 rad/ℓo.

The rotary to linear transmission value of $n\u2113=2\u2009rad/\u2113o$ in this case is relatively small. This case will, therefore, require substantially larger *n _{r}* values, which, in turn, will create larger friction losses.

This clearly benefited PEA, which was the energetically optimal actuator type for all cost functions in this case (Fig. 8). For *C*_{mech}, the PEA hopper was 77% more energetically economical than SEA (0.21 *m _{o}gℓ_{o}* versus 0.46

*m*). For PEA, energy was primarily lost to damping, accounting for 65% of the losses. The frictional losses in the gearbox accounted for 30% of losses. For SEA, the trend was reversed: 75% was lost to gearbox friction, whereas 21% was lost to damping. For

_{o}gℓ_{o}*C*

_{therm}, the PEA hopper was 41% more energetically economical than SEA (0.44

*m*versus 0.67

_{o}gℓ_{o}*m*). For

_{o}gℓ_{o}*C*

_{el}, the PEA hopper was 39% more energetically economical than SEA (0.77

*m*versus 1.1

_{o}gℓ_{o}*m*). For both actuators, the losses were dominated by thermal losses and gear friction (Fig. 8).

_{o}gℓ_{o}In contrast to the previous cases, there was now a significant penalty for SEA to choose a large *n*, as it directly led to a large *n _{r}* and meant that there were significant frictional losses (Eq. (17)). PEA, with its smaller optimal rotary gearbox transmission

*n*therefore suffered far lower frictional losses (Table 4).

_{r}Cost function | Actuator type | k (m/_{o}gℓ)_{o} | n (rad/ℓ)_{o} |
---|---|---|---|

Positive motor work, C_{mech} | PEA | 19.8 | 27.0 |

Positive motor work, C_{mech} | SEA | 41.8 | 306 |

Thermal losses, C_{therm} | PEA | 31.2 | 90 |

Thermal losses, C_{therm} | SEA | 91.8 | 406 |

Positive electrical work, C_{el} | PEA | 34.4 | 63.0 |

Positive electrical work, C_{el} | SEA | 119 | 388 |

Cost function | Actuator type | k (m/_{o}gℓ)_{o} | n (rad/ℓ)_{o} |
---|---|---|---|

Positive motor work, C_{mech} | PEA | 19.8 | 27.0 |

Positive motor work, C_{mech} | SEA | 41.8 | 306 |

Thermal losses, C_{therm} | PEA | 31.2 | 90 |

Thermal losses, C_{therm} | SEA | 91.8 | 406 |

Positive electrical work, C_{el} | PEA | 34.4 | 63.0 |

Positive electrical work, C_{el} | SEA | 119 | 388 |

The optimal motion for both PEA and SEA was now nearly identical for all three cost functions (Fig. 9). For all cost functions, both hoppers used the clamping strategy. There was no longer any oscillatory behavior by the PEA hopper. This lack of oscillatory behavior likely results from the presence of significant rotary gearbox friction. PEA could no longer choose as small of an *n _{r}* value as in the previous two cases, meaning that gear friction induced losses whenever the leg moved. Therefore, the optimizer chose to clutch the leg for all cost functions, avoiding unnecessary leg oscillation and removing any trade-off for

*C*

_{el}.

*C*

_{el}represented primarily the sum of

*C*

_{mech}and

*C*

_{therm}.

## Discussion and Conclusions

In this study, we employed numerical optimization to compare optimal hopping motions of series and parallel elastic actuated hoppers. To make the comparison between the two actuation concepts as fair as possible, we compared the best possible hopper with SEA to the best possible hopper with PEA. To this end, we optimized actuator forces, motion trajectories, and system parameters simultaneously. The analysis was performed for three different cost functions: positive mechanical motor work *C*_{mech}, electrical thermal losses *C*_{therm}, and positive electrical work *C*_{el}. We studied hopping for two representative cases of rotary to linear transmission. In addition, we compared these results to an idealized case without friction in the transmission. Optimizing system parameters such as the rotary gearbox ratio and the stiffness of the springs was crucial. The different actuator types did require very different parameter values to perform optimally.

Given that *C*_{el} is the most complete cost function, we would conclude that for electrically driven hoppers, SEA is the most energetically economic actuator solution when the majority of the transmission ratio is accommodated for by a very high efficiency rotary to linear transmission, such as a prismatic ball screw. For a hopper with a smaller transmission (such as a robot hopping with a bent knee), however, PEA performed better. These optimal choices are primarily a consequence of the fact that SEA requires larger overall transmission ratios than PEA, as could be seen clearly when looking at the frictionless transmission. When the transmission ratios are achieved in an efficient way, SEA is the better actuator type. If they are subject to friction losses, PEA becomes the better actuator type.

Throughout our analysis, it became evident, that each configuration had a unique optimal motion profile and a set of parameters that differed greatly between the two actuation concepts. This clearly illustrated the necessity of our optimization approach. From a conceptual point of view, this simultaneous optimization of motion and morphology extends the paradigm of system-based design to include the future motion of the robotic system already at the design stage. We strongly believe that such a combined and simultaneous optimization of robot and motion will be at the core of future robotic design. This technique is slowly making its way into robotics, for example, in the optimization of an assistive device for sit-to-stand motions [41], to create stable running motions in a human-like robot [18] and to optimize limb lengths for particular tasks [19]. Here, we extended its principle to make a discrete design comparison between two fundamentally different morphologies. In the future, this approach can potentially be applied to more complex models and robots to achieve truly optimal performance during the design phase.

There are further important considerations for the actuator problem that must be explored before implementation in hardware. For PEA, the motor inertia is not isolated from collisions, which can lead to significant energetic losses. Therefore, it is imperative that the foot velocity is near-zero at touchdown. PEA achieves this by utilizing the oscillation strategy and clamping strategy, which let the spring passively retract the foot and minimize foot velocity with minimal energetic effect. Overall, we would argue that these optimal hopping strategies for PEA are more sensitive to timing issues and would be harder to implement in hardware. These strategies require a precise tuning of the spring stiffness, rotor inertia, and foot mass to achieve the correct spring natural frequency. However, small deviations from the optimal parameter values (e.g., caused by modeling errors or disturbances in the motion) can get in the way of this phase matching and might result in large collision losses.

In the present study, to avoid biasing the results toward either actuation configuration, we did not constrain the hopping frequency. The primary constraint on the hopping motion was that the hopper had to reach a height of *y* = 1.3*ℓ*_{0}. This constraint ensured that the hoppers had both a stance phase and a flight phase. Timing was not considered, however, in finding the energetically optimal motion. From the motion plots, it is clear that SEA has significantly shorter optimal hopping periods, suggesting a potential advantage over PEA. In our prior work, the effect of frequency was implicitly included by enforcing an average velocity during two-dimensional motion [27] and was found to have little effect on the energetics at most speeds. Still, future studies could consider exploring hopping at a fixed frequency to see additional sensitivity effects.

There are additional mechanical and sensing complexities that must be taken into account when implementing either actuator type in hardware. For example, in PEA, the gearbox is not isolated from collisions. To prevent unnecessary damage to the gears, PEA would likely need an additional mechanism that softens these collisions. Furthermore, the spring could be disengaged to allow for easier leg retraction during swing. Karssen [16], for example, implements a clutch that disengages the spring during flight, and a small series elastic element that protects the gearbox. SEA would need additional components as well. Since SEA moves the proximal end of the spring, which, in turn, changes the leg length, it would need additional sensors to measure the location of each. More moving parts for SEA further means that additional bearings and mechanical components are necessary in practice.

Despite the simplicity of our models compared to a final hardware implementation, the motion strategies shown here can provide templates for further optimizations with more complex models and gaits. Given that hopping can be considered the archetype of legged locomotion [42,43], one would expect that similar motion patterns will be found in active running and walking gaits of multilegged robots. They can, thus, be used to systematically initialize optimizations or be employed as a tool to check the viability of more complex optimization results. This could enable roboticists to detect and avoid local minima that do not represent the global optimum in their optimizations. Finally, the actuation strategies that we observed in the optimal solutions provide an interesting departure point for further hardware development. For example, the clamping strategy that the optimizer discovered for PEA, where the leg was held fixed during the flight phase, could be implemented more efficiently with an actual clutch rather than with a DC motor.

The results presented in this work have been developed for the particular case of 1D hopping with electrically driven actuators. Certain parameters, such as the mass of the motors, are based off of our previous hardware [33]. The obtained motions, parameter choices, and cost values are specific to this choice of problem. From the work presented here, it is clear that even for this particular problem, the optimal actuation type is case dependent. It is impossible to state conclusively whether in general SEA or PEA is the better actuation concept for legged robots. Naturally, the cases presented in this work only represent a sample of a wide range of possible hopper configurations and focus on energetics. There likely is a gradual transition from hoppers, which would benefit from SEA to hoppers that should use PEA as friction in the transmission becomes more important. That is, ball screws are available with a wide variety of transmissions, knees can be operated at many different angles, and hybrids of the two transmission systems even exist [44]. Each of these transmission types can benefit a particular actuator. Additionally, more complexity has been introduced by robotic prototypes that utilize both SEA and PEA within a single design [16]. Future studies should take into account the vast variety of possible configurations.

The highly complex dependence on actuation type, parameters, and motion profile emphasizes the necessity of optimization to explore design issues. Still, it is important here to state that the methods we use are local, and so we are unable to guarantee their global optimality. In hopes of escaping local minima, all of the results presented in this paper represent optima that have been tested from a variety of different initial conditions.

Given that our evaluation took into account different cost functions that penalized force and work independently, our results can still be considered fairly general. Furthermore, this paper should also be understood as an introduction of a systematic methodology that enables a fair comparison of different robot morphologies. In a sense, the comparison of PEA and SEA is just an application of this method. In the future, this methodology can and should be extended to include other models as well as more complex motions. Moreover, the methodology can be utilized to look at other cost functions to maximize performance metrics other than energetics, such as maximizing hopping height.

Overall, the presented work gives a deeper understanding of the optimal motion for both SEA and PEA driven hopping motions. Understanding this underlying optimal motion can help designers ensure that they are not losing energy by unintentionally fighting the natural dynamics of a particular type of actuator. This understanding also gives insight into the coupling of morphology and motion. By using trajectory optimization, we can see how the choice of transmission ratio and stiffness change for a particular implementation and how these parameters and the resulting motion mutually influence each other. Recognizing these relationships opens the possibility to make robots more efficient, faster, and more robust.

## Funding Data

National Science Foundation (Grant No. 1453346).

This value is twice the maximum continuous torque the motor can provide. It is the maximum value suggested through personal communication with Maxon Motor.