About this Journal Submit a Manuscript Table of Contents
Journal of Robotics
Volume 2011 (2011), Article ID 735407, 15 pages
Research Article

Output-Based Control of Robots with Variable Stiffness Actuation

Dipartimento di Elettronica, Informatica e Sistemistica, Università di Bologna, Viale Risorgimento 2, 40136 Bologna, Italy

Received 12 April 2011; Accepted 25 August 2011

Academic Editor: G. Muscato

Copyright © 2011 Gianluca Palli and Claudio Melchiorri. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


The output-based control of a redundant robotic manipulator with relevant and adjustable joint stiffness is addressed. The proposed controller is configured as a cascade system that allows the decoupling of the actuators dynamics from the arm dynamics and the consequent reduction of the order of the manipulator dynamic model. Moreover, the proposed controller does not require the knowledge of the whole robot state: only the positions of the actuators and of the joints are necessary. This approach represents a significant simplification with respect to previously proposed state feedback techniques. The problem of controlling simultaneously the position trajectory and the desired stiffness in both the joint and work space is investigated, and the relations between the manipulator redundancy and the selection of both the joint and work space stiffness of the manipulator are discussed. The effectiveness of the proposed approach is verified by simulations of a 3 degrees of freedom planar manipulator.

1. Introduction

Physical human-robot interaction currently represents one of the most challenging issues in robotics. To this end, a new class of robots is under development considering safety and dependability as a primary goal. While early applications of the variable stiffness concept were developed for prosthetic purposes due to the bio-inspired ability of these actuation systems of emulating the behavior of human muscles [1], variable stiffness actuators have been recently introduced in robotics to achieve a suitable tradeoff between accuracy in tasks execution and safety while operating in unstructured environments and in particular in the case of physical interaction/cooperation with humans. In [2], it has been shown how, due to the high reduction ratio of the gearboxes usually adopted in the transmission systems of robotic arms, the kinetic energy of the actuation system is predominant with respect to the one of the arm structure. Moreover, while the inertia of the robot arm can be reduced by adopting lightweight structures and composite materials, it is more difficult to reduce the inertia of the actuators and of the gearboxes because of technological limitations. Due to these considerations, Bicchi and Tonietti [2] proposed to adopt variable stiffness actuation to make the robot more rigid when the system velocities (and then kinetic energy) are low and to make the robot more compliant (or soft) during fast motions (that means high kinetic energy) so that the transmission system may store at least part of the kinetic energy of the actuation system in case unpredictable impacts occur. On the other hand, the kinetic energy of the robotic arm structure can be absorbed during impacts by suitable soft covers mounted on the robot [3], emulating in this way the function and the behavior of the human skin. The research effort in this field is motivated by the potential benefits given by the introduction of robots in human everyday activities, such as assistance to elder or handicapped people, homework activities, and entertainment. It is clear that purposely designed robots are needed to accomplish this kind of tasks.

Although several prototypes of single degree-of-freedom (DOF) joints with relevant and adjustable compliance for robotic manipulators have been designed and developed [48] with the aim of verifying the effectiveness of the variable stiffness approach, the implementation of a robotic arm with variable stiffness actuation is still an issue, and only recently, the German Aerospace Center (DLR) has started the development of such a robot. From the mechanical point of view, the main drawbacks of variable stiffness joints are a more complex mechanical design and the necessity of a pair of actuators to adjust simultaneously both the joint position and stiffness. It follows that a more complicated modeling and control approach is necessary, especially when kinematic structures with multiple DOF are considered [9]. In [10], a general approach to the feedback linearization and the simultaneous decoupled control of both the position and stiffness in a serial manipulator actuated by means of variable stiffness joints is proposed, while in [11], a robust controller for the compensation of the model uncertainties in variable stiffness robotic manipulators is presented.

In this paper, the control of robotic manipulators with variable joint stiffness is addressed in the general case of a redundant arm by means of a simplified (with respect to previous ones) output-feedback control approach. The proposed control is configured as a cascade system composed by three independent controllers, one for each part in which the dynamics of the robotic arm with variable stiffness joint can be subdivided. A singular perturbation analysis is then applied to the resulting system to show the separation of the dynamics of the different components, that are the arm, the positioning actuators, and the stiffness actuators, obtaining in this way a reduced order model of the manipulator and then a simplified control structure. The proposed control approach is then validated by means of simulations of a three-link planar manipulator with variable stiffness actuation. In this paper, we refer to a variable stiffness joint as a robotic joint actuated by means of a variable stiffness actuator.

2. System Dynamics

Starting from the general dynamic model of a robot with 𝑛 variable stiffness joints as proposed in [10], a second-order dynamic model is used to describe the stiffness behavior of the joints (or, more in general, of the transmissions), also in view of the nature of the mechanisms used to adjust the joint stiffness [6, 10].

A scheme of the working principle of a variable stiffness joint is shown in Figure 1: the robot link (on the right side) is moved by means of a positioning actuator (on the left side), and these two parts are connected by an elastic element modulated by means of a second actuator, also referred to as the stiffness actuator (in the lower center part). In [10], it has been shown how the conceptual design shown in Figure 1 can be also used to describe the behavior of antagonistic variable stiffness joints, such as the VSA by Bicchi et al. [12], by means of a suitable change of coordinates in the dynamic model of the device. The overall model of a robotic manipulator equipped with this variable stiffness actuators is then composed by the dynamics of an 𝑛 DOF arm driven by 𝑛 positioning actuators through elastic elements and by 𝑛 stiffness actuators that modify the characteristic of the couplings that connect the positioning actuators to the links. It is assumed that each joint has its own positioning and stiffness actuator and that no direct coupling between the actuators of different joints is present. Let 𝑞,𝜃,𝑘𝑛 be, respectively, the generalized joint positions, the positions of the the positioning actuators, and the vector of the joint stiffness. Under the simplifying modeling assumption used in [13] (viz., the kinetic energy of the actuators are due only to their own spinning), the dynamic model of the system can be written as 𝑀(𝑞)̈𝑞+𝑁(𝑞,̇𝑞)+𝐾[𝑞𝜃]+𝜂𝑞=𝜏𝑒,(1)𝐵̈𝜃+𝐾[𝜃𝑞]+𝜂𝜃=𝜏𝜃,(2)𝛾(𝑞,𝜃,𝑘)̈𝑘𝛽(𝑞,𝜃,𝑘)+𝜂𝑘=𝜏𝑘,(3)

Figure 1: Working principle scheme of a variable stiffness joint.

where 𝑀(𝑞) is the inertia matrix of the robotic arm, the vector 𝑁(𝑞,̇𝑞) contains the centrifugal, Coriolis, and gravity forces, 𝐾=diag{𝑘1,,𝑘𝑛}>0 is the joint stiffness matrix, 𝐵=diag{𝑏1,,𝑏𝑛}>0 is the inertia matrix of the actuators, 𝜏𝜃,𝜏𝑘𝑛 are the positioning and the stiffness actuators command torques, respectively, 𝜏𝑒 is the joint torque contribution due to the effect of an external load; that is,𝜏𝑒=𝐽(𝑞)𝑇𝐹𝑒,(4)

where 𝐽(𝑞) is the manipulator Jacobian matrix and 𝐹𝑒 is the external wrench. Equation (3) covers a situation where the stiffness actuation modifies the transmission configuration in order to change its compliance, for example, by precompressing elastic elements or by moving some mechanical parts. As shown in [10], the general model (1)–(3) can be used to describe the characteristics of different variable stiffness joint implementations both in the antagonistic case, such as the VSA proposed by Bicchi et al. [12], and in the nonantagonistic case, such as the device proposed by Wolf and Hirzinger [6], by means of a suitable change of coordinates in the dynamic model of the system. In (3), 𝛾(𝑞,𝜃,𝑘) and 𝛽(𝑞,𝜃,𝑘) are used to model the stiffness variation as a function of the device configuration, and they depend on the particular implementation of the robot transmission system; refer to [10] for a detailed description of the general model (1)–(3). Note that 𝐾=diag{𝑘1,,𝑘𝑛} and 𝑘=[𝑘1,,𝑘𝑛]𝑇 are alternative representations of the joint stiffness. It is important to stress that in this paper, the joint stiffness matrix 𝐾 in (1)-(2) is not constant and that, in general, is a function of time 𝐾=𝐾(𝑡),(5)

as expressed by (3), that models the dynamics of the joint stiffness variation. Moreover, we consider 𝑘𝑖(𝑡)>0,forall𝑡, since it has no physical meaning to consider negative stiffness, while if the stiffness drops to zero, an unactuated system would be obtained.

As a nonrestrictive assumption, all the effects due to frictions, dead zones, nonmodeled dynamics, parameters variations, and so forth, not considered in other terms of the system dynamics (1)–(3), are collected in the functions 𝜂{𝑞,𝜃,𝑘}(𝑡). It is then possible to refer to the nominal system dynamics by considering 𝜂{𝑞,𝜃,𝑘}(𝑡)=0.

As stated before, the overall system is composed by 3𝑛 rigid bodies, and therefore, the state vector of the system can be defined as 𝜒=𝑞𝑇̇𝑞𝑇𝜃𝑇̇𝜃𝑇𝑘𝑇̇𝑘𝑇𝑇6𝑛,(6)

while it is assumed that only the state subset information 𝜁𝑚=𝑞𝑇𝜃𝑇𝑘𝑇𝑇3𝑛(7)

are directly measurable, typically by means of optical encoders. Note that the joint stiffness 𝑘 is usually reconstructed by means of the knowledge of the elastic elements characteristics and by measuring the actuators and joints positions [9].

In any case, the control objective is the simultaneous and independent regulation of the following set of outputs: 𝜁=𝑞𝑘2𝑛,(8)

namely, the joint positions (and thus, through the robot direct kinematics, the end-effector pose) and the joint stiffness, acting on the system control inputs 𝜌=𝜏𝜃𝜏𝑘2𝑛,(9)

by means of an output-feedback controller based on the knowledge of 𝜁𝑚. As general assumption that takes into account the physical limitations of the system, the state variables are considered bounded; that is, 𝜒<𝐗 for some properly defined norm .

3. Control Strategy

In the problem considered here, the main control goal is to achieve the tracking of suitable workspace trajectories assigned to the manipulator end-effector while maintaining a desired level of apparent stiffness. As stated previously, this allows to ensure a harmless behavior of the manipulator in case of collision, in particular with humans.

The term apparent stiffness requires a more detailed discussion. Usually, active stiffness control refers to the case in which the manipulator stiffness is modified by means of the control only and has been introduced in [14]. A more general and notable example of active regulation of both static and dynamic parameters of the manipulator is probably the impedance control [15] that has been successfully used for the control of lightweight robots [16] and for safety evaluation of manipulators [17]. On the other hand, passive stiffness control indicates the possibility of changing the mechanical stiffness of the system. Note that in the case of series-elastic actuators [18, 19], the stiffness is fixed by design and cannot be changed anymore without using an active stiffness control approach. In this paper, the discussion is posed in terms of apparent stiffness, because the desired workspace stiffness is achieved by modulating both the mechanical compliance of the joints, that is, the passive stiffness, and the compliance due to the control. This principle is similar to the combined compliance control [20], and it has been already used also in robotic hands control; see [21].

An important point to note is that in case of an unpredictable impact, only the mechanical stiffness of the manipulator is experienced by the environment (or the human) during the very short time interval in which the energy transfer due to the impact occurs [17]. This is due to the limits in both the sensor and the control bandwidth that prevent an immediate reaction of the control system. Another motivation of this fact is the limited stiffness of the robot structure itself (if compared with the time constant of the force transmission along the structure) and to the fact that the impact may occur in any point along the manipulator, while the torque sensors, when available, are located in the joints. This latter effect, as, in general, the energy exchange from the robot to the environment during impacts, can be attenuated by means of suitable soft covers mounted on the robots [3] that allow early transmission of relatively lower forces along the structure together with some storage and dissipation of the energy due to the impact.

Recalling that variable stiffness actuation has been introduced for safety purposes, and that its main role is to absorb at least part of the kinetic energy of the actuation system (i.e., usually predominant with respect to the one of the arm structure because of the high reduction ratio gearboxes used to connect the actuators to the links) during unpredictable impacts, another important point to note is that the selection of suitable values of the joint stiffness depends mainly on the task to be performed (in particular on the required velocities), on both the manipulator and the actuators inertia, that together with the task specifications determines the overall system kinetic energy, and, the most important, on the criteria used to evaluate the robot safety. To this end, the definition of suitable safety criteria for robotic applications represents a key topic. While in early investigations these criteria have been inherited from automotive industries, it has been successively shown that they are not completely suitable for robotic applications because of the very different scenario [17]. Therefore, great attention has been devoted to the definition of adequate principles for the measurements and the quantification of the robot safety [17]. It is not the topic of this paper to deal with the selection of the joint stiffness: it is here assumed that the joint stiffness trajectory is given by a suitable safety planner. An overview of the complete system is depicted in Figure 2.

Figure 2: Scheme of the overall output-based controller.
3.1. Preliminaries

Consider the generic second-order dynamic system 𝐴(𝑣,𝑤)̈𝑣+𝑓(̇𝑣,𝑣,𝑤)+𝛼=𝜏,(10)

where 𝑣,̇𝑣 are the state variables, 𝑣 is the measurable output, 𝑤 is an exogenous measurable signal, 𝐴(,) is a bounded symmetric positive definite (invertible) matrix, 𝑓(,,) is bounded and Lipschitz in all its arguments, and 𝛼 is a generic function of the time collecting all the disturbances and uncertainties of the system. It is assumed that 𝐴1(𝑣,𝑤)[𝑓(̃𝑣,𝑣,𝑤)𝑓(̇𝑣,𝑣,𝑤)]𝜎̃𝑣̇𝑣,𝜎>0.(11)

Due to the general separation principle introduced in [22] for SISO systems and then extended to MIMO systems in [23], semiglobal stability (i.e., imposing that the domain of attraction of the equilibrium contains a prescribed compact set, [24]) of the desired equilibrium of (10) in nominal conditions (𝛼=0) can be achieved by output feedback, since (1) state feedback linearization is clearly possible under the previously defined hypotheses and (2) the system is uniform complete observable. While the verification of the former point is trivial, the latter can be easily proven by dynamic extension, [25]. It is straightforward to show that (10) is a generic form that can be used to represent (1)–(3) (see Property 5 in the appendix). Then, the following theorems hold.

Theorem 1. The system (10) in nominal conditions (𝛼=0) and under the effect of the output-based dynamic compensator 𝜏=𝐴(𝑣,𝑤)𝜏+Γ3𝑧3𝑧2+𝑓𝑧2,𝑣,𝑤,(12)𝜏=Λ1𝑣𝑣𝑑Λ2𝑧2̇𝑣𝑑+̈𝑣𝑑,(13)̇𝑧1=Γ1𝑧1𝑣+𝑧2,(14)̇𝑧2=Γ2𝑧1𝑣+𝜏,(15)̇𝑧3=𝜏,(16) where 𝑣𝑑 is the desired value of 𝑣, Λ1,Λ2,Γ1,Γ2,andΓ3 that are positive definite (diagonal) matrices and has a unique steady-state exponentially stable equilibrium point with arbitrary fast error dynamics.

Proof. By defining the state vector composed by the system tracking error 𝜀1 plus the compensator estimation error 𝜀2 as 𝜀=𝑣𝑣𝑑𝑇̇𝑣̇𝑣𝑑𝑇|𝑧1𝑣𝑇𝑧2̇𝑣𝑇𝑧3̇𝑣𝑇𝑇=𝜀𝑇1|𝜀𝑇2𝑇,(17) then (10)–(16) can be written as ̇𝜀=𝐻𝜀+Θ̂𝛼,(18) where 𝐻=0𝐼000Λ1Λ20Λ2Γ3Γ300Γ1𝐼000Γ2Γ3Γ3000Γ3Γ3=𝐻11𝐻120𝐻22,(19)Θ=0𝐼0𝐼𝐼𝑇=Θ𝑇1Θ𝑇2𝑇,̂𝛼=𝐴1𝑓(̇𝑣,𝑣,𝑤)𝑓𝑧2,𝑣,𝑤.(20) The arguments of 𝐴(,) are omitted for brevity.
At first, the autonomous linear system obtained by neglecting ̂𝛼 is considered. The eigenvalues of matrix 𝐻, as can be easily seen by looking at the partitioning in (19), can be freely selected by means of a suitable choice of the parameters Λ1,Λ2,Γ1,Γ2, and Γ3.
Due to the linearity of the nominal system under the effect of the proposed controller, the origin is the unique exponentially-stable equilibrium point of the autonomous (nominal) system and the rate of convergence can be made arbitrarily fast.
Now, the implications of the previous hypothesis (11) are considered: on this basis and from (20) it is possible to write 𝐴1𝑓𝑧2,𝑣,𝑤𝑓(̇𝑣,𝑣,𝑤)𝜎𝑧2̇𝑣,𝜎>0.(21) The effects of 𝑓() on (18) can be seen as a vanishing perturbation of the nominal system, then the exponential stability of the equilibrium point can be recovered by means of a suitable choice of the controller parameters. In particular, by considering a symmetric positive definite matrix 𝑄, it is well known that a unique symmetric positive definite matrix 𝑃 exists that satisfies the Lyapunov equation 𝑃𝐻+𝐻𝑇𝑃=𝑄.(22) If and only if the matrix 𝐻 has all negative eigenvalues [26]. The Lyapunov function 𝑉(𝜀)=𝜀𝑇𝑃𝜀 satisfies the relations 𝑐1𝜀2𝑉(𝜀)𝑐2𝜀2,𝜕𝑉𝜕𝜀𝐻𝜀=𝜀𝑇𝑄𝜀𝑐3𝜀2,𝜕𝑉𝜕𝜀=2𝜀𝑇𝑃2𝑃𝜀=𝑐4𝜀,(23) where 𝑐1=𝜆min(𝑃), 𝑐2=𝜆max(𝑃), 𝑐3=𝜆min(𝑄), and 𝑐4=2𝜆max(𝑃). The derivative of 𝑉(𝜀) along the trajectory of the system satisfies ̇𝑉(𝜀)𝑐3𝜀2+𝑐4𝜎𝜀2,(24) meaning that the system is exponentially stable if 𝜎<𝑐3/𝑐4. It is possible to show that this ratio can be maximized by choosing 𝑄=𝐼 [26].

Theorem 2. The system (18) under the effects of nonvanishing bounded disturbances 𝛼̇𝜀=𝐻𝜀+Θ̂𝛼+𝐴1𝛼(25) is ultimate uniformly bounded (UUB), [26].

Proof. Theorem 1 shows that the system (18) has the origin 𝜀=0 as unique exponentially stable equilibrium point. It is here supposed that 𝜀(𝑡)<𝜍,forall𝑡0. By means of a suitable choice of the control parameters, it possible to ensure that the following condition holds: 𝐴1𝛼𝛿<𝑐3𝑐4𝑐1𝑐2𝜗𝜍,(26) for some positive constant 𝜗<1, where 𝛿 is a positive constant representing a suitable bound of the disturbances. Taking into consideration the same Lyapunov function adopted in Theorem 1, the derivative of 𝑉(𝜀) along the trajectory of the perturbed system satisfies ̇𝑉(𝜀)𝑐3𝜀2+𝑐4𝜎𝜀2+𝑐4𝛿𝜀=[1𝜗]𝑐3𝜎𝑐4𝜀2𝜗𝑐3𝜎𝑐4𝜀2+𝑐4𝛿𝜀[1𝜗]𝑐3𝜎𝑐4𝜀2,𝜀𝑐4𝛿𝜗𝑐3𝜎𝑐4.(27) Then, for all 𝜀(𝑡0)<𝑐1/𝑐2𝜍, the trajectory of the system satisfies the relations 𝜀(𝑡)𝑐2𝑐1exp[1𝜗]𝑐3𝜎𝑐4𝑐4𝑡𝑡0𝜀𝑡0,𝑡0𝑡<𝑡0+𝑇,𝜀(𝑡)𝑐4𝑐3𝜎𝑐4𝑐2𝑐1𝛿𝜗,𝑡𝑡0+𝑇,(28) for some finite time interval 𝑇, see [26] for additional details on the proof.

Property 1. Due to the structure of 𝐻 and to the partitioning of 𝜀, the eigenvalues of the dynamic compensator (i.e., of matrix 𝐻22) depend only on Γ1,Γ2 and Γ3, while those of the controlled system (i.e., of 𝐻11) depend on Λ1 and Λ2. This separation property is a consequence of the system linearization obtained by output feedback.

Property 2. Note that (16) represents the internal model of the system, and by assuming constant values of 𝛼, namely, 𝛼ss, the steady-state equilibrium point of the system is 𝜀ss=𝐻1Θ𝛼ss=0000Γ1𝑇3𝑇𝛼ss,(29) thus resulting that only 𝑧3 is affected by the external disturbance.
Generally speaking, by applying a singular perturbation analysis, from Property 1, and by assuming that the eigenvalues of 𝐻22 are fast enough with respect to those of 𝐻11, it is possible to assume that the system ̇𝜀2=𝐻22𝜀2+Θ2̂𝛼+𝐴1𝛼(30) is always working in steady-state conditions, and its equilibrium point is 𝜀2=𝐻122Θ2𝐴1𝛼ss.(31) This last relation is a generalization of (29). It follows that 𝑧1𝑣,𝑧2̇𝑣,𝑧3𝑧2Γ13𝐴1𝛼ss(32) Then, 𝑧3 (together with 𝑧2) provides a filtered estimation of the disturbance. This effect on the system can be seen as a disturbance compensation action (i.e., it compensates for the effects of 𝛼 on the system state), and this compensation action can be annihilated by posing Γ3=0 in (13)–(16). Note that from (19), this is a singular condition for the controller.

Property 3. The control system (10)–(16) can be used also as a regulator by neglecting the time derivative of the reference signal ̇𝑣𝑑 and ̈𝑣𝑑. This property will be exploited for the design of the fast actuators controllers, since, in this case, tracking performance are not necessary, and these controllers are assumed fast enough to work always in steady-state conditions with respect to the slow arm controller.

Property 4. By redefining the control action in (13) as 𝜏=Λ1𝑣𝑣𝑑Λ2𝑧2̇𝑣𝑑+̈𝑣𝑑+𝑢,(33) where 𝑢 is an auxiliary input signal, the system dynamics becomes ̇𝜀=𝐻𝜀+𝐺𝑢+Θ̂𝛼+𝐴1𝛼,𝐺=0𝐼000𝑇.(34) Due to the structure of matrix 𝐺, it is possible to note that the auxiliary input 𝑢 does not influence the compensator error dynamics 𝜀2, and therefore, the resulting dynamics of the tracking error 𝜀1 is ̇𝜀1=0𝐼Λ1Λ2𝜀1+0𝐼𝑢,(35) showing a steady state tracking error (𝑣𝑣𝑑)=Λ11𝑢.

3.2. Actuators Control Design

It is important to note that both (2) and (3) are in the form of (10) (see Property 5 in the appendix). This essentially means that both the positioning and the stiffness actuators dynamics can be decoupled from the arm dynamics (1) by exploiting the properties of the controller (12)–(16). Even if this decoupling can be achieved by means of a unique controller, the use of two separated controllers, one for the positioning actuators and one for the stiffness actuators, is chosen to have a more clear view of the overall system structure. Due to the implications of Theorem 1, the possibility of assigning an arbitrary dynamics to the actuators allows to obtain a two time-scale system, where the actuators represent the fast system and the arm the slow one, and then to apply a singular perturbation approach to reduce the order of the resulting system 𝑀(𝑞)̈𝑞+𝑁(𝑞,̇𝑞)+𝐾[𝑞𝜃]+𝜂𝑞=𝜏𝑒,̇𝜀𝜃=𝐻𝜃𝜀𝜃+Θ𝜃̂𝛼𝜃+𝐵1𝛼𝜃,(36)̇𝜀𝑘=𝐻𝑘𝜀𝑘+Θ𝑘̂𝛼𝑘+𝛾1𝛼𝑘,(37) where 𝜀𝜃=(𝜃𝜃𝑑)𝑇̇𝜃𝑇|(𝑧1𝜃𝜃)𝑇(𝑧2𝜃̇𝜃)𝑇(𝑧3𝜃̇𝜃)𝑇𝑇,𝜀𝑘=(𝑘𝑘𝑑)𝑇̇𝑘𝑇|(𝑧1𝑘𝑘)𝑇(𝑧2𝑘̇𝑘)𝑇(𝑧3𝑘̇𝑘)𝑇𝑇,(38)

and 𝜃𝑑 is the vector of the desired positioning actuators configuration, 𝑘𝑑=[𝑘1𝑑,,𝑘𝑛𝑑] is the vector of the desired joint stiffness, and ̂𝛼{𝜃,𝑘}, 𝛼{𝜃,𝑘} summarize the disturbances acting on the positioning and stiffness actuators, respectively. In the case of the actuators controllers, Property 3 has been exploited, since, with respect to the arm dynamics, these controllers are always working in steady-state conditions, then tracking performance are not necessary but just the regulation of the desired position and stiffness. Note that the robotic arm dynamics remains the same as in (1). On the basis of Theorem 1, for both (36) and (37), it is possible to write 𝐻1̇𝜀=𝜀+𝐻1Θ̂𝛼,(39) where the distinction between the positioning and the stiffness actuators has been omitted for brevity. As long as the eigenvalues of 𝐻 can be freely assigned, it is clear how the time scale of (36) and (37) can be made arbitrarily faster than the time scale of (1); therefore, (39) represents the fast manifold dynamics. With some abuse of notation, as 𝐻10 (remember that this can be achieved by a suitable choice of the controller parameters) the fast manifold dynamics (39) degenerates, and due to Theorem 1, it is possible to state that the fast manifold has a unique root expressed by (29). From (29)–(37), since the disturbance does not affect the coupling between the actuators and the arm, it is then possible to define the reduced-order model of the robotic manipulator with variable joint stiffness as 𝑀(𝑞)̈𝑞+𝑁(𝑞,̇𝑞)+𝐾𝑑𝑞𝜃𝑑+𝜂𝑞=𝜏𝑒,(40) where the actual actuator configuration 𝜃 and 𝐾 have been substituted with the desired ones 𝜃𝑑 and 𝐾𝑑 in the arm dynamic equation.

3.3. Trajectory Planning for the Actuators

In this section, the problem of defining the desired trajectories for the actuators given the manipulator trajectories, both in terms of position and stiffness, is addressed. Note that in case the desired manipulator stiffness is defined in the joint space only, the stiffness actuator reference signals are fully defined by 𝐾𝑑, then the more general case of workspace stiffness reference is here addressed.

It is important to note that once the desired joint stiffness 𝐾𝑑 has been determined, the torque provided to the joints by the actuators can be regulated by acting on 𝜃𝑑, that is, by a suitable selection of the positioning actuators configuration. To this end, (40) can be rewritten as 𝑀(𝑞)̈𝑞+𝑁(𝑞,̇𝑞)+𝜂𝑞=𝜏𝑒+𝜏act,(41)𝜏act=𝐾𝑑𝜃𝑑𝑞,(42)

where 𝜏act is the vector of the torques provided by the arm virtual actuators obtained from the joined action of the stiffness and the positioning actuators. It follows that 𝜃𝑑=𝐾1𝑑𝜏act𝑑+𝑞,(43)

where 𝜏act𝑑 is the desired torque to be applied on the arm joints by the actuation system and computed by the arm controller that will be defined in the following.

3.3.1. Position Trajectory

Let us define the direct kinematic problem and the manipulator Jacobian as 𝑝=Ψ(𝑞)̇𝑝=𝜕Ψ(𝑞)𝜕𝑞̇𝑞=𝐽(𝑞)̇𝑞,(44)

where 𝑝𝑚 is the manipulator end-effector position in an 𝑚-dimensional task space. Differently from the case of state feedback linearization [10], where 𝑞𝑑4 and 𝑘𝑑2, that is, the joint trajectories must be defined up to the fourth time derivative, while the joint stiffness must be defined up to the second time derivative, in this case, due the reduced order of the arm dynamics (40), the desired manipulator trajectory can be specified as usual for rigid robots both in the task space and in the joint space. In other words, the robot trajectory will be defined in terms of desired positions, velocities, and accelerations, while no particular requirements are necessary for the stiffness desired values. In any case, smooth-enough stiffness profiles are preferable to avoid abrupt oscillations of the actuators and of the arm.

If the manipulator trajectory is defined in the task space, the vectors 𝑝𝑑, ̇𝑝𝑑, and ̈𝑝𝑑 indicate, respectively, the desired position, velocity, and acceleration of the end effector. It is assumed that both task- and joint-space trajectories are limited and continuous (at least with a piece-wise continuous acceleration profile). Obviously, in case of task space control, singular configurations of the manipulator must be avoided or explicitly considered in the controller. In that case, it is then possible to define the desired joint-space trajectory as 𝑞𝑑=Ψ1𝑝𝑑,(45)̇𝑞𝑑=𝐽𝑞𝑑#̇𝑝𝑑,̈𝑞𝑑=𝐽𝑞𝑑#̈𝑝𝑑̇𝐽𝑞𝑑̇𝑞𝑑,(46) where (45) indicates the inverse kinematic function. The solution of the manipulator inverse kinematics can be also obtained by means of a closed-loop inverse kinematic algorithm (CLIK), [27], as specified in the following. In case of manipulator redundancy (𝑛>𝑚), the Jacobian pseudoinverse 𝐽(𝑞𝑑)# must be considered. Note that the manipulator redundancy can be exploited by defining secondary tasks [28], for example, to avoid obstacles, to deal with general task constraints, or to optimize the manipulator joint configuration on the basis of the desired task space stiffness, as will be shown in the following.

3.3.2. Stiffness Reference Signals

In the following discussion, the leading superscript 𝑤 will be used to refer to workspace-related variables. In static conditions, the workspace stiffness of the arm 𝑤𝐾𝑑 can be specified as 𝐹𝑒=𝑤𝐾𝑑𝑑𝑝,(47) where the infinitesimal end-effector workspace displacement 𝑑𝑝 and the environmental force 𝐹𝑒 are considered.

Several works can be found in the literature about the workspace stiffness control of both redundant and nonredundant robots. In particular, in [29], the different components that determine the total workspace stiffness are analyzed, and the orthogonal stiffness decomposition control is introduced, while in [30], the stability properties of the workspace stiffness components and a minimal parametrization of the manipulator stiffness and compliance are reported. Since variable stiffness actuation allows the independent control of the joint compliance, an interesting approach for the regulation of the workspace compliance of redundant robots with compliant actuation can be the combined control approach proposed in [20], where the desired manipulator stiffness in the operational space is achieved by combining a decoupled stiffness joint control together with a workspace stiffness controller.

The relation between the workspace stiffness 𝑤𝐾𝑑 and the overall joint stiffness 𝐾𝑎 can be expressed as 𝐾𝑎=𝐽𝑇𝑤𝐾𝑑𝐽,(48) where the dependence of the Jacobian from the joint configuration 𝑞 has been omitted for brevity. Note that the overall joint stiffness 𝐾𝑎 may be the result of several effects, as specified in the following. For this reason, this quantity has been distinguished from the mechanical joint stiffness 𝐾𝑑. In the right-hand side of (48), since the workspace stiffness must be, in general, symmetric and positive definite and due to the transformation through the manipulator Jacobian, only 𝑚[𝑚+1]/2 independent terms are present. At the same time, by considering the case of independent regulation of the joint mechanical stiffness only, matrix 𝐾𝑎 is diagonal (𝐾𝑎=𝐾𝑑), then there are 𝑛 independent terms in the left-hand side of (48). It follows that it is possible to assign arbitrary workspace stiffness by selecting suitable joint mechanical stiffness only if 𝑛=𝑚[𝑚+1]/2, that means a certain redundancy in the manipulator, at least nonconsidering the trivial case 𝑚=1, that is, a single DOF mechanism. If 𝑛>𝑚[𝑚+1]/2, there exist infinite combinations of joint stiffness that allow to achieve the desired workspace stiffness. Then, this redundancy must be managed by suitable projection of the workspace stiffness on the joint space [29], and it can be exploited maybe to minimize (or maximizing depending on the task) the overall joint stiffness. If 𝑛<𝑚[𝑚+1]/2, it is not possible to control all the components of the workspace stiffness be means of the joint stiffness only, and suitable couplings between the components of the workspace stiffness must be introduced.

It is now important to remark that the introduction of variable stiffness joints is useful to reduce the energy transferred from the manipulator to the environment (or humans) during unpredictable impacts by decoupling the actuators inertia from the link during fast robot motions (i.e., by reducing the joint stiffness) or, on the other hand, to increase the precision of the manipulator during slow and limited movements by imposing high joint stiffness values to allow adequate external disturbance force/torque rejection. Moreover, variable stiffness actuators are usually designed to have high compliance in rest conditions and to become rigid if necessary under the action of the control system. In this sense, it is then clear that the role of the joint mechanical stiffness, or passive stiffness, is completely different from the active stiffness, that is a control strategy that to allow to rigid robots to show a “compliant’’ behavior with respect to environmental uncertainties, for example, during assembly operations [14].

By introducing the concept of combined compliance control reported in [20] but with different goals, it is possible to split the overall workspace stiffness into two components, one given by the mechanical joint stiffness 𝐾𝑑 only, the passive stiffness, that acts at all the frequencies but that shows its major effects at high frequency and then in case of impacts, and in a component due to the control, the active stiffness 𝐾𝑐, that is significant only within the control bandwidth. It is also important to remark that due to the cascade structure of the proposed controller, it is possible to design the arm controller with a purposely defined bandwidth to decouple the passive stiffness from the active stiffness.

The contribution to the workspace stiffness due to the mechanical joint stiffness 𝐾𝑑 only, namely, 𝑤𝐾𝑗, is given by 𝑤𝐾𝑗1=𝐽𝐾1𝑑𝐽𝑇.(49)

It is important to note that (49) maps the joint stiffness to the workspace stiffness no matter what is the number of joints or the dimension of the workspace. This point is important in case of manipulator redundancy to deal with suitable selection of the joint stiffness.

On the basis of these considerations, the workspace stiffness 𝑤𝐾𝑑 can be partitioned into two different contributions 𝑤𝐾𝑑=𝑤𝐾𝑗+𝑤𝐾𝑐,(50)

where 𝑤𝐾𝑗 is the contribution due to the passive mechanical joint stiffness and 𝑤𝐾𝑐 is the active stiffness due to the control. Note that since from (49) 𝑤𝐾𝑗 is symmetric positive definite by definition and 𝑤𝐾𝑐 must be symmetric and positive definite to avoid instability of the positioning actuators, also 𝑤𝐾𝑑 is positive definite, thus avoiding system instability: this results in a lower limit of the workspace stiffness given by the effects of the joint mechanical stiffness, since the controller can be used only to increment the overall stiffness. It is important to remark that while in the case of industrial robots the control objective is to make rigid robot softer and then the system compliance is given by the sum of the control compliance and the joint compliance, this last equation shows that the proposed control system is used to the make soft robots (because of the relatively low value of 𝑤𝐾𝑗) stiffer when necessary by means of control (i.e., by acting on 𝑤𝐾𝑐).

From (49) and (50), the induced partitioning of the apparent joint stiffness 𝐾𝑎 is 𝐾𝑎=𝐾𝑑+𝐽𝑇𝑤𝐾𝑐𝐽=𝐾𝑑+𝐾𝑐.(51)

Due to the effect of the joint apparent stiffness, the relation between the infinitesimal joint displacements and the joint torques can be expressed as 𝜏act𝑑=𝐾𝑎𝑑𝑞.(52)

Taking into consideration the relation between the workspace displacement and joint displacement, by inverting (52), it is possible to write 𝑑𝑞=𝐾1𝑎𝜏act𝑑=𝐾1𝑎𝐽𝑇𝐹𝑒=𝐾1𝑎𝐽𝑇𝑤𝐾𝑑𝑑𝑝=𝐾1𝑎𝐽𝑇𝐽𝐾1𝑎𝐽𝑇1𝑑𝑝,(53)

where the relation between the workspace stiffness and the joint apparent stiffness 𝑤𝐾𝑑1=𝐽𝐾1𝑎𝐽𝑇 has been exploited (see Theorem 3 in the appendix). It is possible to demonstrate that the same relation can be written in the form (Theorem 3 in the appendix) [20] 𝑑𝑞=𝐾1𝑑𝐽𝑇𝐽𝐾1𝑑𝐽𝑇1𝑑𝑝=𝐽#𝑑𝑝,(54)

where 𝐽#=𝐾1𝑑𝐽𝑇[𝐽𝐾1𝑑𝐽𝑇]1 is the pseudoinverse Jacobian matrix, weighted by the joint mechanical stiffness 𝐾𝑑. This result is important, because it shows that the relations between workspace and joint displacements are affected only by the mechanical joint stiffness. From the meaning of the pseudoinverse, the inverse kinematic solution obtained from (54) can be considered as providing the manipulator configuration which minimizes the deviation of the joints weighted by the joint stiffness; that is, it minimizes the potential energy stored in the joint compliance, ensuring the minimum energy transfer due to the inertia of the actuators during unexpected collisions. It follows that the desired workspace stiffness has no effect on the inverse kinematics solution of the manipulator. This fact suggests that the regulation of the manipulator configuration 𝑞 must be performed on the basis of the inverse kinematic solution obtained by means of 𝐽# both to stabilize the manipulator dynamics in the Jacobian null space and to achieve minimum energy transfer from the actuators to the arm.

3.4. Arm Control Design

Several robust controllers for robotic arms have been presented in the literature; see for example, [3133] for some well-known examples. Anyway, in case of robots with variable joint stiffness, the application suggests a particular implementation of the controller.

First of all, aiming to clarify the role of the passive and the active stiffness, the static relation between the environmental force 𝐹𝑒 and the deviation of the end effector in the workspace is considered 𝐹𝑒=𝑤𝐾𝑑𝑝𝑑𝑝𝜏act𝑑=𝐽𝑇𝐹𝑒=𝐽𝑇𝑤𝐾𝑑𝑝𝑑𝑝,(55) and for small-enough displacements, 𝜏act𝑑=𝐽𝑇𝑤𝐾𝑑𝐽𝑞𝑑𝑞=𝐾𝑎𝑞𝑑𝑞.(56) From (43) and (51), it follows that 𝜃𝑑=𝐾1𝑑𝐾𝑎𝑞𝑑𝑞+𝑞=𝐾1𝑑𝐾𝑑+𝐾𝑐𝑞𝑑𝑞+𝑞=𝑞𝑑+𝐾1𝑑𝐾𝑐𝑞𝑑𝑞.(57)

By substituting (57) into (41)-(42), it follows that𝑀(𝑞)̈𝑞+𝑁(𝑞,̇𝑞)+𝐾𝑑+𝐾𝑐𝐾𝑎𝑞𝑞𝑑+𝜂𝑞=𝜏𝑒,(58) that means that the robotic arm behaves as a joints-position controlled system with the desired apparent stiffness.

Looking at the right hand side of (57), the first term represents the passive stiffness, since the positioning actuators are commanded just on the basis of the manipulator joint desired position trajectory, and then, no feedback from the actual joint positions is present, while the second term represents the active stiffness, since it implies the joint positions feedback on the positioning actuators. It is important to remark that while the passive stiffness acts at high frequency since it does not require the intervention of any controller, the active stiffness acts within the arm controller bandwidth. Moreover, if the manipulator configuration is obtained by the inverse kinematic solution by using (54), the minimum value for the potential energy stored in the joints as a consequence of the end-effector displacement is achieved.

With the aim of controlling the workspace motion of the manipulator, a specialized version of the controller (12)–(16) that in particular exploits Property 4 is now applied to the arm reduced model (41)–(43) (see Property 6 in the appendix). The torque of the virtual actuators can be then computed as 𝜏act𝑑=𝑀(𝑞)𝐽#𝐹̇𝐽𝑧2+𝜏𝑗+𝑁𝑞,𝑧2𝐽𝑇𝐹𝑒𝑚+𝑀(𝑞)Γ3𝑧3𝑧2,(59)𝐹=𝑤𝑀𝑑1𝑤𝐾𝑑𝑝𝑝𝑑𝑤𝐷𝑑𝐽𝑧2̇𝑝𝑑+𝑤𝑀𝑑̈𝑝𝑑+𝐹𝑒,(60)𝜏𝑗=𝑆𝐾𝑑𝑞𝑞𝑑𝐷𝑑𝑧2̇𝑞𝑑,(61)̇𝑧1=Γ1𝑧1𝑞+𝑧2,(62)̇𝑧2=Γ2𝑧1𝑞+𝐽#𝐹̇𝐽𝑧2+𝜏𝑗,(63)̇𝑧3=𝐽#𝐹̇𝐽𝑧2+𝜏𝑗,(64) where 𝑆=𝐼𝐽#𝐽 is the matrix that projects the joint space control action in the null space of the manipulator Jacobian, 𝐷𝑑=2𝐾𝑑 is the joint space damping, 𝑤𝐷𝑑 and 𝑤𝑀𝑑 represent the desired workspace damping and mass of the manipulator, respectively, and Γ1,Γ2,Γ3 are positive definite (diagonal) matrices. In particular, 𝑞𝑑 and ̇𝑞𝑑 are provided by a closed-loop inverse kinematic solver that exploits (54), while 𝐹𝑒𝑚 and 𝐹𝑒 are the external force measurement and estimation, respectively. On the basis of Property 4, the external force estimation is used as the auxiliary input 𝑢 appearing in (34) in the arm controller. Note that on the basis of (54), the static term 𝐾𝑑[𝑞𝑞𝑑] in the joint space control ensures the minimum norm of the joint space control effort.

It is initially assumed that the force 𝐹𝑒 applied by the end effector to the environment is measurable (i.e., by means of a 6-axis force/torque sensor), then it is possible to impose 𝐹𝑒𝑚=𝐹𝑒=𝐹𝑒. Note that by recalling the boundedness of the arm trajectory and the limited value of the state variables and in particular of both the joint velocity ̇𝑞 and of the controller state 𝑧2, the condition (11) is ensured by the general properties of the robotic arm dynamics [34, 35].

By recalling Property 2, and with a suitable choice of large-enough values of Γ{1,2,3}, the dynamic compensator represented by (62)–(64) can be configured as a third fast system together with the actuators controllers previously defined. Due to the controller structure, the arm response would be slower in two-order than the actuation control system. This fact does not represent a limitation, since (1) the fast control systems are very simple, and then, very fast implementations are possible with current technologies and computational capabilities of modern digital systems, (2) robots with variable stiffness actuation are conceived for the interaction with humans and to operate in a human-like environment, then smoother and bandwidth-limited motions are required with respect to industrial robots, (3) fast high-gain controllers are used in the actuators, while a slow low-gain controller is more suitable for the arm control also for safety reasons. It is then possible to assume that 𝑧1𝑞,𝑧2̇𝑞,𝑧3𝑧2Γ13𝛼,(65) where 𝛼=𝑀(𝑞)1𝜂𝑞 is the effect of the model and parameter uncertainties. Then, by denoting 𝑒=𝑝𝑝𝑑, the reduced-order nominal arm dynamics (41) under the effect of the controller (59)-(60) can be written as 𝑤𝑀𝑑̈𝑒+𝑤𝐷𝑑̇𝑒+𝑤𝐾𝑑𝑒=𝐹𝑒.(66)

In case a direct measure of the force applied by the end effector to the environment is not available, (𝐹𝑒𝑚=0) by assuming the perfect knowledge of the arm dynamics, that is, 𝜂𝑞=0, and by considering the external force 𝐹𝑒 as a disturbance, Property 2 implies that in steady state conditions, 𝑧1𝑞,𝑧2̇𝑞,𝑧3𝑧2Γ13𝑀(𝑞)1𝐽𝑇𝐹𝑒.(67) A suitable estimation of the external force 𝐹𝑒 can be then assumed as 𝐹𝑒=𝐽#𝑇𝑀(𝑞)Γ3𝑧3𝑧2.(68)

This equation defines a virtual force sensor that provides an indirect measure of the external force. Hence, the reduced-order nominal arm dynamics can be written as 𝑤𝑀𝑑̈𝑒+𝑤𝐷𝑑̇𝑒+𝑤𝐾𝑑𝑒=𝐹𝑒+𝑤𝑀𝑑𝐽𝑀(𝑞)1𝐽𝑇𝐹𝑒𝐹𝑒.(69)

The last term in (69) represents the perturbation on the system dynamics due to the use of the force estimation instead of a direct measure. Apart from this spurious term, the robotic manipulator is now configured as an impedance controlled system.

Once the torque 𝜏act𝑑 is determined by means of (59), the reference trajectory for the positioning actuators can be computed through (43).

An important remark is that while the workspace stiffness 𝑤𝐾𝑑 is obtained by combining the mechanical stiffness of the robot transmission together with the control action, only the active part of the desired stiffness implies the feedback of the manipulator joint positions. This is a key feature for safety reasons [2], because no reaction from the controller is needed to achieve the desired compliant behavior at joint level, since the compliance is a structural characteristic of the robot itself.

Note that also the arm controller provides a filtered estimation of the disturbance acting on the joints through (64). These information, also called residuals in the robotic community [36], can be used for fault evaluation and collision detection/reaction [37], and they have recently been successfully used for safe interaction of lightweight robots with human operators [17]. The possibility of exploiting the properties of robots with variable stiffness actuation together with the safety-oriented control techniques related to the use of the residuals justifies the proposed implementation of the controller.

3.5. Remarks

With respect to the control approaches for the system (1)–(3) proposed in [10, 11], the control scheme presented in this paper presents a number of remarkable advantages (i)The complexity of the controller is significantly reduced since neither static nor dynamic inversion of the system is needed. (ii)Redundancy of the manipulator is considered. (iii)The computation of the Lie derivative of the system dynamics up to the fourth order is avoided. (iv)An output-feedback control law has been defined instead of state-feedback controllers. (v)The arm trajectory must now be defined only up to the second order derivative, as in the rigid robot case, while previously, it was necessary to define the position trajectory up to the fourth order and the stiffness trajectory to the second order derivative. (vi)Since the arm control is based on the two signals 𝜃𝑑 and 𝐾𝑑, that are generated online and that act as exogenous reference for the actuators fast system, online compensation of disturbances is permitted. (vii)The joint stiffness is controlled by a specific control system that directly compensates for the effects of disturbances, external forces and couplings due to both positioning actuators and arm dynamics. (viii)The proposed control strategy exploits the mechanical compliance of the transmission to achieve the desired workspace behavior (passive compliance). (ix)In case of perfect knowledge of the arm dynamics, the system is able to estimate and regulate the workspace interaction force without a direct measurement of the force itself. (x)The manipulator residuals are directly provided by the controller, allowing a simple integration with fault and collision detection/reaction techniques. (xi)The extension of the proposed approach to the case of robotic manipulator where the stiffness can be adjusted only for some joints is straightforward, since by applying the proposed positioning and stiffness actuators controller where necessary, the overall robot dynamics is reduced to the rigid case.

4. Simulation of a Three-Link Planar Manipulator

The proposed controller has been validated by means of simulations of a planar three-link robotic arm with variable joint stiffness. Only the Cartesian coordinates 𝑥 and 𝑦 of the workspace end-effector position are considered (not the orientation), and this allows to analyze the effects of the manipulator redundancy in a quite simple manner. The extension of the proposed control approach to the case in which gravity compensation is needed is straightforward. The well-known dynamic model of the three-link planar arm is omitted for brevity; see [35], and only the dynamics of stiffness variation is specified 𝜆2̈𝑘𝑖+𝜆1𝑘2𝑖+𝜆0𝑞𝑖𝜃𝑖2+𝜂𝑘𝑖=𝜏𝑘𝑖.(70)

This stiffness dynamic model [10] is typical of a variable stiffness actuator with antagonistic elastic elements characterized by quadratic force-deformation relation proposed by Migliore et al. in [5]. The parameters of the three-link planar manipulator used in simulation are reported in Table 1, and all the links of the manipulator are considered identical for simplicity. In the simulation scheme, the trajectories are generated using proper filters to compute the position trajectory in the workspace together with its derivatives up to the second order. In this manner, the conditions given by Theorem 1 for the stability of the system are implicitly met. It is assumed that the joint stiffness reference values are provided by a suitable safety controller on the basis of some requirements given by the particular application and contest. Moreover, the matrices that characterize the arm impedance controller have been considered diagonal to simplify the analysis of the simulation results. The controller parameters adopted in the simulations are summarized in Tables 2 and 3. Three case studies have been considered: in the first case, whose results are reported in Figure 3, the manipulator end-effector is driven along several point-to-point trajectories. During this simulation, friction in the actuators is considered together with parameter uncertainties in the arm dynamics. External forces are applied along the 𝑥 and 𝑦 directions at time 𝑡=4[s] and 𝑡=15[s], respectively. From these plots, it is possible to note that as expected on the basis of Property 2, the end-effector position error is mainly affected by the external force, that is considered measurable during this test, and exactly related to the desired workspace stiffness. It is important to note that the joint mechanical stiffness is not significantly affected by the remaining parts of the system, then the independent regulation of the actuation stiffness is achieved. The joint stiffness tracking error is not reported, since it is very limited and not meaningful. In Figure 3, also the joint reference trajectory computed by the inverse kinematic solver is reported (upper left plot): this information is useful for the regulation of the null-space behavior of the system. Finally, also the positioning actuators trajectories are reported.

Table 1: Nominal parameters of the three-link planar manipulator.
Table 2: Parameters of the three-link planar manipulator controller used in the simulations.
Table 3: parameters of the positioning and stiffness actuators fast controllers used in the simulations.
Figure 3: Response of the system during test no. 1.

In Figure 4, the results of the case study no. 2 are reported. During this test, a different set of point-to-point trajectories have been evaluated together with an higher values of the desired workspace stiffness and larger (measurable) external forces. The same limitations on the maximum workspace velocity considered in the previous case are assumed during this test. Also, in this simulation, the external forces are applied along 𝑥 and 𝑦 at time 𝑡=4[s] and 𝑡=15[s], respectively. As in the previous case, the end-effector position error is determined by the external force and by the desired workspace stiffness, while the commanded joint stiffness is the same of the case no. 1, hence showing the independent control of the workspace and of the joints stiffness.

Figure 4: Response of the system during test no. 2.

For the simulation of case no. 3, whose results are reported in Figure 5, the same condition of case no. 2 are considered, but the perfect knowledge of the arm parameters is assumed, and the external force has been considered nonmeasurable: from the comparison with Figure 4, it is possible to note that the proposed virtual force sensor is able to correctly estimate the external force; however, the response of the system in case no. 3 shows larger oscillations during transients with respect to case no. 2, but the response in static conditions is identical to those obtained considering direct measure of the external force, that is, in case no. 2.

Figure 5: Response of the system during test no. 3.

5. Conclusions

In this paper, the problem of output-based regulation of robotic manipulators with variable stiffness actuation has been discussed. A control scheme has been proposed whose main features are the simultaneous and decoupled workspace stiffness-position control with limited error, disturbance compensation and estimation-control of the environment interaction force. Moreover, the case of manipulator redundancy has been explicitly addressed and some nice features and implications related to the control of the both the joint mechanical stiffness and workspace stiffness have been highlighted. The control objective is achieved by means of a multilevel controller that decouples the dynamics of both the positioning and the stiffness actuators from the arm dynamics. This result is obtained by applying a singular perturbation analysis to both the actuators and the arm dynamics, configuring in this way a multilevel control structure that presents fast controllers for the actuators and a slow controller for the arm. Moreover, the proposed controller presents several advantages with respect to previous approaches for this class of devices, as the rigid-robot-like trajectory specification, the reduced complexity and the availability of the system residuals for collision detection/reaction purposes to cite some of them. The extension of the proposed analysis to an hybrid case, in which both rigid and elastic joints are present, with constant or variable stiffness, can be easily achieved. Also, the case of joints with different stiffness variation models can be easily considered. The proposed control approach is validated by means of simulation of a three-link planar manipulator, and the experimental validation of the proposed approach will be considered in the future.


Property 5. The general second-order dynamic system (10) can be used in place of each of (1)–(3) by posing: (i)In case of (1) 𝑣=𝑞,𝑤=[𝜃𝑘],𝛼=𝜂𝑞,𝜏=𝜏𝑒,𝐴(𝑣,𝑤)=𝑀(𝑞),𝑓(̇𝑣,𝑣,𝑤)=𝑁(𝑞,̇𝑞)+𝐾[𝑞𝜃].(A.1)(ii)In case of (2), 𝑣=𝜃,𝑤=[𝑞𝑘],𝛼=𝜂𝜃,𝜏=𝜏𝜃,𝐴(𝑣,𝑤)=𝐵,𝑓(̇𝑣,𝑣,𝑤)=𝐾[𝜃𝑞].(A.2)(iii)In case of (3), 𝑣=𝑘,𝑤=[𝑞𝜃],𝛼=𝜂𝑘,𝜏=𝜏𝑘,𝐴(𝑣,𝑤)=𝛾(𝑞,𝜃,𝑘),𝑓(̇𝑣,𝑣,𝑤)=𝛽(𝑞,𝜃,𝑘).(A.3)

Property 6. The general second-order dynamic system (10) can be used in place of (41)–(43) by posing 𝑣=𝑞,𝑤=[𝜃𝑘],𝛼=𝜂𝑞+𝜏𝑒,𝜏=𝜏act𝑑,𝐴(𝑣,𝑤)=𝑀(𝑞),𝑓(̇𝑣,𝑣,𝑤)=𝑁(𝑞,̇𝑞).(A.4)

Theorem 3. Let 𝐾𝑑, 𝐽, and 𝑤𝐾𝑑 be 𝑛×𝑛 diagonal positive definite, 𝑚×𝑛 and 𝑚×𝑚 symmetric positive definite matrices, respectively. The following definitions are given: 𝑤𝐾𝑗1=𝐽𝐾1𝑑𝐽𝑇,(A.5)𝑤𝐾𝑐=𝑤𝐾𝑑𝑤𝐾𝑗,(A.6)𝐾𝑎=𝐾𝑑+𝐽𝑇𝑤𝐾𝑐𝐽,(A.7) where the properties of symmetry and positive definition of the matrices 𝑤𝐾𝑗, 𝑤𝐾𝑐, and 𝐾𝑎 should hold. Note that (A.5)–(A.7) correspond to (49)–(51). Then, 𝐽𝐾1𝑎𝐽𝑇=𝑤𝐾𝑑1,(A.8)𝐾1𝑎𝐽𝑇𝐽𝐾1𝑎𝐽𝑇1=𝐾1𝑑𝐽𝑇𝐽𝐾1𝑑𝐽𝑇1.(A.9)

Proof. Let consider (A.9) first. By applying the matrix inversion lemma [38] to (A.7) 𝐾1𝑎=𝐾1𝑑𝐾1𝑑𝐽𝑇𝑤𝐾𝑐1+𝑤𝐾𝑗11𝐽𝐾1𝑑.(A.10) It follows that 𝐽𝐾1𝑎𝐽𝑇=𝑤𝐾𝑗1𝑤𝐾𝑗1𝑤𝐾𝑐1+𝑤𝐾𝑗11𝑤𝐾𝑗1.(A.11) By applying the matrix inversion lemma also to (A.6), it results 𝑤𝐾𝑑1=𝑤𝐾𝑗1𝑤𝐾𝑗1𝑤𝐾𝑐1+𝑤𝐾𝑗11𝑤𝐾𝑗1.(A.12) By noticing that (A.11) and (A.12) are equivalent, it follows that the relation (A.8) is satisfied.
Now (A.9) is considered. From (A.10) and (A.8), 𝐾1𝑎𝐽𝑇𝐽𝐾1𝑎𝐽𝑇1=𝐾1𝑑𝐾1𝑑𝐽𝑇𝑤𝐾𝑐1+𝑤𝐾𝑗11𝐽𝐾1𝑑𝐽𝑇𝑤𝐾𝑑.(A.13) By introducing (A.6) 𝐾1𝑎𝐽𝑇𝐽𝐾1𝑎𝐽𝑇1=𝐾1𝑑𝐾1𝑑𝐽𝑇𝑤𝐾𝑐1+𝑤𝐾𝑗11𝐽𝐾1𝑑𝐽𝑇𝑤𝐾𝑐+𝑤𝐾𝑗=𝐾1𝑑𝐽𝑇𝑤𝐾𝑐1+𝑤𝐾𝑗11𝑤𝐾𝑗1𝑤𝐾𝑐+𝑤𝐾𝑗+𝑤𝐾𝑐+𝑤𝐾𝑗=𝐾1𝑑𝐽𝑇𝑤𝐾𝑐1+𝑤𝐾𝑗11𝑤𝐾𝑗1+𝑤𝐾𝑐1𝑤𝐾𝑐+𝑤𝐾𝑐+𝑤𝐾𝑗=𝐾1𝑑𝐽𝑇𝐽𝐾1𝑑𝐽𝑇1,(A.14) where (A.5) has been finally used.


This research has been partially supported by the EC Seventh Framework Programme (FP7) under Grant agreement no. 216239 as part of the IP DEXMART (DEXterous and autonomous dual-arm/hand robotic manipulation with sMART sensory-motor skills: a bridge from natural to artificial cognition) and by the Italian MIUR with the PRIN2007CCRNFA-004/PRIN2008-ZRLPWS projects.


  1. C. E. English and D. Russell, “Mechanics and stiffness limitations of a variable stiffness actuator for use in prosthetic limbs,” Mechanism and Machine Theory, vol. 34, no. 1, pp. 7–25, 1999. View at Scopus
  2. A. Bicchi and G. Tonietti, “Fast and soft arm tactics: dealing with the safety-performance trade-off in robot arms design and control,” IEEE Robotics and Automation Magazine, vol. 11, no. 2, pp. 22–33, 2004.
  3. G. Berselli and G. Vassura, “Differentiated layer design to modify the compliance of soft pads for robotic limbs,” in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '09), pp. 1285–1290, Kobe, Japan, May 2009. View at Publisher · View at Google Scholar · View at Scopus
  4. A. Bicchi, G. Tonietti, M. Bavaro, and M. Piccigallo, “Variable stiffness actuators for fast and safe motion control,” in Proceedings of the Robotics Research: The Eleventh International Symposium (ISRR '05), vol. 15, pp. 527–536, Springer, Italy, 2005.
  5. S. A. Migliore, E. A. Brown, and S. P. DeWeerth, “Novel nonlinear elastic actuators for passively controlling robotic joint compliance,” Journal of Mechanical Design, vol. 129, no. 4, pp. 406–412, 2007. View at Publisher · View at Google Scholar · View at Scopus
  6. S. Wolf and G. Hirzinger, “A new variable stiffness design: matching requirements of the next robot generation,” in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '08), pp. 1741–1746, Pasadena, Calif, USA, May 2008. View at Publisher · View at Google Scholar · View at Scopus
  7. R. Schiavi, G. Grioli, S. Sen, and A. Bicchi, “VSA-II: a novel prototype of variable stiffness actuator for safe and performing robots interacting with humans,” in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '08), pp. 2171–2176, Pasadena, Calif, USA, 2008. View at Publisher · View at Google Scholar
  8. D. Hyun, H. S. Yang, J. Park, and Y. Shim, “Variable stiffness mechanism for human-friendly robots,” Mechanism and Machine Theory, vol. 45, no. 6, pp. 880–897, 2010. View at Publisher · View at Google Scholar · View at Scopus
  9. G. Palli, C. Melchiorri, T. Wimböck, M. Grebenstein, and G. Hirzinger, “Feedback linearization and simultaneous stiffness-position control of robots with antagonistic actuated joints,” in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '07), pp. 4367–4372, Rome, Italy, April 2007. View at Publisher · View at Google Scholar · View at Scopus
  10. G. Palli, C. Melchiorri, and A. de Luca, “On the feedback linearization of robots with variable joint stiffness,” in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '08), pp. 1753–1759, Pasadena, Calif, USA, 2008. View at Publisher · View at Google Scholar · View at Scopus
  11. G. Palli and C. Melchiorri, “Robust control of robots with variable joint stiffness,” in Proceedings of the International Conference on Advanced Robotics (ICAR '09), pp. 1–6, Munich, Germany, June 2009. View at Scopus
  12. A. Bicchi, G. Tonietti, M. Bavaro, and M. Piccigallo, “Variable stiffness actuators for fast and safe motion control,” Springer Tracts in Advanced Robotics, vol. 15, pp. 527–536, 2005. View at Publisher · View at Google Scholar
  13. M. W. Spong, “Modeling and control of elastic joint robots,” Journal of Dynamic Systems, Measurement and Control, vol. 109, no. 4, pp. 310–319, 1987. View at Scopus
  14. J. K. Salisbury, “Active stiffness control of a manipulator in cartesian coordinates,” in Proceedings of the 19th IEEE Conference on Decision and Control including the Symposium on Adaptive Processes, vol. 19, pp. 95–100, December 1980. View at Publisher · View at Google Scholar · View at Scopus
  15. N. Hogan, “Impedance control: an approach to manipulation, parts I–III,” ASME—Journal of Dynamic Systems, Measurement, and Control, vol. 107, pp. 1–24, 1985. View at Scopus
  16. A. Albu-Schäffer, C. Ott, and G. Hirzinger, “A unified passivity-based control framework for position, torque and impedance control of flexible joint robots,” International Journal of Robotics Research, vol. 26, no. 1, pp. 23–39, 2007. View at Publisher · View at Google Scholar · View at Scopus
  17. S. Haddadin, A. Albu-Schffer, and G. Hirzinger, “Safety analysis for a human-friendly manipulator,” International Journal of Social Robotics, vol. 2, no. 3, pp. 235–252, 2010. View at Publisher · View at Google Scholar
  18. D. W. Robinson, J. E. Pratt, D. J. Paluska, and G. A. Pratt, “Series elastic actuator development for a biomimetic walking robot,” in Proceedings of the IEEE/ASME International Conference on Advanced Intelligent Mechatronics (AIM '99), pp. 561–568, Atlanta, Ga, USA, September 1999. View at Scopus
  19. B. T. Knox and J. P. Schmiedeler, “A unidirectional series-elastic actuator design using a spiral torsion spring,” Journal of Mechanical Design, vol. 131, no. 12, pp. 125001–125010, 2009. View at Publisher · View at Google Scholar · View at Scopus
  20. K. Yokoi, H. Maekawa, and K. Tanie, “A method of compliance control for a redundant manipulator,” in Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems, vol. 3, pp. 1927–1934, July 1992.
  21. T. Wimbock, C. Ott, A. Albu-Schaffer, A. Kugi, and G. Hirzinger, “Impedance control for variable stiffness mechanisms with nonlinear joint coupling,” in Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems, pp. 3796–3803, Nice, France, Septemper 2008. View at Publisher · View at Google Scholar
  22. A. Teel and L. Praly, “Global stabilizability and observability imply semi-global stabilizability by output feedback,” Systems and Control Letters, vol. 22, no. 5, pp. 313–325, 1994. View at Scopus
  23. B. Yang and W. Lin, “On semi-global stabilizability of MIMO nonlinear systems by output feedback,” Automatica, vol. 42, no. 6, pp. 1049–1054, 2006. View at Publisher · View at Google Scholar · View at Scopus
  24. A. Isidori, Nonlinear Control Systems II, Springer, Berlin, Germany, 1999.
  25. A. Tornambè, “Output feedback stabilization of a class of non-minimum phase nonlinear systems,” Systems and Control Letters, vol. 19, no. 3, pp. 193–204, 1992. View at Scopus
  26. H. K. Khalil, Nonlinear Systems, Prentice Hall, New York, NY, USA, 3rd edition, 2002.
  27. L. Sciavicco and B. Siciliano, Modeling and Control of Robot Manipulators, McGraw–Hill, New York, NY, USA, 1996.
  28. R. V. Patel and F. Shadpey, Control of Redundant Robot Manipulators: Theory and Experiments, Springer, New York, NY, USA, 2009.
  29. H. Y. Choi, W. K. Chung, and Y. Youm, “Stiffness analysis and control of redundant manipulators,” in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '94), vol. 1, pp. 689–695, San Diego, Calif, USA, May 1994. View at Publisher · View at Google Scholar
  30. M. M. Svinin, S. Hosoe, M. Uchiyama, and Z. Luo, “On the stiffness and stiffness control of redundant manipulators,” in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '02), pp. 2393–2399, Washington, DC, USA, May 2002. View at Publisher · View at Google Scholar
  31. W. Yu and X. Li, “PD control of robot with velocity estimation and uncertainties compensation,” International Journal of Robotics and Automation, vol. 21, no. 1, pp. 1–9, 2006. View at Scopus
  32. E. Shin and K. Lee, “Robust output feedback control of robot manipulators using high-gain observer,” in Proceedings of the IEEE Conference on Control Applications, vol. 1, pp. 881–886, Kohala Coast, Hawaii, USA, August 1999. View at Publisher · View at Google Scholar
  33. P. Tomei, “Robust adaptive friction compensation for tracking control of robots,” in Proceedings of the IEEE Conference on Control Applications, vol. 1, pp. 875–880, Kohala Coast, Hawaii, USA, August 1999.
  34. R. Kelly, V. Santibánez, and A. Loría, Control of Robot Manipulators in Joint Space, Springer, New York, NY, USA, 1st edition, 2005.
  35. B. Siciliano and O. Khatib, Eds., Springer Handbook of Robotics, Springer, New York, NY, USA, 2008.
  36. A. de Luca and R. Mattone, “Sensorless robot collision detection and hybrid force/motion control,” in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '05), pp. 999–1004, Barcelona, Spain, April 2005. View at Publisher · View at Google Scholar
  37. A. de Luca, F. Flacco, R. Schiavi, and A. Bicchi, “Nonlinear decoupled motion-stiffness control and collision detection/reaction for the VSAII variable stiffness device,” in Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems, pp. 5487–5494, Saint-Louis, Mo, USA, October 2009. View at Publisher · View at Google Scholar
  38. T. Anderson, An Introduction to Multivariate Statistical Analysis, Wiley-Interscience, New York, NY, USA, 2003.