Next Article in Journal
Shape Memory Polymer Composite Actuator: Modeling Approach for Preliminary Design and Validation
Previous Article in Journal
A Multistate Friction Model for the Compensation of the Asymmetric Hysteresis in the Mechanical Response of Pneumatic Artificial Muscles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling, Identification, and Control of a Discrete Variable Stiffness Actuator (DVSA)

by
Irfan Hussain
1,*,†,
Ahmad Albalasie
2,†,
Mohammad I. Awad
1,† and
Dongming Gan
1
1
Khalifa University Center for Autonomous Robotic Systems (KUCARS), Khalifa University of Science and Technology, P.O. Box 127788, Abu Dhabi, United Arab Emirates
2
Department of Mechanical and Mechatronics Engineering, Birzeit University, P.O Box 14, Birzeit West Bank, Palestine
*
Author to whom correspondence should be addressed.
These authors have contributed equally to this work.
Actuators 2019, 8(3), 50; https://doi.org/10.3390/act8030050
Submission received: 23 April 2019 / Revised: 25 May 2019 / Accepted: 28 May 2019 / Published: 27 June 2019

Abstract

:
A branch of robotics, variable impedance actuation, along with one of its subfields variable stiffness actuation (VSA) targets the realization of complaint robotic manipulators. In this paper, we present the modeling, identification, and control of a discrete variable stiffness actuator (DVSA), which will be developed for complaint manipulators in the future. The working principle of the actuator depends on the involvement of series and parallel springs. We firstly report the conceptual design of a stiffness varying mechanism, and later the details of the dynamic model, system identification, and control techniques are presented. The dynamic parameters of the system are identified by using the logarithmic decrement algorithm, while the control schemes are based on linear quadratic control (LQR) and computed torque control (CTC), respectively. The numerical simulations are performed for the evaluation of each method, and results showed the good potentialities for the system. Future work includes the implementation of the presented approach on the hardware.

1. Introduction

Safety is the most important element during human–robot interactions in smart manufacturing and domestic scenarios [1,2,3]. Bio-inspired robotics is inspired by soft manipulation with improvement towards human–robot interaction [4,5,6]. Remarkable progress has been made towards complaint robots, and human–robot interaction (HRI) has become the most important feature in the robot design and development. Although far from reaching the performance of their biological counterpart, robots have been offering a wide range of applications in many different fields, from medicine to industry [7]. Human–robot interaction (HRI) aims at developing hardware (actuation) and software (algorithms) to allow more natural and effective communication and interaction between humans and robots [2,8,9]. The complaint actuators are growing in the scientific community because of the performance limitations of traditional stiff actuation approaches in terms of safety, energy efficiency, and the ability to interact with the environment [10].
The possible advantages that can be obtained with complaint actuators are very well known, but most of the robots in industry still lack these features. Towards this direction, many complaint actuation systems have been proposed in the last two decades. The research on complaint manipulators was driven by the need for more safety and robustness in human–robot interactions [1,11]. Robustness and soft interaction are mainly achieved by either actively regulating the compliance of robotic joints [1,2] or tuning the intrinsic softness, acting on the passive characteristics of the robot body ware [12,13,14]. Incorporating complaint elements in actuators comes with many challenges from the mechatronic and control points of view. For this reason, the use of complaint actuators in robotic manipulators is still limited and challenging [11].
Many designs and prototypes of variable stiffness actuators are presented by the research community. Some remarkable examples are VSA-I [15], VSA-II [16], AMASC [17], and the biologically inspired joint stiffness control mechanism [6]. An alternative approach of varying the stiffness was achieved through the principle of a lever mechanism [6,18,19,20,21,22,23]. Another new method of stiffness variation was proposed through the addition and subtraction of elastic element involvement in actuator mechanism. The arrangements of elastic elements can be done in series [24] or in parallel, which are also called series–parallel elastic actuators (SPEAs). The approach of series–parallel arrangement is efficient in terms of energy and peak torque [25,26,27]. Based on the concept of SPEAs, we firstly proposed a passive version of the discrete variable stiffness actuator (DVSA) (see Figure 1) [28]. In our all previous works, we mainly focused on the passive versions of the actuators. In this regard, we proposed various design topologies and optimized their design specifications according to the intended applications ranging from industrial to haptics and medical. Following this motivation, in [24], we presented an approach of realizing the arm exoskeleton (TELEXOS-I) for applications such as haptic-teleoperation and rehabilitation. The key motivation behind our design topology is the need for instantaneous switching between stiffness levels when applied for virtual reality (VR) or remote environment (RE) stiffness mapping applications, in addition to the need of low-energy consumption. Moreover, in [29], we presented a case scenario of using another developed passive discrete variable stiffness joint (pDVSJ-II) in a remote palpation application. Palpation is defined as the application of the fingers with light pressure to the surface of the human body for the purpose of determining the condition of the parts beneath during physical diagnosis. Palpation with a robotic device, or the artificial tactile exploration of soft viscoelastic and nonhomogeneous objects, is an important area of study for various fields including medical, VR, or teleoperation applications.
Now that the ability of passive actuators has been successfully demonstrated to change their stiffness at different levels using the proposed topologies, we are in the process of developing a compact and active version (with a motor). In [30], we introduced the concept of the active binary-controlled variable stiffness actuator (BcVSA), which is an example that falls under the category of DVSA. One of the benefits of our approach is the selection of a desired level of stiffness without the need for succession involvement, which improves the response time in switching the level of stiffness. In order to control the actuator, we presented the model predictive control (MPC) technique. The control design and identification of the dynamic parameters are quite challenging tasks and play a great role in the efficiency of the actuator. In fact, system identification is an important and crucial problem in a control system. Any uncertainty in the dynamic parameters will eventually affect the performance of the system if the controller does not address it. In the state of the art, several algorithms for identifying the dynamic parameters have been introduced to enhance the performance of the system. One of these algorithms is the logarithmic decrement algorithm, and it is introduced with suitable measuring techniques in [31,32]. Other publications have concentrated on identifying dynamic parameters for various types of robots based on the logarithmic decrement algorithm, which is presented in [33,34].
In this paper, we use the same concept of the BcVSA for varying stiffness but propose system identification and different control techniques. Thus, the contribution of this paper lies in identifying the dynamic parameters of our active actuator and implementing new suitable control techniques. The aim is to evaluate the ability of the proposed controllers using the identified dynamic parameters to control our actuator and to test the controllers robustness in different cases while changing the actuator’s stiffness in two cases—offline (when the system is stopped) and online (while the system is running). Both methods are simulated using actuator dynamics under the effect of various combinations of active (engaged) and inactive (disengaged) clutches. In particular, we have used the logarithmic decrement algorithm to perform system identification and have evaluated two control schemes using the identified parameters in the model. The control techniques implemented are the linear quadratic regulator (LQR) [12,13] and the computed torque control (CTC) [14,35,36,37]. We performed numerical simulations to evaluate the effectiveness of our approach. The purpose of the numerical simulations was to choose the proper control scheme for experiments on the hardware of the system in the near future.
The rest of the paper is organized as it follows. The conceptual design and the stiffness modeling are presented in Section 2. Section 3 describes the system identification by using a logarithmic decrement algorithm. Section 4 presents the two control schemes based on linear quadratic control (LQR) and computed torque control (CTC), respectively. The simulation results are reported in Section 5. Moreover, conclusions and future work are detailed in Section 6.

2. Design and Stiffness Model of the Discrete Variable Stiffness Actuator

The motivation for the proposed actuator lies in the series–parallel variable stiffness mechanism, as shown in Figure 2, which shows the conceptual 3D CAD representation and schematic illustration of the proposed discrete variable stiffness actuator. A DC motor is coupled with the sun gear of the first planetary gear set, which is engaged with its plane gears. The torsional springs are serially connected to the inline clutches at one end and, to another planet, gears of the second planet gear set. The sun gear of this planet gear stage drives the load shaft (see Figure 2a). If any of the clutches are active, the actuator torque is transmitted to load through the corresponding spring. Here, a stiffness bit consists of a clutch, a torsional spring, and a planet gear, which transmit the power from actuator to load. Several stiffness bits can be connected in parallel to achieve the overall stiffness of the actuator.
The novelty of our discrete variable stiffness actuator lies in the design topology (see Figure 2a). Three springs with different stiffness values each are serially attached to clutches. The clutches are connected to the planet gears which in turn are attached to the sun gear. The other end of each spring is connected to a planet gear that contributes part of the resultant torque on the other sun gear. The load shaft is attached to the sun gear of the system. If any clutch is disengaged, the corresponding spring will be detached. The achievable levels of stiffness in binary representation are reported in Table 1, where “0” represents inactive stiffness and “1” active stiffness. The clutch reduces the stiffness switching time. Moreover, scalability can be achieved either by changing the value of the lower stiffness (K0) or by adding extra stiffness-bits.
In order to derive the stiffness model, we can use the actuator’s kinematics. If the motor’s shaft is clamped (blocked from motion) and an external load is applied, the external torque will rotate the sun gear. In the scenario where no stiffness bit is active, the spring will rotate freely, while in the scenario where any of the stiffness bits are active, the load torque will be transmitted through the corresponding spring, which deflects and produces a counter torque on the human hand.
The resultant torque can be written as follows:
τ Σ = ( τ 0 + τ 1 + τ 2 )
where τ 0 , τ 1 ,   and   τ 2 are the torque transmitted by each planet gear to the sun gear, which can be further written as
τ n = N n β n ( 2 n ( K 0 ) ( θ ϕ n φ ) ) , n { 0 , 1 , 2 }
ϕ n = θ ( t O N , n ) , n { 0 , 1 , 2 }
β n = { 0 ,                           if   the   clutch   is   off 1 ,                           if     the   clutch   is   on
where:
  • N: the gear ratio;
  • β: binary number;
  • θ: the angular position at present time;
  • ϕ: the angular position at activation time;
  • φ: the backlash angle.
The total torque ( τ Σ ) and the total torsional stiffness ( K Σ ) can be obtained as
τ Σ = 0 n N n β n ( 2 n ( K 0 ) ( θ ϕ n φ ) ) , n = { 0 , 1 , 2 }
K Σ = δ τ Σ δ θ = 0 n N n β n ( 2 n ( K 0 ) ) , n = { 0 , 1 , 2 } .
As can be seen in Equation (5), the involvement of a spring does not depend on other springs, which means that the level of stiffness can be uniquely selected at any position without the need of returning to the equilibrium position. Moreover, the system is scalable.

3. System Identification of Discrete Variable Stiffness Actuator

This section focuses on identifying the unknown dynamical parameters of the discrete variable stiffness actuator model. If we assume that there is no gears backlash, the dynamic equations of the system can be written as
I m θ ¨ m + B m θ ˙ m + N K e q ( θ m θ L ) = T i n
I L θ ¨ L + B L θ ˙ L + N K e q ( θ L θ m ) = T e x t
where:
  • θ m , θ ˙ m ,     and     θ ¨ m are the angular position, the angular velocity, and the angular acceleration of the motor, respectively;
  • θ L , θ ˙ L ,     and   θ ¨ L are the angular position, the angular velocity, and the angular acceleration of the motor, respectively;
  • T i n is the input torque of the motor;
  • T e x t is the torque needed to lift the external load.
Note that some of the dynamic parameters shown in the model are fixed in the design phase. In particular, the stiffness of each spring (hence K e q ) and the gear ratio (N). However, our system may have eight possible values for K e q depending upon the status of clutches, as shown in Table 1. The goal of this section is to identify the remaining unknown model parameters—notably, the load inertia ( I L ) , the motor inertia ( I m ), the load damping ( B L ), and the motor damping ( B m ).
The parameters identification is carried out in two steps. At both steps, the two-Degrees-of-Freedom dynamic model is converted into a single-DOF one by constraining the other. Firstly, we clamp θ L and set T i n equal to zero to identify I m and B m . Thus, the differential equation for the system can be written as
θ ¨ m + B m I m θ ˙ m + N K e q I m ( θ m ) = 0 .
The system is excited by any initial condition (IC), e.g., θ m ( 0 ) = 0.1 rad, and the response of θ m is recorded. The logarithmic decrement algorithm is applied to identify both parameters. For a more detailed description on the logarithmic decrement algorithm and how it works, the reader is encouraged to refer to [31,32,33,34]. The system vibrates under the effect of the IC while the vibrations ( θ m ) are recorded. As a result, the periodic time ( τ d 1 ), the damped natural frequency ( ω d 1 ), and the natural frequency ( ω n 1 ) are calculated from the response of θ m . On the other hand, the damping ratio ( ζ 1 ) of the system is calculated by using the following equation:
δ 1 = 1 m   ln   ( θ 1 θ m + 1 )
where:
  • θ 1 : the peak amplitude corresponding to time t 1 ;
  • θ m + 1 : the peak amplitude corresponding to time t m + 1 ;
  • m: number of complete cycles between θ 1 and θ m + 1 ;
  • δ 1 : the logarithmic decrement for the motor side under the effect of initial condition (IC) θ m ( 0 ) .
Therefore,
ζ 1 = δ 1 ( δ 1 2 + ( 2 π ) 2 )
where ζ1 is the damping ratio.
As a result, I m and B m are calculated using the following equations:
I m = N K e q ω n 1 2
B m = 2 ζ 1 ω n 1 I m .
The whole process is repeated to identify I L and B L , but in this case θ m is clamped while keeping θ L free. The dynamic system as previously becomes a single DOF system, where it is excited by IC ( θ L ( 0 ) = θ L 0 = 0.1 rad). Hence, τ d 2 , ω d 2 , ω n 2 , δ 2 , and ζ 2 are calculated from the response of θ L . Similarly, I L and B L are calculated using the following equations:
I L = N K e q ω n 2 2
B L = 2 ζ 2 ω n 2 I L
where:
  • I L : the inertia load;
  • B L : the damping in the load side;
  • ζ 2 : the damping ratio for the load side under the effect of the IC θ L ( 0 ) ;
  • ω n 2 : the natural frequency under the constraints of the second experiment.

4. Control System Design of Discrete Variable Stiffness Actuator

In this section, two control techniques are implemented on the dynamic model of the DVSA. The first control method is based on a type-one servo system design using linear quadratic regulator (LQR). This controller is classified as an optimal controller because it minimizes a quadratic cost function with respect to the dynamics of the system, and it has the capability to maximize energy saving as well as coping with the tracking of the desired response, as detailed in [12,13]. On the other hand, the computed torque control (CTC) scheme is classified as a model-based control. In the following, we present both methods and numerical simulations to evaluate the performance of each for our DVSA.

4.1. Type-One Servo System Design Based on LQR

The system dynamic model is linear as shown in Equations (7) and (8). The idea is to construct a tracking system using LQR control to track the reference signals θ m and θ L . Thus, differential equations of our system can be written as
θ ¨ m = T i n B m θ ˙ m + N K e q θ L N K e q θ m I m
θ ¨ L = T e x t B L θ ˙ L N K e q θ L + N K e q θ m I L .
The corresponding state space representation for the binary controlled variable stiffness actuator (BcVSA) can be derived as follows:
q 1 = θ m q 2 = θ ˙ m q 3 = θ L q 4 = θ ˙ L
where q = [ q 1 q 2 q 3 q 4 ] T is called the state vector.
Therefore, the state space representation is defined as follows:
q ˙ ( t ) = A q ( t ) + B u ( t )
y 1 ( t ) = C q ( t ) + D u ( t )
where:
  • q = [ q 1 q 2 q 3 q 4 ] T is called the state vector and u = [ T i n T e x t ] T is called the input vector.
    A = [ 0 1 0 0 N K e q I m B m I m N K e q I m 0 0 0 0 1 N K e q I L 0 N K e q I L B L I L ] ,   B = [ 0 0 1 I m 0 0 0 0 1 I L ]
  • C: the output matrix, which is an identity matrix.
  • D: the feedforward matrix, which is a zero matrix.
Consequently, a type-one servo system with an integral action is designed (see Figure 3).
Thus, another state is added to the system to eliminate the steady state error. However, the following equations are used to construct the control system of the BcVSA. More details on this method can be found in [12].
e ˙ = A ^ e + B ^ u e
where:
  • e ( t ) = [ q e ( t ) ξ e ( t ) ] = [ q ( t ) q ( ) ξ ( t ) ξ ( ) ] is the dynamic error on the new state space vector.
  • u e ( t ) = u ( t ) - u ( ) is the dynamic error on the input vector.
    A ^ = [ A 0 C 2 0 ] ,   B ^ = [ B 0 ] ξ ˙ ( t ) = r ( t ) C 2 q ( t )
  • C 2 = [ 1 0 0 0 0 0 1 0 ] is the output-tracking matrix.
  • r ( t ) = [ q 1 q 3 ] represents the desired responses.
    u e ( t ) =   K ^ e =   K q e ( t ) + K I ξ e ( t )
where K ^ is the gain matrix calculated in MATLAB by using the following command:
K ^ = l q r ( A ^ , B ^ , Q ^ , R ^ )
where:
  • K ^ = [ K     K I ] ;
  • Q ^ is a semi-positive definite weighting matrix or positive definite weighting matrix;
  • R ^ is a positive definite weighting matrix.
The main features of this tracking system include its ability to track the reference commands and to regulate the other states. This controller, in reality, is classified as an optimal controller because it is minimizing a quadratic cost function to calculate the optimal response, which is achieved by solving the reduced form for the Riccati equation [12,13].

4.2. Tracking System Based on Computed Torque Control (CTC)

In general, this method is based on cancelling the nonlinear terms from the dynamic equations. However, the model of the BcVSA does not include any nonlinear term. We write the dynamic equation in the following form:
M θ ¨ +   B θ ˙ +   K θ ˙   =   T
where:
  • M = [ I m 0 0 I L ] is the mass matrix.;
  • B = [ B m 0 0 B L ] is the damping matrix;
  • K = N K e q [ 1 1 1 1 ] is the stiffness matrix;
  • T = [ T i n T e x t ] T is the torque input vector;
  • θ = [ θ m θ L ] T is the angular position for the motor and for the load, respectively.
Based on that, the CTC provides the system with an input vector ( T ) to cancel the damping and the stiffness matrices and to convert the model to a unit mass model but with ‘n’ DOF (see Equation (26)). However, this is achieved by identifying the mass, the damping, and the stiffness matrices, respectively.
T =   M ^ T ^ + B ^ θ ˙ + K ^ θ ˙
where M ^ ,   B ^ , and K ^ are the estimated mass, the estimated damping, and the estimated stiffness matrices. In a case of have a perfect estimation process, it is expected that the controller will work perfectly. Otherwise, the performance of the controller will be affected based on the error in the estimation process.
I . θ ¨ = { T ^ , w i t h   p e r f e c t   e s t i m a t e s T ^ + ε ( t ) , w i t h   g o o d   e s t i m a t e s
T ^ = θ ¨ d k d ( θ ˙ θ ˙ d ) k p ( θ ˙ θ d )
where:
  • ε ( t ) is the estimation error;
  • k d is the gain matrix of derivative term;
  • k p is the proportional gain matrix.
Based on that,
e ¨ + k d e ˙ + k p e = { 0 ,       w i t h   p e r f e c t   e s t i m a t e s ε ( t ) ,     w i t h   g o o d   e s t i m a t e s
where e is the error between the actual and the estimated values for the angular positions of the motor and the load, respectively.
In general, it is recommended that high gains for k d and k p are used to reject the errors.

5. Numerical Simulations

In this section, we present the numerical simulations of our model based on the system identification and control methods detailed in previous sections. The proposed results are intended to be used on the hardware platform of our actuator in the near future. This section is divided into three major subsections. In particular, the first subsection focuses on verifying the proposed method for identifying the unknown dynamic parameters, while the second and the third subsections evaluate the performance of both control techniques, i.e., LQR and CTC, respectively. Both control techniques are tested under different possible conditions, i.e., while altering the status of each clutch in different combinations.

5.1. Evaluation of the System Identifications Algorithm

In order to verify the previously introduced parameters identification method, we built a Simulink scheme based on the dynamical equations of the BcVSA to identify the same parameters. The main idea is to assume that the parameter values are the same as previously identified and re-determine the parameter values under different conditions. Hence, the proposed algorithm will be considered correct if the estimated values confirm the assumed values or vice versa. Accordingly, in the first test, θ L is clamped while θ m is free, so the goal is to identify I m and B m . This test is repeated under various conditions, e.g., by choosing different values of K e q through the selection of different combinations of active and inactive clutches. Moreover, we also simulated this test under the effect of various values of ICs (less than 0.3 rad). The resulting parameters are detailed in Table 2.
In the second test, θ L is free while θ m is clamped to identify I L and B L . Again, the test is performed under the effect of various combinations of active and inactive clutches. The resultant relative errors for the estimated parameters are shown in Table 3.
As a result, the logarithmic decrement algorithm method showed good performance in estimating the assumed parameters without any significant relative error, as detailed in Table 2 and Table 3.

5.2. Evaluation of the Servo System Design Based on a Linear Quadratic Regulator (LQR)

A cubic trajectory planning algorithm was used to create the desired trajectories for θ m and θ L , respectively. The design parameters in this scheme are Q and R matrices. Nevertheless, suitable parameters were selected to perform the required task. Under these conditions, we tested our system by selecting different stiffness levels through the status (engaged/disengaged) of the clutches.

5.2.1. All Clutches Are Active ( K e q = 7 K 0 ) the Entire Time

In the first case, the status of all the clutches is active, i.e., K e q = 7 K 0 . The resultant simulated response and the input signals are shown in Figure 4 and Figure 5, respectively. As is clear from the plots, the LQR controller was able to track the desired reference trajectory.

5.2.2. Altering the Clutches’ Status While the System Is Running at Different Time Intervals

In this case, we used the same design parameters as chosen in the previous case, but this time we tested the controller ability to control the actuator by altering the stiffness level while the system is running. The cubic trajectory-planning algorithm was used for three seconds. Different stiffness levels ( K e q ) were selected at different sub-time intervals as follows:
  • In the first-time interval [0 ≤ t < 1] second, Bit 21 is on i.e., ( K e q = 2 K 0 ).
  • Bit 20 and Bit 22 are on i.e., ( K e q = 5 K 0 ) in the next time interval [1 ≤ t < 2] second.
  • In the third time interval [2 ≤ t ≤ 3] second, all clutches are activated ( K e q = 7 K 0 ).
The simulation result for this case are shown in Figure 6. The corresponding input signal is shown in Figure 7. As a result, the controller successfully tracked the reference trajectory.

5.3. Evaluating the Tracking System Based on Computed Torque Control (CTC)

In order to evaluate the performance of the CTC, we repeated the simulations with CTC while using the same status of the clutches as previously chosen for LQR.

5.3.1. All Clutches Are Active ( K e q = 7 K 0 ) the Entire Time

In this test, all clutches are active i.e., K e q = 7 K 0 . The simulation results are reported in Figure 8, while the corresponding input signal is shown in Figure 9. The CTC scheme also showed its ability to track the desired reference trajectory with little error. The major source of these errors are the saturation blocks on the control scheme, because the control signals must not exceed 35 Nm. when θ = 10 ° (which is fixed according to actuator design criteria). However, the controller has the capability of tracking the desired reference signals without any errors in principle as long as the required input signals are lower than the saturated levels.

5.3.2. Altering the Clutches Status at Different Time Intervals While the System Is Running

Similarly, in this case, the same procedure of applying different stiffness levels was achieved by altering the status of clutches at different sub-time intervals. Details are as follows:
  • In the first time interval [0 ≤ t < 1] second, the Bit 21 clutch is on.
  • Bit 20 and Bit 22 are activated, i.e., ( K e q = 5 K 0 ) in the second interval [1 ≤t < 2] second.
  • In the third interval [2 ≤ t ≤ 3] second, all clutches are on ( K e q = 7 K 0 ).
The results of this case are presented in Figure 10 while the applied input signals are shown in Figure 11.
In summary, both control schemes showed good performance in controlling our BcVSA model in a simulation environment. It was observed that the CTC method has comparatively fewer errors but is comparatively less efficient in terms of energy consumption. In addition, the CTC mainly depends on our knowledge about the model itself and on reducing any source of uncertainty on the dynamic parameters. However, the preceding problem can be partially solved by increasing the value of the gain matrices k d and k p . On the other hand, the LQR method showed a good capability to track the desired signals with relatively less energy consumption, in spite of the fact that the presented results were achieved by selecting higher values for the entries of the Q matrix and smaller values for the entries of the R matrix. These choices mean that the priority is focused on decreasing the error between the desired and the simulated responses instead of focusing on minimizing the energy consumption. Overall, it was observed that the CTC method is more reliable for our system, as the error between the desired and the simulated responses are less than the LQR (see Figure 6 and Figure 8). Table 4 shows the quantitative values of power and energy consumption for each case using both types of controllers.
We aim to study the robustness of both control approaches in the simulation and the prototype of our actuator in the future.

6. Conclusions

In this work, the design, modeling, system identification, and control of discrete variable stiffness actuators are presented. We used a logarithmic decrement algorithm for the system identification. In conclusion, the implemented approach for identifying the dynamic parameters is working perfectly as shown in the numerical results, but it is expected that, during the practical phase, the efficiencies of all the gears will cause an increase in errors, i.e., the backlashes may cause this problem. Nevertheless, these errors can be dramatically reduced by taking into account the efficiency of the gears into the dynamic equations. Additionally, this paper presented two control methods based on LQR and CTC along with their simulated results. Both methods were simulated under the effect of various combinations of active (engaged) and inactive (disengaged) clutches. Based on relatively fewer errors and a shorter settling time as compared to the LQR, results of the CTC are more reliable for our model. On the contrary, the LQR method was proven to be more energy-efficient as compared to the CTC, as shown in the simulated results. In future work, we plan to test the presented approaches on our hardware prototype and compare the results with the simulated ones.

Author Contributions

Conceptualization, M.I.A. and I.H.; methodology, A.A., I.H., M.I.A.; software I.H., A.A.; validation, I.H., A.A. and D.G.; writing—original draft preparation, I.H., M.I.A., A.A.; writing—review and editing, I.H., M.I.A. and A.A.; supervision, D.G.; project administration, D.G.; funding acquisition, D.G.

Funding

The research leading to these results was funded by the ADEK Award for Research Excellence (AARE) 2017, and Partially Supported Khalifa University of Science and Technology under Award No. CIRA-2018-55 and RC1-2018-KUCARS

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

References

  1. Tsagarakis, N.G.; Sardellitti, I.; Caldwell, D.G. A new variable stiffness actuator (CompAct-VSA): Design and modelling. In Proceedings of the 2011 IEEE/RSJ International Conference on Intelligent Robots and Systems, San Francisco, CA, USA, 25–30 September 2011; pp. 378–383. [Google Scholar]
  2. Bicchi, A.; Peshkin, M.; Colgate, J.E. Safety for Physical Human–Robot Interaction (Handbook of Robotics 2008); Springer: Berlin, Germany, 2008; pp. 1335–1348. [Google Scholar]
  3. Lasota, P.A.; Fong, T.; Shah, J.A. A survey of methods for safe human–robot interaction. Found. Trends Robot. 2017, 5, 261–349. [Google Scholar] [CrossRef]
  4. Hogan, N. Impedance control: An approach to manipulation: Part III—Applications. J. Dyn. Syst. Meas. Control 1985, 107, 8–16. [Google Scholar] [CrossRef]
  5. Shin, D.; Sardellitti, I.; Park, Y.-L.; Khatib, O.; Cutkosky, M. Design and Control of a Bio-inspired Human-friendly Robot. Int. J. Robot. Res. 2010, 9, 571–584. [Google Scholar] [CrossRef]
  6. Migliore, S.A.; Brown, E.A.; DeWeerth, S.P. Biologically inspired joint stiffness control. In Proceedings of the IEEE International Conference on Robotics and Automation, Barcelona, Spain, 18–22 April 2005; pp. 4508–4513. [Google Scholar]
  7. Bischoff, R.; Kurth, J.; Schreiber, G.; Koeppe, R.; Albu-Schaeffer, A.; Beyer, A.; Eiberger, O.; Haddadin, S.; Stemmer, A.; Grunwald, G.; et al. The KUKA-DLR lightweight robot arm—A new reference platform for robotics research and manufacturing. In Proceedings of the 41st International Symposium on Robotics (ISR 2010) and 6th German Conference on Robotics (ROBOTIK 2010), Munich, Germany, 7–9 June 2010; pp. 1–8. [Google Scholar]
  8. Bianchi, M.; Valenza, G.; Greco, A.; Nardelli, M.; Battaglia, E.; Bicchi, A.; Scilingo, E.P. Towards a novel generation of haptic and robotic interfaces: Integrating affective physiology in human–robot interaction. In Proceedings of the 2016 25th IEEE International Symposium on Robot and Human Interactive Communication (RO-MAN), New York, NY, USA, 26–31 August 2016; pp. 125–131. [Google Scholar]
  9. Albu-Schaffer, A.; Haddadin, S.; Ott, C. The DLR lightweight robot: Design and control concepts for robots in human environments. Ind. Robot Int. J. 2007, 34, 376–385. [Google Scholar] [CrossRef]
  10. Haddadin, S.; Schaffer, A.A.; Hirzinger, G. Safety evaluation of physical human–robot interaction via crash testing. Robot. Sci. Syst. 2007, 3, 217–224. [Google Scholar]
  11. Pratt, G.; Williamson, M. Series Elastic Actuators. In Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems, Pittsburgh, PA, USA, 5–9 August 1995; Volume 1, pp. 399–406. [Google Scholar]
  12. Ogata, K. Modern Control Engineering, 5th ed.; Prentice Hall: Upper Saddle River, NJ, USA, 2009; pp. 793–806. [Google Scholar]
  13. Kirk, D. Optimal Control Theory: An Introduction; Courier Corporation: New York, NY, USA, 2012; pp. 209–227. [Google Scholar]
  14. Piltan, F.; Yarmahmoudi, M.; Shamsodini, M.; Mazlomian, E.; Hosainpour, A. PUMA-560 robot manipulator position computed torque control methods using Matlab/Simulink and their integration into graduate nonlinear control and Matlab courses. Int. J. Robot. Autom. 2012, 3, 167–191. [Google Scholar]
  15. Tonietti, G.; Schiavi, R.; Bicchi, A. Design and Control of a Variable Stiffness Actuator for Safe and Fast Physical Human/Robot Interaction. In Proceedings of the 2005 IEEE International Conference on Robotics and Automation, Barcelona, Spain, 18–22 April 2005; pp. 526–531. [Google Scholar]
  16. Schiavi, R.; Grioli, G.; Sen, S.; Bicchi, A. 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 2008, Pasadena, CA, USA, 19–23 May 2008; pp. 2171–2176. [Google Scholar]
  17. Hurst, J.; Chestnutt, J.; Rizzi, A. The Actuator with Mechanically Adjustable Series Compliance. IEEE Trans. Robot. 2010, 26, 597–606. [Google Scholar] [CrossRef]
  18. Jafari, A.; Tsagarakis, N.; Vanderborght, B.; Caldwell, D. A novel actuator with adjustable stiffness (AwAS). In Proceedings of the 2010 IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS), Taipei, Taiwan, 18–22 October 2010; pp. 4201–4206. [Google Scholar]
  19. Jafari, A.; Tsagarakis, N.; Caldwell, D. AwAS-II: A New Actuator with Adjustable Stiffness based on the Novel Principle of Adaptable Pivot point and Variable Lever ratio. In Proceedings of the 2011 IEEE International Conference on Robotics and Automation, Shanghai, China, 9–13 May 2011; pp. 4638–4643. [Google Scholar]
  20. Visser, L.C.; Carloni, R.; Stramigioli, S. Energy efficient variable stiffness actuators. IEEE Trans. Robot. 2011, 27, 865–875. [Google Scholar] [CrossRef]
  21. Fumagalli, M.; Barrett, E.; Stramigioli, S.; Carloni, R. The mVSA-UT: A miniaturized differential mechanism for a continuous rotational variable stiffness actuator. In Proceedings of the 2012 4th IEEE RAS&EMBS International Conference on Biomedical Robotics and Biomechatronics (BioRob), Rome, Italy, 24–27 June 2012; pp. 1943–1948. [Google Scholar]
  22. Groothuis, S.S.; Rusticelli, G.; Zucchelli, A.; Stramigioli, S.; Carloni, R. The Variable Stiffness Actuator vsaUT-II: Mechanical Design, Modeling, and Identification. IEEE/ASME Trans. Mechatron. 2014, 19, 589–597. [Google Scholar] [CrossRef]
  23. Awad, M.I.; Gan, D.; Cempini, M.; Cortese, M.; Vitiello, N.; Dias, J.; Dario, P.; Seneviratne, L. Modeling, Design and Characterization of A Novel Passive Variable Stiffness Joint (pVSJ). In Proceedings of the IEEE International Conference of Intelligent Robotics & Systems (IROS), Daejeon, Korea, 9–14 October 2016; pp. 323–329. [Google Scholar]
  24. Awad, M.I.; Gan, D.; Thattamparambil, J.; Az-zu’bi, A.; Stefanini, C.; Dias, J.; Seneviratne, L. Novel Passive Discrete Variable Stiffness Joint (pDVSJ): Modeling, Design, and Characterization. In Proceedings of the IEEE International Conference of Robotics and Biomimetics (ROBIO), Qingdao, China, 3–7 December 2016; pp. 1808–1813. [Google Scholar]
  25. Mathijssen, G.; Lefeber, D.; Vanderborght, B. Variable Recruitment of Parallel Elastic Elements: Series–Parallel Elastic Actuators (SPEA) With Dephased Mutilated Gears. IEEE/ASME Trans. Mechatron. 2015, 20, 594–602. [Google Scholar] [CrossRef]
  26. Mathijssen, G.; Furnémont, R.; Beckers, S.; Verstraten, T.; Lefeber, D.; Vanderborght, B. Cylindrical cam mechanism for unlimited subsequent spring recruitment in Series-Parallel Elastic Actuators. In Proceedings of the IEEE International Conference on Robotics and Automation (ICRA), Seattle, WA, USA, 26–30 May 2015; pp. 857–862. [Google Scholar]
  27. Mathijssen, G.; Furnémont, R.; Verstraten, T.; Brackx, B.; Premec, J.; Jiménez, R.; Lefeber, D.; Vanderborght, B. +SPEA introduction: Drastic actuator energy requirement reduction by symbiosis of parallel motors, springs and locking mechanisms. In Proceedings of the IEEE International Conference on Robotics and Automation (ICRA), Stockholm, Sweden, 16–21 May 2016; pp. 676–681. [Google Scholar]
  28. Awad, M.I.; Gan, D.; Hussain, I.; Az-zu’bi, A.; Stefanini, C.; Khalaf, K.; Zweiri, Y.; Taha, T.; Dias, J.; Seneviratne, L. Design of A Novel Passive Binary-Controlled Variable Stiffness Joint (BpVSJ) Towards Passive Haptic Interface Application. IEEE Access 2018, 6, 63045–63057. [Google Scholar] [CrossRef]
  29. Awad, M.I.; Hussain, I.; Gan, D.; Az-zu’bi, A.; Stefanini, C.; Khalaf, K.; Zweiri, Y.; Taha, T.; Dias, J.; Seneviratne, L. Passive Discrete Variable Stiffness Joint (pDVSJ-II): Modeling, Design, Characterization and Testing Towards Passive Haptic Interface. ASME. J. Mech. Robot. 2018, 11, 011005. [Google Scholar] [CrossRef]
  30. Hussain, I.; Albalasie, A.; Awad, M.I.; Gan, D.; Seneviratne, L. Modeling, Multiple Model Predictive Control and Numerical Simulations of a Novel Binary-Controlled Variable Stiffness Actuator (BcVSA). Front. Robot. AI 2018, 5, 68. [Google Scholar] [CrossRef]
  31. Rao, S. Mechanical Vibrations, 5th ed.; Prentice Hall: Upper Saddle River, NJ, USA, 2011; pp. 164–166. [Google Scholar]
  32. Magalas, L.; Malinowski, T. Measurement techniques of the logarithmic decrement. Solid State Phenom. 2003, 89, 247–260. [Google Scholar] [CrossRef]
  33. Aribowo, W.; Yamashita, T.; Terashima, K.; Masui, Y.; Saeki, T.; Kamigaki, T.; Kawamura, H. Vibration control of semiconductor wafer transfer robot by building an integrated tool of parameter identification and input shaping. IFAC Proc. 2011, 44, 14367–14373. [Google Scholar] [CrossRef]
  34. Tsetserukou, D.; Kawakami, N.; Tachi, S. Vibration damping control of robot arm intended for service application in human environment. In Proceedings of the 8th IEEE-RAS International Conference on Humanoid Robots, Daejeon, Korea, 1–3 December 2008; pp. 441–446. [Google Scholar]
  35. Park, J.; Kim, K. Biped robot walking using gravity-compensated inverted pendulum mode and computed torque control. In Proceedings of the 1998 IEEE International Conference on Robotics & Automation, Leuven, Belgium, 20–20 May 1998; Volume 4, pp. 3528–3533. [Google Scholar]
  36. Lammerts, I. Adaptive Computed Reference Computed Torque Control of Flexible Manipulators; Eindhoven University of Technology: Eindhoven, The Netherlands, 1993. [Google Scholar]
  37. Spong, M. Partial feedback linearization of underactuated mechanical systems. In Proceedings of the IEEE/RSJ International Conference on Intelligent Robots and Systems (IROS’94), Munich, Germany, 12–16 September 1994; Volume 1. [Google Scholar]
Figure 1. The passive version of the discrete variable stiffness actuator (DVSA).
Figure 1. The passive version of the discrete variable stiffness actuator (DVSA).
Actuators 08 00050 g001
Figure 2. (a) The CAD representation of discrete variable stiffness actuator. (b) The schematic representation of the working principle of the actuator.
Figure 2. (a) The CAD representation of discrete variable stiffness actuator. (b) The schematic representation of the working principle of the actuator.
Actuators 08 00050 g002
Figure 3. Type-one servo system for controlling the binary controlled variable stiffness actuator (BcVSA).
Figure 3. Type-one servo system for controlling the binary controlled variable stiffness actuator (BcVSA).
Actuators 08 00050 g003
Figure 4. The responses of the BcVSA based on the linear quadratic control (LQR) scheme under the condition when all the clutches are active ( K e q = 7 K 0 ) during the entire simulation time. (Top) the input trajectory response; (Bottom) the output response.
Figure 4. The responses of the BcVSA based on the linear quadratic control (LQR) scheme under the condition when all the clutches are active ( K e q = 7 K 0 ) during the entire simulation time. (Top) the input trajectory response; (Bottom) the output response.
Actuators 08 00050 g004
Figure 5. The input signals of BcVSA based on the LQR scheme under the condition when all the clutches are active ( K e q = 7 K 0 ) during the entire simulation time.
Figure 5. The input signals of BcVSA based on the LQR scheme under the condition when all the clutches are active ( K e q = 7 K 0 ) during the entire simulation time.
Actuators 08 00050 g005
Figure 6. The responses of DVSA based on the LQR scheme in the scenario of altering the status of clutches while the system is running. (Top) the input trajectory response; (Bottom) the output response.
Figure 6. The responses of DVSA based on the LQR scheme in the scenario of altering the status of clutches while the system is running. (Top) the input trajectory response; (Bottom) the output response.
Actuators 08 00050 g006
Figure 7. The input signals for the LQR scheme while altering the clutches’ status at different time intervals while the system is running.
Figure 7. The input signals for the LQR scheme while altering the clutches’ status at different time intervals while the system is running.
Actuators 08 00050 g007
Figure 8. The responses of BcVSA by using the computed torque control (CTC) scheme under the condition when all the clutches are active ( K e q = 7 K 0 ) for the entire simulation time. (Top) the input trajectory response; (Bottom) the output response.
Figure 8. The responses of BcVSA by using the computed torque control (CTC) scheme under the condition when all the clutches are active ( K e q = 7 K 0 ) for the entire simulation time. (Top) the input trajectory response; (Bottom) the output response.
Actuators 08 00050 g008
Figure 9. The input signals for the CTC scheme when all the clutches are active ( K e q = 7 K 0 ) for the entire simulation time.
Figure 9. The input signals for the CTC scheme when all the clutches are active ( K e q = 7 K 0 ) for the entire simulation time.
Actuators 08 00050 g009
Figure 10. The responses of the BcVSA based on the CTC scheme in the scenario of altering the status of clutches while the system is running. (Top) the input trajectory response; (Bottom) the output response.
Figure 10. The responses of the BcVSA based on the CTC scheme in the scenario of altering the status of clutches while the system is running. (Top) the input trajectory response; (Bottom) the output response.
Actuators 08 00050 g010
Figure 11. The input signals for the CTC scheme during altering status of the clutches at different time intervals while the system is running.
Figure 11. The input signals for the CTC scheme during altering status of the clutches at different time intervals while the system is running.
Actuators 08 00050 g011
Table 1. The total possible stiffness levels of the DVSA.
Table 1. The total possible stiffness levels of the DVSA.
Stiffness Bit 22 (4K0)Stiffness Bit 21 (2K0)Stiffness Bit 20 (K0)Stiffness Level
0000
001K0
0102 K0
0113 K0
1004 K0
1015 K0
1106 K0
Table 2. Relative errors for the estimated parameters under the effect of active and inactive clutches in the case of the first test.
Table 2. Relative errors for the estimated parameters under the effect of active and inactive clutches in the case of the first test.
No.Stiffness Bit 22 (4K0)Stiffness Bit 21 (2K0)Stiffness Bit 20 (K0) I m _ e s t   ( % ) B m _ e s t   ( % )
10000.39380.1075
20010.35840.1160
30100.19860.0445
40110.15620.0451
51000.24440.1348
61010.02780.0235
71100.01110.0035
Table 3. Relative errors for the estimated parameters under the effect of active and inactive clutches in the case of the second test.
Table 3. Relative errors for the estimated parameters under the effect of active and inactive clutches in the case of the second test.
No.Stiffness Bit 22 (4K0)Stiffness Bit 21 (2K0)Stiffness Bit 20 (K0) I L _ e s t   ( % ) B L _ e s t   ( % )
10000.00980.0853
20010.55820.2997
30100.12690.1097
40110.28210.1274
51000.46790.2155
61010.30620.1494
71100.21160.1056
Table 4. The quantitative values of power and energy consumption for each case using both types of controllers.
Table 4. The quantitative values of power and energy consumption for each case using both types of controllers.
Clutches Status The Mean of the Consumed Power [W]The Energy Consumption [J]
θ m θ L 2DOFs θ m θ L 2DOF
5.2.1 All Clutches are Active ( K e q = 7 K 0 ) the Entire Time (LQR)0.960.631.592.901.904.8
5.2.2 Altering the Clutches’ Status while the System is Running At Different Time Intervals (LQR)0.820.521.342.471.554.02
5.3.1 All Clutches are Active ( K e q = 7 K 0 ) the Entire Time (CTC)0.770.961.732.322.885.2
5.3.2 Altering the Clutches Status at Different Time Intervals while The System is Running (CTC)0.630.801.431.882.404.28

Share and Cite

MDPI and ACS Style

Hussain, I.; Albalasie, A.; Awad, M.I.; Gan, D. Modeling, Identification, and Control of a Discrete Variable Stiffness Actuator (DVSA). Actuators 2019, 8, 50. https://doi.org/10.3390/act8030050

AMA Style

Hussain I, Albalasie A, Awad MI, Gan D. Modeling, Identification, and Control of a Discrete Variable Stiffness Actuator (DVSA). Actuators. 2019; 8(3):50. https://doi.org/10.3390/act8030050

Chicago/Turabian Style

Hussain, Irfan, Ahmad Albalasie, Mohammad I. Awad, and Dongming Gan. 2019. "Modeling, Identification, and Control of a Discrete Variable Stiffness Actuator (DVSA)" Actuators 8, no. 3: 50. https://doi.org/10.3390/act8030050

APA Style

Hussain, I., Albalasie, A., Awad, M. I., & Gan, D. (2019). Modeling, Identification, and Control of a Discrete Variable Stiffness Actuator (DVSA). Actuators, 8(3), 50. https://doi.org/10.3390/act8030050

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop