Next Article in Journal
On Joan Robinson’s Output Symmetry Theorems with Various Taxations under Third-Degree Price Discrimination: A Generalization
Next Article in Special Issue
Classical Preisach Model Based on Polynomial Approximation and Applied to Micro-Piezoelectric Actuators
Previous Article in Journal
Bayesian Inferential Approaches and Bootstrap for the Reliability and Hazard Rate Functions under Progressive First-Failure Censoring for Coronavirus Data from Asymmetric Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Adaptive Sliding Mode Control Anticipating Proportional Degradation of Actuator Torque in Uncertain Serial Industrial Robots

1
Faculty of Automobile Engineering, Hung Yen University of Technology and Education, Hung Yen 17817, Vietnam
2
School of Electrical and Electronic Engineering, Hanoi University of Science and Technology, Hanoi 11615, Vietnam
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(5), 957; https://doi.org/10.3390/sym14050957
Submission received: 7 April 2022 / Revised: 1 May 2022 / Accepted: 4 May 2022 / Published: 7 May 2022

Abstract

:
The paper focuses on faulty actuator problems in an industrial robot using servomotors, and provides an adaptive sliding mode control law to overcome this circumstance. Because of multifarious reasons, robot actuators can undergo a variety of failures, such as locked or stuck joints, free-swinging joints, and partial or total loss of actuation effectiveness. The robot behavior will become worsen if the system controller has not been designed with adequate faulty tolerance. The proportional degradation of actuator torque at unknown degrees of loss, which is one type of partial loss of actuation effectiveness, is considered in this study to design a suitable controller. The robot model is constructed with uncertain parameters and unknown friction, whereas the controller uses only the approximate parameters. Symmetry and skew-symmetry give important contributions in robot modeling and transformation, as well as in the process of proving the system stability. An adjustable coefficient vector of the proposed controller can adaptively reach the upper bounds of an uncertain parametric vector, which guarantees the criterion of Lyapunov stability. In the numerical simulation stage, the selected industrial robot is a Serpent 1 robot with three degrees of freedom. A quasi-physical model based on MATLAB/Simscape Multibody for the robot is built and used in order to increase the reliability of the simulation performance closer to reality. Simulation results illustrate the efficiency of the proposal control methodology in the presence of the mentioned failure. The controller can still deliver satisfactory responses to the robot system under reasonable levels of actuator torque degradation.

1. Introduction

Robots have long been used in many fields and have been developing to become the core implementation component of Industry 4.0 [1,2,3]. Currently, the growing demand for the reliability of robotic systems in production lines leads researchers and manufactures to keep enhancing the performance of robot in some fault situations. Robotic systems typically respond to a faulty state by shutting down the entire system. The downtime due to the failed parts will reduce productivity and increase costs. In a study about defective industrial robots [4], almost all common failures can be detectable by using a proposed training criteria. Many studies have shown that several certain failure kinds of robotic systems can be overcome by reasonable strategies. One of them is an unexpected decrease in robotic actuator torque, which is investigated in this paper. Generally, some attractive research topics related to drive failure are: (i) detection and location, (ii) analysis and diagnosis, and (iii) failure tolerance control. Many diverse studies related to the topics (i) and (ii) have been published. Robot drive failures are detected in [5,6] by Support-Vector–Machines and a Partial-Least-Squares-based statistical technique, respectively. A non-linear observer in [7] and a failure identification method [8] are developed to recognize faulty states in robotic actuators. In [9], a discrete time framework is established and performed to analyze robot failures by utilizing a decision-making system. Failures in both actuators and sensors of a robot system can be identified by an adaptation-based methodology in [10]. Several analyses regarding manipulator responses in state of a free-moving joint, a locked joint, and actuator torque attenuation are studied in [11,12,13,14], respectively.
Under attraction of the topic (iii), many related investigations have been provided. In [15], an adaptive scheme is designed to counteract the effects of actuator torque degradation on a robot system. In [16], a fault tolerance controller using boundary estimation and an adaptation technique is developed for dealing with partial loss of robot actuators. The result in [17] gives an adaptive control law synthesized with a backstepping scheme for robotic manipulators in case of joint constraints and actuator faults. The sliding mode methodology has been exploited and developed in many various engineering fields, such as angle estimations for salient-pole wound rotor synchronous generators [18], 7-DOF full-car suspension systems [19], position control of an uncertain underwater vehicle [20], fault tolerance control based on an adaptive sliding mode scheme for a quadcopter UAV [21], attitude control of quadrotor UAVs [22], etc. Some diverse interesting applications could be given in more detail in the followings. A novel adaptive PID-type fast terminal sliding mode controller is synthesized for a class of non-linear systems with disturbances [23]. Hyper-chaotic systems can be synchronized by using an integral terminal sliding mode approach [24]. The study in [25] contributes a PID-type terminal sliding surface strategy to enhance the performance of PMSM wind energy conversion systems. For saturated uncertain non-linear systems, a barrier function based adaptive sliding mode scheme is proposed without using the information of disturbance bounds [26]. A sliding mode controller with a leakage-type adaptation law is provided to handle an uncertain Euler–Lagrange system under actuator saturation [27]. The authors in [28] developed a super twisting fractional-order terminal sliding mode method for a type of wind energy conversion systems with non-linearities and disturbances.
Following that strong development trend, the sliding mode methodology is also applied in the robotic control discipline. Robot actuator failures can be rebuilt in a finite time by developing a fast terminal sliding mode observer based fashion [29]. Reference [30] presents a fault tolerance controller using a sliding mode approach to overcome the kind of robot actuator failures: additive time-varying and constant torque. Various studies handle with a mix of robotic joint faults: identification, diagnosis, analysis, and control. Some of them are: [31] for flexible-joint manipulators, [32] for robots subject to free-moving joints, and [33] for locked joints occurring in a specific type of manipulators.
On the other hand, the actuator fault problems in general mechanical systems have been received considerable interest up to now. The study in [34] provides an LMI-based robust adaptive fault-tolerant controller to stabilize a Lipschitz non-linear system. An adaptive controller combined with a sliding mode method is proposed in [35] for uncertain multi-input systems under disturbances and actuation failures. In [36], an overview of direct adaptive compensation control for faults in redundant actuators is given. For a class of uncertain SISO non-linear systems, a partial failed-actuator condition is dealt with a robust adaptive algorithm proposed in [37]. For an unknown MISO non-linear system suffering from actuator faults, some compensation fashions are: adaptive control using a backstepping technique with uncertain degrees of actuator loss [38], mix of a neural-network adaptation strategy and an automatic switching procedure [39], and failure compensation using event-triggering control and an adaptation method with the boundary estimation of failed-actuator parameters [40].
Inspired by the above research, this paper addresses an adaptive sliding mode strategy to anticipate multiple failed actuators for a non-redundant serial industrial robot with uncertain parameters. The investigated kind of failures is unknown proportional degradation of actuator torque. The reaching/switching gains of the proposed sliding mode control method are alterable and regulated by an adaptation law. The effects of actuator torque degradation can be reduced, and system responses are maintained at acceptable qualities under reasonable degrees of loss. The system stability is confirmed by a Lyapunov-based approach. The next parts of this paper are as follows. Section 2 provides the expression of the proportional degradation of the actuator torque in a robot dynamic model. Section 3 contributes to the design of an adaptive sliding mode controller anticipating the mentioned fault. In Section 4, a quasi-physical robot model is constructed as a virtual plant. Section 5 shows the efficiency of the proposed controller through several numerical simulations, and Section 6 gives some important remarks.

2. Proportional Degradation of the Actuator Torque in a Robot Model

Industrial robots considered in this study belong to a class of n-DOF serial manipulators with revolute joints. The matrix form of dynamic equations for that robot kind is:
M ( q ) q ¨ + C ( q , q ˙ ) q ˙ + G ( q ) + τ f = τ a
where q = [ q 1 , q 2 , . . . , q n ] T R n × 1 is the joint angle variable, τ a = [ τ 1 a , τ 2 a , . . . , τ n a ] T R n × 1 is the actual torque produced by actuators, M q R n × n is the general inertia matrix, C q , q ˙ R n × n is the Coriolis/centrifugal matrix, G q R n × 1 is the gravity vector, and τ f R n × 1 is the rotational friction torque. In Equation (1), torque τ a generated by robot actuators suffering the proportional degradation of torque (PDT) can be represented as:
τ a = Λ τ c
where τ c = [ τ 1 c , τ 2 c , . . . , τ n c ] T R n × 1 is the control torque, Λ = diag λ R n × n is the diagonal matrix representing the coefficient degradation of actuator torque, λ = λ 1 , λ 2 , , λ n T , λ i 0 , 1 , and ( i = 1 , 2 , , n ). Actuator i is in the normal operation state if λ i = 1 or in a PDT faulty state if 0 < λ i < 1 . The total loss of actuator torque ( λ i = 0 ) will lead to either the free-swinging joint failure or the locked-joint failure if the fail-safe brake is activated. Therefore, the total torque loss is not considered in this investigation. Since matrix M q is a positive definite and has the symmetry property, it is invertible. Thus, after substituting (2) into (1), Equation (1) can be rewritten in the second-order differential equation as:
q ¨ = M 1 ( q ) [ C ( q , q ˙ ) q ˙ G ( q ) τ f ] + M 1 ( q ) Λ τ c
Compactifying (3) gives:
q ¨ = k + Δ τ c
where:
k : = M 1 ( q ) [ C ( q , q ˙ ) q ˙ G ( q ) τ f ] Δ : = M 1 ( q ) Λ
Because M , C , G , τ f , and Λ are not achieved precisely, k and Δ in (4) and (5) are also not obtained exactly. However, it is possible to interpret k = k ^ + k ˜ and Δ = Δ ^ + Δ ˜ , where k ^ and Δ ^ are the nominal values of k and Δ , respectively. Estimation errors k ˜ and Δ ˜ are assumed to be matched uncertainties. In the following section, under the presence of uncertainties, an adaptive sliding mode control method is designed to make joint variable q follow a given desired joint trajectory q d R n × 1 as well as tracking error e = q d q converge to zero. Besides the symmetry characteristic of matrix M ( q ) , the skew-symmetry property in robot dynamics always exists itself, but its expression depends on the structure of the Coriolis/centrifugal matrix. By exploiting the results provided in [41], the generalized inertia matrix M q can be obtained by:
M ( q ) = i = 1 n m i ( J T i 0 ) T J T i 0 + J R i T I i J R i
and the Coriolis/centrifugal matrix C q , q ˙ satisfying the mentioned skew-symmetry property is formulated as:
C ( q , q ˙ ) = 1 2 M q 1 n q ˙ + M q q ˙ 1 n M q q ˙ 1 n T
The vector of the gravitational term G q is:
G ( q ) = P q T
where m i is the mass of robotic body i, J R i R 3 × n is the rotational Jacobian matrix, J T i 0 R 3 × n is the translational Jacobian matrix, symmetric matrix I i R 3 × 3 is the inertia tensor of body i with respect to the frame attached at the centroid of body i and parallel to the attached-body frame i, 1 n is the n-by-n identity matrix, ⊗ signifies Kronecker product operator, and P is the potential energy of the whole robot body. Hence, the parametric vectors and matrices in robot model (1) can be calculated by (6)–(8).

3. Adaptive Sliding Mode Controller

Usually, the robust ability of a sliding mode controller will be gained by an appropriate choice of reaching/switching gains if the uncertain system parameters are constrained by known certain bounds. However, the boundaries of system uncertainties are mostly unknown in reality. Therefore, an adaptive sliding mode controller (ASMC) could be designed to fine-tune the controller parameters. The sliding variable vector s = s 1 , s 2 , , s n T is selected as:
s = Be + e ˙
where B = diag b is the n-by-n diagonal matrix with b = b 1 , b 2 , , b n T , b i > 0 ( i = 1 , 2 , . . . , n ) is selected such that the dynamics of sliding manifold s i = b i e i + e ˙ i = 0 is stable. Taking the time derivative for both sides of (9) yields:
s ˙ = B e ˙ + q ¨ d q ¨ = B e ˙ + q ¨ d ( k + Δ τ c ) = B e ˙ + q ¨ d ( k ^ + k ˜ + ( Δ ^ + Δ ˜ ) τ c ) = B e ˙ + q ¨ d ( k ^ + Δ ^ τ c ) ( k ˜ + Δ ˜ τ c ) = B e ˙ + q ¨ d ( k ^ + Δ ^ τ c ) w
where w : = k ˜ + Δ ˜ τ c = w 1 , w 2 , , w n T is the vector of uncertainty combination. Coefficient λ i > 0 ; therefore, Λ , Δ , and Δ ^ are all invertible. The adaptive sliding mode controller is designed in the form as follows:
τ s l = Δ ^ 1 ( B e ˙ + q ¨ d k ^ ) τ s w = Δ ^ 1 ( Z ^ + Σ ) sgn ( s ) τ c = τ s l + τ s w = Δ ^ 1 [ B e ˙ + q ¨ d k ^ + ( Z ^ + Σ ) sgn ( s ) ]
where τ s l is the sliding control term, τ s w is the switching control term using adaptive gains, Z ^ = diag z ^ R n × n is the diagonal matrix with z ^ = z ^ 1 , z ^ 2 , , z ^ n T , z ^ i is the adjustable coefficient, Σ = diag σ R n × n , σ = σ 1 , σ 2 , , σ n T , σ i > 0 is the arbitrary additive positive parameter, and sgn s : = sgn s 1 , sgn s 2 , , sgn s n T , sgn . is the signum function. Substituting (11) into (10) gives:
s ˙ = B e ˙ + q ¨ d { k ^ + Δ ^ Δ ^ 1 [ B e ˙ + q ¨ d k ^ + ( Z ^ + Σ ) sgn s ] } w = ( Z ^ + Σ ) sgn s w
Assumption 1.
Vector w is changeable but bounded by an unknown constant vector z = z 1 , z 2 , , z n T such that w i z i , i = 1 , 2 , , n .
Under Assumption 1, although z i is undetermined, the criterion of sliding mode can be guaranteed if an adaptation law can fine-tune z ^ i convergent to z i . The adaptation law is defined as:
z ^ ˙ i = ψ i s i = ψ i s i sgn ( s i )
where ψ i > 0 ( i = 1 , 2 , . . . , n ) is the selected positive factor which specifies the adaptation speed. The larger ψ i , the faster matching convergence. By denoting Z = diag z R n × n , the adaptation error matrix and vector are Z ˜ = Z Z ^ and z ˜ = z z ^ , respectively. In the Lyapunov stability theory, there are some types of symmetric, positive definite parameter matrices often used in a chosen Lyapunov function that supports the proof of the Lyapunov criterion. Let us give a Lyapunov function candidate as:
V = 1 2 s T s + 1 2 z ˜ T Ψ 1 z ˜
where Ψ = diag ψ R n × n are the diagonal matrices with ψ = ψ 1 , ψ 2 , , ψ n T , Ψ 1 = diag ψ 1 1 , ψ 2 1 , , ψ n 1 T > 0 . Taking the time derivative of V in (14) we obtain:
V ˙ = s T s ˙ + z ˜ T Ψ 1 z ˜ ˙
After inserting (12) into (15) we obtain:
V ˙ = s T [ ( Z ^ + Σ ) sgn ( s ) w ] + z ˜ T Ψ 1 z ˜ ˙
Vector z is constant according to Assumption 1. Therefore, z ˜ ˙ = z ^ ˙ , and substituting it into (16) leads to:
V ˙ = s T ( Z ^ + Σ ) sgn ( s ) s T w ( z z ^ ) T Ψ 1 z ^ ˙
From the adaptive law (13), the vector type of the adaptive law is:
z ^ ˙ = Ψ s = Ψ S sgn ( s )
where S = diag s R n × n . Substituting (18) into (17) yields:
V ˙ = s T ( Z ^ + Σ ) sgn ( s ) s T w ( z z ^ ) T Ψ 1 Ψ S sgn ( s ) = s T Z ^ sgn ( s ) s T Σ sgn ( s ) s T w z T S sgn ( s ) + z ^ T S sgn ( s )
Due to matrices S , Z , and Z ^ being diagonal, the commutative property of the matrix multiplication is satisfied. Hence:
z ^ T S = i = 1 n z ^ i s i = s T Z ^
z T S = i = 1 n z i s i = s T Z
Substituting (20) and (21) into (19) gives:
V ˙ = s T Σ sgn ( s ) s T Z sgn ( s ) s T w = i = 1 n σ i s i sgn ( s i ) z i s i sgn ( s i ) s i w i = i = 1 n σ i s i z i s i s i w i i = 1 n σ i s i z i s i + s i w i = i = 1 n σ i z i + w i s i
Applying Assumption 1 to (22) gives:
V ˙ i = 1 n σ i s i 0
Equation (23) expresses that the global stability of the robot system using the proposal adaptive sliding mode controller (11) is guaranteed by the Lyapunov criterion. The convergence to zero of sliding variable s i within finite time is proven. Consequently, tracking error e i converges to zero.

4. Quasi-Physical Model of the Serpent 1 Robot

The effectiveness of the designed ASMC can be verified through simulation by using a conventional mathematical robot model as a virtual plant. However, in order to increase the reliability of the results, a quasi-physical robot model will be used instead of the mathematical robot model. In this section, a Serpent 1 robot is selected as an illustrative object, and its quasi-physical model will be built. The SCARA Serpent 1 robot has three rotary joints using servomotors. The robot end-effector can be moved vertically by a pneumatic cylinder. This study concentrates on the problem of actuator torque loss; therefore, the pneumatic cylinder can be inactive. In fact, this neglect does not affect the position and orientation of the end-effector (EE) projected on the working plane. The initial configuration of the Serpent 1 robot with attached-body frames is depicted in Figure 1, and the D-H parameters are described in Table 1.
In Figure 1, besides the base frame, which is O 0 x 0 y 0 z 0 , another fixed frame O w x w y w z w (called the world frame) is also depicted with a 0 = 0.15 and d 0 = 0.3725 . Based on the three-dimensional CAD assembly files of the Serpent 1 robot depicted in Figure 2, the approximated values of mass m i (kg), body centroid r C i (m) with respect to the attached-body coordinate system, and inertia tensor I i (kg m2) can be obtained by any professional mechanical design software. In this way, the nominal values of robot parameters have been achieved and are shown in (24):
Base : m 0 = 7.645 Link 1 : m 1 = 3 . 529 , r C 1 T = 85.088 , 0 , 69 . 331 × 10 3 , I 1 = 57 . 219 0 20.2 0 88 . 702 0 20.2 0 33 . 877 × 10 3 Link 2 : m 2 = 1.397 , r C 2 T = 80 . 403 , 0 , 3.111 × 10 3 , I 2 = 2 . 220 0 0.381 0 6 . 817 0 0.381 0 5 . 433 × 10 3 Link 3 : m 3 = 1.315 , r C 3 T = 8.101 , 0 , 244 . 765 × 10 3 , I 3 = 48 . 294 0 1.831 0 48 . 480 0 1.831 0 0 . 423 × 10 3
The gravity acceleration is chosen as 9.807 m/s2. Because the robot parameters are not perfectly exact, the controller uses nominal values (24), whereas the virtual robot model (quasi-physical model) is set up with 10% deviation from its nominal values. Figure 3 depicts a quasi-physical model of the Serpent 1 robot based on MATLAB/Simscape Multibody similar to the approach presented in [42]. Every joint in the quasi-physical model is affected by rotational friction with parameters described in Table 2.

5. Simulation Results

The initial robot configuration in joint-space is q 0 = q ˙ 0 = 0 , 0 , 0 T , and in the base frame O 0 x 0 y 0 z 0 is point O = 0.4 , 0 , d 3 T with three orientation vectors ( n = 1 , 0 , 0 T , o = 0 , 1 , 0 T , and a = 0 , 0 , 1 T ). Vectors n , o , and a are the unit vectors of axes x 3 , y 3 , and z 3 with respect to the base frame, respectively. In the base frame, consider two points A and B with attached coordinates 0.36 , 0 , d 3 T and 0 , 0.36 , d 3 T , respectively. The given task of the robot is to move its end-effector from point A to point B along the straight line AB (Figure 4). The motion trajectory is divided into two phases. Firstly, the robot end-effector slides along the x 0 -axis from point O to point A (phase 1), and subsequently moves from point A to point B (phase 2). According to robot forward kinematics, position vector p and three orientation vectors of the end-effector are calculated in (25):
p : = [ p x , p y , p z ] T = [ a 1 c 1 + a 2 c 12 , a 1 s 1 + a 2 s 12 , d 3 ] T n : = [ n x , n y , n z ] T = [ c 123 , s 123 , 0 ] T o : = [ o x , o y , o z ] T = [ s 123 , c 123 , 0 ] T = [ n y , n x , 0 ] T a : = [ a x , a y , a z ] T = [ 0 , 0 , 1 ] T
where s 1 : = sin ( q 1 ) , s 12 : = sin ( q 12 ) = sin ( q 1 + q 2 ) , s 123 : = sin ( q 123 ) = sin ( q 1 + q 2 + q 3 ) , c 1 : = cos ( q 1 ) , c 12 : = cos ( q 12 ) = cos ( q 1 + q 2 ) , and c 123 : = cos q 123 = cos q 1 + q 2 + q 3 .
The requirement for the end-effector orientation during the movement on trajectory AB is that vector n is kept perpendicular to line AB. This implies that the angle q 123 described in Figure 4 (i.e., the angle between vector n and the x 0 -axis) is maintained at π / 4 rad through phase 2. Let us call the angle q 123 the “angle of vector n ”, and the “angle of vector o ” is named for the angle ( π / 2 + q 123 ). In each phase, the speed and acceleration along the desired trajectory in the base frame are mainly converted to velocity v x and acceleration a c c x along the x 0 -axis. The velocity and acceleration components along the y 0 -axis are calculated as v y = v x and a c c y = a c c x , respectively. The z 0 -components v z = 0 and a c c z = 0 , since there is no movement along the z 0 -axis. The motion trajectory along the x 0 -axis is designed by utilizing trapezoidal acceleration profiles (Figure 5 and Figure 6). The design of corresponding orientation also uses that technique with the same time-frame (Figure 7 and Figure 8). Then, the desired position trajectories combined with the corresponding orientation can be converted to the desired joint variable trajectories.
Remark 1.
A sliding surface (9) is selected with larger coefficient b i , then error e i more rapidly converges to zero. Saturation function s a t s i / h i can be implemented instead of the sign function in (11) and (13) for chattering reduction. The greater the value of h i is, the lesser the amplitude of the chattering phenomenon is; however, this will prolong the time to reach the sliding surface. The larger the coefficient ψ i in (13), the faster the convergence of z ^ i to z i ; however, actuators will be required to provide more torque in the beginning of operation. The increase of factor σ i not only improves the convergence speed of sliding variable s i to zero but also raises the risk of the chattering phenomenon.
The value range of the design parameters is quite wide and depends on specific plants. After some simulations have been tested, the controller parameters are set as an illustration as follows: sliding surface (9) with b 1 = b 2 = b 3 = 15 ; saturation function sat s i / h i with h 1 = h 2 = h 3 = 0 . 1 ; adaptive speed parameters in (13) are ψ 1 = ψ 2 = ψ 3 = 60 ; additive positive parameters in (11) are σ 1 = σ 2 = σ 3 = 1 . Simulations are performed in three cases:
Case 1: Parameters in (24) are assumed to be ideal and set in the quasi-physical model of the robot (Figure 3). In this case, rotational friction at joints is not considered in the quasi-physical model. All actuators work in normal operation λ = 1 , 1 , 1 T .
Case 2: Assume that the deviation of parameters in (24) is of 10 % , and consider that the rotational friction components (Table 2) exist in the respective joints. These parameters are set in the quasi-physical model of the robot (Figure 3). Every robot actuator suffers 10 % loss in PDT, corresponding to λ = 0 . 9 , 0 . 9 , 0 . 9 T .
Case 3: The conditions of the quasi-physical robot model are initiated as in Case 2. However, the level of torque loss increases to 50 % corresponding to λ = 0 . 5 , 0 . 5 , 0 . 5 T .
To evaluate the performance qualities of the end-effector position and orientation in the work-space, we judge the magnitude of the position error vector e p = e p x , e p y , e p z T and the absolute error e n between the desired and response angle of vector n that are defined in (26) and (27), respectively:
e p = e p x 2 + e p y 2 + e p z 2 = ( p x d p x ) 2 + ( p y d p y ) 2 + ( p z d p z ) 2
e n = atan 2 ( n y d , n x d ) atan 2 ( n y , n x )
where p d = [ p x d , p y d , p z d ] T is the desired vector of p and n d = [ n x d , n y d , 0 ] T is the desired vector of n . Simulation results for aforementioned cases are shown below.

5.1. Position Trajectory Response Analysis

When the model is correct and all robot actuators work in the normal operating range (Figure 9), the position response of the end-effector gives excellent quality (Figure 10) with a maximum error e p of 1.1 × 1 0 5 m. In the remaining two cases, the end-effector finally follows the desired trajectory (Figure 11 and Figure 12), despite the existence of model parameter deviations, joint friction, and actuator loss in PDT being 10% (Figure 13), 50% (Figure 14) when implementing the ASMC controller. In case 2, the largest position deviation of the end-effector e p is 6.6 × 1 0 4 m from the initial time and 4.6 × 1 0 4 m counts from 1 s. For case 3, these maximum position deviations are 2.8 × 1 0 3 m and 4.2 × 1 0 4 m. The moving trajectories along the x 0 , y 0 , and z 0 axes of the end-effector are p x , p y , and p z reflected in Figure 15, Figure 16 and Figure 17, respectively, for cases 1, 2, and 3. The deviation in the z 0 -axis is zero because the end-effector of the Serpent 1 robot moves only in the O 0 x 0 y 0 plane. The tracking errors e x and e y along the axes x 0 and y 0 are tiny for all three cases. Case 1 gives the best quality, and the other two cases also have small errors within ± 5 × 1 0 4 m after 1 s and the errors are always ensured to converge to zero.

5.2. Translational Velocity and Acceleration Response Analysis

Case 1 (Figure 18) shows that the translational velocity along the x 0 and y 0 axes of the end-effector has met the requirements for the motion stages with tiny errors (below 1 × 1 0 4 m/s). In Cases 2 and 3, the translational velocity still follows the design velocity, but the deviations gradually appear more obvious, especially during acceleration and deceleration (Figure 19 and Figure 20). Moreover, the uncertainties, friction, and actuator torque loss also have certain influences on the uniform motion phase. The y 0 axial velocity largely fluctuates during the first 1 s of acceleration (ranges approximately ± 5 × 1 0 3 m/s case 1, and approximately ± 0 . 02 m/s in case 2). However, these deviations occur only in a short time and are quickly eliminated to zero. The translational acceleration responses are also shown in Figure 21, Figure 22 and Figure 23 for cases 1, 2, and 3, respectively. Actuator loss in torque has a major influence on acceleration. The greater the degree of attenuation, the greater the effect. Similar to velocity, the acceleration also fluctuates during the first 1 s, and overshoot happens at several transitions between the different profiles of the acceleration set (curved and straight profiles). It can be seen that the ASMC controller still shows reliable results in limiting these effects and making the acceleration follow the desired acceleration profile.

5.3. Orientation Responses Analysis

The desired orientation of the end-effector in the base coordinate is represented in Figure 24. The Serpent 1 robot is considered with motion in the horizontal plane; hence, vector a is always parallel and opposite to axis z 0 . The orientation responses of the end-effector in the three cases are shown in Figure 25, Figure 26 and Figure 27, respectively. The projections on plane O 0 x 0 y 0 partly illustrate that the orientation of the end-effector has met the requirements in the process of moving along the trajectory. However, for a more detailed assessment of the orientation response quality, Figure 28, Figure 29 and Figure 30, which reflect the angle of vector n (i.e., q 123 ) and the angle of vector o (i.e., q 123 π / 2 ), are also considered.
In case 1, we can see that angular errors of vectors n , o are ensured to be zero when the robot move uniformly on OA and AB (Figure 28). The greatest error e n is 4 . 1 × 1 0 5 rad only exists on OA when accelerating (in the first 1 s) and decelerating (from time 2 s to time 3 s). Furthermore, the direction requirement of vector n is maintained at q 123 = π / 4 rad.
In case 2 (Figure 29), in about the first 1 s, maximum overshoot of tracking error e n is witnessed at the amplitude of 0.095 rad = 5 . 4 0 and then, this error converges to zero. In the motion along AB (from time 4 s to time 10 s), e n is kept in a vicinity of zero with the greatest deviation of e n is 0.0026 rad = 0 . 1 5 0 and decreases while approaching B.
In case 3, the orientation response is relatively similar to case 2 but with worse control quality due to higher torque loss in actuators. (Figure 30). Vector n takes 1.2 s to track the reference angular trajectory. In the transient time period of 1.2 s, the amplitude of error e n has the greatest value of 0.106 rad = 6 . 0 7 0 . In the motion along AB, error e n slowly varies around the origin with the greatest amplitude of 0.0032 rad = 0 . 1 8 0 and converges to zero while approaching B.

6. Conclusions

This paper proposed a fault-tolerant controller for a class of serial industrial robots with uncertainties and proportional torque loss in actuators. The proposed adaptive sliding mode controller (ASMC) is synthesized to have the fault-tolerant ability and reduce the effect of this situation. The parameters of the control algorithm are chosen based on the Lyapunov criterion to guarantee the global stability of the robot system. The effectiveness of the proposed ASMC is verified by applying on a 3-DOF Serpent 1 robot. Simulation results show that the robot is completely able to track the reference trajectory with infinitesimal errors when a 10 % torque loss exists in actuators. Requirements for positions, velocities, accelerations, and orientations are satisfied. It is undeniable that the response of the system worsens as the decline degree of actuator torque increases. Nevertheless, when the torque loss is up to 50 % , the controller still ensures the stability of the robot system and brings an acceptable response. Theoretically, the proposal control algorithm is proven to be able to handle the study problem, and simulation results have shown its effectiveness. However, when it is applied for a real robot system, there are some arising problems to face, such as adjustment of control parameters to minimize the chattering magnitude, measurement noises, signal processing techniques, effects of sampling time and transient time of inner control loops embedded inside motor-drive boards, etc. Additionally, in this study, the reduction in the chattering phenomenon is just the usage of saturation functions instead of signum functions, and external disturbances have not been taken into account. This will be our focus in the next experimental phase.

Author Contributions

Conceptualization, L.N.T., L.A.V., T.V.T., B.T.T. and T.L.N.; data curation, L.A.V. and T.V.T.; formal analysis, T.V.T. and T.L.N.; funding acquisition, L.A.V. and T.V.T.; investigation, L.N.T.; methodology, L.N.T.; project administration, L.A.V. and B.T.T.; resources, T.V.T.; software, T.L.N.; supervision, B.T.T., L.A.V.; validation, B.T.T., L.N.T.; visualization, L.A.V., T.V.T. and T.L.N.; writing—original draft, L.N.T.; writing—review and editing, L.N.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Ministry of Education and Training of Vietnam under grant number B2020-SKH-01.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This research is supported by the Ministry of Education and Training of Vietnam in project code B2020-SKH-01. The representative organization in charge of the project is Hung Yen University of Technology and Education.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bi, M. Control of Robot Arm Motion Using Trapezoid Fuzzy Two-Degree-of-Freedom PID Algorithm. Symmetry 2020, 12, 665. [Google Scholar] [CrossRef] [Green Version]
  2. Elawady, W.M.; Bouteraa, Y.; Elmogy, A. An Adaptive Second Order Sliding Mode Inverse Kinematics Approach for Serial Kinematic Chain Robot Manipulators. Robotics 2020, 9, 4. [Google Scholar] [CrossRef] [Green Version]
  3. Mourtzis, D.; Angelopoulos, J.; Panopoulos, N. Closed-Loop Robotic Arm Manipulation Based on Mixed Reality. Appl. Sci. 2022, 12, 2972. [Google Scholar] [CrossRef]
  4. Sathish, V.; Ramaswamy, S.; Butail, S. Training data selection criteria for detecting failures in industrial robots. IFAC-PapersOnLine 2016, 49, 385–390. [Google Scholar] [CrossRef]
  5. Caccavale, F.; Cilibrizzi, P.; Pierri, F.; Villani, L. Actuators fault diagnosis for robot manipulators with uncertain model. Control Eng. Pract. 2009, 17, 146–157. [Google Scholar] [CrossRef]
  6. Muradore, R.; Fiorini, P. A PLS-based statistical approach for fault detection and isolation of robotic manipulators. IEEE Trans. Ind. Electron. 2012, 59, 3167–3175. [Google Scholar] [CrossRef]
  7. McIntyre, M.; Dixon, W.; Dawson, D.; Walker, I. Fault identification for robot manipulators. IEEE Trans. Robot. 2005, 21, 1028–1034. [Google Scholar] [CrossRef]
  8. De Luca, A.; Mattone, R. An identification scheme for robot actuator faults. In Proceedings of the 2005 IEEE/RSJ International Conference on Intelligent Robots and Systems, IROS, Edmonton, AB, Canada, 2–6 August 2005; pp. 2613–2617. [Google Scholar] [CrossRef] [Green Version]
  9. Caccavale, F.; Marino, A.; Muscio, G.; Pierri, F. Discrete-time framework for fault diagnosis in robotic manipulators. IEEE Trans. Control. Syst. Technol. 2013, 21, 1858–1873. [Google Scholar] [CrossRef]
  10. Zeng, Y.; Xing, Y.R.; Ma, H.J.; Yang, G.H. Adaptive fault diagnosis for robot manipulators with multiple actuator and sensor faults. In Proceedings of the 27th Chinese Control and Decision Conference (CCDC), Qingdao, China, 23–25 May 2015; pp. 6569–6574. [Google Scholar] [CrossRef]
  11. Domski, W.; Mazur, A. Emergency control of a space 3R manipulator in case of one joint failure. In Proceedings of the 22nd International Conference on Methods and Models in Automation and Robotics (MMAR), Miedzyzdroje, Poland, 28–31 August 2017; pp. 384–389. [Google Scholar] [CrossRef]
  12. Mu, Z.; Han, L.; Xu, W.; Li, B.; Liang, B. Kinematic analysis and fault-tolerant trajectory planning of space manipulator under a single joint failure. Robot. Biomimetics 2016, 3, 16. [Google Scholar] [CrossRef] [Green Version]
  13. Truc, L.N.; Nghia, N.T.; Thanh, N.T.; Thien, N.M.; Nguyen, T.L. Effect of Actuator Torque Degradation on Behavior of a 6-DOF Industrial Robot. Univers. J. Mech. Eng. 2020, 8, 114–128. [Google Scholar] [CrossRef]
  14. Truc, L.N.; Quang, N.P.; Nguyen, H.Q. Impact Analysis of Actuator Torque Degradation on the IRB 120 Robot Performance using Simscape-Based Model. Int. J. Electr. Comput. Eng. (IJECE) 2021, 11, 4850–4864. [Google Scholar] [CrossRef]
  15. Liu, G. Control of robot manipulators with consideration of actuator performance degradation and failures. In Proceedings of the 2001 IEEE International Conference on Robotics & Automation, Seoul, Korea, 21–26 May 2001; pp. 2566–2571. [Google Scholar] [CrossRef]
  16. Lei, R.h.; Chen, L. Adaptive fault-tolerant control based on boundary estimation for space robot under joint actuator faults and uncertain parameters. Def. Technol. 2019, 15, 964–971. [Google Scholar] [CrossRef]
  17. Jin, X. Adaptive finite-time tracking control for joint position constrained robot manipulators with actuator faults. In Proceedings of the 2016 American Control Conference (ACC), Boston, MA, USA, 6–8 July 2016; pp. 6018–6023. [Google Scholar] [CrossRef]
  18. Lopac, N.; Bulic, N.; Vrkic, N. Sliding Mode Observer-Based Load Angle Estimation for Salient-Pole Wound Rotor Synchronous Generators. Energies 2019, 12, 1609. [Google Scholar] [CrossRef] [Green Version]
  19. Aljarbouh, A.; Fayaz, M.; Qureshi, M.S.; Boujoudar, Y. Hybrid Sliding Mode Control of Full-Car Semi-Active Suspension Systems. Symmetry 2021, 13, 2442. [Google Scholar] [CrossRef]
  20. Vu, M.T.; Le, T.H.; Thanh, H.L.N.N.; Huynh, T.T.; Van, M.; Hoang, Q.D.; Do, T.D. Robust Position Control of an Over-actuated Underwater Vehicle under Model Uncertainties and Ocean Current Effects Using Dynamic Sliding Mode Surface and Optimal Allocation Control. Sensors 2021, 21, 747. [Google Scholar] [CrossRef]
  21. Nguyen, N.; Hong, S. Fault Diagnosis and Fault-Tolerant Control Scheme for Quadcopter UAVs with a Total Loss of Actuator. Energies 2019, 12, 1139. [Google Scholar] [CrossRef] [Green Version]
  22. Zheng, B.; Wu, Y.; Li, H.; Chen, Z. Adaptive Sliding Mode Attitude Control of Quadrotor UAVs Based on the Delta Operator Framework. Symmetry 2022, 14, 498. [Google Scholar] [CrossRef]
  23. Goel, A.; Mobayen, S. Adaptive nonsingular proportional–integral–derivative-type terminal sliding mode tracker based on rapid reaching law for nonlinear systems. J. Vib. Control 2021, 27, 2669–2685. [Google Scholar] [CrossRef]
  24. Alattas, K.A.; Mostafaee, J.; Sambas, A.; Alanazi, A.K.; Mobayen, S.; Vu, M.T.; Zhilenkov, A. Nonsingular Integral-Type Dynamic Finite-Time Synchronization for Hyper-Chaotic Systems. Mathematics 2021, 10, 115. [Google Scholar] [CrossRef]
  25. Nasiri, M.; Mobayen, S.; Member, S.; Arzani, A. PID-type terminal sliding mode control for permanent magnet synchronous generator based enhanced wind energy conversion systems. CSEE J. Power Energy Syst. 2021, 1–10. [Google Scholar] [CrossRef]
  26. Shao, K.; Zheng, J.; Tang, R.; Li, X.; Man, Z.; Liang, B. Barrier Function Based Adaptive Sliding Mode Control for Uncertain Systems With Input Saturation. IEEE/ASME Trans. Mechatronics 2022, 1–11. [Google Scholar] [CrossRef]
  27. Shao, K.; Tang, R.; Xu, F.; Wang, X.; Zheng, J. Adaptive sliding mode control for uncertain Euler–Lagrange systems with input saturation. J. Frankl. Inst. 2021, 358, 8356–8376. [Google Scholar] [CrossRef]
  28. Sami, I.; Ullah, S.; Ali, Z.; Ullah, N.; Ro, J.S. A Super Twisting Fractional Order Terminal Sliding Mode Control for DFIG-Based Wind Energy Conversion System. Energies 2020, 13, 2158. [Google Scholar] [CrossRef]
  29. Bing, X.; Yin, S. An fast reconstruction approach for actuator fault in robot manipulators. In Proceedings of the 14th International Workshop on Variable Structure Systems (VSS), Nanjing, China, 1–4 June 2016; pp. 414–419. [Google Scholar] [CrossRef]
  30. Freddi, A.; Longhi, S.; Monteriù, A.; Ortenzi, D.; Proietti Pagnotta, D. Fault Tolerant Control Scheme for Robotic Manipulators Affected by Torque Faults. IFAC-PapersOnLine 2018, 51, 886–893. [Google Scholar] [CrossRef]
  31. Yoo, S.J. Actuator fault detection and adaptive accommodation control of flexible-joint robots. IET Control. Theory Appl. 2012, 6, 1497–1507. [Google Scholar] [CrossRef]
  32. Shin, J.H.; Lee, J.J. Fault detection and robust fault recovery control for robot manipulators with actuator failures. In Proceedings of the 1999 IEEE International Conference on Robotics and Automation, Detroit, MI, USA, 10–15 May 1999; pp. 861–866. [Google Scholar] [CrossRef]
  33. She, Y.; Xu, W.; Su, H.; Liang, B.; Shi, H. Fault-tolerant analysis and control of SSRMS-type manipulators with single-joint failure. Acta Astronaut. 2016, 120, 270–286. [Google Scholar] [CrossRef]
  34. Azmi, H.; Khosrowjerdi, M.J. Robust adaptive fault tolerant control for a class of Lipschitz nonlinear systems with actuator failure and disturbances. Proc. Inst. Mech. Eng. Part I J. Syst. Control. Eng. 2015, 230, 13–22. [Google Scholar] [CrossRef]
  35. Mirkin, B.; Gutman, P.O.; Shtessel, Y. Adaptive continuous control with sliding mode for plants under nonlinear perturbations, external disturbances and actuator failures. In Proceedings of the 11th International Workshop on Variable Structure Systems, Mexico City, Mexico, 26–28 June 2010; pp. 250–255. [Google Scholar] [CrossRef]
  36. Tao, G. Direct adaptive actuator failure compensation control: A tutorial. J. Control. Decis. 2014, 1, 75–101. [Google Scholar] [CrossRef]
  37. Zhang, J.X.; Yang, G.H. Robust Adaptive Fault-Tolerant Control for a Class of Unknown Nonlinear Systems. IEEE Trans. Ind. Electron. 2017, 64, 585–594. [Google Scholar] [CrossRef]
  38. Wang, W.; Wen, C. Adaptive actuator failure compensation control of uncertain nonlinear systems with guaranteed transient performance. Automatica 2010, 46, 2082–2091. [Google Scholar] [CrossRef]
  39. Yang, Q.; Ge, S.S.; Sun, Y. Adaptive actuator fault tolerant control for uncertain nonlinear systems with multiple actuators. Automatica 2015, 60, 92–99. [Google Scholar] [CrossRef]
  40. Xing, L.; Wen, C.; Liu, Z.; Su, H.; Cai, J. Adaptive compensation for actuator failures with event-triggered input. Automatica 2017, 85, 129–136. [Google Scholar] [CrossRef]
  41. Truc, L.N.; Quyen, N.; Quang, N.P. Dynamic model with a new formulation of Coriolis/centrifugal matrix for robot manipulators. J. Comput. Sci. Cybern. 2020, 36, 89–104. [Google Scholar] [CrossRef]
  42. Truc, L.N.; Lam, N.T. Quasi-physical modeling of robot IRB 120 using Simscape Multibody for dynamic and control simulation. Turk. J. Electr. Eng. Comput. Sci. 2020, 28, 1949–1964. [Google Scholar] [CrossRef]
Figure 1. The Serpent 1 robot with all attached frames.
Figure 1. The Serpent 1 robot with all attached frames.
Symmetry 14 00957 g001
Figure 2. The three–dimensional CAD assembly model of the Serpent 1 robot.
Figure 2. The three–dimensional CAD assembly model of the Serpent 1 robot.
Symmetry 14 00957 g002
Figure 3. The quasi–physical model of the Serpent 1 robot based on Simscape Multibody.
Figure 3. The quasi–physical model of the Serpent 1 robot based on Simscape Multibody.
Symmetry 14 00957 g003
Figure 4. Path (a) and orientation (b) planning for the Serpent 1 robot.
Figure 4. Path (a) and orientation (b) planning for the Serpent 1 robot.
Symmetry 14 00957 g004
Figure 5. Phase 1: The desired motion trajectory of the end–effector (vector p ) from O to A.
Figure 5. Phase 1: The desired motion trajectory of the end–effector (vector p ) from O to A.
Symmetry 14 00957 g005
Figure 6. Phase 2: The desired motion trajectory of the end–effector (vector p ) from A to B.
Figure 6. Phase 2: The desired motion trajectory of the end–effector (vector p ) from A to B.
Symmetry 14 00957 g006
Figure 7. The desired orientation of the end–effector (angle of vector n ) in phase 1.
Figure 7. The desired orientation of the end–effector (angle of vector n ) in phase 1.
Symmetry 14 00957 g007
Figure 8. The desired orientation of the end–effector (angle of vector n ) in phase 2.
Figure 8. The desired orientation of the end–effector (angle of vector n ) in phase 2.
Symmetry 14 00957 g008
Figure 9. Control torque (a) and actuator torque (b) under normal operation (case 1).
Figure 9. Control torque (a) and actuator torque (b) under normal operation (case 1).
Symmetry 14 00957 g009
Figure 10. The end–effector (EE) trajectory response when projecting onto (a) O 0 x 0 y 0 , (b) O 0 y 0 z 0 , (c) O 0 z 0 x 0 , and (d) in the work–space (case 1).
Figure 10. The end–effector (EE) trajectory response when projecting onto (a) O 0 x 0 y 0 , (b) O 0 y 0 z 0 , (c) O 0 z 0 x 0 , and (d) in the work–space (case 1).
Symmetry 14 00957 g010
Figure 11. The end–effector (EE) trajectory response when projecting onto (a) O 0 x 0 y 0 , (b) O 0 y 0 z 0 , (c) O 0 z 0 x 0 , and (d) in the work–space (case 2).
Figure 11. The end–effector (EE) trajectory response when projecting onto (a) O 0 x 0 y 0 , (b) O 0 y 0 z 0 , (c) O 0 z 0 x 0 , and (d) in the work–space (case 2).
Symmetry 14 00957 g011
Figure 12. The end–effector (EE) trajectory response when projecting onto (a) O 0 x 0 y 0 , (b) O 0 y 0 z 0 , (c) O 0 z 0 x 0 , and (d) in the work–space (case 3)).
Figure 12. The end–effector (EE) trajectory response when projecting onto (a) O 0 x 0 y 0 , (b) O 0 y 0 z 0 , (c) O 0 z 0 x 0 , and (d) in the work–space (case 3)).
Symmetry 14 00957 g012
Figure 13. Control torque (a) and actuator torque (b) under 10% attenuation of torque (case 2).
Figure 13. Control torque (a) and actuator torque (b) under 10% attenuation of torque (case 2).
Symmetry 14 00957 g013
Figure 14. Control torque (a) and actuator torque (b) under 50% attenuation of torque (case 3).
Figure 14. Control torque (a) and actuator torque (b) under 50% attenuation of torque (case 3).
Symmetry 14 00957 g014
Figure 15. The p x , p y , and p z components of the working point (a), and the corresponding tracking errors (b) (case 1).
Figure 15. The p x , p y , and p z components of the working point (a), and the corresponding tracking errors (b) (case 1).
Symmetry 14 00957 g015
Figure 16. The p x , p y , and p z components of the working point (a), and the corresponding tracking errors (b) (case 2).
Figure 16. The p x , p y , and p z components of the working point (a), and the corresponding tracking errors (b) (case 2).
Symmetry 14 00957 g016
Figure 17. The p x , p y , and p z components of the working point (a), and the corresponding tracking errors (b) (case 3).
Figure 17. The p x , p y , and p z components of the working point (a), and the corresponding tracking errors (b) (case 3).
Symmetry 14 00957 g017
Figure 18. Speed responses along the x 0 , y 0 , and z 0 axes of the end–effector (a) and the corresponding tracking errors (b) (case 1).
Figure 18. Speed responses along the x 0 , y 0 , and z 0 axes of the end–effector (a) and the corresponding tracking errors (b) (case 1).
Symmetry 14 00957 g018
Figure 19. Speed responses along the x 0 , y 0 , and z 0 axes of the end–effector (a) and the corresponding tracking errors (b) (case 2).
Figure 19. Speed responses along the x 0 , y 0 , and z 0 axes of the end–effector (a) and the corresponding tracking errors (b) (case 2).
Symmetry 14 00957 g019
Figure 20. Speed responses along the x 0 , y 0 , and z 0 axes of the end–effector (a) and the corresponding tracking errors (b) (case 3).
Figure 20. Speed responses along the x 0 , y 0 , and z 0 axes of the end–effector (a) and the corresponding tracking errors (b) (case 3).
Symmetry 14 00957 g020
Figure 21. Acceleration responses along the x 0 , y 0 , and z 0 axes of the end–effector (a) and the corresponding tracking errors (b) (case 1).
Figure 21. Acceleration responses along the x 0 , y 0 , and z 0 axes of the end–effector (a) and the corresponding tracking errors (b) (case 1).
Symmetry 14 00957 g021
Figure 22. Acceleration responses along the x 0 , y 0 , and z 0 axes of the end–effector (a) and the corresponding tracking errors (b) (case 2).
Figure 22. Acceleration responses along the x 0 , y 0 , and z 0 axes of the end–effector (a) and the corresponding tracking errors (b) (case 2).
Symmetry 14 00957 g022
Figure 23. Acceleration responses along the x 0 , y 0 , and z 0 axes of the end–effector (a) and the corresponding tracking errors (b) (case 3).
Figure 23. Acceleration responses along the x 0 , y 0 , and z 0 axes of the end–effector (a) and the corresponding tracking errors (b) (case 3).
Symmetry 14 00957 g023
Figure 24. Desired directional vectors of the end–effector.
Figure 24. Desired directional vectors of the end–effector.
Symmetry 14 00957 g024
Figure 25. Directional vector n (a) and vector o (b) of the end–effector when the robot moves along the trajectory in plane O 0 x 0 y 0 (case 1).
Figure 25. Directional vector n (a) and vector o (b) of the end–effector when the robot moves along the trajectory in plane O 0 x 0 y 0 (case 1).
Symmetry 14 00957 g025
Figure 26. Directional vector n (a) and vector o (b) of the end–effector when the robot moves along the trajectory in plane O 0 x 0 y 0 (case 2).
Figure 26. Directional vector n (a) and vector o (b) of the end–effector when the robot moves along the trajectory in plane O 0 x 0 y 0 (case 2).
Symmetry 14 00957 g026
Figure 27. Directional vector n (a) and vector o (b) of the end–effector when the robot moves along the trajectory in plane O 0 x 0 y 0 (case 3).
Figure 27. Directional vector n (a) and vector o (b) of the end–effector when the robot moves along the trajectory in plane O 0 x 0 y 0 (case 3).
Symmetry 14 00957 g027
Figure 28. Angular responses of directional vectors n , o of the end–effector (a) and the respective tracking errors (b) (case 1).
Figure 28. Angular responses of directional vectors n , o of the end–effector (a) and the respective tracking errors (b) (case 1).
Symmetry 14 00957 g028
Figure 29. Angular responses of directional vectors n , o of the end–effector (a) and the respective tracking errors (b) (case 2).
Figure 29. Angular responses of directional vectors n , o of the end–effector (a) and the respective tracking errors (b) (case 2).
Symmetry 14 00957 g029
Figure 30. Angular responses of directional vectors n , o of the end–effector (a) and the respective tracking errors (b) (case 3).
Figure 30. Angular responses of directional vectors n , o of the end–effector (a) and the respective tracking errors (b) (case 3).
Symmetry 14 00957 g030
Table 1. The D-H parameters of the Serpent 1 robot.
Table 1. The D-H parameters of the Serpent 1 robot.
Joint i θ i ( rad ) d i ( m ) a i ( m ) α i ( rad )
1 q 1 d 1 = 0 a 1 = 0.25 α 1 = 0
2 q 2 d 2 = 0 a 2 = 0.15 α 2 = 0
3 q 3 d 3 = 0 a 3 = 0.070 α 3 = + π
Table 2. Rotational friction parameters (Reprinted/adapted with permission from [42]. 2020, Le Ngoc Truc. See [42] for detail descriptions of notations).
Table 2. Rotational friction parameters (Reprinted/adapted with permission from [42]. 2020, Le Ngoc Truc. See [42] for detail descriptions of notations).
Joint i τ brki ( Nm ) ω brki ( rad / s ) τ Ci ( Nm ) k vi ( Nm / ( rad / s ) )
1 30 × 10 3 0.01 29 × 10 3 0.001
2 23 × 10 3 0.01 20 × 10 3 0.001
3 10 × 10 3 0.01 9 × 10 3 0.001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Truc, L.N.; Vu, L.A.; Thoan, T.V.; Thanh, B.T.; Nguyen, T.L. Adaptive Sliding Mode Control Anticipating Proportional Degradation of Actuator Torque in Uncertain Serial Industrial Robots. Symmetry 2022, 14, 957. https://doi.org/10.3390/sym14050957

AMA Style

Truc LN, Vu LA, Thoan TV, Thanh BT, Nguyen TL. Adaptive Sliding Mode Control Anticipating Proportional Degradation of Actuator Torque in Uncertain Serial Industrial Robots. Symmetry. 2022; 14(5):957. https://doi.org/10.3390/sym14050957

Chicago/Turabian Style

Truc, Le Ngoc, Le Anh Vu, Tran Van Thoan, Bui Trung Thanh, and Tung Lam Nguyen. 2022. "Adaptive Sliding Mode Control Anticipating Proportional Degradation of Actuator Torque in Uncertain Serial Industrial Robots" Symmetry 14, no. 5: 957. https://doi.org/10.3390/sym14050957

APA Style

Truc, L. N., Vu, L. A., Thoan, T. V., Thanh, B. T., & Nguyen, T. L. (2022). Adaptive Sliding Mode Control Anticipating Proportional Degradation of Actuator Torque in Uncertain Serial Industrial Robots. Symmetry, 14(5), 957. https://doi.org/10.3390/sym14050957

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