Next Article in Journal
Fundamentals and Applications of Hybrid LWFA-PWFA
Next Article in Special Issue
Adaptive Fuzzy Backstepping Sliding Mode Control for a 3-DOF Hydraulic Manipulator with Nonlinear Disturbance Observer for Large Payload Variation
Previous Article in Journal
External Unilateral Fixator of Own Design for the Treatment of Selected Mandibular Fractures in Horses
Previous Article in Special Issue
Adaptive Super-Twisting Control for Mobile Wheeled Inverted Pendulum Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling and Discrete-Time Terminal Sliding Mode Control of a DEAP Actuator with Rate-Dependent Hysteresis Nonlinearity

School of Automation, Beijing Institute of Technology, Beijing 100081, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(13), 2625; https://doi.org/10.3390/app9132625
Submission received: 28 April 2019 / Revised: 11 June 2019 / Accepted: 25 June 2019 / Published: 28 June 2019
(This article belongs to the Special Issue The Application of Sliding Mode Control in Robots)

Abstract

:
Dielectric electro-active polymer (DEAP) materials, also called artificial muscle, are a kind of EAP smart materials with extraordinary strains up to 30% at a high driving voltage. However, the asymmetric rate-dependent hysteresis is a barrier for trajectory tracking control of DEAP actuators. To overcome the barrier, in this paper, a Hammerstein model is established for the asymmetric rate-dependent hysteresis of a DEAP actuator first, in which a modified Prandtl-Ishlinskii (MPI) model is used to represent the static hysteresis nonlinear part, and an autoregressive with exogenous inputs (ARX) model is used to represent the linear dynamic part. Applying Levenberg-Marquardt (LM) algorithm identifies the parameters of the Hammerstein model. Then, based on the MPI model, an inverse hysteresis compensator is obtained to compensate the hysteresis behavior. Finally, a compound controller consisting of the hysteresis compensator and a novel discrete-time terminal sliding mode controller (DTSMC) without state observer is proposed to achieve the high-precision trajectory tracking control. Stability analysis of the closed-loop system is verified by using Lyapunov stability theorem. Experimental results based on a DEAP actuator show that the proposed controller has better tracking control performance compared with a conventional discrete-time sliding mode controller (DSMC).

1. Introduction

Dielectric electro-active polymer (DEAP) materials are a promising actuation technology due to the merits of lightweight, large strain, high energy density, low-noise, and fast response [1,2]. DEAP materials, as a kind of soft electroactive material, are composed of a polymer film sandwiched between two silver metal soft electrodes layers [3]. The working principle of DEAP materials is shown in Figure 1. It can be simply divided into three steps [3]: (1) a large driving voltage is applied to the silver metal soft electrodes layers to produce a high electric field; (2) the silver metal soft electrodes layers can be regarded as a flat capacitor, a Maxwell stress will be produced by electrostatic forces; (3) the polymer film is extruded under the Maxwell stress to result in actuation. Based on this basic working principle, DEAP materials have been widely applied in high-precision trajectory tracking control as actuators and achieved great potential for applications in bioinspired soft robots filed, including an arm wrestling robot [4], a flying robot [5], and a swimming robot [6]. However, like other smart materials such as shape memory alloys (SMA), piezoelectrics and magnetostrictives, the inherent asymmetric rate-dependent hysteresis nonlinearity is ubiquitous in DEAP actuators [7], which causes tracking inaccuracy, leads to a decrease in the performance of the closed-loop system, or may even result in oscillations. Therefore, the asymmetric rate-dependent hysteresis modeling and tracking control of such DEAP actuators have been important issues.
The asymmetric rate-dependent hysteresis models of the smart materials can be roughly classified into three categories. The first category is physics-based hysteresis models. For instance, Ref. [8] presented a dynamic Jiles–Atherton model to describe the asymmetric rate-dependent hysteresis of giant magnetostrictive actuators (GMAs). The second category is phenomenological hysteresis models. For example, Ref. [9] proposed an MPI model to characterize the asymmetric rate-dependent hysteresis of dielectric elastomer actuators (DEAs), where dynamical envelope functions were used to describe the rate-dependent hysteresis, and polynomial operators were introduced to characterize the asymmetric . Ref. [10] proposed an MPI model to describe the asymmetric rate-dependent hysteresis of magnetostrictive actuators, in which a rate-dependent PI model was utilized to describe the rate-dependent hysteresis, and dead-zone operators were introduced to characterize the asymmetric. The third category is the intelligent methods-based hysteresis models. Ref. [11] developed a nonlinear ARX recurrent neural network hysteresis model to characterize the asymmetric rate-dependent hysteresis of SMA actuators. Refs. [7,12] developed a Preisach-type dynamic nonlinear ARX fuzzy hysteresis model to describe the asymmetric rate-dependent hysteresis of DEAP actuators. It should be pointed out that the unmodeled hysteresis nonlinearity, residual dynamics, parameter uncertainties and external disturbances will reduce the trajectory tracking accuracy. Therefore, it is necessary to design a robust controller to deal with the nonlinear disturbances.
Various control strategies have been developed for trajectory tracking control of DEAP actuators such as PID control [7,13], internal model control [14] and adaptive sliding mode control (SMC) [15]. Among those control strategies, SMC scheme has been widely used to control the hysteresis nonlinearity systems, because SMC scheme can guarantee the global robustness of the systems, and can simply and effectively cope with parameters variation and external disturbances [16]. However, the conventional SMC scheme can only make the system states converge to the equilibrium point asymptotically with infinite settling time. To overcome the shortcoming, terminal sliding mode controller (TSMC) scheme was proposed [17], in which a nonlinear switching surface was employed to achieve finite time convergence [18]. In the literature, Ref. [19] developed a TSMC scheme for precision tracking control and finite time convergence of a piezoelectric-driven micropositioning system. Ref. [20] proposed an adaptive-gain fast nonsingular TSMC scheme for a piezo positing stage system which can quickly converge in finite time and reduce chattering. The aforementioned SMC and TSMC schemes are all applied in continuous-time domain; however, with the wide use of computers in control systems, the DSMC theory [21] has attracted considerable attention in recent years [22,23,24]. It must be pointed out that the continuous-time SMC cannot directly convert into its discrete counterpart through simple equivalence. DSMC and DTSMC schemes have been successfully applied to other smart actuators. For instance, Ref. [25] proposed a novel DSMC scheme for precision tracking control of a piezoelectric actuator. In addition, for making the system state converge to the equilibrium point asymptotically within finite steps, Refs. [26,27] proposed DTSMC schemes for precision tracking control of piezoelectric actuators, and the stability of the control systems and the effectiveness of the proposed control schemes were demonstrated through experimental studies. However, to the knowledge of the author, DTSMC scheme has not been applied for precision trajectory tracking of DEAP actuators in the literature.
This paper focuses on the asymmetric rate-dependent hysteresis modeling and trajectory tracking control of a DEAP actuator under different input frequencies. A Hammerstein model with hysteresis nonlinearity is established to describe the asymmetric rate-dependent hysteresis of the DEAP actuator. To estimate the parameters of the model, LM algorithm [28,29] is utilized. Based on the MPI model, an inverse hysteresis compensator is obtained to suppress the asymmetric hysteresis. In order to achieve precision tracking control and robustness of the closed-loop system, a DTSMC with a new discrete-time terminal sliding function is proposed. The main contributions of this paper are summarized as follows:
(1) A Hammerstein model with hysteresis nonlinearity is established to account the asymmetric rate-dependent hysteresis of the DEAP actuator and the lumped nonlinear disturbance term is well formulated. It provides with a detailed model for the closed-loop system.
(2) A DTSMC scheme based on a new discrete-time terminal sliding function is first proposed for the DEAP actuator. The proposed control scheme not only improves the trajectory tracking control accuracy and error convergence but also suppresses the chattering phenomenon. In addition, the proposed control scheme does not require a state observer.
(3) Compared with other DTSMC schemes, a terminal-sliding-mode-type switching control law is defined in the proposed DTSMC scheme, which improves the robustness feature of the closed-loop system in the presence of unmodeled hysteresis nonlinearity, residual dynamics, parameter uncertainties and external disturbances.
(4) Finite steps convergence and stability of the closed-loop system are theoretically verified and an experimental investigation is conducted to validate the effectiveness of the proposed control scheme to the DEAP actuator. Then, a comparison between the proposed control scheme and the DSMC scheme is undertaken.
The remainder of the paper is organized as follows. In Section 2, the analytical formulation of the Hammerstein model and the inverse MPI model are established first, and then the parameters of the Hammerstein model are identified by LM algorithm. In Section 3, a compound controller is proposed, and the system stability is proved. Experimental results of the proposed controller to the DEAP actuator are presented in Section 4. Conclusions are given in Section 5.

2. Hammerstein Model for DEAP Actuator

In this section, in order to characterize asymmetric hysteresis loops of the DEAP actuator, a Hammerstein model is established, as shown in Figure 2. It consists of two parts: nonlinear static part and linear dynamic part. The nonlinear static part is represented by an MPI model to describe the asymmetric hysteresis, and the linear dynamic part is represented by an ARX model to describe the rate-dependent hysteresis. In the Figure 2, v ( k ) and y ( k ) represent the input and output of the DEAP actuator, respectively, u ( k ) is the corresponding unmeasurable internal variable, and d 0 ( k ) is the external disturbance.

2.1. MPI Model

In this subsection, an MPI model is established to characterize the asymmetric hysteresis of the DEAP actuator. The MPI model is formulated as the cascade of the PI model and the dead-zone (DZ) model [30], in which the PI model is a weighted superposition of play operators, and the DZ model is a weighted superposition of dead-zone operators. The play operator and the dead-zone operator are shown in Figure 3 and Figure 4, respectively.
The play operator is expressed as [31]
P r h i ( k ) = m a x m i n v ( k ) + r h i , P r h i ( k 1 ) , v ( k ) r h i
while the initial condition is defined as follows:
P r h i ( 0 ) = m a x m i n v ( 0 ) + r h i , P r h i 0 , v ( 0 ) r h i
The output of the PI model can be formulated as
W ( k ) = q v ( k ) + [ ϑ 1 , ϑ 2 , , ϑ n ] · P r h 1 ( k ) , P r h 2 ( k ) , , P r h n ( k ) T = q v ( k ) + ϑ T · H r h [ v ( k ) , P r h 0 ]
where q is a given positive constant; v and W represent the input and output of the PI model, respectively; ϑ = [ ϑ 1 , ϑ 2 , , ϑ n ] T is the vector of weights; r h = [ r h 1 , , r h n ] T is the vector of thresholds with 0 < r h 1 < < r h n < + ; H r h denotes the vector of play operators; P r h 0 = [ P r 10 , P r 20 , , P r n 0 ] is the vector of initial states; n is the number of play operators.
The formulation of the dead-zone operator is expressed as
S r d [ W ] ( k ) = m a x W ( k ) r d , 0 , r d > 0 W ( k ) , r d = 0 m i n W ( k ) r d , 0 , r d < 0
The output of the DZ model can be formulated as
Z ( k ) = [ θ 0 , θ 1 , θ 2 , , θ m ] · S r d 0 [ W ] ( k ) , S r d 1 [ W ] ( k ) , S r d 2 [ W ] ( k ) , , S r d m [ W ] ( k ) T = θ T · S r d [ W ( k ) ]
where W and Z denote the input and output of the DZ model, respectively; θ = [ θ 0 , θ 1 , θ 2 , , θ m ] T is the vector of weights; r d = [ r d 0 , r d 1 , r d 2 , , r d n ] T is the vector of thresholds with 0 = r d 0 < r d 1 < < r d m < + ; m is the number of dead-zone operators.
Substituting Equations (3) into (5), the mathematical formulation of the MPI model can be derived as
u ( k ) = [ θ 0 , θ 1 , θ 2 , , θ m ] · [ S r d 0 [ W ] ( k ) , S r d 1 [ W ] ( k ) , S r d 2 [ W ] ( k ) , , S r d m [ W ] ( k ) ] T = θ T · S r d q v ( k ) + ϑ T · H r h [ v ( k ) , P r h 0 ] ( k )
The form of the inverse MPI model is the same as the MPI model. From Equation (6), the inverse MPI model can be expressed as [31]
u 1 [ W ] ( k ) = q v ( k ) + ϑ T P r h θ T · S r d [ W ] , P r h 0 ( k )
where q is a positive constant; P r h is the vector of inverse play operators; P r h 0 is the intinal state; ϑ = [ ϑ 1 , ϑ 2 , , ϑ n ] T and r h = [ r h 1 , r h 2 , , r h n ] T denote the weighting vector and the threshold vector of inverse play operators, respectively; S r d is the vector of inverse dead-zone operators; θ = [ θ 0 , θ 1 , , θ m ] T and r d = [ r d 0 , r d 1 , , r d m ] T denote the weighting vector and the threshold vector of inverse dead-zone operators, respectively.
Here
q = 1 q , ϑ i = ϑ i ( q + j = 1 i ϑ j ) ( q + j = 1 i 1 ϑ j ) , i = 1 , 2 , , n r h i = q r h i + j = 1 i 1 ϑ j ( r h i r h j ) , i = 1 , 2 , , n P r i 0 = ( q + j = 1 i ϑ j P r i 0 ) + j = i + 1 n ϑ j P r j 0 , i = 1 , 2 , , n
and
θ 0 = 1 θ 0 , θ l = θ l ( θ 0 + k = 1 l θ k ) ( θ 0 + k = 1 l 1 θ k ) , l = 1 , 2 , , m r d l = k = 0 l θ k ( r h l r h k ) , l = 0 , 1 , 2 , , m
The parameters of the inverse MPI model can be obtained by using Equations (8) and (9).

2.2. ARX Model

In this subsection, a linear dynamic ARX model is established to characterize the rate-dependent hysteresis of the DEAP actuator, which can be defined as
A ( z ) y ( k ) = B ( z ) u ( k ) + d 0 ( k )
where A ( z ) : = 1 + a 1 z 1 + a 2 z 1 + + a n a z n a and B ( z ) : = b 1 z 1 + b 2 z 1 + + b n b z n b are polynomials in the unit backward shift operator z 1 [ z 1 y ( k ) = y ( k 1 ) ] with the known orders n a and n b ; y ( k ) is the system output sequence; u ( k ) is the input of the ARX model and also the output of the MPI model; d 0 ( k ) is the external disturbance.
Define the parameter vectors ζ a : = [ a 1 , a 2 , , a n a ] T n a , ζ b : = [ b 1 , b 2 , , b n b ] T n b and ζ : = [ ζ a , ζ b ] T n a + n b . Define the information vectors φ a : = [ y ( k 1 ) , y ( k 2 ) , , y ( k n a ) ] T n a , φ b : = [ u ( k 1 ) , u ( k 2 ) , , u ( k n b ) ] T n b and ϕ : = [ φ a , φ b ] T n a + n b .
Equation (10) can be turned into
y ( k ) = i = 1 n a a i y ( k i ) + i = 1 n b b i u ( k i ) + d 0 ( k ) = φ a T ζ a + φ b T ζ b + d 0 ( k ) = ϕ T ζ + d 0 ( k )

2.3. Parameters Identification

The goal of this subsection is to obtain the parameters values of the MPI model and the ARX model. Specific identification steps are as follows:
Step (1): The asymmetric hysteresis characterization of the DEAP actuator is modeled by using the MPI model. A low frequency input signal v of 0 . 1 Hz is given to the DEAP actuator, and N pairs of discrete input and output sequences ( v , y ) are measured. Then, based on the input data v and the output data y, LM algorithm is employed to identify the MPI model parameters.
Step (2): After the MPI model is determined, a hybrid frequency input signal in range of 0 . 1 to 0 . 5 Hz is given to the DEAP actuator, and the discrete input and output sequences ( v , y ) are measured. Based on the identified MPI model and v, the intermediate unmeasurable variable u can be obtained. Then, based on the intermediate variable u and the output data y, the parameters of the ARX model are identified by LM algorithm.
The accuracy of the identified model can be expressed by root-mean-square error (RMSE) and mean absolute error (MAE), which are defined as follows:
R M S E = i = 1 N y ( i ) y ^ ( i ) 2 / N
M A E = 1 N i = 1 N ( | y ( i ) y ^ ( i ) | )
where y ( i ) is the output of the DEAP actuator, y ^ ( i ) is the output of the Hammerstein model, and N is the amount of experimental data. The identified parameters of the PI model, DZ model and ARX model are given in Table 1, Table 2 and Table 3, respectively. Figure 5 shows the asymmetric rate-dependent hysteresis loops of the DEAP actuator in the range of 0 . 1 to 0 . 5 Hz. It illustrates that the width of the hysteresis loops will gradually increase with increasing frequency of the input signal [7,12].

2.4. Model Validation

To validate the Hammerstein model for the DEAP actuator, input signals at a frequency of 0 . 1 Hz, 0 . 2 Hz, 0 . 3 Hz, 0 . 4 Hz, 0 . 5 Hz, and a hybrid frequency in range of 0 . 1 to 0 . 5 Hz are applied to the DEAP actuator, respectively. Figure 6 shows the input-output relationship between the Hammerstein model and the DEAP actuator under different input frequencies. Table 4 shows the modeling errors. It can be seen that the maximal RMSE and maximal MAE are 0 . 0271 mm and 0 . 0182 mm, respectively. Therefore, the identified parameters of the Hammerstein model are correct, and the Hammerstein model is effective in characterizing asymmetric rate-dependent hysteresis of the DEAP actuator.

3. Compound Controller Design

In this section, a compound controller is developed for the high-precision trajectory tracking control of the DEAP actuator, as shown in Figure 7. It consists of a feedforward hysteresis compensator and a feedback controller. The inverse MPI model is used to be the hysteresis compensator to offset the hysteresis behavior. A novel DTSMC is proposed to be the feedback controller to improve the trajectory tracking precision and maintain the system stability in the presence of unmodeled hysteresis nonlinearity, inverse compensation errors, residual dynamics, parameter uncertainties and external disturbances.

3.1. Design of the Novel DTSMC

In this subsection, the nonlinear disturbance term d ( k ) is considered, which describes the lumped effect of unmodeled hysteresis nonlinearity, inverse compensation errors, residual dynamics, parameter uncertainties and external disturbances. In order to design the proposed controller, Equation (11) is rewritten as follows:
y ( k ) = i = 1 n a a i y ( k i ) + i = 1 n b b i u ( k i ) + d ( k )
Assumption 1.
The nonlinear disturbance term d ( k ) is assumed to be unknown but bounded, which can be expressed as follows:
| d ( k ) | δ 0
where δ 0 is a known positive constant.
Assumption 2.
The change rate of the nonlinear disturbance term d ( k ) is assumed to be unknown but bounded, which can be expressed as follows:
| d ( k ) d ( k 1 ) | T δ
where T is the sampling period, k N , and δ is a known positive constant.
As reported in [32], the nonlinear disturbance term d ( k ) possesses the following property.
Property 1.
d ( k ) = O ( T ) , d ( k ) d ( k 1 ) = O ( T 2 ) , and d ( k ) 2 d ( k 1 ) + d ( k 1 ) = O ( T 3 ) .
In order to complete the controller design and its stability analysis, a controllable canonical representation can be obtained according to the discrete-time plant model (14)
X ( k + 1 ) = A X ( k ) + R u ( k ) + B U ( k ) + D ( k ) y ( k ) = C X ( k )
in which
U ( k ) = u ( k n b + 1 ) , u ( k n b ) , , u ( k 1 ) T ( n b 1 ) × 1 ,
X ( k ) = x ( k n a + 1 ) , , x ( k 1 ) , x ( k ) T n a × 1 ,
A = 0 1 0 0 0 1 a n a a n a 1 a 1 n a × n a , B = 0 0 0 0 0 0 b n b b n b 1 b 2 n a × ( n b 1 ) ,
C = [ 0 , , 0 , 1 ] 1 × n a , R = [ 0 , , 0 , b 1 ] T n a × 1 , D ( k ) = 0 , , 0 , d ( k ) T n a × 1 .
The control objective is to make the system output y ( k ) track a desired displacement in the presence of the nonlinear disturbance D ( k ) . To achieve this objective, a novel DTSMC is presented in this subsection. First, the nonlinear disturbance D ( k ) is estimated. Then, the DTSMC is designed based on a new discrete-time terminal sliding function and a terminal-sliding-model-type reaching law.
Remark 1.
In order to obtain the nonlinear disturbance term D ( k ) , the technique of perturbation estimation [33] is employed under the Assumption 1, the nonlinear disturbance term is estimated by its one-step delayed value, which can be expressed as
D ^ ( k ) = D ( k 1 ) = X ( k ) A X ( k 1 ) R u ( k 1 ) B U ( k 1 )
The state-space model (17) can be rewritten as follows:
X ( k + 1 ) = A X ( k ) + R u ( k ) + B U ( k ) + D ^ ( k ) D ˜ ( k ) y ( k ) = C X ( k )
where D ˜ ( k ) = D ^ ( k ) D ( k ) is the nonlinear disturbance estimation error, it can be further expressed as
D ˜ ( k ) = A X ( k ) X ( k 1 ) + R u ( k ) u ( k 1 ) + B U ( k ) U ( k 1 ) X ( k + 1 ) X ( k )
Define
e ( k ) = y ( k ) y * ( k )
as the trajectory tracking error of the DEAP actuator, where y * ( k ) is the desired displacement, y ( k ) is the actual displacement. Based on the trajectory tracking error, a discrete-time terminal sliding function is proposed as
s ( k ) = κ 1 e ( k ) + κ 2 e ( k 1 ) α
where κ 1 > 0 and κ 2 > 0 are positive constants. The notation e ( k 1 ) α | e ( k 1 ) | α s g n ( e ( k 1 ) ) , α denotes a ratio of two odd integers p and q, where 0 < p < q . The design of 0 < α < 1 will produce a steeper slope of the sliding surface when the system states are close to the origin, and it results in a faster convergence speed [34].
In this work, the DTSMC is designed based on equivalent control, the following reaching law is used
Δ s ( k ) = s ( k + 1 ) s ( k ) = 0
The discrete-time terminal sliding function (22) can be deduced as follows:
s ( k ) = κ 1 e ( k + 1 ) + κ 2 e ( k ) α
s ( k ) = κ 1 y ( k + 1 ) y * ( k + 1 ) + κ 2 e ( k ) α
s ( k ) = κ 1 C ( A X ( k ) + R u ( k ) + B U ( k ) + D ( k ) ) y * ( k + 1 ) + κ 2 e ( k ) α
Based on the nominal model, the equivalent control input u e q ( k ) can be derived as follows:
u e q ( k ) = ( CR ) 1 y * ( k + 1 ) + κ 1 1 s ( k ) κ 1 1 κ 2 e ( k ) α CAX ( k ) CBU ( k ) C D ^ ( k )
It must be pointed out that the equivalent control law u e q ( k ) is derived based on the assumptions that the initial state is on the sliding surface and the estimation error D ˜ ( k ) = 0 . However, in practice, the above assumptions are almost non-existent. Therefore, in order to deal with the problem, defining a terminal-sliding-model-type switching control law as
u s w ( k ) = ( CR ) 1 ξ 1 s ( k ) + ξ 2 s ( k ) β
where ξ 1 > 0 , ξ 2 > 0 , s ( k ) β | s ( k ) | β s g n ( s ( k ) ) with 0 < β < 1 .
Combining the equivalent control law (27) and the switching control law (28) can obtain the overall control law
u d ( k ) = u e q ( k ) + u s w ( k ) = ( CR ) 1 [ y * ( k + 1 ) + κ 1 1 s ( k ) κ 1 1 κ 2 e ( k ) α CAX ( k ) CBU ( k ) C D ^ ( k ) ξ 1 s ( k ) ξ 2 s ( k ) β ]
Substituting Equations (29) into (21) yields
e ( k + 1 ) = C D ( k ) D ^ ( k ) + κ 1 1 s ( k ) κ 1 1 κ 2 e ( k ) α ξ 1 s ( k ) ξ 2 s ( k ) β
Based on the proposed overall control law (29), the expression of s ( k + 1 ) becomes
s ( k + 1 ) = κ 1 e ( k + 1 ) + κ 2 e ( k ) α = s ( k ) κ 1 ξ 1 s ( k ) κ 1 ξ 2 s ( k ) β + κ 1 C D ( k ) D ^ ( k )
Remark 2.
In order to alleviate the chattering phenomenon which may be induced by the s g n ( · ) , we employ the smoother saturation function s a t ( · ) to replace the s g n ( · ) in this work,
s a t ( s ( k ) ) = s g n ( s ( k ) ) , i f | s ( k ) | > σ s ( k ) / σ , i f | s ( k ) | σ
where the parameter σ ensures the band width 2 σ of the sliding mode function s ( k ) . In practice, the trial-and-error method is usually employed to obtain the parameter σ and make a compromise between the chattering phenomenon and the tracking accuracy.
Remark 3.
The switching control law u s w ( k ) can drive the sliding mode dynamics into the DTSM domain with the band width of 2 σ and enhance the robustness feature in the presence of the estimation error D ˜ ( k ) and other perturbations. In addition, D ˜ ( k ) and other perturbations can be tolerated by selecting suitable control gains ξ 1 and ξ 2 .
Remark 4.
The proposed DTSM control law u d ( k ) can accelerate the convergence rate of tracking error by adding the non-linear term e ( k 1 ) α to the discrete-time terminal sliding function s ( k ) . In addition, it can be observed from Equation (29) that the designed DTSMC only needs the system output information y ( k ) ( i . e . , y ( k ) = C X ( k ) ) , that is to say, the state observer is not needed.

3.2. Stability Analysis

In order to analyze the stability of the closed-loop system with the proposed DTSMC, the following assumption condition about the discrete-time terminal sliding surface is required.
Assumption 3.
There is a positive constant ϵ such that | κ 1 C δ | ϵ = O ( T 2 ) and k N .
During the stability analysis, the following proposition is needed.
Proposition 1.
ϵ κ 1 ξ 2 Φ β Φ ( 1 κ 1 ξ 1 ) .
The proof of the proposition can be seen in the Appendix A.
Theorem 1.
Considering the state-space model (17) along with the discrete-time terminal sliding surface (22) and Assumption 3, if we employ the DTSM control law (29), the following sliding mode dynamics can be guaranteed.
(1) The discrete-time sliding mode dynamics from any initial state will enter the region Δ within K * steps, where
Δ = s ( k ) : | s ( k ) | σ = ψ ( β ) · m a x ϵ κ 1 ξ 2 1 β , κ 1 ξ 2 1 κ 1 ξ 1 1 1 β
K * = k * + 1 , k * = s 2 ( 0 ) σ 2 μ 2 , μ = κ 1 ξ 1 σ + [ ψ β ( β ) 1 ] ϵ , ψ ( β ) = 1 + β β 1 β β 1 1 β , Φ = m a x ϵ κ 1 ξ 2 1 β , κ 1 ξ 2 1 κ 1 ξ 1 1 1 β .
where the operator · gives the nearest integer that is less than or equal to the original value.
For any 0 < β < 1 and 0 Φ 1 , the following relationships are met [35]:
1 < ψ ( β ) < 2 , ψ ( β ) Φ [ ψ ( β ) Φ ] β + ψ ( β ) 1
(2) Once the sliding mode dynamics enter the region Δ, it will never escape in the subsequent time.
Proof of Theorem 1.
(1) Define a discrete-time Lyapunov function V ( k ) = s 2 ( k ) , it can be obtained that
Δ V ( k ) = V ( k + 1 ) V ( k ) = s 2 ( k + 1 ) s 2 ( k ) = [ s ( k ) s ( k + 1 ) ] [ s ( k ) + s ( k + 1 ) ] = κ 1 ξ 1 s ( k ) + κ 1 ξ 2 s ( k ) β κ 1 C D ( k ) D ^ ( k ) × 2 s ( k ) κ 1 ξ 1 s ( k ) κ 1 ξ 2 s ( k ) β + κ 1 C D ( k ) D ^ ( k )
In this work, we consider the following two cases if s ( k ) σ .
Case 1: Consider s ( k ) > σ = ψ ( β ) · Φ > 0 , we can obtain that
s ( k ) s ( k + 1 ) = κ 1 ξ 1 s ( k ) + κ 1 ξ 2 s ( k ) β κ 1 C D ( k ) D ^ ( k ) κ 1 ξ 1 s ( k ) + κ 1 ξ 2 ψ β ( β ) Φ β κ 1 C D ( k ) D ^ ( k )
By Proposition 1 and Assumption 3 , Equation (36) can be deduced as
s ( k ) s ( k + 1 ) κ 1 ξ 1 s ( k ) + κ 1 ξ 2 ψ β ( β ) Φ β κ 1 C D ( k ) D ^ ( k ) κ 1 ξ 1 σ + ψ β ( β ) 1 ϵ : = μ
Referring to s ( k ) > 0 , the one-step forward sliding function s ( k + 1 ) > 0 , it can be obtained that
s ( k ) > κ 1 ξ 1 s ( k ) + κ 1 ξ 2 s ( k ) β κ 1 C D ( k ) D ^ ( k )
then
s ( k ) + s ( k + 1 ) = 2 s ( k ) κ 1 ξ 1 s ( k ) κ 1 ξ 2 s ( k ) β + κ 1 C D ( k ) D ^ ( k ) κ 1 ξ 1 s ( k ) + κ 1 ξ 2 ψ β ( β ) Φ β κ 1 C D ( k ) D ^ ( k ) μ
Based on the Equations (39) and (37) , Equation (35) can be further calculation as
Δ V ( k ) = [ s ( k ) s ( k + 1 ) ] [ s ( k ) + s ( k + 1 ) ] = κ 1 ξ 1 s ( k ) + κ 1 ξ 2 s ( k ) β κ 1 C D ( k ) D ^ ( k ) × 2 s ( k ) κ 1 ξ 1 s ( k ) κ 1 ξ 2 s ( k ) β + κ 1 C D ( k ) D ^ ( k ) μ 2
Hence, if s ( k ) > 0 , Δ V ( k ) μ 2 holds.
Case 2: Consider s ( k ) < σ = ψ ( β ) · Φ < 0 with Assumption 3 and Proposition 1, the proof procedure of Δ V ( k ) μ 2 is similar to that of the Case 1.
According to the analysis results of Case 1 and Case 2 when s ( k ) σ , it can be derived from Equation (40) that s 2 ( k ) s 2 ( 0 ) k μ 2 . We can directly solve s 2 ( 0 ) k μ 2 = σ 2 to obtain the solution k * defined in Equation (33). Hence, existing K * = k * + 1 N results in | s ( K * ) | σ .
(2) Next, we will divide two cases to prove that once s ( k ) Δ the sliding mode state s ( k + 1 ) Δ . Define 0 ς 1 , we can obtain | s ( k ) | = ς σ = ς ψ ( β ) Φ .
By Equation (31), it can be obtained that
s ( k + 1 ) = s ( k ) κ 1 ξ 1 s ( k ) κ 1 ξ 2 s ( k ) β + κ 1 C D ( k ) D ^ ( k ) = ( 1 κ 1 ξ 1 ) s ( k ) κ 1 ξ 2 | s ( k ) | β s g n ( s ( k ) ) + κ 1 C D ( k ) D ^ ( k )
Case 1: The proof of σ s ( k + 1 ) σ with s ( k ) 0 .
When s ( k ) 0 , Equation (41) becomes
s ( k + 1 ) = ( 1 κ 1 ξ 1 ) ς ψ ( β ) Φ κ 1 ξ 2 [ ς ψ ( β ) Φ ] β s g n ( ς ψ ( β ) Φ ) + κ 1 C D ( k ) D ^ ( k )
If 0 s ( k ) < 1 , by Proposition 1 and Assumption 3, it can be derived from Equation (42) as
s ( k + 1 ) ( 1 κ 1 ξ 1 ) ς ψ ( β ) Φ [ ς ψ ( β ) ] β ϵ + ϵ ( 1 κ 1 ξ 1 ) ς ψ ( β ) Φ + 1 [ ς ψ ( β ) ] β ( 1 κ 1 ξ 1 ) Φ ( 1 κ 1 ξ 1 ) Φ ς ψ ( β ) ς ψ ( β ) β + 1 ( 1 κ 1 ξ 1 ) Φ < ψ ( β ) Φ = σ
Otherwise, when s ( k ) 1 , utilizing Proposition 1 and Assumption 3, Equation (42) can be deduced as
s ( k + 1 ) ( 1 κ 1 ξ 1 ) ς ψ ( β ) Φ + 1 [ ς ψ ( β ) ] β ϵ ( 1 κ 1 ξ 1 ) ς ψ ( β ) Φ ς ψ ( β ) Φ σ
Next, the proof of σ s ( k + 1 ) with s ( k ) 0 .
Similarly, if 0 s ( k ) < 1 , by Proposition 1 and Assumption 3, it can be derived from Equation (42) as
s ( k + 1 ) ( 1 κ 1 ξ 1 ) ς ψ ( β ) Φ κ 1 ξ 2 [ ς ψ ( β ) Φ ] β s g n ( ς ψ ( β ) Φ ) ε ( 1 κ 1 ξ 1 ) ς ψ ( β ) Φ κ 1 ξ 2 [ ς ψ ( β ) Φ ] β ε ( 1 κ 1 ξ 1 ) Φ ς ψ ( β ) [ ς ψ ( β ) ] β 1 ( 1 κ 1 ξ 1 ) ψ ( β ) Φ σ
Otherwise, when s ( k ) 1 , utilizing Proposition 1 and Assumption 3, Equation (42) can be deduced as
s ( k + 1 ) ( 1 κ 1 ξ 1 ) ς ψ ( β ) Φ κ 1 ξ 2 [ ς ψ ( β ) Φ ] β s g n ( ς ψ ( β ) Φ ) ε ( 1 κ 1 ξ 1 ) ς ψ ( β ) Φ κ 1 ξ 2 [ ς ψ ( β ) Φ ] β ε ϵ ( 1 κ 1 ξ 1 ) Φ > σ
Until now, we have proved that σ s ( k + 1 ) σ with s ( k ) 0 .
Case 2: The proof of σ s ( k + 1 ) σ with s ( k ) < 0 .
When s ( k ) < 0 , Equation (41) becomes
s ( k + 1 ) = ( 1 κ 1 ξ 1 ) ς ψ ( β ) Φ + κ 1 ξ 2 [ ς ψ ( β ) Φ ] β s g n ( ς ψ ( β ) Φ ) + κ 1 C D ( k ) D ^ ( k )
If 1 < s ( k ) < 0 , by Proposition 1 and Assumption 3, it can be derived from Equation (47) as
s ( k + 1 ) ( 1 κ 1 ξ 1 ) ς ψ ( β ) Φ + [ ς ψ ( β ) ] β ( 1 κ 1 ξ 1 ) Φ + ϵ ( 1 κ 1 ξ 1 ) ς ψ ( β ) Φ + [ ς ψ ( β ) ] β ( 1 κ 1 ξ 1 ) Φ + ( 1 κ 1 ξ 1 ) Φ ( 1 κ 1 ξ 1 ) Φ ψ ( β ) ψ ( β ) Φ = σ
Otherwise, when s ( k ) 1 , utilizing Proposition 1 and Assumption 3, Equation (47) can be deduced as
s ( k + 1 ) ( 1 κ 1 ξ 1 ) Φ ς ψ ( β ) + [ ς ψ ( β ) ] β + ϵ ϵ ( 1 κ 1 ξ 1 ) Φ ψ ( β ) Φ = σ
Next, the proof of σ s ( k + 1 ) with s ( k ) < 0 .
Similarly, if 1 < s ( k ) < 0 , by Proposition 1 and Assumption 3, it can be derived from Equation (47) as
s ( k + 1 ) ς ψ ( β ) Φ + κ 1 ξ 1 ς ψ ( β ) Φ + κ 1 ξ 2 [ ς ψ ( β ) Φ ] β ϵ ς ψ ( β ) Φ + κ 1 ξ 1 ς ψ ( β ) Φ + ϵ [ ς ψ ( β ) ] β ϵ ( 1 κ 1 ξ 1 ) ς ψ ( β ) Φ + [ ς ψ ( β ) ] β 1 ϵ ϵ ς ψ ( β ) [ ς ψ ( β ) ] β + 1 ( 1 κ 1 ξ 1 ) Φ > σ
Otherwise, when s ( k ) 1 , utilizing Proposition 1 and Assumption 3, Equation (47) can be deduced as
s ( k + 1 ) ( 1 κ 1 ξ 1 ) Φ ς ψ ( β ) [ ς ψ ( β ) ] β + 1 ( 1 κ 1 ξ 1 ) ψ ( β ) Φ ψ ( β ) Φ = σ
To this end, we have completed the proof of σ s ( k + 1 ) σ with s ( k ) Δ . That is to say, the sliding mode state s ( k ) will enter the region Δ within K * steps and stay in the DTSM domain thereafter by using the proposed DTSM control law. □

4. Experimental Results and Discussion

4.1. Experimental Setup

In order to verify the effectiveness of the proposed DTSMC, the controller is deployed on the DEAP actuator for the high-precision trajectory tracking control. The block diagram and the experimental platform of the DEAP actuator are shown in Figure 8 and Figure 9, respectively. The actuator is manufactured by Danfoss PolyPower A/S Inc. Because DEAP materials are flexible materials, the DEAP actuator needs to be suspended vertically in the experimental platform to generate a certain pre-tension force. Considering the real time executions and the simplicity of the controllers development, a personal computer running Windows 10 operating system and Matlab 2017a software is employed. A high resolution split LVDT displacement sensor is utilized to measure the output displacement of the DEAP actuator. Considering the driving voltage of DEAP materials is 0–2500 V, a high voltage source S15-3P is employed. For realizing the data interaction between MATLAB 2017a and the DEAP actuator, a data acquisition card PCI-1710 is used. In this work, the sampling time T is set to 0 . 01 s, and a white Gaussian noise is added to the controller as the external disturbance.

4.2. Experimental Results

In this subsection, the proposed DTSMC is tested for the high-precision trajectory tracking control of the DEAP actuator. First, we select the optimal control parameters of the proposed controller. If α is fixed, a lager κ 2 can guarantee a better performance. In addition, if choosing β = 1 / 2 which implies 1 / β = 1 / ( 1 β ) = 2 and σ = O ( T 2 ) , we can obtain the smallest boundary of the final sliding mode motion. Disturbances can be well suppressed with the increasing of ξ 1 and ξ 2 ; however, the chattering phenomenon will be also increased in the control signal. Through trial and error, the selected optimal control parameters are presented in Table 5.
This work employs the real experimental data to verify the effectiveness of the proposed controller. For comparison, the conventional DSMC with the hysteresis compensator is employed. The control parameters of the DSMC are the same as those of the proposed controller, which are presented in Table 5.
During experimental studies, different frequencies’ input signals are generated by S15-3P and applied to drive the DEAP actuator. The corresponding tracking results of the DEAP actuator are shown in Figure 10, Figure 11, Figure 12, Figure 13, Figure 14 and Figure 15.
The desired displacement has a frequency of 0.1 Hz, the tracking results of the DSMC and the proposed DTSMC are depicted in Figure 10. The RMSE of the DSMC is 0 . 0565 mm, and the MAE is 0 . 0207 mm. In contrast, for the proposed DTSMC, the RMSE is 0 . 0536 mm, and the MAE is 0 . 0122 mm. Compared with the DSMC, the proposed DTSMC reduces the RMSE and MAE by 5 % and 41 % . In addition, from Figure 10c,d, one can see that the hysteresis phenomenon has been substantially compensated by the DTSMC, and the DSMC can also basically compensate the hysteresis phenomenon at a low frequency of 0 . 1 Hz.
In Figure 11, the desired displacement has a frequency of 0.2 Hz. The RMSE of the DSMC is 0 . 0942 mm, and the MAE is 0 . 0504 mm. In contrast, for the proposed DTSMC, the RMSE is 0 . 0612 mm, and the MAE is 0 . 0215 mm. Compared with the DSMC, the proposed DTSMC reduces the RMSE and MAE by 35 % and 57 % . In addition, compared with the desired displacement of 0 . 1 Hz, the proposed DTSMC increases the RMSE and MAE by 12 % and 43 % . However, the DSMC increases the RMSE and MAE by 40 % and 59 % , which are all obviously larger than those of the proposed DTSMC. Furthermore, from the Figure 12, Figure 13 and Figure 14, it can be seen that the hysteresis compensations of the proposed DTSMC gradually exceeds those of the DSMC with increasing frequency of the input signal. For more detail trajectory tracking-error results, see the computational results in Table 6. In addition, the overshoots of the DSMC at the turning points are significantly larger than those of the proposed DTSMC. Therefore, it can be concluded that the proposed DTSMC can improve the tracking control accuracy of the DEAP actuator with increasing frequency to 0 . 5 Hz.
In order to further demonstrate the effectiveness of the proposed DTSMC, a hybrid frequency input signal in range of 0 . 1 to 0 . 5 Hz is used to be the desired displacement, the tracking results of the two controllers are shown in Figure 15. Table 6 gives the computational results of two indices. It can be seen that the RMSE and MAE of the proposed DTSMC is 0 . 0500 mm and 0 . 0242 mm, respectively. In contrast, the RMSE and MAE of the DSMC is 0 . 1536 mm and 0 . 0990 mm, respectively. As compared with the DSMC, the proposed DTSMC considerably reduces the RMSE and MAE by 67 % and 76 % . From Figure 15c,d, one can see that the hysteresis compensation of the proposed DTSMC is more obvious than that of the DSMC. Therefore, the proposed DTSMC has better performance than the DSMC in following a hybrid frequency input signal.

4.3. Discussion

In this paper, we propose a compound controller for the high-precision trajectory tracking control of the DEAP actuator. It consists of a feedforward hysteresis compensator (i.e., the inverse MPI model) and a feedback controller (i.e., the novel DTSMC). In order to achieve the high-precision trajectory tracking control, the inverse MPI model is first applied to compensate the hysteresis of the DEAP actuator, and then the proposed DTSMC is employed to overcome the estimation error D ˜ ( k ) and other perturbations. From Figure 10 and Table 6, it can be seen that the DSMC and the proposed DTSMC can achieve good tracking control when the input signal frequency is 0 . 1 Hz, that is to say, the hysteresis characteristic is well compensated, and the estimation error D ˜ ( k ) and other perturbations are well overcome. From Figure 11, Figure 12, Figure 13, Figure 14 and Figure 15, and Table 6, one can see that, whether the input signal is a high frequency or a hybrid frequency, the proposed DTSMC can always achieve excellent tracking control performances. In contrast, the tracking errors of the DSMC will be significantly increased with increasing input signal frequency, and it has large overshoots at the turning points. To summarize, the proposed DTSMC is effective for the high-precision trajectory tracking control of the DEAP actuator. It provides faster error convergence, more precise trajectory tracking control, less chattering and better robustness than the DSMC. However, it should be pointed out that the nonlinear disturbance term upper bound and its change rate need to be known in the design process of the DTSMC, which are difficult to achieve in practice. In future work, we will try to apply the adaptive technology to completely compensate the uncertainty of unknown upper bound. Another possible research topic of this work may be the design of an adaptive discrete-time terminal sliding mode control without a inverse hysteresis compensator to simplify the design process of the controller and optimize the controller performance.

5. Conclusions

The inherit asymmetric rate-dependent hysteresis characterization in the DEAP actuator is a barrier to achieving high-precision trajectory tracking control. In this paper, we establish a Hammerstein model to characterize the asymmetric rate-dependent hysteresis and employ LM algorithm to identify the parameters of the MPI model and the ARX model. In order to suppress the asymmetric hysteresis behavior, the inverse MPI model is used to be the hysteresis feedforward compensator. Furthermore, to overcome the challenge of the estimation error D ˜ ( k ) and other perturbations, based on the hysteresis feedforward compensator, a novel DTSMC is proposed for the high-precision tracking control of the DEAP actuator, and the system stability is analyzed. A conventional DSMC with the hysteresis compensator is implemented for comparison.
The above two control schemes are tested under different frequencies’ input signals, including 0 . 1 Hz, 0 . 2 Hz, 0 . 3 Hz, 0 . 4 Hz, 0 . 5 Hz and a hybrid frequency. Noting that, in the DSMC and DTSMC schemes, the switching control may make the chattering occur. To eliminate the chattering, we employ the smoother saturation function to replace the signum function in this work. When the input signal frequency is increased, or the input signal is a hybrid frequency, compared with the DSMC with hysteresis compensator, experimental results show that the proposed DTSMC with the hysteresis compensator can improve the tracking accuracy, robustness, and error convergence of the DEAP actuator. The hysteresis compensation effect of the proposed control scheme is significantly higher than that of the DSMC scheme.

Author Contributions

M.L. proposed the controller and wrote the manuscript. M.L. and Z.J. designed experimental and analyzed the results. Y.L. and Q.W. provided some improvement strategies to improve the accuracy of the model and helped verify the effectiveness of the control strategy. M.L. and Y.L. revised and edited the paper.

Funding

This research was funded by the National Science Foundation of China (61472037, 61375100).

Conflicts of Interest

All authors declare that there is no conflict of interest

Appendix A

The proof of the Proposition 1:
If σ = ψ ( β ) ϵ κ 1 ξ 2 1 β , then Φ = ϵ κ 1 ξ 2 1 β , ϵ = κ 1 ξ 2 Φ β . Because ϵ κ 1 ξ 2 1 β κ 1 ξ 2 1 κ 1 ξ 1 1 1 β , ϵ 1 β κ 1 ξ 2 ( 1 κ 1 ξ 1 ) β , ϵ β 1 ( 1 κ 1 ξ 1 ) β κ 1 ξ 2 , ϵ β Φ ( 1 κ 1 ξ 1 ) β , ϵ Φ ( 1 κ 1 ξ 1 ) . It can be obtained that κ 1 ξ 2 Φ β Φ ( 1 κ 1 ξ 1 ) .
If σ = ψ ( β ) κ 1 ξ 2 1 κ 1 ξ 1 1 1 β , then Φ = κ 1 ξ 2 1 κ 1 ξ 1 1 1 β , Φ 1 β = κ 1 ξ 2 1 κ 1 ξ 1 , Φ Φ β = κ 1 ξ 2 1 κ 1 ξ 1 , κ 1 ξ 2 Φ β = Φ ( 1 κ 1 ξ 1 ) , Φ β = 1 κ 1 ξ 1 κ 1 ξ 2 Φ = 1 κ 1 ξ 1 κ 1 ξ 2 κ 1 ξ 2 1 κ 1 ξ 1 1 1 β , Φ β = κ 1 ξ 2 1 κ 1 ξ 1 1 1 β 1 = κ 1 ξ 2 1 κ 1 ξ 1 β 1 β , ϵ κ 1 ξ 2 κ 1 ξ 2 1 κ 1 ξ 1 β 1 β , ϵ κ 1 ξ 2 κ 1 ξ 2 1 κ 1 ξ 1 1 1 β β . It can be obtained that ϵ κ 1 ξ 2 Φ β .
Hence, based on the above formulas, we can obtain that ϵ κ 1 ξ 2 Φ β Φ ( 1 κ 1 ξ 1 ) .

References

  1. Benslimane, M.Y.; Kiil, H.E.; Tryson, M.J. Dielectric electro-active polymer push actuators: Performance and challenges. Polym. Int. 2010, 59, 415–421. [Google Scholar] [CrossRef]
  2. Gu, G.Y.; Zhu, J.; Zhu, L.M.; Zhu, X. A survey on dielectric elastomer actuators for soft robots. Bioinspir. Biomim. 2017, 12, 011003. [Google Scholar] [CrossRef] [PubMed]
  3. Iskandarani, Y.; Karimi, H.R. Dynamic characterization for the dielectric electroactive polymer fundamental sheet. Int. J. Adv. Manuf. Technol. 2013, 66, 1457–1466. [Google Scholar] [CrossRef]
  4. Kovacs, G.; Lochmatter, P.; Wissler, M. An arm wrestling robot driven by dielectric elastomer actuators. Smart Mater. Struct. 2007, 16, S306. [Google Scholar] [CrossRef]
  5. Shintake, J.; Rosset, S.; Schubert, B.E.; Floreano, D.; Shea, H.R. A foldable antagonistic actuator. IEEE/ASME Trans. Mechatron. 2015, 20, 1997–2008. [Google Scholar] [CrossRef]
  6. Godaba, H.; Li, J.; Wang, Y.; Zhu, J. A soft jellyfish robot driven by a dielectric elastomer actuator. IEEE Robot. Autom. Lett. 2016, 1, 624–631. [Google Scholar] [CrossRef]
  7. Truong, B.N.M.; Ahn, K.K. Inverse modeling and control of a dielectric electro-active polymer smart actuator. Sens. Actuators A Phys. 2015, 229, 118–127. [Google Scholar] [CrossRef]
  8. Liu, Y.; Gao, X.; Li, Y. Giant magnetostrictive actuator nonlinear dynamic Jiles–Atherton model. Sens. Actuators A Phys. 2016, 250, 7–14. [Google Scholar] [CrossRef]
  9. Zou, J.; Gu, G. Modeling the viscoelastic hysteresis of dielectric elastomer actuators with a modified rate-dependent prandtl–ishlinskii model. Polymers 2018, 10, 525. [Google Scholar] [CrossRef]
  10. Al Janaideh, M.; Aljanaideh, O. Further results on open-loop compensation of rate-dependent hysteresis in a magnetostrictive actuator with the Prandtl-Ishlinskii model. Mech. Syst. Signal Process. 2018, 104, 835–850. [Google Scholar] [CrossRef]
  11. Wang, H.; Song, G. Innovative NARX recurrent neural network model for ultra-thin shape memory alloy wire. Neurocomputing 2014, 134, 289–295. [Google Scholar] [CrossRef]
  12. Truong, B.N.M.; Ahn, K.K. Modeling, control and experimental investigation of time-average flow rate of a DEAP actuator based diaphragm pump. Int. J. Precis. Eng. Manuf. 2017, 18, 1119–1129. [Google Scholar] [CrossRef]
  13. Rizzello, G.; Naso, D.; York, A.; Seelecke, S. Modeling, identification, and control of a dielectric electro-active polymer positioning system. IEEE Trans. Control Syst. Technol. 2015, 23, 632–643. [Google Scholar] [CrossRef]
  14. Jones, R.W.; Sarban, R. Inverse grey-box model-based control of a dielectric elastomer actuator. Smart Mater. Struct. 2012, 21, 075019. [Google Scholar] [CrossRef]
  15. Hoffstadt, T.; Maas, J. Adaptive sliding-mode position control for dielectric elastomer actuators. IEEE/ASME Trans. Mechatron. 2017, 22, 2241–2251. [Google Scholar] [CrossRef]
  16. Qian, J.; Ji, C.; Pan, N.; Wu, J. Improved sliding mode control for permanent magnet synchronous motor speed regulation system. Appl. Sci. 2018, 8, 2491. [Google Scholar] [CrossRef]
  17. Zhihong, M.; Yu, X.H. Terminal sliding mode control of MIMO linear systems. IEEE Conf. Decis. Control 1996, 4, 4619–4624. [Google Scholar]
  18. Vo, A.T.; Kang, H.J. An adaptive neural non-singular fast-terminal sliding-mode control for industrial robotic manipulators. Appl. Sci. 2018, 8, 2562. [Google Scholar] [CrossRef]
  19. Wang, G.; Xu, Q. Adaptive terminal sliding mode control for motion tracking of a micropositioning system. Asian J. Control 2018, 20, 1241–1252. [Google Scholar] [CrossRef]
  20. Dinh, T.X.; Ahn, K.K. Adaptive-gain fast nonsingular terminal sliding mode for position control of a piezo positioning stage. Proc. Inst. Mech. Eng. Part I J. Syst. Control Eng. 2018, 232, 994–1014. [Google Scholar] [CrossRef]
  21. Furuta, K. Sliding mode control of a discrete system. Syst. Control Lett. 1990, 14, 145–152. [Google Scholar] [CrossRef]
  22. Chakrabarty, S.; Bandyopadhyay, B. A generalized reaching law for discrete time sliding mode control. Automatica 2015, 52, 83–86. [Google Scholar] [CrossRef]
  23. Nguyen, M.L.; Chen, X.; Yang, F. Discrete-time quasi-sliding-mode control with prescribed performance function and its application to piezo-actuated positioning systems. IEEE Trans. Ind. Electron. 2018, 65, 942–950. [Google Scholar] [CrossRef]
  24. Huber, O.; Acary, V.; Brogliato, B. Lyapunov stability and performance analysis of the implicit discrete sliding mode control. IEEE Trans. Autom. Control 2016, 61, 3016–3030. [Google Scholar] [CrossRef]
  25. Ma, H.; Li, Y.; Xiong, Z. Discrete-Time Sliding-Mode control with enhanced power reaching law. IEEE Trans. Ind. Electron. 2019, 66, 4629–4638. [Google Scholar] [CrossRef]
  26. Xu, Q. Digital integral terminal sliding mode predictive control of piezoelectric-driven motion system. IEEE Trans. Ind. Electron. 2016, 63, 3976–3984. [Google Scholar] [CrossRef]
  27. Xu, Q.; Cao, Z. Piezoelectric positioning control with output-based discrete-time terminal sliding mode control. IET Control Theory Appl. 2017, 11, 694–702. [Google Scholar] [CrossRef]
  28. Levenberg, K. A method for the solution of certain non-linear problems in least squares. Quart. Appl. Math. 1944, 2, 164–168. [Google Scholar] [CrossRef] [Green Version]
  29. Marquardt, D.W. An algorithm for least-squares estimation of nonlinear parameters. J. Soc. Ind. Appl. Math. 1963, 11, 431–441. [Google Scholar] [CrossRef]
  30. Krejci, P.; Kuhnen, K. Inverse control of systems with hysteresis and creep. IEE Proc.-Control Theory Appl. 2001, 148, 185–192. [Google Scholar] [CrossRef] [Green Version]
  31. Kuhnen, K. Modeling, identification and compensation of complex hysteretic nonlinearities: A modified Prandtl-Ishlinskii approach. Eur. J. Control 2003, 9, 407–418. [Google Scholar] [CrossRef]
  32. Abidi, K.; Jian-Xin, X. A discrete-time integral sliding mode control approach for output tracking with state estimation. IFAC Proc. Vol. 2008, 41, 14199–14204. [Google Scholar] [CrossRef]
  33. Su, W.C.; Drakunov, S.V.; Ozguner, U. An O (T2) boundary layer in sliding mode for sampled-data systems. IEEE Trans. Autom. Control 2000, 45, 482–485. [Google Scholar]
  34. Abidi, K.; Xu, J.X.; She, J.H. A discrete-time terminal sliding-mode control approach applied to a motion control problem. IEEE Trans. Ind. Electron. 2009, 56, 3619–3627. [Google Scholar] [CrossRef]
  35. Li, S.; Du, H.; Yu, X. Discrete-time terminal sliding mode control systems based on Euler’s discretization. IEEE Trans. Autom. Control 2014, 59, 546–552. [Google Scholar] [CrossRef]
Figure 1. The working principle of dielectric electro-active polymer (DEAP) materials. (a) Power cut: Initial state of DEAP materials. (b) Power supply: Actuation direction of DEAP materials.
Figure 1. The working principle of dielectric electro-active polymer (DEAP) materials. (a) Power cut: Initial state of DEAP materials. (b) Power supply: Actuation direction of DEAP materials.
Applsci 09 02625 g001
Figure 2. Hammerstein model for DEAP actuator.
Figure 2. Hammerstein model for DEAP actuator.
Applsci 09 02625 g002
Figure 3. The play operator.
Figure 3. The play operator.
Applsci 09 02625 g003
Figure 4. The one side dead-zone operator.
Figure 4. The one side dead-zone operator.
Applsci 09 02625 g004
Figure 5. The rate-dependent hysteresis loops of the DEAP actuator.
Figure 5. The rate-dependent hysteresis loops of the DEAP actuator.
Applsci 09 02625 g005
Figure 6. The hysteresis loops of the Hammerstein model under different frequencies. (a) 0.1 Hz. (b) 0.2 Hz. (c) 0.3 Hz. (d) 0.4 Hz. (e) 0.5 Hz. (f) 0.1/0.2/0.3/0.4/0.5 Hz. Red dotted line: output of the Hammerstein model. Blue solid line: output of the DEAP actuator.
Figure 6. The hysteresis loops of the Hammerstein model under different frequencies. (a) 0.1 Hz. (b) 0.2 Hz. (c) 0.3 Hz. (d) 0.4 Hz. (e) 0.5 Hz. (f) 0.1/0.2/0.3/0.4/0.5 Hz. Red dotted line: output of the Hammerstein model. Blue solid line: output of the DEAP actuator.
Applsci 09 02625 g006
Figure 7. The structure of the compound controller.
Figure 7. The structure of the compound controller.
Applsci 09 02625 g007
Figure 8. Block diagram of the DEAP actuator.
Figure 8. Block diagram of the DEAP actuator.
Applsci 09 02625 g008
Figure 9. Experimental platform of the DEAP actuator.
Figure 9. Experimental platform of the DEAP actuator.
Applsci 09 02625 g009
Figure 10. The tracking results under 0.1 Hz input signal. (a) Displacement output. (b) Tracking error. (c) Desired displacement versus actual displacement of the DSMC. (d) Desired displacement versus actual displacement of the DTSMC.
Figure 10. The tracking results under 0.1 Hz input signal. (a) Displacement output. (b) Tracking error. (c) Desired displacement versus actual displacement of the DSMC. (d) Desired displacement versus actual displacement of the DTSMC.
Applsci 09 02625 g010
Figure 11. The tracking results under 0.2 Hz input signal. (a) Displacement output. (b) Tracking error. (c) Desired displacement versus actual displacement of the DSMC. (d) Desired displacement versus actual displacement of the DTSMC.
Figure 11. The tracking results under 0.2 Hz input signal. (a) Displacement output. (b) Tracking error. (c) Desired displacement versus actual displacement of the DSMC. (d) Desired displacement versus actual displacement of the DTSMC.
Applsci 09 02625 g011
Figure 12. The tracking results under 0.3 Hz input signal. (a) Displacement output. (b) Tracking error. (c) Desired displacement versus actual displacement of the DSMC. (d) Desired displacement versus actual displacement of the DTSMC.
Figure 12. The tracking results under 0.3 Hz input signal. (a) Displacement output. (b) Tracking error. (c) Desired displacement versus actual displacement of the DSMC. (d) Desired displacement versus actual displacement of the DTSMC.
Applsci 09 02625 g012
Figure 13. The tracking results under 0.4 Hz input signal. (a) Displacement output. (b) Tracking error. (c) Desired displacement versus actual displacement of the DSMC. (d) Desired displacement versus actual displacement of the DTSMC.
Figure 13. The tracking results under 0.4 Hz input signal. (a) Displacement output. (b) Tracking error. (c) Desired displacement versus actual displacement of the DSMC. (d) Desired displacement versus actual displacement of the DTSMC.
Applsci 09 02625 g013
Figure 14. The tracking results under 0.5 Hz input signal. (a) Displacement output. (b) Tracking error. (c) Desired displacement versus actual displacement of the DSMC. (d) Desired displacement versus actual displacement of the DTSMC.
Figure 14. The tracking results under 0.5 Hz input signal. (a) Displacement output. (b) Tracking error. (c) Desired displacement versus actual displacement of the DSMC. (d) Desired displacement versus actual displacement of the DTSMC.
Applsci 09 02625 g014
Figure 15. The tracking results under the hybrid frequency input signal. (a) Displacement output. (b) Tracking error. (c) Desired displacement versus actual displacement of the DSMC. (d) Desired displacement versus actual displacement of the DTSMC.
Figure 15. The tracking results under the hybrid frequency input signal. (a) Displacement output. (b) Tracking error. (c) Desired displacement versus actual displacement of the DSMC. (d) Desired displacement versus actual displacement of the DTSMC.
Applsci 09 02625 g015
Table 1. The identified weighting vectors of the Prandtl-Ishlinskii (PI) model.
Table 1. The identified weighting vectors of the Prandtl-Ishlinskii (PI) model.
Index n12345678
ϑ n −1.31881.9548−0.93920.84911.36190.31450.26480.2226
r h n 0.00070.00580.01960.04660.09090.15710.24950.3724
Index n9101112131415
ϑ n 0.04560.25841.5619−0.2742−0.3569−0.3978−0.4230
r h n 0.53030.72740.96821.25701.59821.99612.4551
Table 2. The identified weighting vectors of the dead-zone (DZ) model.
Table 2. The identified weighting vectors of the dead-zone (DZ) model.
Index m12345678
θ m −17.5517−15.7191−7.802013.524558.4796−30.1681−0.99780.4027
r d m 0.00010.00130.00640.02010.04910.10190.18880.3220
Index m9101112131415
θ m −0.05770.16850.26330.18060.28440.26200.5739
r d m 0.51580.78621.15101.63022.24533.02013.9799
Table 3. The identified parameters of the autoregressive with exogenous inputs (ARX) model.
Table 3. The identified parameters of the autoregressive with exogenous inputs (ARX) model.
a 1 a 2 b 1 b 2
−1.95600.9545−0.30950.3086
Table 4. The modeling errors of the Hammerstein model under different frequencies.
Table 4. The modeling errors of the Hammerstein model under different frequencies.
Frequency/Hz0.10.20.30.40.50.1/0.2/0.3/0.4/0.5
RMSE (mm)0.00520.00980.01450.01900.02310.0271
MAE (mm)0.00400.00760.01130.01490.01820.0083
Table 5. The control parameters of different controllers.
Table 5. The control parameters of different controllers.
Controller κ 1 κ 2 α ξ 1 ξ 2 β σ
DSMC3.50-7/110.250.501/20.01
DTSMC3.500.557/110.250.501/20.01
Table 6. The tracking-error results of the DSMC and DTSMC controllers at different frequencies.
Table 6. The tracking-error results of the DSMC and DTSMC controllers at different frequencies.
Frequency/Hz0.10.20.30.40.50.1/0.2/0.3/0.4/0.5
RMSE (mm)DSMC0.05650.09420.10540.15140.22260.1536
DTSMC0.05360.06120.07500.07610.08470.0500
MAE (mm)DSMC0.02070.05040.07880.12550.18340.0990
DTSMC0.01220.02150.03140.03890.04560.0242

Share and Cite

MDPI and ACS Style

Li, M.; Wang, Q.; Li, Y.; Jiang, Z. Modeling and Discrete-Time Terminal Sliding Mode Control of a DEAP Actuator with Rate-Dependent Hysteresis Nonlinearity. Appl. Sci. 2019, 9, 2625. https://doi.org/10.3390/app9132625

AMA Style

Li M, Wang Q, Li Y, Jiang Z. Modeling and Discrete-Time Terminal Sliding Mode Control of a DEAP Actuator with Rate-Dependent Hysteresis Nonlinearity. Applied Sciences. 2019; 9(13):2625. https://doi.org/10.3390/app9132625

Chicago/Turabian Style

Li, Mengmeng, Qinglin Wang, Yuan Li, and Zhaoguo Jiang. 2019. "Modeling and Discrete-Time Terminal Sliding Mode Control of a DEAP Actuator with Rate-Dependent Hysteresis Nonlinearity" Applied Sciences 9, no. 13: 2625. https://doi.org/10.3390/app9132625

APA Style

Li, M., Wang, Q., Li, Y., & Jiang, Z. (2019). Modeling and Discrete-Time Terminal Sliding Mode Control of a DEAP Actuator with Rate-Dependent Hysteresis Nonlinearity. Applied Sciences, 9(13), 2625. https://doi.org/10.3390/app9132625

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