# A Configurable Architecture for Two Degree-of-Freedom Variable Stiffness Actuators to Match the Compliant Behavior of Human Joints

^{1}Centro “E. Piaggio” and Dipartimento di Ingegneria dell’Informazione, University of Pisa, Pisa, Italy^{2}Soft Robotics for Human Cooperation and Rehabilitation, Fondazione Istituto Italiano di Tecnologia, Genoa, Italy

Living beings modulate the impedance of their joints to interact proficiently, robustly, and safely with the environment. These observations inspired the design of soft articulated robots with the development of Variable Impedance and Variable Stiffness Actuators. However, designing them remains a challenging task due to their mechanical complexity, encumbrance, and weight, but also due to the different specifications that the wide range of applications requires. For instance, as prostheses or parts of humanoid systems, there is currently a need for multi-degree-of-freedom joints that have abilities similar to those of human articulations. Toward this goal, we propose a new compact and configurable design for a two-degree-of-freedom variable stiffness joint that can match the passive behavior of a human wrist and ankle. Using only three motors, this joint can control its equilibrium orientation around two perpendicular axes and its overall stiffness as a one-dimensional parameter, like the co-contraction of human muscles. The kinematic architecture builds upon a state-of-the-art rigid parallel mechanism with the addition of nonlinear elastic elements to allow the control of the stiffness. The mechanical parameters of the proposed system can be optimized to match desired passive compliant behaviors and to fit various applications (e.g., prosthetic wrists or ankles, artificial wrists, etc.). After describing the joint structure, we detail the kinetostatic analysis to derive the compliant behavior as a function of the design parameters and to prove the variable stiffness ability of the system. Besides, we provide sets of design parameters to match the passive compliance of either a human wrist or ankle. Moreover, to show the versatility of the proposed joint architecture and as guidelines for the future designer, we describe the influence of the main design parameters on the system stiffness characteristic and show the potential of the design for more complex applications.

## 1 Introduction

In the last years, articulated soft robots inspired by the musculoskeletal system of vertebrate animals received increased attention from researchers, since they represent promising solutions to enhance the interactions of robots with unknown and dynamic environments, i.e., the real world (Albu-Schaeffer et al., 2008). In line with this, Variable Impedance Actuators and the subgroup of Variable Stiffness Actuators (VSAs) were widely investigated recently (Vanderborght et al., 2013).

Variable Stiffness Actuators usually rely on two main approaches to modulate the stiffness. In the first one, the stiffness is controlled using software together with a mechanism with fixed impedance properties. In the other one, the stiffness is modulated through a mechanical reconfiguration of the system (Vanderborght et al., 2013). Some literature refers to the first approach as *active VSA* and the second one as *passive VSA* (Vanderborght et al., 2013). However, to avoid confusion possibly induced by the notion of passivity, we will refer to the first approach as *software-controlled VSA* and to the second one as *physically compliant VSA*. Software-controlled VSAs allow the design of lightweight devices that can in theory simulate any desired stiffness behavior. However, this apparent stiffness emulated by the system relies on an accurate sensing strategy and control computations, and it has been shown that even if the impacts are detected timely, the motors could not be able to react fast enough by solely an impedance control and that the system should therefore be considered as stiff during the impacts (Haddadin et al., 2007). To address these limitations, physically compliant VSAs have been developed with an inherent compliance. They can present several advantages such as shock absorptions, better performances in cyclic or explosive tasks (Albu-Schaeffer et al., 2008; Wolf et al., 2016), and a possible embodiment of specific behavior to improve the control strategies (Visser et al., 2011).

As VSAs are inspired by musculoskeletal systems, their applications as prostheses or part of humanoid devices would seem straightforward. Even though there are some advantages in using VSAs compared to stiff actuators, finding concrete use-cases of VSAs is an on-going research topic and requires attention (Sensinger and ff. Weir, 2008; Blank et al., 2014; Stillfried et al., 2018).

To test hypotheses in this objective, we need multi-degree-of-freedom (DoF) joints that show abilities similar to those of a human joint, such as the same functional range of motion, a variable stiffness mechanism, the overall shape and mass, and so on. Integrating an inherent compliance inside these systems enables us to test both the variable stiffness ability but also the capabilities due to the compliance, such as the possibility of exploiting the natural dynamic of the system (e.g., the resonance frequency) (Stillfried et al., 2018). Besides, artificial wrists can benefit from a compliant architecture to enhance their manipulation abilities in tight or cluttered spaces (Negrello et al., 2019). Moreover, the availability of various levels of stiffness can enhance the performances in activities of daily living in prosthetic applications (Kanitz et al., 2018). However, the mechanical complexity of such mechanisms increases and constitutes a challenge.

Up to now, most of the proposed physically compliant VSAs have only 1 DoF in position (Vanderborght et al., 2009; Vanderborght et al., 2013; Petit et al., 2015; Lemerle et al., 2019). The classical approach to get multi DoF physically compliant VSAs is to put in series several 1-DoF VSAs, such as the solutions proposed in the DLR Hand Arm System (Grebenstein et al., 2011; Romano et al., 2014; Savin et al., 2019). This approach uses the design of serial manipulators (SMs) that are easier to model and control compared to the parallel manipulators (PMs). The latter, on the other hand, allow more compact designs and have better performances in terms of output torques than SMs for the same motor sizes. Nevertheless, the performances of PMs are very sensitive to their geometric parameters (Siciliano and Khatib, 2007). Alternative works explored the independent setup approach, consisting mainly of a serial arrangement of position motor units with several DoF in position (using either SMs or PMs) and an additional motor unit dedicated solely to the control of the stiffness (Weckx et al., 2014; Barrett et al., 2017; Tan et al., 2017; Rodriguez-Cianca et al., 2019). Finally, other works draw inspiration from the antagonistic arrangements of the musculoskeletal system (Koganezawa, 2017; Stoeffler et al., 2018; Malosio et al., 2019). In the last category, several antagonistic motors are used jointly to control the stiffness and the position of the joint. Theses architectures are solely based on PMs. For more details on the classification and principle of work of VSAs, the reader is invited to refer to (Vanderborght et al., 2013; Wolf et al., 2016).

One limitation of physically compliant VSAs is that one mechanical implementation corresponds to one embedded passive behavior. Although there are devices that are human-like in terms of shapes and dimensions, and are able to perform activities of daily living (ADL) (Grebenstein et al., 2011; Weckx et al., 2014; Koganezawa, 2017; Stoeffler et al., 2018; Rodriguez-Cianca et al., 2019), there is, for now, no multi-DoF VSA embedding a compliant behavior that matches human joints. It is worth mentioning that there are alternatives to design compact 2 DoF joints with inherent compliance (not necessarily with variable stiffness in the following examples). A first concept is to design an underactuated mechanism (Casini et al., 2017). Its main drawback could be its limited dexterity to do all the desired tasks. Another approach is to have a switchable stiffness to get lightweight and compact devices (Montagnani et al., 2013). This type of solution could be suitable when lightweight systems are needed. However, they do not allow for a continuous and smooth control of the stiffness. Therefore, they lose some advantages of such abilities, e.g., during explosive or cyclic tasks, when a VSA performs better than a fixed compliant mechanism (Garabini et al., 2011). An additional idea is to use tendon driven mechanisms with compliant pneumatic actuation (Toedtheide and Haddadin, 2018). Then, this system could be lightweight and somehow compact if the remote actuation units are excluded. Yet, this is not fully satisfying in terms of compactness if we consider them, especially when they are pneumatic actuation units. Therefore, all of these approaches do not fully satisfy the criteria of compactness and variable stiffness ability as well as doing various types of tasks. We address this limitation in the following work.

In general, the qualitative representation of the mechanical impedance of human joints is made using stiffness (or compliance) ellipsoids (Mussa-Ivaldi et al., 1985). To modulate the stiffness ellipsoids of the end effector during tasks, it appears that the limb geometry and posture have a large impact on its shape and orientation, whereas the voluntary co-contraction of the participating muscles affects mainly the size of the ellipsoids (Milner, 2002; Perreault et al., 2002). To quantify the mechanical impedance properties of multi-DoF human joints, the traditional approach consists in applying position or force perturbations to the joint of interest and measuring the displacement force or position responses (Mussa-Ivaldi et al., 1985). Yet, this process is complex and long. Investigations are made to simplify and speed up the procedure (Masia et al., 2012; Masia and Squeri, 2015). The measurements usually concern the endpoint impedance (Mussa-Ivaldi et al., 1985; Masia et al., 2012). Yet, there is mechanical impedance data of specific human joints, such as the wrist (Formica et al., 2012; Pando et al., 2014) or the ankle (Kearney et al., 1997; Lee et al., 2014). There is, for now, no data on the mechanical impedance of more complex joints such as the shoulder or hip. Yet, efforts are made toward this direction with the design of low inertia shoulder exoskeleton to measure neuromuscular properties (Hunt and Lee, 2019).

The observations on how human beings modulate the stiffness ellipsoids were successfully exploited to implement a teleimpedance controller to generate human-like motions and desired task space impedance (Ajoudani et al., 2018). Similarly, these observations could be used to reduce the mechanical complexity of human-like VSAs. Indeed, if the human-like compliant behavior is directly implemented in the mechanism, only one additional actuator is ideally required to modulate the overall size of the stiffness ellipsoid, in analogy to the voluntary co-contraction of human muscles. As a matter of fact, from the study of tendon mechanisms (e.g., the Salisbury Hand), we know that to control N DoF in position only N+1 actuators are required. There is then an internal tension that can be used to modulate the stiffness. This was shown for a grasping hand in Mason and Salisbury (1985).

In this work, we propose the concept and modeling of a configurable architecture of a physically compliant 2 DoF VSA that can match the passive behavior of either a human wrist or a human ankle. This joint can control its orientations around 2 perpendicular axes and its overall stiffness as a one-dimensional parameter similar to the co-contraction of human muscles. It is based on a state-of-the-art parallel manipulator that we modified to implement an inherent compliance with a compact architecture. Hence, the variable stiffness (VS) ability of the system relies on the antagonistic arrangement of nonlinear elastic actuators (Vanderborght et al., 2013; Wolf et al., 2016). The compact architecture of the system is based on the PM architecture with the minimum number of required motors to get a VSA for 2 DoF in position (i.e., 3). Thanks to its kinematics, the proposed system could be used as a wrist or ankle joint. It can also be part of more complex joints such as the hip or shoulder joint. Indeed, these joints are represented with 3 DoF as ball-and-socket joints and one of their three rotations is around a longitudinal axis of the body segment attached to the joint (Wilhelms and Gelder, 2001). Therefore, we can design them with a hybrid architecture of a PM with 2 DoF in series with a 1 DoF rotator, like for instance the humeral rotator proposed in Mazzotti et al., (2015), which can be integrated with a 2 DoF shoulder joint presented in Troncossi et al., (2009). A similar idea could be proposed for a hip joint with a femoral rotator. Therefore, our proposed system could also be used as part of these more complex joints. Figure 1 shows the general principles of the system and some of its potential applications.

**FIGURE 1**. The proposed 2-DOF variable stiffness joint allows shaping the stiffness (or compliance) ellipsoid at different postures in the design phase, while the volume can be varied in real-time to make the overall behavior stiffer or softer. The kinematics of the 2-DOF VS joint can be used for the design of wrists and ankles of human-like bionic systems (in orange). In principle, this joint can also be used as part of more complex designs for a shoulder or hip human-like joint (in yellow).

Besides the general description of the proposed system in Section 2.1, the main contributions of this paper are: (i) to prove the human-like VS ability of the system with a compact architecture (Section 4.1), (ii) to study the effect of the main design parameters of the system on its passive compliance to provide general guidelines (Section 4.2), and (iii) to provide sets of design parameters to match the passive compliance of a human wrist or ankle based on an optimization process (Section 4.3). Other contributions of this paper include improvements on the kinematic analysis of the rigid structure (used as a basis of our mechanical architecture) done in Sofka et al., (2006a) and the kinetostatic analysis of our proposed compliant structure (Section 3).

The paper is organized as in the following: Section 2 details the design concept of the proposed system and the methodology of its study. It includes, in particular, a description of its general architecture and its main design parameters. Section 3 provides the model and the kinetostatic analysis of the system. Section 4 gives the main results of the current study that are then discussed in Section 5. Finally, Section 6 summarizes the main contributions and concludes the article.

## 2 Design Concept and Methodology

### 2.1 General Architecture of the System

Based on the state of the art of compact 2 DoF in orientation joints based on PMs (refer to (Bajaj et al., 2019) for a quite extended list), we selected the kinematic architecture of the Omni-Wrist III, designed by Rosheim and Sauter (2002) and studied by Sofka et al., (2006a) and Sofka et al., (2006b). This design has several advantages such as 180° hemispherical singularity free movement (Rosheim and Sauter, 2002). In theory, only two legs are needed to actuate the system but the main version has 4 legs. Its general architecture is shown in Figure 2A. In our design, we propose to use a 3–legs configuration with nonlinear elastic elements added in the serial arrangement of each leg. Figure 2B presents its general architecture. The system is a PM composed of a fixed frame (in dark gray on the scheme) and a moving platform (in light gray), i.e., the end effector, linked together with three independent and identical kinematic serial chains, called respectively leg A (in blue), leg B (in green) and leg C (in red). Each leg is composed of the serial arrangement of one motor *α*, which is the angular deviation of the link 2 (in red), and by a characteristic length *d* (Sofka et al., 2006a). Figure 2C shows the main structure of one KC with its parameters.

**FIGURE 2**. **(A)** is the classic architecture of the Omni-Wrist III (Rosheim and Sauter, 2002) is composed of 4 legs linking a fixed frame and a moving platform. Each actuated leg is composed of the serial arrangement of one motor *KC*_{*}) **(B)** Represents the architecture of the proposed system. We use a configuration with 3 legs only. In each leg, we add elastic elements **(C)** shows the details of one kinematic chain (each KCs are identical). It is composed of 3 linkages (1, 2 and 3) and 4 noncoplanar revolute joints (I, II, III, IV). The joint IV is linked to the end effector and the joint I is linked either to the output of the motor in the Omni-Wrist III or to the output of the elastic element in our proposed architecture.

### 2.2 Variable Stiffness Ability

As it is common in VSA with antagonistic architecture, to get the variable stiffness ability, the stiffness characteristic of the elastic element has to be nonlinear. There are many ways to design nonlinear springs, to cite a few, the physical implementation can rely on tendons and linear springs (Catalano et al., 2011), on a combination of belts and linear springs (Lemerle et al., 2019), on a combination of torsional spring with guide shafts with varying radius (Koganezawa et al., 2005; Koganezawa, 2017), or on hydraulic or pneumatic solutions (Stoeffler et al., 2018).

Nevertheless, nonlinear stiffness alone is not sufficient to grant the variable stiffness control of the system, which is equivalent to prove that we can control independently the equilibrium orientations of the system and its overall stiffness along one dimension. To prove this analytically, the kinetostatic analysis of the system is first presented in Section 3. Then, the final arguments to show the VS control of the proposed system are given in Section 4.1.

In addition, we would like to show the effects of the modulation of the stiffness on the overall compliance of the system. To represent them graphically, we need to derive a model of the system (done in Section 3) and to model the elastic elements. In this study, we make the deliberate choice of adopting a hyperbolic sine spring, whose mechanical characteristics have the form

where

For the sake of clarity, the detailed graphical representation methodology is described after the complete analysis of the system in Section 3.8. The effects of the modulation of the stiffness are shown in Section 4.1.

### 2.3 General Design Guidelines

As previously described, there are various ways to implement the nonlinear elastic elements. In a first approach, we leave this choice to the designer to focus on the design of the KCs and their arrangement. Therefore, as general design guidelines, we study the effects of the purely geometric parameters (simply referred to as geometric parameters) on the compliance of the proposed system. To do this study, the compliance of the system should be modeled as a function of these geometric parameters. This is done in Section 3. The first step in the modeling is to define the geometric parameters. They correspond to the parameters of the KCs *α* and *d* (Sofka et al., 2006a), and the angular arrangement of the first joint of each leg around the basis, namely

In the case of a fixed distribution of the legs around the basis, *α* and *d* are not enough to modulate the compliant behavior of the system (refer to Section 4.2 for more details). That is why we need to extend the modeling done in Sofka et al. (2006a) by adding the angular positions of the legs.

Moreover, to study separately the effects on the compliance due to the kinematic architecture from the effects due to the internal variation of the compliance of the elastic elements, the compliance characteristic of each leg is modeled as a constant. Concretely, to study the impact of the angular position of the legs, two types of configuration are simulated. Firstly the leg A is considered as fixed with *π*}. To study the impact of the angular deviation of the middle link of each leg *α*, this angle is varying in the set {*d* of each leg, this length is varying in the set {25 mm, 50 mm, 75 mm, and 100 mm}. All the simulations and computations are performed using MATLAB. The results are shown in Section 4.2.

### 2.4 Design Parameters to Match Compliance of Human Wrist and Ankle

One of the main goals of this study is to show that we can implement passive compliant behaviors of human joints. In this work, we focus on two different human joints: the wrist and the ankle. To find the appropriate sets of design parameters to match their passive compliant behavior, we derived an optimization process. In the first step in Section 2.4.1, we need to define the behaviors we would like to embed based on the mechanical impedance data of human joints. Then, in Section 2.4.2 we define the design parameters of the proposed architecture. Finally, to derive the cost function of the minimization algorithm, we need to characterize the compliance of the system as a function of the design parameters. This is done in Section 3 and the cost function is defined in Eq. 58. Results of the optimization process are shown in Section 4.3.

##### 2.4.1 Parameters of Human Wrist and Human Ankle

Based on the state-of-the-art, we selected the biomechanical data of a human wrist presented in Pando et al. (2014) and of a human ankle presented in Lee et al. (2014).

In Pando et al. (2014), Pando et al. studied, among others, the passive wrist stiffness in the two dimensions that are the flexion/extension and the radial/ulnar deviation. They extracted from their data the stiffness ellipses around the neutral position of the wrist in case of no voluntary co-contraction of the subjects to evaluate thus the passive stiffness of the articulation. They found that the human wrist stiffness is anisotropic and the direction of the greater magnitude of the average stiffness ellipse is around 12.1° from the pure radial deviation in the counter-clockwise direction toward a pure flexion. The ratio between the magnitude of the two axes of the stiffness ellipse is about 2.69 and using the area of the ellipse [5.61 (Nm/rad)^{2}], we can extract the lengths of the two semi-axes of the ellipse, which are then 2.19 and 0.81 Nm/rad.

In Lee et al. (2014), Lee et al. studied the mechanical impedance of a human ankle with relaxed muscles, i.e., in passive conditions, both in seated and standing positions. They found that for low-frequency motion the stiffness dominates the mechanical impedance of the ankle, which is our application case as we consider the static stiffness in our design. Its characteristics in standing positions are: a tilt angle of 10°, a magnitude in the major principal axis direction of 23.86 Nm/rad, and a magnitude in the orthogonal direction of 8.14 Nm/rad (the ratio is 2.99).

Considering these data, we can derive the average Cartesian compliance ellipse and the associated compliance matrix for both a human wrist and ankle, using similar parametrization as the one proposed in Section 3.7.

##### 2.4.2 Design Parameters Definition

In a second step, we should define all the design parameters of the proposed architecture. They can be divided into geometric parameters (previously described), elastic parameters, and internal modulation parameters. As a reminder, there are 5 geometric parameters: *α*, *d*,

Each nonlinear spring can be characterized by nonlinear functions. Their explicit formulations depend on mechanical implementation. They can have various forms such as polynomial functions of degree higher than 2 (e.g., quadratic springs), trigonometric functions, exponential or hyperbolic functions to name but a few (English and Russell, 1999; Koganezawa et al., 2005; Vanderborght et al., 2009; Catalano et al., 2011; Petit et al., 2015; Lemerle et al., 2019). Therefore the number of their associated parameters depends on the complexity of the model. These parameters are referred to as elastic parameters. Noting

The internal modulation parameter, called *λ*, stands for the one-dimensional quantity we can use to control the stiffness in the proposed system. Unlike the previous parameters, *λ* is not fixed to one mechanical implementation and can therefore be different if we try to match several compliant behaviors. Therefore, there is one additional parameter for each desired behavior.

Therefore, the number of design parameters, noted

where

For this study, we model the nonlinear elastic elements as hyperbolic sine springs for similar reasons as previously described in Section 2.2. Their mechanical characteristics have the form defined in Eq. 1. Therefore, there are 2 elastic parameters for each leg: (i) the spring constant

## 3 Model and Analysis

In this section, we describe the kinetostatic analysis of the proposed system with two goals in mind: (i) to derive the compliance matrix as a function of the design parameters and (ii) to show the variable stiffness ability of the system.

Firstly the kinematic analysis done by Sofka et al. (2006a) needs to be extended to include the angular arrangement of the legs. Then, we describe the inverse and differential kinematics of the system in a closed form, showing that it is a 2 DoF joint. We provide an explicit and complete formulation of the general mapping of the generalized pose of the end effector as a function of a minimal parametrization. This result is then used to analyze the static equilibrium of the system and describe its static stiffness and Cartesian compliance, required for the optimization algorithm. A list of the main symbols used in the analysis can be found in Table 2.

### 3.1 Parametrization of the System

Define the reference coordinate frame attached to the fixed base of the mechanism

where

A possible choice of parameters that easily parametrizes the generalized pose of the end effector is

where the first three components stand for the Euler angles of the end effector, defining, therefore, its angular orientation about respectively the axis **Z**_{0}, namely (roll, yaw, pitch) and the last three components are the position of its center

Denoting with *d* along the axis

where

i.e.,

where, *k* in

Considering now each leg (*A*, *B*, and *C*) of the mechanism, it is possible to parametrize the configuration of *d* and *α* are fixed parameters that define the shared legs geometry, while *η* is different for each leg.

**FIGURE 4**. Kinematic diagram of a generic leg. It is composed of 3 linkages (1, 2, and 3), and 4 noncoplanar joints. It is parameterized by its angular deviation *α*, its characteristic length *d* and its angular orientation *η*. The coordinates frames are also represented. Its associated DH table is given in Table 3.

Calling *η* the joint angles and the base angle of the generic leg, the global transformation matrix going from the system base to the end effector through the leg

where all the intermediate matrices are defined based on the parameters in Table 3. After calculations, we obtain:

where all the non-trivial terms are defined in Supplementary Appendix A.

### 3.2 Closing the Kinematic Chain

When the joint is mounted, i.e., the legs are assembled, the assembly constrains the final frame associated to each leg (and the generic frame

As it is well-known, the solution of Eq. 8, is fundamental to yield both the direct and inverse kinematics of the parallel system. An important aid to facilitate this comes from Sofka et al., (2006a), where it was found that a feasible mounting of each leg holds for

Therefore, after simplifications, Eq. 7 can be written as

where is a family of functions defined as

and

The families

As expected from Sofka et al. (2006a), and as we will confirm with our derivation, the assembled system has 2 DoF. Thus, it is possible to parametrize the position and orientation of the reachable end effector configurations with only two variables. For convenience we choose

To fully characterize the set of feasible end effector configurations it is necessary to specify the other 4 variables in Eq. 3.

To derive the angular orientation

and then, by using Eq. 5 and rearranging the terms, it is possible to extract

a result in agreement with (Sofka et al., 2006a).

Then, the analysis can be extended by deriving an explicit formulation of the position of the center of the end effector

from which it can be shown [based on (10)] that

which gives the distance square of the point *d* and *α* which are two parameters inherent of the mechanical design of the system and could also be found based on a geometric closed form of the system.

Defining such distance as

and using the previous equations, it can be shown that

According to the convention defined previously, we have

therefore, using Eq. 19a, we obtain

Then, the system defined by Eqs. 19b and 19c has 4 pairs

To conclude, the proposed system has 2 DoF in orientation and the explicit formulation of the general mapping *T* is as

This mapping is always defined for all

### 3.3 Inverse Kinematics Model

This subsection aims to find the position of all the joints of each leg of the system as a function of the pose of the end effector,

Based on the chosen parametrization, calculations for each leg are very similar. Therefore, only generic calculations are presented in the following section. Moreover, using the reduction of parametrization, presented in Eq. 9, the aim is equivalent to find only

Firstly, to obtain

From Eq. 23, it is possible to extract

Besides that, it can be shown that

where

As

The complete solution for a generic leg is then

It can be easily adapted for each leg, by replacing the generic parameters and variables with the proper ones. For instance, by replacing in the previous equations, *η* with *A*, namely *B* and *C*.

### 3.4 Differential Kinematics

The explicit inverse kinematics model is given as

where

the function

and

where the matrix

where for any function *f* and variable *f* with respect to

### 3.5 Static Equilibrium

Using the kinetostatic duality, we obtain

where

For a given external wrench **f**, acting on the end effector, the static equilibrium gives us, using Eq. 34, the relationship between the external wrench **f** and the torque of the actuated joints as

As we assume to be working outside of the singularities, the rank of the matrix

where *λ* is a scalar real and

The reader is invited to refer to Section 4.1 for details on the calculation of

### 3.6 Static Stiffness

To define the static stiffness, we first need to define

Define

where

Therefore,

where

Let’s define

Using the chain rule we have that

Based on Eq. 34, we have

and based on Eqs. 29 and 33, we have

Moreover, we have

where

Therefore, we obtain

### 3.7 Cartesian Compliance and Cost Function

To represent the stiffness characteristics of the system, we decide to represent in the space its Cartesian compliance, the quantity denoted by

where *p* in *f* represents only the forces applied on p in the directions *x*, *y* and *z*. This decision is motivated by the fact that, as shown in the previous section, the end-effector of the system has only 2 Dof, which means that the system will not allow motions in one of the three Cartesian directions of space, thus in that same direction the stiffness would be infinite. Therefore its inverse (i.e., the compliance) will be more manageable mathematically, and easier to represent graphically.

According to the minimal parametrization defined in Eq. 13, *p* is a function of

where the matrix

where *u* and *p*, from which we obtain in the static equilibrium

Moreover, using the definition of the static stiffness of the end effector in Eq. 42, we obtain

The matrix

with

So, by injecting Eq. 52 in Eq. 54, we obtain

which gives us, based on Eq. 50, the expression of the Cartesian compliance

Knowing the compliance matrix, it is possible to derive the cost function used in the optimization algorithm such as

where

### 3.8 Graphical Representation Methodology

To get a visual representation of the Cartesian compliance, we plot the ellipsoid associated with the matrix *c*. This ellipsoid is defined using the singular value decomposition of *c*. Its axes are directed by the eigenvectors of the matrix and their half lengths are defined by their respective associated eigenvalue. The ellipsoid is centered in *p*. Hence, the ellipsoid represents the possible deflection of the system in the space for a normalized external wrench, and thus the greater the ellipsoid is, the more compliant is the system.

Figure 5 shows an example of the various representations and views of the compliance ellipsoids for a specific configuration of the system. One way to represent the compliance ellipsoids is to draw them in 3 dimensions. They can be viewed either from a 3D perspective, such as in Figure 5A, or from a top view of the system [i.e., on a projection on the plane (*y*, *z*)] such as in Figure 5B. A second option is a representation in 2D, where the same ellipsoids of the Cartesian compliance are plotted in a plane *flattened representation*.

**FIGURE 5**. Cartesian compliance ellipsoids of the system **(A)** 3D view, **(B)** top view and **(C)** flattened (see text for detail) representations. Simulations are made at the center of the end effector, with

The analogous of meridians and parallels have been plotted also to help the understanding of the corresponding ellipses in the example figures. In each representation, the locus of all the possible points reachable of *p* is also plotted as a colored area even if the ellipsoids are displayed only in some discrete points. This area is referred to as the workspace of the system. The flattened representation has the advantage not to distort the ellipsoids and not to be sensitive to the Cartesian workspace which may be affected by the design parameters of the mechanism. So, this representation will mainly be used thereafter. Finally, the angular orientation of the axes of the first joint of each leg is plotted as small black lines in the different graphs.

Besides, one thing to notice is that the ellipsoids are only ellipses, tangent to the workspace. This is in agreement with the consideration made at the beginning of this subsection, which motivated the choice of representing the system compliance instead of its stiffness.

## 4 Results

### 4.1 Variable Stiffness Ability of the System

In this section, we prove that it is possible to have a decoupled control of the position and the stiffness of the system, i.e., that it is possible to modify the spring load without moving the end effector.

Let us consider first the static stiffness given in Eq. 48 and show that it can be modulated. For a given position

Similarly, the Cartesian compliance of the system defined in Eq. 57 can be controlled if the elastic elements are nonlinear by modifying the positions of the motors

Let us prove now that the compliance can be modulated at a fixed position. The first step is to derive a possible expression of

where

Therefore, *X* in

This system of equations has a space of solution of dimension 1 if and only if

Therefore, the Solutions of Eq. 38 are described as

A possible definition of

where

Therefore, along the space described by

The effect of the internal modulation is studied while varying the parameter *λ* defined in Eq. 37, with the model of the elastic elements defined in Eq. 1. It has to be noticed that with this definition of *λ*, as *λ*, using Eq. 37. Figure 6 represents the Cartesian compliance ellipses for various values^{1} of *λ*. We can see how, acting on *λ*, the ellipses scale homogeneously in all directions, without appreciable changes in shape and orientation like the co-contraction of human muscles (Milner, 2002; Perreault et al., 2002).

**FIGURE 6**. System Compliance in flattened representation, as a consequence of the internal modulation while *λ* is varying in the set {0 Nm, 2 Nm, 5 Nm, 10 Nm}. Simulations are made at the center of the end effector, with

Moreover, Figure 7 shows the evolution of the legs torque and stiffness, and the output force at three different positions of the end effector, highlighting how the actuation of *λ* changes the stiffness but does not change the external equilibrium. Indeed, the external wrench balanced at the end effector

**FIGURE 7**. Evolution of the torque of each leg *α*_{y} and *α*_{z}) for *λ* varying between -10 Nm and 10 Nm and for various equilibrium positions **(A)** corresponds to **(B)** to **(C)** to

As a final note, it is possible to show that the stiffness of the system is increasing when the external wrench is increasing. In addition, although the stiffness increases slightly in all directions, the increase is larger in the direction of the applied external force. For sake of space, the associated figures are not displayed here as this is an expected behavior for a VSA (Vanderborght et al., 2013).

### 4.2 General Design Guidelines

In this section, we give the results of the effects of the geometric parameters on the compliance of the system according to the methodology described in Section 2.3. To show their effects, the Cartesian compliance ellipsoids are plotted at the center of the end effector (i.e., *p* in the previous analysis is

Figure 8 represents the different figures of the simulations when the various design parameters are varying. Figures 8A,B show the results when the angular positions of the legs are varying in a symmetric or asymmetric way. We can see that in both cases, depending on the angular position, it is possible to stretch the ellipses in one of its two main directions by modifying the value of the angle and to rotate them. Figures 8C–E show the effects of *α* on the Cartesian compliance ellipses and the workspace. We can notice that the compliance ellipses are globally increasing when *α* is increasing. Figures 8F–H show the effects of *d* on the Cartesian compliance ellipses and the workspace. We can notice that the compliance ellipses are globally increasing when *d* is increasing.

**FIGURE 8**. Effects of the geometric parameters on the system compliance. **(A)** represents the system compliance in flattened representation, as the angular positions of the two legs *c* = (1 rad/Nm 1 rad/Nm 1 rad/Nm). **(B)** represents the system compliance in flattened representation, when the angular position of one leg *π*}. Simulations are made at the center of the end effector, with **(C)** and **(D)** represent the workspace of the system in 3D view and side view when *α* is varying in the set {**(E)** is the associated flattened representation of the compliance ellipsoids. Simulations are made at the center of the end effector, with **(F)** and **(G)** represent the workspace of the system in 3D view and side view when *d* is varying in the set {25 mm, 50 mm, 75 mm, and 100 mm}. **(H)** is the associated flattened representation of the compliance ellipsoids. Simulations are made at the center of the end effector, with

As it is shown in the different figures, all the design parameters do not have the same effect nor the same intensity on the compliance characteristic of the system. For instance, the effects of *α* are more important in terms of stretching the ellipses, especially at the borders of the workspace, whereas *d* has more a global influence everywhere that can be compensating by using different elastic elements. Moreover, it is important to notice that both previous results go together with a stretch of the surface of the feasible positions of the end effector, shown in Figures 8C,D,F,G. This explains that the increase in compliance is also an effect of the increase of the lever arm of the forces. Combining the previous two effects, it is conceivable to change the shape of the ellipses while maintaining the feasible surface and ellipses area roughly constant by changing alpha and d simultaneously in opposite directions. However, it appears that *α* and *d* are not enough to modulate the ratio of the ellipse at (0,0) position where we observe always a circle for an even distribution of the legs as shown in Figures 8E,H. This shows a limitation of using only these two parameters to modify the embedded behavior in the system and justify the use of the angular positions of the legs as geometric parameters.

### 4.3 Design Parameters for Human-Like Joints

In this section, we derive sets of parameters to match the passive compliant behavior of a human wrist (Pando et al., 2014) and a human ankle (Lee et al., 2014). In both cases, the idea is to find the best set of parameters to match the targeted ellipse at the neutral position, defined in Section 2.4.1. Regarding the input of the algorithm, we would like to match 3 characteristics. Therefore, with 12 design parameters the solutions are not unique. We can then set additional constraints besides the realistic values of the design parameters. As we are considering passive behavior, we set *λ* equal to 0. This means that in absence of external load, no motor torque is required to get the passive behavior of the system. Moreover, we consider an even distribution of the legs around the basis. Figure 9 represents solutions for a wrist and an ankle considering theses constraints. Their associated design parameters are reported in Table 4. Therefore, it is possible to get feasible systems based on the proposed architecture with a passive compliant behavior similar to the human wrist or a human ankle.

**FIGURE 9**. **(A)** Compliance ellipses at the neutral position in absence of external load for a wrist joint. The dashed black line stands for the targeted ellipse extracted from human biomechanical data in Pando et al. (2014). The blue line represents compliant behavior of the proposed system with an even distribution of the legs around the basis. **(B)** Compliance ellipses at the neutral position in absence of external load for an ankle joint. The dashed black line stands for the targeted ellipse extracted from human biomechanical data in Lee et al. (2014)). The blue line represents compliant behaviour of the proposed system with an even distribution of the legs around the basis.

## 5 Discussion

In this work, we proposed and studied a new concept to design a 2 DoF VSA. To reduce the complexity of the mechanical implementation and increase the compactness of the system, we use only 3 motors. This is the minimum number of motors to be able to have a 2 DoF mechanism fully actuated in positions and partially actuated in stiffness. Therefore, our solution is more compact than most of the solutions existing in the current state-of-the-art (refer for instance to the works of (Grebenstein et al., 2011; Romano et al., 2014; Weckx et al., 2014; Koganezawa, 2017; Tan et al., 2017)). To the best knowledge of the authors, only the work in (Stoeffler et al., 2018) follows the same approach of using only 3 motors.

Obviously, a limitation of our proposed solution relies in the loss of DoF in the control of stiffness of the system with respect to the existing 2 DoF VSAs. Indeed, our system can control the overall size of its endpoint stiffness ellipse but not its stiffness independently in the two directions of motion. Yet this approach is motivated by the fact that human beings modulate their stiffness through both their posture (exploiting the geometry and properties of their limbs) and co-contraction (Milner, 2002; Perreault et al., 2002). More interestingly, their effects are quite distinct as the limb geometry and posture have a large impact on the shape and orientation of the ellipse and the co-contraction on the overall size (Perreault et al., 2002). Therefore, by implementing physically the human joint geometry and allowing the control of the overall size of the stiffness, we can design human-like joints with the potential of natural control of both the position and stiffness. To the best knowledge of the authors, there is no other work following this approach.

However, the challenge of the mechanical implementation of the passive compliance of human joints remains. That is why we derived a model of our system and studied the effect of its main design parameters on the characteristics of the compliance ellipses. Additionally to the guidelines given in Section 4.2, it appears that the influences of the angular positions of the leg are interesting as they allow us to stretch or rotate the compliance ellipses in certain directions. Combining all of the effects, it will be possible to embed some specific behavior in the design for a specified application. And indeed, we derived sets of design parameters to match the passive behavior of a human wrist and a human ankle.

For now, biomechanical data exists only at the neutral position of the joints. However, in the future, additional data could be available on various positions of the workspace of the joints. In this case, we may want to match several ellipses. In the proposed system, it is possible to increase the number of desired ellipses. Figures 10A,B present solutions for 2 desired ellipses. The case 1.1 corresponds to a solution with an additional constraint of the vector *λ* equal to 0. As it is shown, this solution may be insufficient. Case 1.2 stands for a solution without any additional constraints. In this case, the two ellipses are well matched. It shows that the system could match additional passive behaviors, for instance by relaxing the constraint of passivity, previously explained (*λ* equal to 0). The associated design parameters of these 2 cases are given in Table 5.

**FIGURE 10**. **(A)** (case 1.1) represents the compliance ellipses at the neutral position in absence of external load with *λ* null. **(B)** (case 1.2) represents the compliance ellipses at the neutral position in absence of external load without any additional constraints. **(C)** (case 2.1) represents the compliance ellipses in a large portion of the workspace. In all figures the dashed black lines stand for the targeted ellipses and the blue lines stand for the solutions.

To go further, we could also apply the proposed concept to more complex joints (such as the shoulder or the hip). Unfortunately, such biomechanical data are not yet available. We are more than aware that this is not an easy task and we can only encourage the biomechanical scientific community to pursue their current efforts in this direction, such as the work proposed in Hunt and Lee (2019).

As another further step, we could imagine having a lot of desired ellipses. The case 2.1 for instance corresponds to an isotropic compliant behavior of the system in a large portion of the workspace. Such a system could be conceived for applications using teleoperated robots where we would like to have the same behavior over the workspace. Concretely, the required behavior is now a compliance circle over 25 positions over a grid of 5 × 5 positions between −60° and 60° in both directions. In this case, we have thus 75 characteristics to match. Similar results are obtained if the vector^{2}*λ* is free and if *λ* is set to 0, ensuring thus the passivity of the compliant behavior. Therefore, only the results when *λ* is null, noted **0**^{25}, are displayed here. Figure 10C represents the result of the optimization algorithm and the parameters are given in Table 5. As it is shown in Figure 10C, the behavior is not exactly matched but it seems acceptable over a smaller range of the workspace.

Obviously, in the case of several desired ellipses, there are some mechanical limitations but we can refer to Section 4.2 to get the general trends and to design useful desired behavior. And it seems possible to design artificial joints that match the passive compliant behavior of human articulations in various positions of the workspace, to investigate the potential of such features in robotic systems.

Another limitation of this study is that we did not explore the space of all the possible spring characteristics. Indeed, we used a model (referred to as hyperbolic sine springs) that we already implemented in previous works (Catalano et al., 2011; Lemerle et al., 2019), as similar elastic behaviors can reproduce some attributes of an antagonistic pair of muscles driving a joint (Garabini et al., 2017). However, to improve the human-likeliness of the output stiffness characteristic, we could study more extensively on the best nonlinear spring characteristic that should be implemented. This was out of the intended scope of this work, but future studies could investigate further on this topic using additional biomechanical data to implement human-like stiffness behaviors.

## 6 Conclusion

This paper presents the concept of a new configurable 2 DoF variable stiffness joint. The system is based on redundant parallel nonlinear elastic actuation, using the antagonistic approach inspired by the human musculoskeletal system. The kinematic structure of the system is based on the architecture of the Omni-Wrist III, designed by Rosheim et al. and we propose to add in series nonlinear elastic systems to obtain a VSA. To get a compact design, only three motors are used and we prove that this allows us to control both the position and the overall stiffness of system. A one-dimensional quantity linked to the stiffness can thus be controlled like the voluntary co-contraction of human muscles. Moreover, we outline the impact of the main design parameters on the compliant behavior of the system as general guidelines. We describe an optimization methodology to fine-tune the mechanical implementation of this type of system to match the specific passive behavior of a human wrist and ankle.

We believe that the proposed design could find many applications, including industrial manipulators and prosthetic devices for various joints, due to its versatility. The sets of parameters obtained to match the passive behavior of a human wrist are particularly relevant as they open the possibility to design a compact system with an embedded passive behavior close to the human one. A new prototype of a compact physically compliant variable stiffness wrist is currently under development.

## Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

## Author Contributions

SL, MC, and GG designed and analyzed the proposed system. All authors contributed to the definition of the possible applications of the system. SL and AB contributed to defining the optimization problem of matching a target stiffness behavior. All authors contributed to the analysis of the state of the art. SL and GG wrote all sections of the manuscript. AB contributed expertise and advice. All authors contributed to manuscript revision, read, and approved the submitted version.

## Funding

This research has received funding from the European Union’s ERC program under the Grant Agreement No. 810346 (Natural Bionics). The content of this publication is the sole responsibility of the authors. The European Commission or its services cannot be held responsible for any use that may be made of the information it contains.

## Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Acknowledgments

The authors would like to thank Giuseppe Averta for his advice on the applications of the proposed system.

## Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/frobt.2021.614145/full#supplementary-material.

## Footnotes

^{1}*λ* is only varying in the positive values as the stiffness of each elastic element is an even function in *λ* in case of a hyperbolic sine spring.

^{2}In this example, the vector *λ* is in

^{3}As a preliminary remark, we should notice that

## References

Ajoudani, A., Fang, C., Tsagarakis, N., and Bicchi, A. (2018). Reduced-complexity representation of the human arm active endpoint stiffness for supervisory control of remote manipulation. *Int. J. Robot Res.* 37, 155–167. doi:10.1177/0278364917744035

Albu-Schaeffer, A., Eiberger, O., Grebenstein, M., Haddadin, S., Ott, C., Wimbock, T., et al. (2008). Soft robotics. *IEEE Robot. Autom. Mag.* 15, 20–30. doi:10.1109/MRA.2008.927979

Bajaj, N. M., Spiers, A. J., and Dollar, A. M. (2019). State of the art in artificial wrists: a review of prosthetic and robotic wrist designs. *IEEE Trans. Robot.* 35, 261–277. doi:10.1109/TRO.2018.2865890

Barrett, E., Reiling, M., Barbieri, G., Fumagalli, M., and Carloni, R. (2017). Mechatronic design of a variable stiffness robotic arm. *IEEE/RSJ Int. Conf. Intell. Robots Syst. (IROS)* 7, 4582–4588. doi:10.1109/iros.2017.8206327

Blank, A. A., Okamura, A. M., and Whitcomb, L. L. (2014). Task-dependent impedance and implications for upper-limb prosthesis control. *Int. J. Robot Res.* 33, 827–846. doi:10.1177/0278364913517728

Casini, S., Tincani, V., Averta, G., Poggiani, M., Della Santina, C., Battaglia, E., et al. (2017). Design of an under-actuated wrist based on adaptive synergies. *IEEE Int. Conf. Robot. Autom. (ICRA)* 12, 6679–6686. doi:10.1109/icra.2017.7989789

Catalano, M. G., Grioli, G., Garabini, M., Bonomo, F., Mancini, M., Tsagarakis, N., et al. (2011). Vsa-cubebot: a modular variable stiffness platform for multiple degrees of freedom robots. *IEEE Int. Conf. Robot. Autom.* 3, 5090–5095. doi:10.1109/icra.2011.5980457

English, C., and Russell, D. (1999). Mechanics and stiffness limitations of a variable stiffness actuator for use in prosthetic limbs. *Mech. Mach. Theor.* 34, 7–25. doi:10.1016/S0094-114X(98)00026-3

Formica, D., Charles, S. K., Zollo, L., Guglielmelli, E., Hogan, N., and Krebs, H. I. (2012). The passive stiffness of the wrist and forearm. *J. Neurophysiol.* 108, 1158–1166. doi:10.1152/jn.01014.2011

Garabini, M., Passaglia, A., Belo, F., Salaris, P., and Bicchi, A. (2011). Optimality principles in variable stiffness control: the vsa hammer. *IEEE/RSJ Int. Conf. Intell. Robots Syst.* 23, 3770–3775. doi:10.1109/iros.2011.6048555

Garabini, M., Santina, C. D., Bianchi, M., Catalano, M., Grioli, G., and Bicchi, A. (2017). *Soft robots that mimic the neuromusculoskeletal system converging clinical and engineering research on neurorehabilitation II*. Berlin: Springer International Publishing, 259–263.

Grebenstein, M., Albu-Schäffer, A., Bahls, T., Chalon, M., Eiberger, O., Friedl, W., et al. (2011). The dlr hand arm system. *IEEE Int. Conf. Robot. Autom.* 7, 3175–3182.

Haddadin, S., Albu-Schäffer, A., and Hirzinger, G. (2007). Safety evaluation of physical human–robot interaction via crash-testing. *Robot. Sci. Syst.* 21, 33–39. doi:10.7551/mitpress/7830.003.0029

Hunt, J., and Lee, H. (2019). Development of a low inertia parallel actuated shoulder exoskeleton robot for the characterization of neuromuscular property during static posture and dynamic movement. *Int. Conf. Robot. Automat. (ICRA)* 33, 556–562. doi:10.1109/icra.2019.8794181

Kanitz, G., Montagnani, F., Controzzi, M., and Cipriani, C. (2018). Compliant prosthetic wrists entail more natural use than stiff wrists during reaching, not (necessarily) during manipulation. *IEEE Trans. Neural Syst. Rehabil. Eng.* 26, 1407–1413. doi:10.1109/TNSRE.2018.2847565

Kearney, R. E., Stein, R. B., and Parameswaran, L. (1997). Identification of intrinsic and reflex contributions to human ankle stiffness dynamics. *IEEE (Inst. Electr. Electron. Eng.) Trans. Biomed. Eng.* 44, 493–504. doi:10.1109/10.581944

Koganezawa, K., Inomata, H., and Nakazawa, T. (2005). Actuator with non-linear elastic system and its application to 3 dof wrist joint. *Int. Conf. Mechatron. Autom.* 19, 1253–1260. doi:10.1109/icma.2005.1626733

Koganezawa, K. (2017). Mechanical stiffness control of a three dof wrist joint that mimics musculo-skeletal system. *IEEE Int. Conf. Adv. Intell. Mechatron. (AIM)* 15, 46–51. doi:10.1109/aim.2017.8013993

Lee, H., Krebs, H. I., and Hogan, N. (2014). Multivariable dynamic ankle mechanical impedance with relaxed muscles. *IEEE Trans. Neural Syst. Rehabil. Eng.* 22, 1104–1114. doi:10.1109/TNSRE.2014.2313838

Lemerle, S., Grioli, G., Bicchi, A., and Catalano, M. G. (2019). A variable stiffness elbow joint for upper limb prosthesis. *IEEE/RSJ Int. Conf. Intell. Robots Syst. (IROS)* 14, 7327–7334. doi:10.1109/iros40897.2019.8970475

Malosio, M., Corbetta, F., Reyes, F. R., Giberti, H., Legnani, G., and Tosatti, L. M. (2019). On a two-dof parallel and orthogonal variable-stiffness actuator: an innovative kinematic architecture. *Robotics* 8, 133. doi:10.3390/robotics8020039

Masia, L., and Squeri, V. (2015). A modular mechatronic device for arm stiffness estimation in human–robot interaction. *IEEE ASME Trans. Mechatron.* 20, 2053–2066. doi:10.1109/TMECH.2014.2361925

Masia, L., Squeri, V., Sandini, G., and Morasso, P. (2012). Measuring end-point stiffness by means of a modular mechatronic system. *IEEE Int. Conf. Robot. Autom.* 21, 2471–2478. doi:10.1109/icra.2012.6225140

Mason, M. T., and Salisbury, J. K. (1985). *Robot hands and the mechanics of manipulation*. Cambridge, MA: MIT Press.

Mazzotti, C., Troncossi, M., and Castelli, V. P. (2015). “*A new powered humeral rotator for upper limb myoelectric prostheses*”, in 14th World congress in mechanism and machine science, IFTOMM, 307–312.

Milner, T. (2002). Contribution of geometry and joint stiffness to mechanical stability of the human arm. *Exp. Brain Res.* 143, 515–519. doi:10.1007/s00221-002-1049-1

Montagnani, F., Controzzi, M., and Cipriani, C. (2013). Preliminary design and development of a two degrees of freedom passive compliant prosthetic wrist with switchable stiffness. *IEEE Int. Conf. Robot. Biomimetr. (ROBIO)* 7, 310–315. doi:10.1109/robio.2013.6739477

Mussa-Ivaldi, F., Hogan, N., and Bizzi, E. (1985). Neural, mechanical, and geometric factors subserving arm posture in humans. *J. Neurosci.* 5, 2732–2743. doi:10.1523/JNEUROSCI.05-10-02732.1985

Negrello, F., Mghames, S., Grioli, G., Garabini, M., and Catalano, M. G. (2019). A compact soft articulated parallel wrist for grasping in narrow spaces. *IEEE Robotics and Automation Letters* 4, 3161–3168. doi:10.1109/LRA.2019.2925304

Pando, A. L., Lee, H., Drake, W. B., Hogan, N., and Charles, S. K. (2014). Position-dependent characterization of passive wrist stiffness. *IEEE (Inst. Electr. Electron. Eng.) Trans. Biomed. Eng.* 61, 2235–2244. doi:10.1109/TBME.2014.2313532

Perreault, E. J., Kirsch, R. F., and Crago, P. E. (2002). Voluntary control of static endpoint stiffness during force regulation tasks. *J. Neurophysiol.* 87, 2808–2816. doi:10.1152/jn.2002.87.6.2808

Petit, F., Friedl, W., Hoppner, H., and Grebenstein, M. (2015). Analysis and synthesis of the bidirectional antagonistic variable stiffness mechanism. *IEEE Transactions on Mechatronics* 20, 684–695. doi:10.1109/TMECH.2014.2321428

Rodriguez-Cianca, D., Weckx, M., Jimenez-Fabian, R., Torricelli, D., Gonzalez-Vargas, J., Sanchez-Villamañan, M., et al. (2019). A variable stiffness actuator module with favorable mass distribution for a bio-inspired biped robot. *Front. Neurorob.* 13, 20. doi:10.3389/fnbot.2019.00020

Romano, F., Fiorio, L., Sandini, G., and Nori, F. (2014). Control of a two-dof manipulator equipped with a pnr-variable stiffness actuator. *IEEE Int. Symp. Intell. Control* 14, 1354–1359. doi:10.1109/isic.2014.6967620

Rosheim, M. E., and Sauter, G. F. (2002). “New high-angulation omnidirectional sensor mount,” in *Free-space laser communication and laser imaging II*. Editors J. C. Ricklin, and D. G. Voelz (New York, NY: SPIE), 163–174.

Savin, S., Golousov, S., Khusainov, R., Balakhnov, O., and Klimchik, A. (2019). Control system design for two link robot arm with maccepa 2.0 variable stiffness actuators. *IEEE/ASME Int. Conf. Adv. Intell. Mechatron. (AIM)* 4, 624–628. doi:10.1109/aim.2019.8868490

Sensinger, J. W., and ff. Weir, R. F. (2008). User-modulated impedance control of a prosthetic elbow in unconstrained, perturbed motion. *IEEE (Inst. Electr. Electron. Eng.) Trans. Biomed. Eng.* 55, 1043–1055. doi:10.1109/TBME.2007.905385

Siciliano, B., and Khatib, O. (2007). *Springer handbook of robotics*. Berlin, Heidelberg: Springer-Verlag.

Sofka, J., Skormin, V., Nikulin, V., and Nicholson, D. J. (2006a). Omni-wrist iii—a new generation of pointing devices. part i. laser beam steering devices—mathematical modeling. *IEEE Trans. Aerospace Electron. Syst.* 42, 718–725. doi:10.1109/TAES.2006.1642584

Sofka, J., , Skormin, V., , Nikulin, V., , and Nicholson, D. J., (2006b). Omni-wrist iii—a new generation of pointing devices. part ii. gimbal systems—controls. *IEEE Trans. Aerospace Electron. Syst.* 42, 726–734. doi:10.1109/TAES.2006.1642585

Stillfried, G., Stepper, J., Neppl, H., Vogel, J., and Höppner, H. (2018). Elastic elements in a wrist prosthesis for drumming reduce muscular effort, but increase imprecision and perceived stress. *Front. Neurorob.* 12, 141. doi:10.3389/fnbot.2018.00009

Stoeffler, C., Kumar, S., Peters, H., Bruls, O., Muller, A., and Kirchner, F. (2018). “*Conceptual design of a variable stiffness mechanism in a humanoid ankle using parallel redundant actuation*,” in IEEE-RAS 18th international conference on humanoids robots, (Humanoids), 462–468.

Tan, D. J., Brouwer, D. M., Fumagalli, M., and Carloni, R. (2017). A 2-dof joint with coupled variable output stiffness. *IEEE Robot. Autom. Lett.* 2, 366–372. doi:10.1109/LRA.2016.2631730

Toedtheide, A., and Haddadin, S. (2018). Cpa-wrist: compliant pneumatic actuation for antagonistic tendon driven wrists. *IEEE Robot. Autom. Lett.* 3, 3537–3544. doi:10.1109/LRA.2018.2853802

Troncossi, M., Gruppioni, E., Chiossi, M., Cutti, A., Davalli, A., and Castelli, V. P. (2009). A novel electromechanical shoulder articulation for upper-limb prostheses: from the design to the first clinical application. *J. Prosth. Orth.* 21, 79–90. doi:10.1097/JPO.0b013e31819f6aed

Vanderborght, B., Albu-Schaeffer, A., Bicchi, A., Burdet, E., Caldwell, D., Carloni, R., et al. (2013). Variable impedance actuators: a review. *Robot. Autonom. Syst.* 61, 1601–1614. doi:10.1016/j.robot.2013.06.009

Vanderborght, B., Tsagarakis, N. G., Semini, C., Van Ham, R., and Caldwell, D. G. (2009). Maccepa 2.0: adjustable compliant actuator with stiffening characteristic for energy efficient hopping. *IEEE Int. Conf. Robot. Autom. (IEEE)* 121, 544–549. doi:10.1109/robot.2009.5152204

Visser, L. C., Stramigioli, S., and Bicchi, A. (2011). Embodying desired behavior in variable stiffness actuators. *IFAC Proc.* 44, 121–124. doi:10.3182/20110828-6-it-1002.02202

Weckx, M., Mathijssen, G, Benali, I. S. M., Furnemont, R, Ham, R. V., Lefeber, D., et al. (2014). A two-degree of freedom variable stiffness actuator based on the maccepa concept. *Actuators* 3, 20–40. doi:10.3390/act3020020

Wilhelms, J., and Gelder, A. V. (2001). Fast and easy reach-cone joint limits. *J. Graph. Tool.* 6, 27–41. doi:10.1080/10867651.2001.10487539

Keywords: artificial joints, articulated soft robotics, variable stiffness, humanoids, prostheses

Citation: Lemerle S, Catalano MG, Bicchi A and Grioli G (2021) A Configurable Architecture for Two Degree-of-Freedom Variable Stiffness Actuators to Match the Compliant Behavior of Human Joints. *Front. Robot. AI* 8:614145. doi: 10.3389/frobt.2021.614145

Received: 05 October 2020; Accepted: 12 January 2021;

Published: 12 March 2021.

Edited by:

Nadia Magnenat Thalmann, Université de Genève, SwitzerlandReviewed by:

Larisa Dunai Dunai, Universitat Politècnica de València, SpainMaja Borry Zorjan, Université de Versailles Saint-Quentin-en-Yvelines, France

Copyright © 2021 Lemerle, Catalano, Bicchi and Grioli. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Simon Lemerle, simon.lemerle@iit.it