Finite-time control strategy for swarm planar underactuated robots via motion planning and intelligent algorithm

Abstract

The swarm planar underactuated robots have the ability to organize each robot to complete task in a finite time and the characteristics of ignoring gravity, energy saving, light weight, and so on. In this paper, we propose control strategy for such robot via motion planning and intelligent algorithm. First, we establish a unified dynamic model and analysis its underactuated characteristics. In order to enable all links to reach the target position smoothly in a finite time, a suitable trajectory is planned and parameters are optimized by the differential evolution algorithm (DEA). Then, the controllers are designed for each planar underactuated robot. Finally, the simulation results show that the proposed strategy is effective.

Keywords

Finite-time control swarm planar underactuated robots motion planning differential evolution algorithm

Introduction

In recent years, the research on swarm robots inspired by biological colonies behavior has become a hot topic in the control field. Each robot is responsible for different tasks, but as a member of the colony, it obeys the self-organizing cooperative control and completes complex tasks together in finite time. At present, the swarm robots have been widely used in exploration,¹ transportation,² medical rehabilitation,³ and so on. When the joint of robot encounters failure, the original robot become underactuated robot.

The underactuated robot^4,5 is an important research direction in the field of nonlinear control, whose the degree of freedoms are more than the number of inputs.^6,7 The nonlinear coupling relationship makes its control full of challenges. The surface vehicle,⁸ tower crane,⁹ and planar underactuated robot^10–12 etc. are all high-frequency research objects in underactuated system.

The planar underactuated robot does not consider gravity, and any point in the plane is equilibrium. But the linear approximate model of the arbitrary equilibrium point is uncontrollable, which makes the robot more difficult to control.¹³ According to the number of links and the position of passive joint in robot, the underactuated characteristics are different.¹⁴ The planar PA robot has angle constraints among links, and it belongs to holonomic system. Wang et al.¹⁵ fully consider the coupling relationship and use Lyapunov function to design the controller to realize the rapid angle control of each link. The planar AP robot has angular acceleration constraints, and it belongs to second-order nonholonomic system. Based on Fourier transform, Wu et al.¹⁶ design a time-dependent controller for the planar AP robot corresponding to the frequency of harmonic term to realize the control target of moving from the initial to target position.

In order to meet the diversified functions and complex operation, the research on the underactuated control strategy of multi-link robot is necessary.^17–19 The planar PAA robot belongs to first-order nonholonomic system with angular velocity constraints. Zhang et al.²⁰ adapt a two-stages control method to solve the problem that the asymptotic stability of such robot at the target position. In stage 1, the controller is designed to stabilize one active link at the target angle and the others at the transition angles; In stage 2, the other links will be moved to the target states simultaneously due to the angle constraint. The passive middle or last joint of the planar underactuated robot belong to the second-order nonholonomic system. Huang et al.²¹ propose a bidirectional trajectory planning and tracking control method based on time scale for planar APA robot, calculate the time scale factor by using intelligent algorithm, and realize the position control through a single trajectory tracking controller.

For the planar AAP robot, the position control of its first active link is realized firstly, and then the fuzzy control method is used to realize the stability of other links.²² Huang et al.²³ propose an unified method for the planar three-link robot. Consider the underactuated states and constraints, he designs trajectory for active link that can be optimized by intelligent algorithm and uses the designed tracking controller to realize the control objectives. Huang et al.²⁴ reduce the planar multi-link robot with a last passive joint to the planar AP robot, and make the state variables converge rapidly based on nilpotent approximation and iterative theory. De Luca and Oriolo²⁵ transform the planar $A^{n} P$ (Active-⋯-Active-Passive) robot into a second-order canonical chained form, calculate the impact coordinates and obtain expected trajectory of the passive link, and realize the stability control of the active link by designing a feedback controller.

The last joint of the planar swarm robot is very fragile because it is used frequently. Once damaged, it will cause huge property losses. Meanwhile, there is no effective control strategy for the swarm robot with a last passive joint in a finite time. Therefore, based on motion planning and intelligent algorithm, we study the control strategy for such system. First, an unified dynamic model is established for the planar $A^{n} P$ robot, and its underactuated constraints are analyzed. Then, a trajectory with adjustable parameters is designed for each active link. Meanwhile, the DEA is used to optimize the parameters so that all links are moved and maintained to the target states in a finite time. Finally, the sliding mode variable structure controllers are designed to achieve the objective. Simulation results prove the effectiveness of the proposed strategy.

Preparations

Model

The unified structure of the planar $A^{n} P$ robot is shown in Figure 1. For the $r$ -th $(r = 1, \dots, n + 1)$ link, $m_{r}$ , $q_{r}$ , $L_{r}$ and $J_{r}$ are the mass, angle, length, and the moment of inertia of the robot, respectively. $l_{r}$ is the distance from $r$ -th joint to its center of the mass, and $τ_{r}$ denote the inputs applied to active joints.

Figure 1.

The unified structure of the swarm planar $A^{n} P$ robot.

The dynamic equations is

M (q) \overset{\cdot\cdot}{q} + H (q, \overset{\cdot}{q}) = τ

(1)

where the vector of $q$ , $\overset{\cdot}{q}$ , $\overset{\cdot\cdot}{q}$ , $τ$ and the matrix of $M (q)$ , $H (q, \overset{\cdot}{q})$ are detailed expression in Yoshikawa et al.²⁶

Underactuated constraint analysis

The underactuated portion in (1) is

M_{(n + 1) 1}^{k} {\overset{\cdot\cdot}{q}}^{k}_{1} + \dots + M_{(n + 1) (n + 1)}^{k} {\overset{\cdot\cdot}{q}}^{k}_{(n + 1)} + H_{n + 1}^{k} = 0

(2)

where $k = 1, 2, \dots, 6$ .

According to (2), the underactuated constraints among the underactuated and active link are obtained.

{\begin{matrix} {\overset{\cdot}{q}}^{k}_{(n + 1)} = \\ - \int_{0}^{T} (\frac{\sum_{i = 1}^{n} M_{(n + 1) i}^{k} {\overset{\cdot\cdot}{q}}^{k}_{i} + H_{(n + 1)}^{k}}{M_{(n + 1) (n + 1)}^{k}}) d t - {\overset{\cdot}{q}}^{k}_{(n + 1) 0} \\ q_{(n + 1)}^{k} = \int_{0}^{T} {\overset{\cdot}{q}}^{k}_{n + 1} d t - q_{(n + 1) 0}^{k} \end{matrix}

(3)

where $q_{(n + 1) 0}^{k}$ and ${\overset{\cdot}{q}}_{(n + 1) 0}^{k}$ are the initial states of the underactuated link, respectively.

Trajectory planning

In order to quickly realized the stability control objective of each robot in a finite time,²⁷ we fully consider the initial and target angles of the system and the termination time $t_{f}$ of the following designed trajectory $T^{k} (t)$ for the active link.

\begin{matrix} T^{k} (t) = q_{n 0}^{k} + (q_{nd}^{k} - q_{n 0}^{k}) (\frac{t}{t_{f}} - \frac{1}{2 π} \sin (\frac{2 π t}{t_{f}})) \\ + A_{n 1}^{k} sech (α_{1}^{k}) + A_{n 2}^{k} sech (α_{2}^{k}) \end{matrix}

(4)

where

{\begin{matrix} α_{1}^{k} = \frac{9 (2 t - t_{m}^{k})}{t_{m}^{k}} \\ α_{2}^{k} = \frac{9 (2 t - t_{f} - t_{m}^{k})}{t_{f} - t_{m}^{k}} \end{matrix}

(5)

The derivatives of $T^{k} (t)$ is

\begin{matrix} {\overset{\cdot}{T}}^{k} (t) = \frac{q_{nd}^{k} - q_{n 0}^{k}}{t_{f}} (1 - \cos (\frac{2 π t}{t_{f}})) \\ - \frac{18 A_{n 1}^{k} sech (α_{1}^{k}) \tanh (α_{1}^{k})}{t_{m}^{k}} \\ - \frac{18 A_{n 2}^{k} sech (α_{2}^{k}) \tanh (α_{2}^{k})}{t_{f} - t_{m}^{k}} \end{matrix}

(6)

\begin{matrix} {\overset{\cdot\cdot}{T}}^{k} (t) = \frac{q_{nd}^{k} - q_{n 0}^{k}}{t_{f}^{2}} (2 π \sin (\frac{2 π t}{t_{f}})) \\ + \frac{324 A_{n 1}^{k} sech (α_{1}^{k}) {(\tanh (α_{1}^{k}))}^{2}}{{(t_{m}^{k})}^{2}} \\ - \frac{324 A_{n 1}^{k} {(sech (α_{1}^{k}))}^{3}}{{(t_{m}^{k})}^{2}} \\ + \frac{324 A_{n 2}^{k} sech (α_{2}^{k}) {(\tanh (α_{2}^{k}))}^{2}}{{(t_{f} - t_{m}^{k})}^{2}} \\ - \frac{324 A_{n 2}^{k} {(sech (α_{2}^{k}))}^{3}}{{(t_{f} - t_{m}^{k})}^{2}} \end{matrix}

(7)

when $t = t_{f}$ , all active links reach the target states, but the underactuated link will rotate freely due to the nonholonomic constraints. In order to indirectly stabilize the passive part to the objective by controlling the active links, the parameters $A_{n 1}^{k}$ , $A_{n 2}^{k}$ , and $t_{m}^{k}$ need to be optimized by the DEA.²⁸

The evaluation function is

h = | q_{(n + 1)}^{k} (t_{f}^{k}) - q_{(n + 1) d}^{k} | + | {\overset{\cdot}{q}}_{(n + 1)}^{k} (t_{f}^{k}) |

(8)

The optimization procedure is as follows:

① The initialization population is randomly generated, and the parameters of the trajectory are assumed to be ${\tilde{A}}_{n 1}^{k}$ , ${\tilde{A}}_{n 2}^{k}$ , and ${\tilde{t}}_{m}^{k}$ .

② According to (4) to (7) and ${\tilde{A}}_{n 1}^{k}$ , ${\tilde{A}}_{n 2}^{k}$ , ${\tilde{t}}_{m}^{k}$ calculate $q_{n + 1} (t_{f})$ and ${\overset{\cdot}{q}}_{n + 1} (t_{f})$ , and substitute them into (8).

③ When $h \leq ε_{1}$ ; $ε_{1}$ is a small positive constant, $A_{n 1}^{k} = {\tilde{A}}_{n 1}^{k}$ , $A_{n 2}^{k} = {\tilde{A}}_{n 2}^{k}$ , $t_{m}^{k} = {\tilde{t}}_{m}^{k}$ . Otherwise, the program is to ④.

④ After mutation, crossover and selection, update ${\tilde{A}}_{n 1}^{k}$ , ${\tilde{A}}_{n 2}^{k}$ , ${\tilde{t}}_{m}^{k}$ , go to ②.

The underactuated link can reach the objectives at the same time when the active links return to the target at $t = t_{f}$ .

Controller design

The sliding mode variable structure controller is designed for active links tracking.

Let $x = [x_{1}, \dots, x_{(n + 1)}, x_{(n + 2)}, \dots, x_{(2 n + 2)}]^{T}$ = $[q_{1}, \dots, q_{(n + 1)}, {\overset{\cdot}{q}}_{1}, \dots, {\overset{\cdot}{q}}_{(n + 1)}]^{T}$ , the state space equation of the planar $A^{n} P$ robot is

\overset{\cdot}{x} = f (x) + g (x) τ

(9)

where

f (x) = [f_{1} f_{2} \dots f_{2 n + 2}] = - M^{- 1} (q) H (q, \overset{\cdot}{q})

(10)

and

g (x) = [\begin{matrix} g_{11} (x) & \dots & g_{1 (n + 1)} (x) \\ g_{21} (x) & \dots & g_{2 (n + 1)} (x) \\ ⋮ & ⋱ & ⋮ \\ g_{(n + 1) 1} (x) & \dots & g_{(n + 1) (n + 1)} (x) \end{matrix}] = M^{- 1} (q)

(11)

The sliding mode surface is

S_{i}^{k} = μ_{k} e_{i}^{k} + {\overset{\cdot}{e}}_{i}^{k} (i = 1 \dots n)

(12)

and

e_{i}^{k} = q_{i}^{k} - T^{k} (t), {\overset{\cdot}{e}}_{i}^{k} = {\overset{\cdot}{q}}_{i}^{k} - {\overset{\cdot}{T}}^{k} (t)

(13)

The derivative of $S_{i}^{k}$ is

{\overset{\cdot}{S}}^{k}_{i} = μ_{k} {\overset{\cdot}{e}}^{k}_{i} + {\overset{\cdot\cdot}{e}}^{k}_{i} = μ_{k} {\overset{\cdot}{e}}^{k}_{i} + f_{i} + \sum_{j = 1}^{n} g_{ij} τ_{j}^{k} - {\overset{\cdot\cdot}{T}}^{k} (t)

(14)

Let

{\overset{\cdot}{S}}_{i}^{k} = - φ_{k} S_{i}^{k} - ε_{k} sgn (S_{i}^{k})

(15)

where $[μ_{k}, φ_{k}, ε_{k}] > [0, 0, 0]$ , and $sgn (S_{i}^{k})$ denote the sign of $S_{i}^{k}$ .

The controllers are designed as

\begin{matrix} τ_{i}^{k} = (- φ_{k} S_{i}^{k} - ε_{k} sgn (S_{i}^{k}) - μ_{k} {\overset{\cdot}{e}}^{k}_{i} - f_{i} \\ - \sum_{j = 1, i \neq j}^{n} g_{ij} τ_{j}^{k} + {\overset{\cdot\cdot}{T}}^{k} (t)) g_{ij, i = j}^{- 1} \end{matrix}

(16)

The Lyapunov function is constructed as:

V = \frac{1}{2} (S_{i}^{k})^{2}

(17)

then

\begin{matrix} \overset{\cdot}{V} = S_{i}^{k} {\overset{\cdot}{S}}^{k}_{i} = - φ_{k} {(S_{i}^{k})}^{2} - ε_{k} | S_{i}^{k} | \leq 0 \end{matrix}

(18)

When $\overset{\cdot}{V} \equiv 0$ , $S_{i}^{k} \equiv 0$ . Thus, according to the LaSalle’s theorem,²⁹ $e_{i}^{k} \to 0$ and ${\overset{\cdot}{e}}^{k}_{i} \to 0$ . the controller (16) realized the control objective of the planar $A^{n} P$ robot.

Simulations

The feasibility and superiority of the proposed strategy in the stability control of the swarm planar $A^{n} P$ $(n = 2, 3)$ robot in finite time are verified by Matlab simulation. The parameters of the DEA and controller (16) are $p_{c} = 0.7$ , $p_{m} = 0.3$ , $Maxgen = 200$ , $e_{i}^{k} = 0.005$ , $μ_{j} = 1.2$ , $φ_{j} = 1$ , and $ε_{j} = 0.3$ .

The structure parameters of a single planar AP robot¹⁶ in the swarm planar underactuated robot are shown in Table 1.

Table 1.

Structure parameters of the planar AP robot.

link $r$	$m_{r}$ (kg)	$L_{r}$ (m)	$l_{r}$ (m)	$J_{r}$ (kg · m ²)
1	1.0	1.0	0.5	0.3333
2	1.0	1.0	0.5	0.3333

We choose the termination time is $t_{f} = 15$ s, the target angle is same as Wu et al.¹⁶ and the different initial states of each link are shown in Table 2.

Table 2.

The initial states of the swarm planar AP robot.

No. $k$	$q_{10}^{k}$ (rad)	$q_{20}^{k}$ (rad)
1	0	0
2	0.1	0.2
3	−0.1	−0.3
4	0.2	−0.1
5	−0.4	0.1
6	0.3	0.1

The parameters of the designed trajectory are optimized by the DEA, and the results are shown in Table 3.

Table 3.

Trajectory parameter optimization results.

No. $k$	$A_{11}^{k}$	$A_{12}^{k}$	$t_{m}^{k}$
1	−0.2723	0.0578	10.0000
2	−1.9322	−0.2704	9.0033
3	1.4172	−0.1777	8.7240
4	−0.3537	1.3247	4.7127
5	−0.2063	−2.2304	3.1014
6	0.2803	−0.7620	6.8077

The simulation results (Figure 2) show that the swarm planar AP robot realize the stable control objective. The angles and angular velocities converge to their target values smoothly. Compared with Wu et al.¹⁶ the control torque of our proposed strategy is smaller.

Figure 2.

Simulation results of swarm planar AP robot: (a) angles of the first link, (b) angles of the second link, (c) angular velocities of the first link, (d) angular velocities of the second link, and (e) control torques .

The structure parameters of a single planar AAP robot²² are shown in Table 4.

Table 4.

Structure parameters of the planar AAP robot.

link $r$	$m_{r}$ (kg)	$L_{r}$ (m)	$l_{r}$ (m)	$J_{r}$ (kg · m ²)
1	0.3	0.3	0.15	0.0023
2	0.3	0.3	0.15	0.0023
3	0.3	0.3	0.15	0.0023

We choose the termination time is $t_{f} = 20$ s, the target angle is same as Liu et al.²² and the different initial states of each link are shown in Table 5.

Table 5.

The initial states of the swarm planar AAP robot.

NO. $k$	$q_{10}^{k}$	$q_{20}^{k}$	$q_{30}^{k}$
1	0	0	0
2	0.2	0.4	0.1
3	0.1	0.3	−0.4
4	0.3	−0.2	0.1
5	−0.4	0.2	0.3
6	−0.3	−0.2	0.1

The optimization results of trajectory parameters are shown in Table 6.

Table 6.

Trajectory parameter optimization results.

NO. $k$	$A_{11}^{k}$	$A_{12}^{k}$	$A_{21}^{k}$	$A_{22}^{k}$	$t_{m}^{k}$
1	−0.1460	0.1540	−0.1028	0.6697	8.4052
2	−0.9901	0.0296	1.7867	1.4031	10.4286
3	−0.3888	−0.0563	1.1684	0.6518	13.5295
4	−0.7338	−0.0704	1.3556	0.6821	12.9062
5	−0.0548	−0.0051	−0.4186	0.7995	9.4230
6	−0.0431	0.8455	−0.2331	−0.3947	5.2324

The simulation results (see Figures 3 –5) show that the swarm planar AAP robot can also achieve stable control in finite time, all links gradually converges to the target states along a continuous and smooth trajectory, and the control torques are maintained within $\pm 0.30$ N· m, which is relatively small. Compared with Liu et al.²² the proposed strategy does not need piecewise control.

Figure 3.

Angles of each link: (a) angles of the first link, (b) angles of the second link, and (c) angles of the third link.

Figure 4.

Angular velocities of each link: (a) angular velocities of the first link, (b) angular velocities of the second link, and (c) angular velocities of the third link.

Figure 5.

Control torques of each active link: (a) control torques of the first active link, (b)control torques of the second active link.

Conclusions

This paper presents an advanced trajectory planning and tracking control strategy for the swarm planar $A^{n} P$ robot, and analyzes two typical robots through simulation. Now the main conclusions are as follows.

A unified dynamic model for swarm planar $A^{n} P$ robot is established, and the underactuated constraints are analyzed based on the above model.

The suitable trajectory with parameters for active links are planned and the parameters are optimized by the DEA.

The controllers are designed for active links to track the optimized trajectory.

The control strategy can be extended to the same type of typical systems with multiple links. Therefore, we will further investigate this kind of control strategy to the other swarm underactuated robot with different structures in the future research.

Footnotes

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work are supported by the Science and Technology Research Project of Hubei Provincial Education Department (D20211506), and the Scientific Research Foundation of Wuhan Institute of Technology (K2021027), and the Hubei Key Laboratory of Intelligent Robot (Wuhan Institute of Technology) (HBIRL202105), and the Graduate Innovative Fund of Wuhan Institute of Technology (CX2021074, CX2021103).

ORCID iD

Zixin Huang

References

St-Onge

Kaufmann

Panerati

, et al. Planetary exploration with robot teams: Implementing higher autonomy with swarm intelligence. IEEE Robot Autom Mag 2020; 27(2): 159–168.

Fricke

Tsuno

, et al. A bio-inspired transportation network for scalable swarm foraging. In: IEEE International conference on robotics and automation (ICRA), Paris, France, 2020, pp.6120–6126. New York, NY: IEEE.

Podevijn

O’Grady

Mathews

, et al. Investigating the effect of increasing robot group sizes on the human psychophysiological state in the context of human–swarm interaction. Swarm Intell 2016; 10: 193–210.

Yang

Sun

Fang

Neuroadaptive control for complicated underactuated systems with simultaneous output and velocity constraints exerted on both actuated and unactuated states. IEEE Trans Neural Netw Learn Syst 2021; 1–11. DOI: 10.1109/TNNLS.2021.3115960

Yang

Chen

Sun

, et al. Adaptive neural network output feedback control of uncertain underactuated systems with actuated and unactuated state constraints. IEEE Trans Syst Man Cybern Syst 2021. DOI: 10.1109/TSMC.2021.3131843

Huang

Qin

Wang

Position control of second order nonholonomic planar underactuated mechanical system. J Wuhan Inst Technol 2021; 43(5): 567–572.

Wang

Chen

, et al. Drift suppression control based on online intelligent optimization for planar underactuated manipulator with passive middle joint. IEEE Access 2021; 9: 38611–38619.

Peng

Jiang

Wang

Event-triggered dynamic surface control of an underactuated autonomous surface vehicle for target enclosing. IEEE Trans Ind Electron 2021; 68(4): 3402–3412.

Sun

Fang

Nonlinear tracking control of underactuated cranes with load transferring and lowering: Theory and experimentation. Automatica 2014; 50(9): 2350–2357.

10.

Wang

Yang

Zhang

Iterative convergence control method for planar underactuated manipulator based on support vector regression model. Nonlinear Dyn 2020; 102: 2711–2724.

11.

Wang

, et al. A general position control method for planar underactuated manipulators with second-order nonholonomic constraints. IEEE Trans Cybern 2021; 51(9): 4733–4742.

12.

Huang

Lai

Zhang

, et al. Trajectory planning and tracking control for positioning of planar three-link underactuated manipulator. In: IEEE Asian control conference (ASCC), Kitakyushu, Japan, 2019, pp. 1243–1247. New York, NY: IEEE.

13.

Huang

Lai

Zhang

, et al. Virtual model reduction-based control strategy of planar three-link underactuated manipulator with middle passive joint. Int J Control Autom Syst 2021; 19(1): 29–39.

14.

Huang

Wang

Review of control method of planar underactuated mechanical system. J Wuhan Inst Technol 2021; 43(4): 448–454.

15.

Wang

Lai

Rapid position control approach based on variable design parameter for planar Acrobot. Electric Machines and Control 2017; 21(9): 110–118.

16.

Wang

, et al. Control strategy based on Fourier transformation and intelligent optimization for planar Pendubot. Inf Sci 2019; 491: 279–288.

17.

Wang

Lai

Zhang

, et al. A new control method for planar four-link underactuated manipulator based on intelligence optimization. Nonlinear Dyn 2019; 96: 573–583.

18.

Zhang

Lai

Wang

, et al. A quick position control strategy based on optimization algorithm for a class of first-order nonholonomic system. Inf Sci 2018; 460–461: 264–278.

19.

Xiong

Lai

A stable control for second-order nonholonomic planar underactuated mechanical system: energy attenuation approach. Int J Control 2018; 91(7): 1630–1639.

20.

Zhang

Lai

Wang

, et al. A novel position-posture control method using intelligent optimization for planar underactuated mechanical systems. Mech Mach Theory 2019; 140: 258–273.

21.

Huang

Lai

Wang

, et al. Position control of planar three-link underactuated manipulator based on trajectory planning. Control Decis 2020; 35(2): 382–388.

22.

Liu

A new method for position control of planar 3-Dof underactuated robots. In: IEEE world congress on intelligent control and automation. Chongqing, China, 2008, pp.9060–9065. New York, NY: IEEE.

23.

Huang

Lai

Zhang

, et al. A general control strategy for planar 3-DoF underactuated manipulators with one passive joint. Inf Sci 2020; 534: 139–153.

24.

Huang

Qin

Wang

A general control strategy of planar multi-link underactuated manipulator with passive last joint based on nilpotent approximation and intelligent optimization. In: IEEE Chinese automation congress (CAC). Shanghai, China, 2020, pp.4029–4031. New York, NY: IEEE.

25.

De Luca

Oriolo

. Trajectory planning and control for planar robots with passive last joint. Int J Rob Res 2002; 21(5–6): 575–590.

26.

Yoshikawa

Kobayashi

Watanabe

Design of a desirable trajectory and convergent control for 3-DOF manipulator with a nonholonomic constraint. In: IEEE international conference on robotics and automation, San Francisco, CA, 2000, pp.1805–1810. New York, NY: IEEE.

27.

Jia

Shan

Finite-time trajectory tracking control of space manipulator under actuator saturation. IEEE Trans Ind Electron 2020; 67(3): 2086–2096.

28.

Yang

Cai

, et al. Differential evolution with auto-enhanced population diversity. IEEE Trans Cybern 2015; 45(2): 302–315.

29.

LaSalle

JP.

Stability theory for ordinary differential equations. J Differ Equ 1968; 4(1): 57–65.