Next Article in Journal
A Method for Enterprise Architecture Model Slicing
Next Article in Special Issue
Image-Based Radical Identification in Chinese Characters
Previous Article in Journal
The Tolerance, Absorption, and Transport Characteristics of Macleaya cordata in Relation to Lead, Zinc, Cadmium, and Copper under Hydroponic Conditions
Previous Article in Special Issue
Transformer-Based Multilingual Speech Emotion Recognition Using Data Augmentation and Feature Fusion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Simulation of Adaptive Radial Basis NN-Based Non-Singular Fast Terminal Sliding Mode Control with Time Delay Estimator for Precise Control of Dual-Axis Manipulator

1
Department of Electrical Engineering, National Central University, Taoyuan City 320317, Taiwan
2
Department of Electronic Engineering, Lunghwa University of Science and Technology, Taoyuan City 333326, Taiwan
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(19), 9605; https://doi.org/10.3390/app12199605
Submission received: 7 August 2022 / Revised: 15 September 2022 / Accepted: 19 September 2022 / Published: 24 September 2022
(This article belongs to the Special Issue Machine Learning for Language and Signal Processing)

Abstract

:

Featured Application

The proposed method is suitable for application in robotic arm systems.

Abstract

Robotic manipulators can reduce the cost of production and improve productivity; however, controlling a manipulator to follow a desired trajectory is a thorny problem. In this study, we introduced various forms of interference to facilitate the modeling of a dual-axis manipulator. The interference associated with the payload is handled by an adaptive radial basis neural network (ARBNN) controller, while other interference is estimated by a time delay estimator (TDE). The control signal is output by a non-singular fast terminal sliding mode controller (NFTSMC) to minimize further interference. Since the proposed controller can deal with the payload, system uncertainties, external disturbances, friction, and backlash, compared with conventional control methods, it has better tracking accuracy and stability.

1. Introduction

The application of industrial manipulators can largely eliminate human error, boost product quality, and increase productivity [1,2]. Numerous researchers have sought to develop control methods capable of ensuring that the manipulator follows a desired trajectory. Su et al. [3] proposed a saturated PID controller to deal with motor failure in mechanical arms resulting from excessive input torque. Martinez et al. [4] proposed a static neural network controller that uses a novel Gaussian basis function to facilitate adaptive tracking control. Their controller switches between a neural network controller and a robust proportional derivative (PD) controller. Wai et al. [5] developed an SMC-based position tracking method based on a high-order dynamics model followed by a fuzzy neural network to simulate the original sliding mode control for stability analysis. Zhu et al. [6] developed a novel composite control scheme that introduces uncertainties into the dynamic modeling of a manipulator and then uses a novel dead zone compensation method to attenuate uncertainties and the corresponding nonlinearities. He et al. [7] proposed an adaptive neural network controller to ensure control over the tracking of uncertain n-links manipulators with full-state constraints. Van et al. [8] achieved trajectory tracking control in a system that combines the advantages of a sliding mode controller with those of a backstepping controller. Liu et al. [9] proposed an adaptive bias radial basis neural network controller with local and global bias to improve approximation accuracy in dynamic models with significant bias. When the input deviates from the approximation domain due to excessive payload, a neural network controller degenerates into a PID controller to return the system to its original state. Note, however, that most previous approaches to manipulator control systems considered a limited number of interference sources and dealt with them as a whole.
Our objective in the current paper was to deal with the above issues from two perspectives. First, we established an enhanced second-order dynamic model of the manipulator, which considers payload-related issues as well as those pertaining to model uncertainties, friction, external disturbances, and backlash. We then developed a composite controller based on NFTSMC, which deals with the two classes of interference (payload and other), respectively, using ARBNN and TDE. The proposed control strategy enhances accuracy while reducing the influence of variations in interferences. It also reduces dependence on parameter adjustment and does not require an interference upper bound.
The remainder of this paper is organized as follows: Section 2 introduces the process of designing our enhanced dynamic model of the dual-axis manipulator. Section 3 outlines our advanced control strategy and stability analysis. Section 4 outlines the simulations used to assess the performance of the proposed algorithm. Conclusions are presented in Section 5.

2. Enhanced Dynamic Model

In this section, we examine a dynamic model of a dual-axis manipulator, which can be expressed as follows:
M q q ¨ + C q , q ˙ q ˙ + G q = U
where M θ refers to the inertia matrix, C θ , θ ˙ is the centripetal and Coriolis force vector, G θ is the gravity vector, U is control torque, and q ,   q ˙ , q ¨ respectively indicate the angle, angular velocity, and angular acceleration of the joint. We first consider friction [10], which can be divided into static friction, viscous friction, and Coulomb friction.
In Figure 1, B , Q C , and Q s , respectively, describe the slope of viscous friction, static friction, and Coulomb friction. Note that when the velocity of the manipulator is increased, static friction decreases rapidly, such that viscous friction dominates. Thus, we consider only viscous friction and Coulomb friction in the current study, the formulas of which can be expressed as follow
f q ˙ = B q ˙ + Q C
and
Q C =     0             q ˙ = 0   Q C +         q ˙ > 0 Q C         q ˙ < 0
Thus, the dynamic model of the manipulator can be rewritten as:
M q q ¨ + C q , q ˙ q ˙ + G q + f q ˙ = U
Next, we consider the external disturbance d t to which the manipulator is most susceptible in practical applications. There are numerous forms of external disturbance, and they can be very difficult to model; therefore, we used white noise with an upper bound to represent external disturbances acting on the two joints of the manipulator, which can be expressed as follows:
d t d ¯ t
where d ¯ t is an unknown positive constant vector. Thus, the dynamic model of the manipulator can be rewritten as
M q q ¨ + C q , q ˙ q ˙ + G q + f q ˙ + d t = U
We also considered the effects of payload properties on the manipulator (a common issue in practical applications). Consider the example in Figure 2, illustrating the process of grasping the payload, moving it to the destination, and putting it down. During this process, the payload is converted into forces and torques acting on the joints and links, such that the dynamic model of the manipulator can be rewritten as follows:
M q q ¨ + C q , q ˙ q ˙ + G q + f q ˙ + d t + U p = U
Next, we considered the uncertainties in the manipulator model associated with system parameters. The basic dynamic model outlined above includes three main parameters: inertia matrix M q , the centripetal force/Coriolis force matrix C q , q ˙ and the gravity vector   G q . In the process of establishing a dynamic model, these parameters are generally derived empirically. Thus, we decompose the system parameters into two classes (known items and model uncertainties), which can be expressed as follows:
M q = M 0 q + Δ M q
C q , q ˙ = C 0 q , q ˙ + Δ C q , q ˙
G q = G 0 q + Δ G q
Thus, the dynamic model of the manipulator can be expressed as follows:
M 0 q q ¨ + C 0 q , q ˙ q ˙ + G 0 q + Δ M q q ¨ + Δ C q , q ˙ q ˙ + Δ G q + f q ˙ + d t + U p = U
Finally, we consider backlash [11,12], which refers to the combined time delay effects of occlusions among gears, chains, and belts used in the transmission of power from the motor to other components (see Figure 3). This phenomenon can be described using a class of nonlinear equations as follows [13]:
V U = Υ U b 2   ,   i f   U ˙ > 0   a n d   U U t 1 + b / 2 Υ U + b 2   ,   i f   U ˙ < 0   a n d   U U t 1 b / 2   V U t   ,       o t h e r w i s e
where Υ refers to the slope, b indicates the size of the dead zone, U indicates the control torque, and V U indicates the control torque after backlash. This relationship is illustrated in Figure 4.
We can rewrite Equation (12) as follows:
V U = Υ U + Γ U
where Γ U is expressed as
Γ U = Υ b 2 ,             i f   U ˙ > 0   a n d   U U t 1 + b / 2 Υ b 2 ,               i f   U ˙ < 0   a n d   U U t 1 b / 2 V U t Υ U   ,           o t h e r w i s e
Our enhanced dynamic model of a dual-axis manipulator under the effects of multiple forms of interference can be expressed as follows:
M 0 q ¨ + C 0 q ˙ + G 0 + Ξ q , q ˙ , q ¨ , t + U P = V U = Υ U + Γ U      
where, for convenience, we specify that
Ξ q , q ˙ , q ¨ , t = Δ M q ¨ + Δ C q ˙ + Δ G + f q ˙ + d t

3. Controller Design

3.1. Preliminaries

As shown in Figure 5a, the radial basis neural network (RBNN) provides a simple control structure and the ability to approximate any nonlinear function with arbitrary accuracy without the need for lengthy calculations [14,15]. RBNN employs an input layer, a hidden layer, and an output layer:
  • The input layer consists of an input vector X i , X i = x i 1   ,   x i 2 x i m T R m × 1 , i = 1 , , n ,   z = 1 ,   , m , where i indicates the number of neural networks and z indicates the number of inputs for each neural network.
  • Note that there is no weight connection between the hidden and input layers, and the hidden layer is composed of neurons. The activation function in a neuron comprises the radial basis function, which is generally a Gaussian function (Figure 5b) and can be expressed as follows:
    μ i j = exp x i z b i j 2 σ i j 2 ,   j = 1 ,   ,   k
    where 𝑗 indicates the number of neurons in the hidden layer of each neural network, while b i j R n × k and σ i j R n × k refer to the mean vector and the width of the Gaussian function of the j t h node in the hidden layer in the i t h neural network.
  • The output layer describes the control signal of RBNN as follows:
    N ^ X i = f 1 , , f n T = W T
    where W refers to the weight between the hidden layer and the output layer and μ is the output of the hidden layer, as follows:
    W = w 11 w 1 k 0 0 0 0 0 0 w n 1 w n 2 w n k T R n k × n
    μ = μ 11 μ 1 k μ n 1 μ n k T R n k × 1
The simple structure of the RBF neural network described above allows it to be used online to estimate expected values over time, the accuracy of which is proportional to the number of neurons in the hidden layer.
We also employed a sliding mode control (SMC) [16,17,18] nonlinear control method. As shown in Figure 6, SMC uses a discontinuous control signal to slide the nonlinear system between two states, forcing the system to finally enter a desired control result through the sliding surface. The system to be controlled is first written out in the form of a state-space equation, an example of which is presented as follows:
x ˙ 1 = x 2     x ˙ 2 = z      
where x i , i = 1 , 2 refers to system states and z is the control input. The next step involves designing sliding surface S , which is a function related to the system state satisfying the following:
S x 1 , x 2 = 0
The sliding mode control method can be divided into two states: approaching and sliding. The approaching state indicates a situation in which the system has not yet reached the sliding surface (i.e., S 0 ). At this point, the following conditions must be met to ensure that the system reaches the sliding surface within a set period of time:
S S ˙ < K S ,     K > 0
The sliding state indicates a situation in which the system has reached the sliding surface ( S = 0 ). At this point, the system rapidly switches between the upper and lower sides of the sliding surface until converging to the point of equilibrium, at which point the system achieves the required control performance. The conditions required for this process are as follows:
lim s 0 S S ˙ < 0   1   S > 0 ,   S ˙ < 0 2   S < 0 ,   S ˙ > 0

3.2. Proposed Controller Design

Figure 7 is the proposed controller, which simultaneously combines the adaptive radial basis neural network (ARBNN) [19], time delay estimator (TDE) [20], and non-singular fast terminal sliding mode control (NFTSMC) [21]. The ARBNN function is to estimate the payload part, using adaptive control to adjust the weights of neurons in the hidden layer. The TDE is used to deal with system uncertainties, external disturbances, friction, and backlash. Finally, the NFTSMC is used to ensure that the entire system remains stable. The proposed controller is to control the output torque of the motor of each joint based on the predefined desired trajectory in the joint space. The predefined trajectory is based on the dynamic model of the manipulator in Section 2. The force and position relationship could be found in [22].

3.2.1. Adaptive Radial Basis NN (ARBNN)

The dynamic model of the dual-axis manipulator is established as follows:
M 0 q q ¨ + C 0 q , q ˙ q ˙ + G 0 q + Ξ q , q ˙ , q ¨ , t + U p = V U D M
where Υ is set to a constant value of 1.0 for the sake of convenience and U D M refers to the control torque of dual-axis manipulator. We can rewrite Equation (25) in state-space form divided into four parts, as follows:
q ¨ = M 0 1 U D M + M 0 1 C 0 q G 0 + M 0 1 U p + M 0 1 Ξ + Γ U D M
where M 0 1 U D M refers to the control signal, M 0 1 C 0 q G 0 is the known term, M 0 1 U p is the unknown payload term, and M 0 1 Ξ + Γ U D M is the unknown interferences and uncertainties beyond the payload term. Equation (26) can then be replaced by the following:
q ¨ = M 0 1 U D M + F + N + Δ q , q ˙ , U , t
where we define the following for convenience: F = M 0 1 C 0 q G 0 , N = M 0 1 U p , and Δ q , q ˙ , U , t = M 0 1 Ξ + Γ U D M .
The trajectory tracking error is defined as e q q d and the actual and desired trajectories are indicated by q and q d , respectively.
We select an NFTSMC surface that combines the advantages of the conventional sliding film and the terminal sliding mode, allowing the error to converge quickly to zero in a finite time [23]. The NFTSMC surface is as follows:
s = e + c 1 e α + c 2 e ˙ p l = s 1     s 2 T
where c 1 , c 2 R n × n are diagonal positive definite matrixes, and 𝑝 and 𝑙 are positive odd numbers satisfying the relation   1 < p / l   <   2   and   α   >   p / l , such that Equation (28) can be derived as
s ˙ = e ˙ + c 1 α e α 1 e ˙ + c 2 p l e ˙ p l 1 e ¨ = s ˙ 1     s ˙ 2 T
where · · means that A B = a 1 b 1     a 2 b 2   a n b n T if A , B R n × 1 . For the sake of simplicity, we stipulate that O 1 c 1 α e α 1 R n × 1 , O 2 c 2 p l e ˙ p l 1 R n × 1 , which are both positive. In accordance with the design procedures for sliding mode control, let s ˙ = 0 and substitute Equation (27) into Equation (29) to obtain the following:
s ˙ = e ˙ + O 1 e ˙ + O 2 M 0 1 U D M + F + N + Δ q , q ˙ , U , t q ¨ d = 0
Thus, if the control strategy based on NFTSMC satisfies Δ + N η R n × 1 (where η is a small positive vector), then the equivalent control signal U e q of this system can be derived as follows:
U e q = O 2 1 e ˙ O 1 e ˙ + M 0 q ¨ d F = u e q 1   u e q 2 T
At the same time, the uncertainties and disturbances are dealt with using switching control signal U r , which ensures that the system continuously moves along the sliding surface:
U r = η + δ sign s = u r 1       u r 2 T
where δ     R n × 1   is a small vector of positive constants. Thus, the overall control input can be expressed as follows:
U D M 1 = M 0 U e q + M 0 U r
Based on the design of the control input signal, we selected a candidate Lyapunov function as V D M 1 = 1 2 s T s , which, when differentiated, gives us V ˙ D M 1 = s T s ˙ . When we combine Equations (30)–(33), the resulting function satisfies the inequality V ˙ D M 1 O 2 T δ s . Furthermore, both O 2 and δ are positive, while s 0 , such that O 2 T δ s 0 , thereby indicating that the entire system is Lyapunov stable.
The NFTSMC-based controller is able to provide the basic control required for nonlinear systems; however, it cannot ensure stability and robustness in the face of multiple complex disturbances and uncertainties. On the one hand, switching control input U r requires a known maximum upper bound on overall interferences; however, this is seldom satisfied in practical scenarios. Note also that this control strategy depends on the design of the sliding surface parameters c 1 and c 2 , the debugging of which would be a tedious process.
We sought to overcome the shortcomings of this controller by employing an additional controller to deal with other (multiple) forms of disturbance and uncertainty while reducing dependence on the NFTSM control strategy. As mentioned previously, the dual-axis manipulator addressed in this paper includes payload U p , external disturbances d t , system uncertainties Δ M q q ¨ ,     Δ C q , q ˙ q ˙ ,     Δ G q , friction force f q ˙ , and backlash V U . We assigned the payload to an estimator and the other factors to the controller.
In dealing with the payload, ARBNN was used for estimation and compensation tasks. Manipulators commonly encounter problems related to the payload when grasping heavy objects or using end-effectors. We previously defined N = M 0 1 U p as the portion of the payload to be estimated by the ARBNN, which can be expressed as follows:
N = N * + Δ N = W * T μ + Δ N
where N * is the optimal estimate achievable by the neural network, Δ N is the neural network reconstruction error, and W * is the optimal weight. We can express the control input U N N estimate N ^ of the neural network as follows:
U N N = N ^ = W ^ T μ = u N N 1     u N N 2 T
where W ^ is the estimated weight, such that estimation error N ˜ can be derived in a straightforward manner, as follows:
N ˜ = N N ^ = N * + Δ N N ^ = W ˜ T μ + Δ N
W ˜ = W ^ W *

3.2.2. Time Delay Estimator (TDE)

The time delay estimator (TDE) is used to estimate the magnitudes of the other (i.e., non-payload-related) disturbances and uncertainties, which is a type of disturbance estimator. TDE is a simple, highly capable estimation method for the estimation of uncertainty terms using dynamic model equations in accordance with known state variables and their derivatives. In other words, TDE modifies the control signal directly using previous observations of system responses and control inputs rather than by adjusting the controller gain. It follows, therefore, that estimation accuracy is inversely proportional to time delay Δ t .
For the dual-axis manipulator system dealt with in this paper, Equation (27) stipulates that the target term Δ q , q ˙ , U , t can be expressed as follows:
Δ q , q ˙ , U , t = q ¨ M 0 1 U + F + N
In estimating this term, we assume that its value at the current time (i.e., t 0 ) is very close to the value at t 0 Δ t , if time delay Δ t is sufficiently small, such that
Δ q , q ˙ , U , t | t = t 0 Δ q , q ˙ , U , t | t = t 0 Δ t
Thus, we define the output value of the TDE as
Δ ^ q , q ˙ , U , t | t = t 0 Δ q , q ˙ , U , t | t = t 0 Δ t = q ¨ M 0 1 U F N | t = t 0 Δ t
The estimation process generates estimation error ε satisfying the following equation:
Δ q , q ˙ , U , t | t = t 0 = Δ ^ q , q ˙ , U , t | t = t 0 + ε
The output value of the time delay estimator is
U T D E = Δ ^ q , q ˙ , U , t | t = t 0 = u T D E 1     u T D E 2 T
Accordingly, we can use ARBNN ( U N N ) to estimate the payload portion without a known upper bound in conjunction with TDE ( U T D E ) to estimate factors other than payload. The equivalent control ( U e q ) and switching control ( U r ) of NFTSMC can then be combined to output the final control signal as long as the estimation error of ARBNN ( Δ N ) and estimation error of TDE ( ε ) satisfy the following inequalities:
Δ N φ R   n × 1
ε ϱ R n × 1
where φ and ϱ are both positive vectors. Accordingly, the control input of the proposed controller can be redesigned as
U D M 2 = M 0 U e q M 0 U N N M 0 U T D E M 0 U r
where U e q = O 2 1 e ˙ O 1 e ˙ F ϖ s + q ¨ d , ϖ R n × n is a small positive definite diagonal matrix, and U r = δ + φ + ϱ sign s .

3.3. Stability Analysis

We define a candidate Lyapunov function V D M 2 (i.e., a positive definite function) as follows:
V D M 2 = 1 2 s T s + 1 2 t r W ˜ T Λ 1 W ˜
where Λ R n k × n k is a positive definite matrix. Differentiating Equation (46) with respect to time, we obtain
V ˙ D M 2 = s T s ˙ + t r W ˜ T Λ 1 W ˜ ˙
According to Equation (30), the first term in Equation (47) can be rewritten as follows:
s T s ˙ = s T e ˙ + O 1 e ˙ + O 2 e ¨ = s T e ˙ + O 1 e ˙ + O 2 M 0 1 U D M 2 + F + N + Δ q , q ˙ , U , t q ¨ d
Substituting Equation (45) into Equation (48), we obtain
s T s ˙ = s T O 2 N ˜ ϖ s + ε δ + φ + ϱ sign s
Then, substituting Equation (49) into Equation (47), we obtain
V ˙ D M 2 = s T O 2 N ˜ ϖ s + ε δ + φ + ϱ sign s + t r W ˜ T Λ 1 W ˜ ˙ = s T O 2 W ˜ T μ + Δ N ϖ s + ε δ + φ + ϱ sign s + t r W ˜ T Λ 1 W ˜ ˙
Based on Equation (37), we know that W ˜ ˙ = W ^ ˙ W ˙ * = W ^ ˙ , such that we obtain the following:
V ˙ D M 2 = s T O 2 W ˜ T μ + Δ N ϖ s + ε δ + φ + ϱ sign s + t r W ˜ T Λ 1 W ^ ˙
The adaptive law used to adjust the RBNN weights is as follows:
W ^ ˙ = Λ μ ( O 2 s ) T ϕ W ^
where ϕ R n k × n k is a positive definite matrix. Substituting Equation (52) into Equation (51) yields
V ˙ D M 2 = s T O 2 ϖ s Ψ t r W ˜ T ϕ W ^
where Ψ = s T O 2 δ + φ + ε sign s Δ N ε .
Thus, we obtain the following:
W ˜ T ϕ W ^ 1 2 W * T ϕ W * W T ϕ W ˜
Accordingly, we can conclude that
V ˙ D M 2 s T O 2 ϖ s Ψ + 1 2 t r W * T ϕ W * W ˜ T ϕ W ˜ ρ V D M 2 + 1 2 t r W * T ϕ W *
where ρ = min λ m i n O 2 ϖ ,   λ m i n ϕ Λ . Based on Equation (54), we can conclude that the error will eventually be bounded.

4. Simulations

In this section, we outline simulations of the dual-axis manipulator system (implemented in MATLAB/Simulink) aimed at assessing the feasibility of the proposed controller. The simulation results were then compared with those of existing controllers.

4.1. Simulation Settings

Simulations were performed using the dynamic parameters of a dual-axis manipulator [24] (see Table 1) under the following uncertainties and disturbances:
Δ M = 0.2 M ,   Δ C = 0.2 C ,   Δ G = 0.2 G ,   B = 0.3 ,   d 1 = d 2 = 0.01 ,   b = 0.1
Q C = 0                           q ˙ = 0 0.285                   q ˙ > 0 0.285               q ˙ < 0 ,   U p = 0.5 sin t 0.3 sin t
The proposed controller was compared against four existing controllers, including a Proportional–Integration–Derivative (PID) controller, NFTSMC, and NFTSMC in conjunction with TDE. All controller parameters were obtained by trial and error, as follows: PID ( P = 100 ,   I = 50 , D = 100 ); NFTSMC ( p = 9 ,   l = 7 ,   α = 1.5 , c 1 = d i a g 1120    320 , c 2 = d i a g 25   15 ); NFTSMC in conjunction with TDE ( Δ t = 10 4 s); proposed controller ( Λ = I , ϕ = I , ϖ = d i a g 1   1 ).
As shown as Figure 8, the performance of the controllers in the simulations was evaluated while following periodic trajectories ( q d = 0.5 sin t , 0.5 cos t T ) and aperiodic trajectories
q d 1 = 0.5 ,                               t = 0 0.5 sin π t π t ,                 t 0   and   q d 2 = 0.3 ,                               t = 0 0.3 sin π t π t ,                 t 0
where t 0 ,   t f and t f = 10   s .

4.2. Simulation Results

4.2.1. Tracking Periodic Trajectories

In the simulations, we assumed that the robotic manipulator would be performing repetitive tasks; i.e., tracking periodic trajectories. Figure 9a presents the position tracking performance for joint 1. Note that in Figure 9a, the PID controller clearly underperformed the other controllers. The reason is that the PID cannot handle the nonlinear model well. Figure 9b presents a partial enlargement of Figure 9a near t = 6.59   s . We can see that NFTSMC and NFTSMC+TDE poorly dealt with disturbances that occurred at the points of position tracking, where the corresponding disturbances were mainly caused by backlash and friction. In contrast, the proposed controller handled these two disturbances very well; thus, the result clearly outperformed the other controllers. Figure 10a presents the position tracking performance for joint 2. Note that we obtain similar results in the use of different controllers. Figure 10b presents a partial enlargement of Figure 10a near t = 4.92   s , from which we can see that the proposed controller is still relatively the best among the four controllers.
Since this simulation uses the differential method to obtain the velocity value of the system, we consider that this method will amplify the noise. Therefore, a simple first-order low-pass filter is used with a time constant of 0.2 for all controllers’ velocity signals in order to obtain more realistic velocity tracking results. Figure 11 and Figure 12 present the velocity tracking performance for joints 1 and 2. Overall, the proposed system was far less effective in controlling velocity than in controlling position. The reason for this phenomenon is that the velocity signal is delayed after the low-pass filter, but it can filter most of the noise caused by the differential methods so that we can compare the tracking results of each controller more clearly. Nonetheless, we can see in Figure 11 and Figure 12 that the proposed controller outperformed the other controllers in controlling the velocity of joint 1.
In order to show the tracking situation of the four controllers more intuitively, the root mean square error (RMSE) of the four controllers in tracking periodic trajectories is calculated, and the results are shown in Table 2. We can see that the RMSE results of the proposed controller are the best. The values of position tracking RMSE for the proposed controller joints 1 and 2 are 0.080 × 10 4 rad and 0.214 × 10 4 rad, respectively. In terms of velocity tracking, although affected by noise, the proposed controller still has pretty good performances, proving the proposed controller’s effectiveness and superiority in tracking periodic trajectories.

4.2.2. Tracking Aperiodic Trajectories

We also examined the performance in controlling aperiodic trajectories. Figure 13a,b present the position tracking results for joint 1 and Figure 14a,b present the position tracking results for joint 2. In the process of tracking aperiodic trajectories, the position tracking results of joint 1 show that the PID controller performs relatively worst, and the proposed controller performs the best. In addition, the position tracking results of joint 2 show that the proposed controller is also the best. The results validate that the proposed controller can also track aperiodic position trajectories well.
Figure 15a,b present the velocity tracking results for joint 1 and Figure 16a,b present the velocity tracking results for joint 2. The results show that the NFTSMC controller performs best in tracking the velocity of aperiodic trajectories, which is better than the proposed controller and NFTSMC+TDE. The reason for this situation is that the TDE is based on a dynamic model to estimate the disturbance term. This process requires acceleration q ¨ ; however, the differentiation used will further amplify the noise and affect the estimation accuracy of the TDE. Therefore, the accuracy of the velocity tracking is reduced. However, we can see the results of Figure 15b and Figure 16b that under the cases of both using TDE, the velocity tracking accuracy of the proposed controller is better than that of NFTSMC+TDE.
Table 3 shows the RMSE of the proposed controller and other controllers when tracking aperiodic trajectories. Compared with tracking periodic trajectories, we can find that the tracking aperiodic trajectories vary more within a short period of time, and the noise changes are more complicated, affecting the tracking accuracy of each controller. However, the proposed controller’s position tracking is still relatively the best, with 1.2 × 10 3 rad and 3.8 × 10 3 rad, respectively. In terms of velocity tracking, because of the influence of noise in the tracking process of TDE, the RMSE results of the proposed controller and NFTSMC+TDE are not as good as NFTSMC. However, comparing the proposed controller with NFTSMC+TDE, we can find that the proposed controller performs better, thus proving that the proposed controller has better tracking accuracy and stability.
We further discussed the weight variances for the periodic and aperiodic trajectories. Figure 17 presents the weight change; the results can be seen in Figure 17a. Since the periodic trajectory is easier to track, thus, the weight changes are minor. The magnitude of the weight changes is about around 1 × 10 5 , which is presented in Figure 17b. Conversely, the aperiodic trajectory has a more significant vibration in the weight changes of Figure 17a,b because it is not easy to track. The results of weight changes can help us understand why the tracking performances of the two trajectories are different.

5. Conclusions

This study established an enhanced second-order dynamic model of manipulators, which considers payload-related issues as well as model uncertainties, friction, external disturbances, and backlash. We then developed a composite controller based on NFTSMC, which deals with the two classes of interference (payload and other), respectively, using ARBNN and TDE. In simulations, the proposed system outperformed state-of-the-art control systems when applied to periodic as well as aperiodic trajectories.

Author Contributions

Conceptualization, J.-W.W. and W.-S.C.; formal analysis, J.-W.W.; funding acquisition, J.-W.W.; methodology, J.-W.W.; W.-S.C. and C.-C.H.; software, W.-S.C.; validation, J.-W.W. and W.-S.C.; writing—original draft preparation, J.-W.W. and W.-S.C.; writing—review and editing, J.-W.W. and C.-C.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Science and Technology of the Republic of China under Grant MOST-111-2221-E-008-099 and Grant MOST-111-2622-E-008-015.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ming, J.S.; Huan, H.X. Review of development and application of industrial robot technology. J. Mech. Electr. Eng. 2015, 32, 1–13. [Google Scholar]
  2. Zhou, L.M.; Jian, F.; Hai, Z.J. The development and the application of the industrial robot technology. Mach. Build. Autom. 2014, 44, 1–4. [Google Scholar]
  3. Su, Y.X.; Müller, P.C.; Zheng, C. Global Asymptotic Saturated PID Control for Robot Manipulators. IEEE Trans. Control Syst. Technol. 2010, 18, 1280–1288. [Google Scholar] [CrossRef]
  4. Martínez, J.I.M. Robust grbf static neuro controller with switch logic for control of robot manipulators. IEEE Trans. Neural Netw. Learn. Syst. 2012, 23, 1053–1064. [Google Scholar] [CrossRef]
  5. Wai, R.-J.; Muthusamy, R. Fuzzy-Neural-Network Inherited Sliding-Mode Control for Robot Manipulator Including Actuator Dynamics. IEEE Trans. Neural Netw. Learn. Syst. 2013, 24, 274–287. [Google Scholar] [PubMed]
  6. Zhu, Y.K.; Qiao, J.; Zhang, Y.; Guo, L. High-precision trajectory tracking control for space manipulator with neutral uncertainty and dead zone nonlinearity. IEEE Trans. Control. Syst. Technol. 2019, 27, 2254–2262. [Google Scholar] [CrossRef]
  7. He, W.; Chen, Y.H.; Yin, Z. Adaptive Neural Network Control of an Uncertain Robot With Full-State Constraints. IEEE Trans. Cybern. 2016, 46, 620–629. [Google Scholar] [CrossRef]
  8. Van, M.; Mavrovouniotis, M.; Ge, S.S. An Adaptive Backstepping Nonsingular Fast Terminal Sliding Mode Control for Robust Fault-Tolerant Control of Robot Manipulators. IEEE Trans. Syst. Man Cybern. Syst. 2019, 49, 1448–1458. [Google Scholar] [CrossRef]
  9. Liu, Q.; Li, D.Y.; Ge, S.S.; Ji, R.H.; Ouyang, Z.; Tee, K.P. Adaptive bias RBF neural network control for a robotic manipulator. Neurocomputing 2021, 447, 213–223. [Google Scholar] [CrossRef]
  10. Jonker, J.B.; Waiboer, R.R.; Aarts, R. Modelling of Joint Friction in Robotic Manipulators with Gear Transmissions. Multibody Dyn. Comput. Methods Appl. 2007, 4, 221–243. [Google Scholar]
  11. Nordin, M.; Gutmanb, P.O. Controlling mechanical systems with backlash—A survey. Automatica 2002, 38, 1633–1649. [Google Scholar] [CrossRef]
  12. Vörös, J. Modeling and identification of systems with backlash. Automatica 2010, 46, 369–374. [Google Scholar] [CrossRef]
  13. Lai, G.Y.; Liu, Z.; Zhang, Y.; Chen, C.L.P. Adaptive Fuzzy Tracking Control of Nonlinear Systems With Asymmetric Actuator Backlash Based on a New Smooth Inverse. IEEE Trans. Cybern. 2016, 46, 1250–1262. [Google Scholar] [CrossRef] [PubMed]
  14. Dash, C.S.K.; Behera, A.K.; Dehuri, S.; Cho, S.-B. Radial basis function neural networks: A topical state-of-the-art survey. Open Comput. Sci. 2016, 6, 33–63. [Google Scholar] [CrossRef]
  15. Strumiłło, P.; Kamiński, W. Radial Basis Function Neural Networks: Theory and Applications. Neural Netw. Soft Comput. 2003, 19, 107–119. [Google Scholar]
  16. Hung, J.Y.; Gao, W.B.; Hung, J.C. Variable structure control: A survey. IEEE Trans. Ind. Electron. 1993, 40, 2–22. [Google Scholar] [CrossRef]
  17. Wang, C.X.; Tang, J.W.; Jiang, B.P.; Wu, Z.T. Sliding-mode variable structure control for complex automatic systems: A survey. Math. Biosci. Eng. 2022, 19, 2616–2640. [Google Scholar] [CrossRef]
  18. Pisano, A.; Usai, E. Sliding mode control: A survey with applications in math. Math. Comput. Simul. 2011, 81, 954–979. [Google Scholar] [CrossRef]
  19. He, W.; Huang, B.; Dong, Y.T.; Li, Z.J.; Su, C.-Y. Adaptive Neural Network Control for Robotic Manipulators With Unknown Deadzone. IEEE Trans. Cybern. 2018, 48, 2670–2682. [Google Scholar] [CrossRef]
  20. Ahmed, S.; Wang, H.P.; Tian, Y. Adaptive High-Order Terminal Sliding Mode Control Based on Time Delay Estimation for the Robotic Manipulators With Backlash Hysteresis. IEEE Trans. Syst. Man Cybern. Syst. 2021, 51, 1128–1137. [Google Scholar] [CrossRef]
  21. Van, M.; Ge, S.S.; Ren, H.L. Finite time fault tolerant control for robot manipulators using time delay estimation and con-tinuous nonsingular fast terminal sliding mode control. IEEE Trans. Cybern. 2017, 47, 1681–1693. [Google Scholar] [CrossRef] [PubMed]
  22. Hernandez, C.B.; Petit, D.; Alpizar, I.R.; Nishi, T.; Kikuchi, S.; Matsubara, T.; Harada, K. Learning force control for contact-rich manipulation tasks with rigid position-controlled robots IEEE Robot. Autom. Lett. 2020, 5, 5709–5716. [Google Scholar] [CrossRef]
  23. Yang, L.; Yang, J. Nonsingular fast terminal sliding-mode control for nonlinear dynamical systems. Int. J. Robust Nonlinear Control 2010, 21, 1865–1879. [Google Scholar] [CrossRef]
  24. Anđelić, N.; Lorencin, I.; Mrzljak, V.; Carl, Z. Friction Modeling of Robot Manipulator Joints. In Proceedings of the International Conference on Innovative Technologies, Belgrade, Serbia, 11–13 September 2019. [Google Scholar]
Figure 1. Relationship between the velocity of the manipulator and corresponding friction in the joint.
Figure 1. Relationship between the velocity of the manipulator and corresponding friction in the joint.
Applsci 12 09605 g001
Figure 2. Diagram of example payload manipulation process.
Figure 2. Diagram of example payload manipulation process.
Applsci 12 09605 g002
Figure 3. Schematic illustration of backlash associated with the links between input torque and output torque.
Figure 3. Schematic illustration of backlash associated with the links between input torque and output torque.
Applsci 12 09605 g003
Figure 4. Relationship between U and V U .
Figure 4. Relationship between U and V U .
Applsci 12 09605 g004
Figure 5. (a) Structure of RBNN. (b) Corresponding Gaussian function.
Figure 5. (a) Structure of RBNN. (b) Corresponding Gaussian function.
Applsci 12 09605 g005
Figure 6. Structure of sliding surface.
Figure 6. Structure of sliding surface.
Applsci 12 09605 g006
Figure 7. Schematic diagram of proposed controller.
Figure 7. Schematic diagram of proposed controller.
Applsci 12 09605 g007
Figure 8. (a) Position and velocity following the trajectory of joint 1. (b) Position and velocity following the trajectory of joint 2.
Figure 8. (a) Position and velocity following the trajectory of joint 1. (b) Position and velocity following the trajectory of joint 2.
Applsci 12 09605 g008
Figure 9. (a) Position tracking of joint 1. (b) Partial enlargement of (a).
Figure 9. (a) Position tracking of joint 1. (b) Partial enlargement of (a).
Applsci 12 09605 g009
Figure 10. (a) Position tracking of joint 2. (b) Partial enlargement of (a).
Figure 10. (a) Position tracking of joint 2. (b) Partial enlargement of (a).
Applsci 12 09605 g010
Figure 11. (a) Velocity tracking of joint 1. (b) Partial enlargement of (a).
Figure 11. (a) Velocity tracking of joint 1. (b) Partial enlargement of (a).
Applsci 12 09605 g011
Figure 12. (a) Velocity tracking of joint 2. (b) Partial enlargement of (a).
Figure 12. (a) Velocity tracking of joint 2. (b) Partial enlargement of (a).
Applsci 12 09605 g012
Figure 13. (a) Position tracking of joint 1. (b) Partial enlargement of (a).
Figure 13. (a) Position tracking of joint 1. (b) Partial enlargement of (a).
Applsci 12 09605 g013
Figure 14. (a) Position tracking of joint 2. (b) Partial enlargement of (a).
Figure 14. (a) Position tracking of joint 2. (b) Partial enlargement of (a).
Applsci 12 09605 g014
Figure 15. (a) Velocity tracking of joint 1. (b) Partial enlargement of (a).
Figure 15. (a) Velocity tracking of joint 1. (b) Partial enlargement of (a).
Applsci 12 09605 g015
Figure 16. (a) Velocity tracking of joint 2. (b) Partial enlargement of (a).
Figure 16. (a) Velocity tracking of joint 2. (b) Partial enlargement of (a).
Applsci 12 09605 g016
Figure 17. (a) Weight change when tracking periodic and aperiodic trajectories. (b) Partial enlargement of (a).
Figure 17. (a) Weight change when tracking periodic and aperiodic trajectories. (b) Partial enlargement of (a).
Applsci 12 09605 g017
Table 1. Parameter values of basic manipulator model.
Table 1. Parameter values of basic manipulator model.
ParametersDescriptionsValues
g Gravity9.81 m / s 2
m 1 Mass of link 12.0 kg
m 2 Mass of link 20.85 kg
l 1 Length of link 10.35 m
l 2 Length of link 20.31 m
l c 1 Half Length of link 10.175 m
l c 2 Half Length of link 20.155 m
I 1 Inertia of link 1 61.25 × 10 3   kgm 2
I 2 Inertia of link 2 20.42 × 10 3   kgm 2
Table 2. RMSE of tracking results for periodic trajectories.
Table 2. RMSE of tracking results for periodic trajectories.
ControllerJoint 1Joint 2
Position
1 × 10 4 ( rad )
Velocity
1 × 10 4 ( rad / s )
Position
1 × 10 4 ( rad )
Velocity
1 × 10 4 ( rad / s )
PID282.0001143.00059.000730.000
NFTSMC15.000694.00065.000717.000
NFTSMC+TDE2.716676.0003.892693.000
Proposed0.080674.0000.214692.000
Table 3. RMSE of tracking results for aperiodic trajectories.
Table 3. RMSE of tracking results for aperiodic trajectories.
ControllerJoint 1Joint 2
Position
1 × 10 3 ( rad )
Velocity
1 × 10 3 ( rad / s )
Position
1 × 10 3 ( rad )
Velocity
1 × 10 3 ( rad / s )
PID26.600110.2006.40052.700
NFTSMC1.90081.1007.80054.300
NFTSMC+TDE2.200107.90015.50056.800
Proposed1.20082.6003.80053.800
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wu, J.-W.; Cen, W.-S.; Ho, C.-C. Numerical Simulation of Adaptive Radial Basis NN-Based Non-Singular Fast Terminal Sliding Mode Control with Time Delay Estimator for Precise Control of Dual-Axis Manipulator. Appl. Sci. 2022, 12, 9605. https://doi.org/10.3390/app12199605

AMA Style

Wu J-W, Cen W-S, Ho C-C. Numerical Simulation of Adaptive Radial Basis NN-Based Non-Singular Fast Terminal Sliding Mode Control with Time Delay Estimator for Precise Control of Dual-Axis Manipulator. Applied Sciences. 2022; 12(19):9605. https://doi.org/10.3390/app12199605

Chicago/Turabian Style

Wu, Jim-Wei, Wen-Shan Cen, and Cheng-Chang Ho. 2022. "Numerical Simulation of Adaptive Radial Basis NN-Based Non-Singular Fast Terminal Sliding Mode Control with Time Delay Estimator for Precise Control of Dual-Axis Manipulator" Applied Sciences 12, no. 19: 9605. https://doi.org/10.3390/app12199605

APA Style

Wu, J. -W., Cen, W. -S., & Ho, C. -C. (2022). Numerical Simulation of Adaptive Radial Basis NN-Based Non-Singular Fast Terminal Sliding Mode Control with Time Delay Estimator for Precise Control of Dual-Axis Manipulator. Applied Sciences, 12(19), 9605. https://doi.org/10.3390/app12199605

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