Optimization of energy efficiency of walking bipedal robots by use of elastic couplings in the form of mechanical springs
Authors: Fabian Bauer, Ulrich Römer, Alexander Fidlin, Wolfgang Seemann
Institute of Engineering Mechanics Karlsruhe Institute of Technology Kaiserstraße 10 76131 Karlsruhe Germany
Contact: [email protected]
Version: accepted manuscript Accepted: 13 September 2015 Published: 12 October 2015 (online)
DOI: 10.1007/s11071-015-2402-9
The final publication is available at link.springer.com Nonlinear Dyn manuscript No. (will be inserted by the editor)
Optimization of energy efficiency of walking bipedal robots by use of elastic couplings in the form of mechanical springs
Fabian Bauer · Ulrich Römer · Alexander Fidlin · Wolfgang Seemann
Received: date / Accepted: date
Abstract This paper presents a method to optimize the en- 1 Introduction ergy efficiency of walking bipedal robots by more than 50 % in a speed range from 0.3 to 2.3m/s using elastic couplings There are many and varied application scenarios of bipedal – mechanical springs with movement speed independent pa- robots. The most impressive ones are the humanoid as surro- rameters. The considered robot consists of a trunk, two stiff gate of workers in disaster response saving human lives and legs and two actuators in the hip joints. It is modeled as un- the exoskeleton as enhancement of disabled people improv- deractuated system to make use of its natural dynamics and ing human lives. These applications share inherently high feedback controlled with input-output linearization. A nu- mobility requirements, which prohibit an external power sup- merical optimization of the joint angle trajectories as well as ply and therefore directly demand for high energy efficiency. the elastic couplings is performed to minimize the average The newly developed robot ATLAS from Boston Dynamics energy expenditure over the whole speed range. The elastic is employed in a disaster response scenario in the Robotics couplings increase the swing leg motion’s natural frequency Challenge of the DARPA Program Maximum Mobility and thus making smaller steps more efficient which reduce the Manipulation (M3). Its 23 kg lithium-ion-battery provides impact loss at the touchdown of the swing leg. The pro- energy for no more than 20 min of running – the program cess of energy turnover is investigated for the robot with and consistently demands to improve the energy efficiency by a without elastic couplings. Furthermore, the influence of the factor of 20 [33]. In order to evaluate the efficiency of lo- elastic couplings’ topology, its degree of nonlinearity, the comotion, the dimensionless specific energetic cost of trans- mass distribution, the joint friction, the coefficient of static port cT = energy input/(weight × distance traveled) is ap- friction and the selected actuator is analyzed. It is shown that plied. For example, the highly developed humanoid ASIMO the optimization of the robot’s motion and elastic coupling from Honda (cT = 3.2) has a 16 times higher specific cost of towards energy efficiency leads to a slightly slower conver- transport than a human (cT = 0.2) at a low walking speed of gence rate of the controller, yet no loss of stability and a v = 0.4m/s [12]. Due to this performance gap between ar- lower sensitivity with respect to disturbances. The optimal tificial bipedal systems and their biological counterpart, en- elastic coupling discovered by the numerical optimization is ergy efficiency can be justifiably stated as one of the major a linear torsion spring between the legs. challenges in the development of walking bipedal robots. There are essentially two reasons for the bad energy ef- Keywords bipedal robot · dynamic walking · nonlinear ficiency of todays bipedal robots. First the conservative con- feedback control · optimization · energy efficiency · elastic trol strategy: the robot is operated close to static equilibrium, coupling · mechanical spring fights against gravitation and tries to suppress any natural dynamics. Humans in contrast walk dynamically, they con- stantly fall over their feet, swinging their legs like pendu- F. Bauer · U. Römer ( ) · A. Fidlin · W. Seemann lums during walking and therefore exploit their natural dy- Institute of Engineering Mechanics, Karlsruhe Institute of Technology, Kaiserstraße 10, 76131 Karlsruhe, Germany namics instead of struggling against them [30]. Second the Tel.: +49-721-60846823 design: the robot consists of rigid links interconnected only Fax: +49-721-60846070 by actuated joints. Humans in contrast consist of both rigid E-mail: [email protected] skeletons and elastic structures serving as elastic boundary 2 Fabian Bauer et al. layer to absorb shocks, as pogo-stick to redirect center of sive robot model consisting of trunk, two rigid legs and two mass motion and as return spring to buffer energy during rotational springs between trunk and legs [22]. However, the the redirection of the swing leg motion [3]. gait without energy costs holds only for one specific speed, Elastic boundary layers acting as shock absorbers are is unstable and comes along with a heavy oscillation of the used in walking robots to reduce the load of the robots’ trunk. On an incline this model has shown a much more sta- mechanical parts, especially the actuator gear boxes [24,25, ble gait than a comparable model without trunk [11]. 28]. Because of their high stiffness compared to other elastic A big step towards an energy efficient and versatile robot structures, they have a minor influence on the overall system was made by enhancing this model with actuation in the hips dynamics. [18,31]. Stiffness and resting length of hip springs were op- Elastic couplings acting as pogo-stick have been used for timized together with the robot’s motion. A significant re- a long time in hopping and running robots [2,36] and their duction of actuation could be shown for one constant speed. influence on the dynamics has been analyzed in detail with However, both studies have three major drawbacks, which the model of the spring loaded inverted pendulum [6,9,21, prevent their practical application: (1) the geometrical and 37]. inertial parameters of the investigated robots were far away Elastic couplings acting as return springs have first been from humanoid robots, (2) the objective function of the op- used solely in passive dynamic running [29,38]. Later an timization – integral of square of actuator torques over one elastic coupling in the form of a hip spring was also used in step – has no physical relationship with energy used and spatial passive dynamic walking to adjust the step frequency does not account for variable step length and (3) the elastic to the frequency of the toddling motion for stabilization [26]. couplings are optimized only for one specific speed whereas However, this is not a necessary measure [1]. The idea of ad- a real robot has to be energy efficient over a range of speeds. justing the step frequency by a hip spring was also applied In order to enable the practical applicability of elastic to prevent stumbling by shortening the swing phase of pla- couplings in bipedal robots, these shortcomings are address- nar passive dynamic walking robots with trunk in simulation ed. A process is developed to optimize the motion and the [41] and experiment [40]. elastic couplings simultaneously over a broad range of speeds By enhancing the simplest walking model [14,20] with with respect to energy efficiency [5]. It allows for a thorough an elastic coupling in form of a hip spring, the preferred investigation of the effects of elastic couplings to develop an speed–step length relationship of humans could be predicted energy efficient robot [4]. [27]. With the same model it could be shown that at a given This paper is organized as follows. Section2 introduces target speed the elastic coupling reduces the step length and the investigated robot model. It consists of the mechanics therewith the impact loss at touchdown of the swing leg model (2.1), the actuator model (2.2) and the control strat- [17]. The robot Cornell Ranger used this effect to walk a egy (2.3). Section3 describes the process of generating the marathon at one battery charge [7]. motion (3.1) as well as selecting the optimal elastic coupling Elastic couplings acting as return springs have also been (3.2) of the robot by simultaneous numerical optimization. used in more complex walking models with segmented legs. Section4 presents the effect of elastic couplings on the in- For a passive model comprising of trunk, two thighs and two vestigated robot. The reduction of specific cost of transport shanks it could be demonstrated that an elastic coupling be- by elastic couplings is explained (4.1) and the influence of tween trunk and thigh is necessary for walking with upright the elastic couplings’ topology, the robot’s mass distribu- trunk [8]. Using an elastic coupling over multiple joints, a tion, joint friction, the coefficient of static friction and the trunkless model with impulsive push off force showed an selected actuator (4.2) as well as stability and sensitivity of increase in speed at constant energy input [16]. the motion (4.3) is displayed. Section5 concludes the in- In addition there was the attempt to rebuild the anatomy vestigation, the used methods and the resulting insights and of the biological example and its elastic couplings in detail. gives an outlook on future research aims. AppendixA pro- The robot FastRunner was inspired by an ostrich and uses a vides a more detailed discussion of the results from Sec. 4.2. network of nonlinear elastic couplings with different modes for the swing phase and the stance phase of each leg. This allows for actuation of the six segmented leg by only one hip actuator [13]. In order to approach human locomotion 2 Robot model capabilities, the robot BioBiped was designed [35]. How- ever, its elastic couplings have only been investigated for a The model of the investigated robot consists of the mechan- hopping motion with given joint torque trajectories without ics model (2.1), the actuator model (2.2) and the control any optimization [34]. strategy (2.3). The approach of modeling the robot as under- On the quest for an energy efficient bipedal robot a col- actuated system and control it with input-output lineariza- lision free gait on level ground was found with a planar pas- tion was taken from [39]. Optimization of energy efficiency of walking bipedal robots by use of elastic couplings 3
2.1 Mechanics model powered by the torque due to the weight rCoMx mg about foot 1 as depicted in Fig.1(b). The resulting motion can be The robot is modeled as planar rigid body system. It is com- described by the angular momentum theorem posed of two rigid kneeless legs which are separately con- ˙ nected to a rigid trunk via two ideal revolute joints in the hip L1 = −rCoMx (qJ,qA)mg (2) as depicted in Fig.1(a). The stance leg foot 1 is modeled a priori as ideal revolute joint and the swing leg foot 2 is as- with the angular momentum L1 of the total system about sumed to move freely. It will be checked a posteriori, if the foot 1. Since the trajectories of the joint angles qJ are pre- unilateral contact and stiction conditions are fulfilled. The scribed by the control input u, Eq. (2) corresponds to the re- ground penetration of foot 2 ("scuffing") is not considered, maining dynamics. In order to reflect this situation in the de- since the planar model without knees cannot lift its foot. The scribing equations and to facilitate the control design, Eq. (1) T is transformed into the mixed partial feedback linearized configuration (q = [q1, q2, q3] ) of the robot is described by T normal form. The mass matrix M does not depend on the two joint angles (qJ = [q1, q2] ) and one absolute orienta- absolute angle q and the torques u do not directly act on it. tion (qA = q3). Two actuators between trunk and legs (cf. A T Therefore, the equation of motion (1) can be subdivided into Fig.3(a)) apply torques ( u = [u1, u2] ) along the respective joint angles (qJ) and control their evolution. The dynamics of the multibody system can be described by M (q ) M (q )q¨ Q (q,q˙) B JJ J JA J J + J = J u . (3) M(q)q¨ + Q(q,q˙) = Bu (1) MAJ (qJ) MAA (qJ) q¨A QA (q,q˙) 0 with the mass matrix M(q), the vector of generalized forces The angular acceleration of the absolute orientation Q(q,q˙) and the input matrix B. The elastic couplings are −1 introduced in Sec. 2.1.3, their restoring torques are included q¨A = −MAA (qJ)(MAJ (qJ)q¨ J + QA (q,q˙)) (4) in the vector of generalized forces. can be calculated with the second row of Eq. (3) and elimi- The later prescription of the joint angle trajectories (qJ) by the control (cf. Sec. 2.3) can be regarded as introduction nated in the first row. Inverse dynamics determine the joint of virtual constraints. No torques are applied about the point torques foot 1 or more precisely about the joint towards ground. u = B−1 M(q )v + Q(q,q˙) , (5) Therefore, the evolution of the absolute orientation is sub- J J −1 ject to the remaining dynamics. The bipedal gait consists of M(qJ) = MJJ (qJ) − MJA (qJ)MAA (qJ)MAJ (qJ) , (6) two distinct phases which are distinguished by the number −1 Q(q,q˙) = QJ (q,q˙) − MJA (qJ)MAA (qJ)QA (q,q˙) (7) of feet on the ground. Adding this discrete state variable to the continuous state variables converts the system to a hy- as a function of the angular acceleration of the joint angles brid system. v = q¨ J, which becomes the new control input. Furthermore, the absolute orientation of the trunk qA is substituted by the absolute orientation of the stance leg (cf. Fig.1(c))
θ = q1 + q3 = cθ q . (8)
The canonical change of coordinates
q I 0 q q˜ = J = J = Hq (9) θ cθ qA
transforms the mass matrix to
T M˜ (q ) = H−1 M(q )H−1 . (10) Fig. 1 Mechanics model of the robot (a) with single support phase as J J rotation of the total system about foot 1 (b) and double support phase as instantaneous impact of foot 2 (c) The total angular momentum with respect to foot 1
q˙ L = M˜ (q ) M˜ (q ) J (11) 1 AJ J AA J θ˙ 2.1.1 Single support phase is the generalized momentum conjugate to θ and can be During the single support phase solely foot 1 is on the ground. determined directly with the mass matrix M˜ and the angu- It corresponds to a rotation of the total robot about foot 1 lar velocities q˜˙. Hence, the system equation for the robot’s 4 Fabian Bauer et al. stance phase in first order form can be displayed as does not change over the impact, whereby +(−) indicates the respective quantity directly before (after) the impact. The re- q q˙ 0 J J sulting motion can be described by the momentum theorem d θ M˜ −1 (q )L − M˜ (q )q˙ 0 x˙ = = AA J 1 AJ J J + v in integral form dt q˙ J 0 I L 0 Z t+ 1 −rCoMx (qJ,θ)mg − + − − ˆ Ma qa q˙ a −Ma qa q˙ a = F2,a = lim F2,a(t)dt (19) =: f(x) + g(x)v (12) t+→t− t− with the generalized impulse ˜ ˜ −1 ˜ with MAJ (qJ) = MAA (qJ)MAJ (qJ). − T The ground contact force acting in foot 1 can be derived ∂r2 (qa ) Fˆ 2,a = Fˆ 2 (20) from the free model in Fig.1(c) and the introduction of the ∂qa constraint equation of the augmented system. The condition for plastic impact r = r˙ = r¨ = 0 . (13) 1 1 1 ∂r (q−) r˙ q+ = 2 a q˙ + = 0 (21) T 2 a a T T ∂qa The configuration vector qa = q ,r1 is augmented by the position vector r1 of foot 1. The equation of motion of is stated on kinematic level. After combining Eq. (19)-(21) the free model is given by into a system of linear equations they can be solved for the unknowns. The angular velocity of the enhanced system af- Ma (qa)q¨ a + Qa (qa,q˙ a) = Bau + F + F (14) 1,a 2,a ter impact as well as the impulse acting on foot 2 with the mass matrix Ma (qa), the vector of the generalized q˙ + ∆ (q−) forces Q (q ,q˙ ), the input matrix B and the vector of the a q˙ a − a a a a ˆ = − q˙ (22) F2 ∆ Fˆ (q ) generalized external forces F1,a, F2,a of the augmented sys- 2 tem. Since the swing leg’s foot does not interact with the result as linear mapping of the angular velocity before the ground during the single support phase, F2,a = 0. The gen- impact using the identities and abbreviations eralized external force I I T − − − − ∂r1 0 qa = q , q˙ a = q˙ , (23) F1,a = F1 = (15) 0 0 ∂qa F1 T !−1 ∂r2 −1 ∂r2 ∂r2 I can be expressed as projection of the ground contact force ∆ Fˆ = − Ma , (24) 2 ∂qa ∂qa ∂qa 0 acting on foot 1 with the position vector r1 onto the config- − qa=qa uration space. Writing Eq. (14) as T −1 ∂r2 I ∆ q˙ = M ∆ ˆ + . (25) a a F2 0 Mqq (qa) Mqr (qa) q¨ Qq (qa,q˙ a) B 0 ∂qa − + = u + qa=qa Mrq (qa) Mrr (qa) r¨1 Qr (qa,q˙ a) 0 F1 (16) After the impact the role of the legs is switched, the former stance leg becomes the new swing leg and vice versa. Left and substituting Eq. (13) into the second line results in and right leg are not distinguished, since only symmetrical −1 gaits are investigated. Using the switching matrix F1 (q,q˙,u) = Mrq (q)M (q)(Bu − Q(q,q˙)) + Qr (q,q˙) , (17) 0 1 0 R = 1 0 0 (26) after elimination of q¨ with Eq. (1). 0 0 1 2.1.2 Double support phase the mapping of the multibody state variables over the in- stance of impact can be stated as During the instantaneous double support phase both feet are + − on the ground. It is modeled as plastic impact of foot 2 onto q ∆ qq + = − − , (27) the ground with simultaneous contact opening in foot 1 (no q˙ ∆ q˙ (q )q˙ interaction: F1,a = 0), as depicted in Fig.1(c). After the ∆ q = R , (28) touchdown, foot 2 is at rest. It neither slips horizontally nor − − ∆ q˙ q = [R0]∆ q˙ q . (29) rebounds vertically. The impact duration is assumed to be a infinitesimally short, therefore the configuration of the sys- The angular momentum about foot 2 does not change during tem impact
+ − + − qa = qa (18) L2 = L2 . (30) Optimization of energy efficiency of walking bipedal robots by use of elastic couplings 5
By changing the reference point of the angular momentum from foot 1− to 2− just before the impact and switching the role of the legs from 2− to 1+ just after the impact the an- gular momentum of the total system with respect to foot 1 is