Download PDF
Research Article  |  Open Access  |  28 Jun 2023

Fixed-time integral sliding mode tracking control of a wheeled mobile robot

Views: 630 |  Downloads: 141 |  Cited:   0
Complex Eng Syst 2023;3:10.
10.20517/ces.2023.14 |  © The Author(s) 2023.
Author Information
Article Notes
Cite This Article

Abstract

This paper presents a fixed-time integral sliding mode control scheme for a nonholonomic wheeled mobile robot (WMR). To achieve the trajectory tracking mission, the dynamic model of a WMR is first transformed into a second-order attitude subsystem and a third-order position subsystem. Two novel continuous fixed-time disturbance observers are proposed to estimate the external disturbances of the two subsystems, respectively. Then, trajectory tracking controllers are designed for two subsystems by utilizing the reconstructed information obtained from the disturbance observers. Additionally, an auxiliary variable that incorporates the Gaussian error function is introduced to address the chattering problem of the control system. Finally, the proposed control scheme is validated by a wheeled mobile robotic experimental platform.

Keywords

Wheeled mobile robot, trajectory tracking, disturbance observer, fixed-time stability, integral sliding mode control

1. INTRODUCTION

During the past decades, the wheeled mobile robot (WMR) has attracted extensive attention as it is widely used in various fields. The research on the WMR mainly includes robot positioning, motion planning, and motion control, among which the motion control is a fundamental problem. There are three main parts of the motion control, including point stabilization, path planning, and trajectory tracking [1]. The trajectory tracking control is a significant field in motion control, which has been studied extensively in recent years [2]. In practical engineering applications, a WMR is a highly coupled system with nonholonomic constraints and external disturbances. Hence, it is significant to design an anti-interference trajectory tracking control scheme with superior performance. At present, the design of the tracking controller of a WMR is mainly based on two types: one is to consider only the kinematic model [3], while the other is to design on the basis of kinematic and dynamic models [4]. The kinematic model-based control only considers the linear velocity and angular velocity as the control inputs. Compared with the kinematic model, the introduction of dynamic models can solve the external disturbance problem and the crucial nonholonomic constraint problem [5].

In [6], the system with nonhonolomic constraints was transformed into an extended chain system by coordinate transformation. On this basis, some scholars have designed the trajectory tracking control schemes by transforming the kinematic model of a WMR into a chain structure [7]. In practice, there is a problem called "excellent velocity tracking" [8] when designing a trajectory tracking controller only based on a kinematic system. Thus, it is more reasonable to take the force or torque as inputs of the control system instead of the speed. Meanwhile, external disturbances can be further taken into account. Nevertheless, the design process of the controller that simultaneously incorporates both the kinematic and dynamic models is complicated. The work of Zhai and Song [9] transformed the dynamic error system into second-order and third-order subsystems. And an intermediate variable related to the position error is introduced to tackle the problem of constructing a control method for a third-order system using the terminal sliding mode control. However, the aforementioned control schemes can only achieve finite time stability. It is noteworthy that the upper limit of the convergence time is unknown and dependent on the initial states of the control system. To overcome this problem, fixed-time stable control methods are proposed [10]. In reference [11], a new integral sliding mode-based control (ISMC) scheme was developed and applied on the dynamic model of the WMR to enable the WMR to track the desired trajectory in a fixed time. However, there exists the singularity problem, making the WMR unable to track the arbitrary trajectories and limiting its practical application when the desired angular velocity is zero.

In the practical motion environment, there are external disturbances and uncertainties that can deteriorate the performance of the control system. To cope with the problem, an observer-based control scheme is an efficient method with disturbance-rejection performance [12]. The traditional observers can only achieve asymptotic stability of the observation errors, whereas the finite time disturbance observers were designed to improve the performance of the observer [13]. On this basis, the fault-tolerant attitude control problem of spacecraft under external disturbances was solved by the introduction of a continuous finite-time observer [14], which also restrains the chattering phenomenon. Zhang et al. put forward a novel continuous practical fixed-time disturbance observer and applied it on a WMR, which can not only avoid the chattering problem but also improve the ability to attenuate disturbance [15]. Different from the work of Zhang, the Gaussian error function, which is sometimes called probability integral [16], can also be used to develop a control scheme that improves the chattering problem [17].

Motivated by the above discussions, an integral sliding mode-based fixed-time trajectory tracking control scheme is proposed by combining the kinematic model with the dynamic model of a WMR in this paper. (1) A continuous fixed-time disturbance observer using the Gaussian error function is proposed, which avoids the chattering problem and estimates the external disturbance of a WMR accurately. (2) An auxiliary variable incorporating variable exponential coefficients is introduced to simplify the design process of the controller for the third-order subsystem and avoid the singularity problem simultaneously. (3) The reliability and effectiveness of the designed control scheme are verified by a comparative experiment conducted on a wheeled mobile experimental platform.

2. PRELIMINARIES AND PROBLEM STATEMENT

2.1. Preliminaries

Lemma 1[18] Consider the following system as

$$ \begin{equation} \dot{\mathit{\boldsymbol{x}}} = f(\mathit{\boldsymbol{x}})\mbox{, } \;\mathit{\boldsymbol{x}}(0) = \mathit{\boldsymbol{x_0}}\mbox{, } \;{\mathit{\boldsymbol{x}}} \in \boldsymbol{R} \end{equation} $$

If there exists a positive definite Lyapunov function $$ V({\mathit{\boldsymbol{x}}}) $$, which satisfies $$ \dot{V}({\mathit{\boldsymbol{x}}}) \leq-m_1 V^{a}({\mathit{\boldsymbol{x}}})-n_1 V^{b}({\mathit{\boldsymbol{x}}})+\varrho $$, where $$ m_1, n_1 $$, and $$ \varrho $$ are all positive constants. $$ 0<a<1, b>1 $$ are real numbers. Then the origin of the system (2) is fixed-time stable, and the settling time is bounded by $$ t_1\leq\frac{1}{m_1\vartheta (1-a)}+\frac{1}{n_1\vartheta (b-1)} $$ with $$ 0<\vartheta<1 $$.

Lemma 2[16] The Gaussian error function is defined as follows:

$$ \begin{equation} \operatorname{erf}(x) = \frac{2}{\sqrt{\pi}}\int_{0}^{x}e^{-2t^2}\rm{d}t \end{equation} $$

where $$ e $$ is the natural constant. If $$ 0\leq x<1 $$, then the Gaussian error function will satisfy $$ \frac{1}{2}x\leq \operatorname{erf}(x)\leq 2x $$.

Lemma 3[19] For $$ x \in \mathit{\boldsymbol{{\rm{R}}}} $$ and $$ \mu>0 $$, one gets the following chain of inequalities: $$ x\tanh(\frac{x}{\mu})<x\operatorname{erf}(\frac{x}{\mu})<|x| $$.

Lemma 4[20] The following inequality will hold $$ |x|-\frac{\varepsilon}{\kappa} \leq x \tanh (\kappa x) $$ for any $$ \kappa>0 $$ and for any $$ \varepsilon\in \mathit{\boldsymbol{{\rm{R}}}} $$, where $$ \varepsilon = e^{-(\varepsilon+1)} $$. Then, $$ \varepsilon = 0.2785 $$ can be obtained.

2.2. Dynamic model of WMR

A nonholonomic WMR system is shown in Figure 1. It consists of two balance wheels and two driving wheels, and the line between the balance wheels is perpendicular to the line between the driving wheels. The distance between the driving wheel and the barycentric coordinate is $$ R $$, and $$ r $$ is the radius of the driving wheel. The position and attitude control is achieved by independent direct current motors, which provide the appropriate torques to the driving wheels. One assumes that the center of mass of the WMR coincides with the geometric center. Then, the dynamic model of the WMR is expressed in the form of [21]

$$ \begin{equation} \left\{\begin{array}{c} \begin{aligned} &\dot{x} = v \cos{\theta}\\ &\dot{y} = v \sin{\theta}\\ &\dot{\theta} = \omega\\ &J \dot{\omega} = u_{1}+d_{1}\\ &m \dot{v} = u_{2}+d_{2} \end{aligned} \end{array}\right.\\ \end{equation} $$

Fixed-time integral sliding mode tracking control of a wheeled mobile robot

Figure 1. Physical model of a WMR. WMR: wheeled mobile robot.

with $$ u_{1} = \frac{R}{r}(\tau_{1}-\tau_{2}) $$, and $$ u_{2} = \frac{1}{r}(\tau_{1}+\tau_{2}) $$. $$ \tau_{1} $$ and $$ \tau_{2} $$ present the control torques. $$ v $$ and $$ \omega $$ are the linear and angular velocities of the WMR, respectively. $$ m $$ denotes the mass, and $$ J $$ the moment of inertia. $$ (x, y) $$ is the actual coordinates. $$ \theta $$ is the orientation of the vehicle counterclockwise from the positive direction of the $$ X $$ axis. $$ (\dot{x}, \dot{y}, \dot{\theta}) $$ denotes the motion of the WMR. $$ d_{1} $$ and $$ d_{2} $$ represent the external disturbances.

The reference trajectory is defined as

$$ \begin{equation} \left\{\begin{array}{c} \begin{aligned} &\dot{x}_{r} = v_r \cos{\theta_{r}}\\ &\dot{y}_{r} = v_r \sin{\theta_{r}}\\ &\dot{\theta}_{r} = \omega_{r}\\ \end{aligned} \end{array}\right.\\ \end{equation} $$

where $$ x_{r} $$, $$ y_{r} $$, and $$ \theta_{r} $$ denote the position and attitude states of the virtual WMR, respectively.

Assumption 1: Suppose $$ \omega_{r} $$, $$ \dot{\omega}_{r} $$, $$ v_{r} $$, and $$ \dot{v}_{r} $$ are satisfied with $$ |\omega_{r}|\leq\omega_{r_{\rm{max}}} $$, $$ |\dot{\omega}_{r}|\leq\omega_{1_{\rm{max}}} $$, $$ |v_r|\leq v_{r_{\rm{max}}} $$, And $$ |\dot{v}_r|\leq v_{1_{\rm{max}}} $$, where $$ \omega_{r_{\rm{max}}} $$, $$ \omega_{1_{\rm{max}}} $$, $$ v_{r_{\rm{max}}} $$, and $$ v_{1_{\rm{max}}} $$ are positive constants.

Assumption 2: Suppose the $$ d_{1} $$, $$ d_{2} $$, and their derivatives exist with bounds, which is given by $$ |d_1|\leq k_{1m}, |d_2|\leq k_{2m} $$, where $$ k_{1m} $$ and $$ k_{2m} $$ are all positive constants.

Then, the tracking errors of the WMR are expressed as

$$ \begin{equation} \left[\begin{array}{c} {x_{e}}\\ {y_{e}}\\ {\theta_{e}}\end{array}\right] = \left[\begin{array}{ccc} {\cos\theta} & {\sin\theta} & {0} \\ {-\sin\theta} & {\cos\theta} & {0}\\ {0} & {0} & {1} \end{array}\right] \left[\begin{array}{c} {x-x_{r}} \\ {y-y_{r}} \\ {\theta-\theta_{r}}\end{array}\right]. \end{equation} $$

Furthermore, the error dynamics system could be transformed in the form of

$$ \begin{equation} \left\{\begin{array}{c} \begin{aligned} &\dot{x}_{e} = \omega y_{e}-v+v_{r} \cos{\theta_{e}}\\ &\dot{y}_{e} = v_{r} \sin{\theta_{e}}-\omega x_{e}\\ &\dot{\theta}_{e} = \omega-\omega_{r}\\ &J \dot{\omega} = u_{1}+d_{1}\\ &m \dot{v} = u_{2}+d_{2} \end{aligned} \end{array}\right.\\ \end{equation} $$

To simplify the whole design process, the system (6) can be divided into two subsystems, which contain a second-order subsystem:

$$ \begin{equation} \left\{\begin{array}{c} \begin{aligned} &\dot{\theta}_{e} = \omega-\omega_{r}\\ &J \dot{\omega} = u_{1}+d_{1}\\ \end{aligned} \end{array}\right. \end{equation} $$

and a third-order subsystem:

$$ \begin{equation} \left\{\begin{array}{c} \begin{aligned} &\dot{x}_{e} = \omega y_{e}-v+v_{r} \cos{\theta_{e}}\\ &\dot{y}_{e} = v_{r} \sin{\theta_{e}}-\omega x_{e}\\ &m \dot{v} = u_{2}+d_{2} \end{aligned} \end{array}\right.\\ \end{equation} $$

3. FIXED-TIME TRAJECTORY CONTROL

In this section, a fixed-time sliding mode control scheme is developed to realize the fast and high-accuracy trajectory tracking control of a WMR under external disturbances. Firstly, a fixed-time disturbance observer and a new fixed-time sliding mode surface are proposed for the second-order subsystem (7). On this basis, a fixed-time controller is constructed to make the error variables, $$ \theta_{e} $$ and $$ \dot{\theta}_{e} $$, converge into a small region around the origin. Then, a fixed-time controller is developed for the third-order subsystem (8), which guarantees that the system state variables, $$ x_e, y_e, $$ and $$ v $$, are all uniformly ultimately bounded, and the tracking errors, $$ x_e $$ and $$ y_e $$, can converge into a small region around the origin in a fixed time.

3.1. Tracking control laws design for the second-order subsystem

Fixed-time disturbance observer

Firstly, for the attitude error subsystem (7), define an auxiliary variable as

$$ \begin{equation} \varsigma_1 = \omega-\varpi_1 \end{equation} $$

where $$ \varpi_1 $$ satisfies

$$ \begin{equation} \dot{\varpi}_1 = \frac{1}{J}u_1+l_{11}\operatorname{erf}(\varsigma_1)+l_{12}|\varsigma_1|^{\gamma_1}\operatorname{erf}(\varsigma_1) \end{equation} $$

The parameters $$ l_{11} $$ and $$ l_{12} $$ are positive constants with $$ l_{11}>k_{1m}/J $$. Let variable exponential coefficient $$ \gamma_1 = \frac{\lambda_0\varsigma_1^2}{1+\mu_0\varsigma_1^2} $$ with $$ \lambda_0 $$ and $$ \mu_0 $$ satisfying $$ 0<\mu_0<1 $$ and $$ 1+\mu_0<\lambda_0 $$.

Theorem 1For the second-order subsystem (7), if the disturbance observer is constructed as

$$ \begin{equation} \hat{d}_1 = J\Big(l_{11}\operatorname{erf}(\varsigma_1)+l_{12}|\varsigma_1|^{\gamma_1}\operatorname{erf}(\varsigma_1)\Big) \end{equation} $$

then it can estimate $$ d_1 $$ accurately in a fixed time. That is to say, the observation error $$ \tilde{d}_1 = d_1-\hat{d}_1 $$ can converge into a small region within a fixed time.

Proof of Theorem 1 Select a Lyapunov function as $$ V_2 = \varsigma^2_1 $$, differentiating it, one has

$$ \begin{equation} \begin{aligned} \dot{V}_2 & = 2\varsigma_1(\dot{\omega}-\dot{\varpi}_1)\\ & = 2\varsigma_1\Big(\frac{1}{J}(u_1+d_1)-\big(\frac{1}{J}u_1+l_{11}\operatorname{erf}(\varsigma_1)+l_{12}|\varsigma_1|^{\gamma_1}\operatorname{erf}(\varsigma_1)\big) \Big)\\ & = 2\varsigma_1\Big(-l_{11}\operatorname{erf}(\varsigma_1)-l_{12}|\varsigma_1|^{\gamma_1}\operatorname{erf}(\varsigma_1) +\frac{1}{J}d_1\Big)\\ & = -2\Big(l_{11}\varsigma_1\operatorname{erf}(\varsigma_1)-\frac{1}{J}\varsigma_1d_1+l_{12}\varsigma_1|\varsigma_1|^{\gamma_1}\operatorname{erf}(\varsigma_1)\Big)\\ &\leq -2\Big(l_{11}\varsigma_1{\rm{tanh}(\varsigma_1)} - \frac{1}{J}\varsigma_1d_1+l_{12}\varsigma_1|\varsigma_1|^{\gamma_1}\operatorname{erf}(\varsigma_1)\Big)\\ &\leq -2\Big(l_{11}|\varsigma_1|-l_{11}\epsilon_1-\frac{k_{1m}}{J}|\varsigma_1|+l_{12}\varsigma_1|\varsigma_1|^{\gamma_1}\operatorname{erf}(\varsigma_1)\Big) \\ &\leq -2l_{12}|\varsigma_1||\varsigma_1|^{\gamma_1}\operatorname{erf}(|\varsigma_1|)+2l_{11}\epsilon_1 \\ & = -2l_{12}|\varsigma_1|^{\gamma_1+1}\operatorname{erf}(|\varsigma_1|)+2l_{11}\epsilon_1 \end{aligned} \end{equation} $$

where $$ \epsilon_1 $$ is a positive constant.

Case 1 When $$ V_2 > 1 $$ and $$ |\varsigma_1|> 1 $$, one has $$ \operatorname{erf}(|\varsigma_1|)>\operatorname{erf}(1) $$ and $$ \frac{\lambda_0 \varsigma_1^2}{1+\mu_0 \varsigma_1^2}\geq \frac{\lambda_0}{1+\mu_0}>1 $$. Then (12) can be rewritten as

$$ \large \begin{equation} \begin{aligned} \dot{V}_2 &\leq -2(l_{12}\operatorname{erf}(1)-l_{11}\epsilon_1)|\varsigma_1|^{\frac{\lambda_0}{1+\mu_0}+1} \\ & \leq -2(l_{12}\operatorname{erf}(1)-l_{11}\epsilon_1)V^{\frac{\lambda_0+\mu_0+1}{2(1+\mu_0)}}_2 \end{aligned} \end{equation} $$

As $$ l_{12}\operatorname{erf}(1)-l_{11}\epsilon_1>0 $$ and $$ \bar{\gamma}_1 = \frac{\lambda_0+\mu_0+1}{2(1+\mu_0)}>1 $$, then all the solutions of $$ \{V_2> 1\} $$ will reach the set $$ \{V_2 \leq 1\} $$ within a fixed time $$ t_{f1} \leq \frac{1}{2(l_{12}\operatorname{erf}(1)-l_{11}\epsilon_1)(\bar{\gamma}_1-1)} $$.

Case 2 In the converse case $$ V_2 \leq 1 $$, one has

$$ \begin{equation} \begin{aligned} \dot{V}_2 &\leq -2l_{12}|\varsigma_1||\varsigma_1|^{\gamma_1}\operatorname{erf}(|\varsigma_1|)+2l_{11}\epsilon_1\\ \end{aligned} \end{equation} $$

As $$ 1+\mu_0 \varsigma^2_1 \geq 1 $$ and $$ |\varsigma_1|\leq 1 $$, it can be obtained that $$ \min\big(|\varsigma_1|^{\gamma_1}\big) \geq \min\big(|\varsigma_1|^{\lambda_0 \varsigma^2_1}\big) = e^{\frac{-\lambda_0}{2e}} $$. Considering the Lemma Lemma3, then (14) is converted into the following form

$$ \large \begin{equation} \begin{aligned} \dot{V}_2 &\leq -2l_{12}|\varsigma_1||\varsigma_1|^{\gamma_1}\operatorname{tanh}(|\varsigma_1|)+2l_{11}\epsilon_1\\ &\leq 2l_{12}|\varsigma_1||\varsigma_1|^{\gamma_1}+2l_{12}|\varsigma_1|^{\gamma_1}\epsilon_2+2l_{11}\epsilon_1\\ &\leq -2l_{12}|\varsigma_1||\varsigma_1|^{\gamma_1}+2l_{12} \epsilon_2+2l_{11}\epsilon_1\\ & \leq -2l_{12}|\varsigma_1|^{\gamma_1+1}+2l_{11}\epsilon_1+2l_{12}\epsilon_2\\ & \leq -2l_{12}e^{\frac{-\lambda_0}{2e}}|\varsigma_1|+2l_{11}\epsilon_1+l_{12}\epsilon_2\\ & \leq -b_1V^{\frac{1}{2}}_2+\tilde{\epsilon}\\ & \leq -l_{13} b_1V^{\frac{1}{2}}_2-(1-l_{13})b_1V^{\frac{1}{2}}_2+\tilde{\epsilon} \end{aligned} \end{equation} $$

with $$ b_1 = 2l_{12}e^{\frac{-\lambda_0}{2e}} $$, and $$ \tilde{\epsilon} = 2l_{11}\epsilon_1+2l_{12}\epsilon_2 $$. When $$ 0<l_{13}<1 $$, and $$ \tilde{\epsilon}-(1-l_{13})b_1V^{\frac{1}{2}}_2 \geq 0 $$, (15) can be simplified as $$ \dot{V}_2 \leq -l_{13} b_1V^{\frac{1}{2}}_2 $$. Then, the solution of $$ V_2 $$ will reach a small set $$ \Delta_1 $$, which is defined as $$ \Delta_1 = \big\{\varsigma_1|V_1(\varsigma_1)\leq(\frac{\tilde{\epsilon}}{b_1(1-l_{13})})^2\big\} $$ within a settling time $$ t_{f2}\leq \frac{2}{b_1l_{13}} $$.

In view of the above two cases, the auxiliary variable $$ \varsigma_1 $$ will converge into a small set $$ \Delta_1 = \big\{\varsigma_1|V_1(\varsigma_1)\leq (\frac{\tilde{\epsilon}}{b_1(1-l_{13})})^2\big\} $$ within settling time $$ t_f = t_{f1}+t_{f2} $$.

Then, the disturbance observation error

$$ \begin{equation} \begin{aligned} \tilde{d}_1& = d_1-\hat{d}_1\\ & = d_1- J\Big(l_{11}\operatorname{erf}(\varsigma_1)+l_{12}|\varsigma_1|^{\gamma_1}\operatorname{erf}(\varsigma_1)\Big)\\ \end{aligned} \end{equation} $$

The disturbance $$ d_1 $$ is bounded according to Assumption 1. Thus, the disturbance observer (11) can estimate $$ d_1 $$ accurately, and the observation error $$ \tilde{d}_1 $$ can remain in a small set

$$ \Delta_2 = \Big\{\varsigma_1\big| |\varsigma_1|\leq k_{1m}+J\big(l_{11}\operatorname{erf}(\Delta_1)+l_{12}|\Delta_1|^{\bar{\gamma}_1}\operatorname{erf}(\Delta_1)\big)\Big\} $$
after a fixed time, where $$ \bar{\gamma}_1 = \frac{\lambda_0 \Delta_1^2}{1+\mu_0\Delta_1^2} $$.

3.1.2. Fixed-time sliding mode controller

For the subsystem (7), define $$ \omega_{e} = \omega-\omega_{r} $$. A fixed-time integral sliding mode surface is introduced as follows [22]

$$ \begin{equation} s_{1} = \omega_{e}+\int_{0}^{t} \big(k_{11}(\lceil{\theta_e}\rfloor^{p_1} + \lceil{\theta_e}\rfloor^{q_1})+k_{12}(\lceil{\omega_e}\rfloor^{p_2} + \lceil{\omega_e}\rfloor^{q_2})\big)\rm{d} \tau \end{equation} $$

with $$ 0<p_i<1 $$, $$ q_i>1 $$, and $$ (i = 1, 2) $$. For any $$ x \in \rm{R} $$, $$ \alpha\in \rm{R}^+ $$, the notation is defined as $$ \lceil{x}\rfloor^{\alpha} = |x|^{\alpha}{\rm{sign}}(x) $$. Based on the sliding mode surface as (17), the fixed-time controller is designed as follows:

$$ \begin{equation} u_1 = -J\Big(k_{11}(\lceil{\theta_e}\rfloor^{p_1} + \lceil{\theta_e}\rfloor^{q_1})+k_{12}(\lceil{\omega_e}\rfloor^{p_2} + \lceil{\omega_e}\rfloor^{q_2})+\alpha_1\lceil{s_1}\rfloor^{p_3} + \alpha_2\lceil{s_1}\rfloor^{q_3}+\alpha_3\operatorname{erf}(s_1)-\dot{\omega}_r\Big)-\hat{d}_1 \end{equation} $$

where $$ \alpha_i, k_{1i}, (i = 1, 2) $$ are positive constants, $$ \alpha_3 $$ satisfies $$ \alpha_3\geq\frac{k_{1m}}{J} $$. In addition, $$ p_i, q_i, (i = 1, 2, 3) $$ are all positive odd integers with $$ 0<p_i<1, q_i>1 $$.

Theorem 2For the second-order system (7), if the fixed-time controller is constructed in the form of (18), then the real sliding mode variable will converge into a small set within a fixed time.

Proof of Theorem 2 Choose a Lyapunov function as $$ V_3 = \frac{1}{2}s^2_1 $$ and refer to Lemma 2 to Lemma 4, the time derivative of $$ V_3 $$ is

$$ \begin{equation} \begin{aligned} \dot{V}_3 & = s_1\Big(\dot{\omega}_e+k_{11}(\lceil{\theta_e}\rfloor^{p_1} + \lceil{\theta_e}\rfloor^{q_1})+k_{12}(\lceil{\omega_e}\rfloor^{p_2} + \lceil{\omega_e}\rfloor^{q_2}\Big)\\ & = s_1\Big(\frac{1}{J}(u_1+d_1)-\dot{\omega}_r+ k_{11}(\lceil{\theta_e}\rfloor^{p_1} + \lceil{\theta_e}\rfloor^{q_1})+k_{12}(\lceil{\omega_e}\rfloor^{p_2} + \lceil{\omega_e}\rfloor^{q_2})\Big)\\ & = s_1\Big(-\alpha_1\lceil{s_1}\rfloor^{p_3} -\alpha_2 \lceil{s_1}\rfloor^{q_3}-\alpha_3\operatorname{erf}(s_1)+\frac{1}{J}\tilde{d}_1\Big)\\ &\leq -\alpha_1|s_1|^{p_3+1}-\alpha_2|s_1|^{q_3+1}-\alpha_3 s_1 \tanh(s_1)+\frac{1}{J}|s_1||\tilde{d}_1|\\ &\leq -2^{\bar{p}_3} \alpha_1V^{\bar{p}_3 }_2 -2^{\bar{q}_3} \alpha_{2}V^{\bar{q}_3}_2+\bar{\vartheta}_1 \\ \end{aligned} \end{equation} $$

where $$ \bar{p}_3 = \frac{p_3+1}{2}, \bar{q}_3 = \frac{q_3+1}{2}, \bar{\vartheta}_1 = \alpha_3\vartheta_1 $$ with $$ \vartheta_1 $$ being a positive constant. By using Lemma Lemma4, the second-order system (7) is fixed-time stable. The sliding mode surface $$ s_1 $$ will converge into a small region $$ \Delta_3 = \Big\{s_1|V(s_1)\leq {\rm{min}}\big\{\big(\frac{c_2}{\alpha_1 2^{\bar{p}_3}}\big)^{\frac{1}{\bar{p}_3}}, \big(\frac{c_2}{\alpha_2 2^{\bar{q}_3}}\big)^{\frac{1}{\bar{p}_3}}\big\}\Big\} $$ around the origin in a fixed time $$ t_{s_1} $$, which is determined by $$ t_{s_1}\leq \frac{1}{\alpha_12^{\bar{p}_3}\phi_1(1-\bar{p}_3)}+\frac{1}{\alpha_2 2^{\bar{q}_3}\phi_1(\bar{q}_3-1)} $$. Then, one can obtain that variables $$ \theta_{e} $$ and $$ \omega_{e} $$ converge to zero along the real sliding mode in a fixed time[23].

3.2. Tracking control laws design for the third-order subsystem

After the angular error $$ \theta_{e} $$ converges to zero according to Theorem 2, one can obtain that $$ \sin{\theta_e} $$ equals zero, and $$ \cos{\theta_e} $$ equals 1. The system (8) can be simplified as

$$ \begin{equation} \left\{\begin{array}{c} \begin{aligned} &\dot{x}_{e} = \omega_r y_{e}-v+v_{r}\\ &\dot{y}_{e} = -\omega_r x_{e}\\ &m\dot{v} = u_{2}+d_{2} \end{aligned} \end{array}\right.\\ \end{equation} $$

3.2.1. Fixed-time disturbance observer

Introduce the following auxiliary variable for the simplified third-order subsystem (20)

$$ \begin{equation} \varsigma_2 = v-\varpi_2 \end{equation} $$

where $$ \varpi_2 $$ satisfies

$$ \begin{equation} \dot{\varpi}_2 = \frac{1}{m}u_2+l_{21}{\operatorname{erf}(\varsigma_2)}+l_{22}|\varsigma_2|^{\gamma_2}\operatorname{erf}(\varsigma_2) \end{equation} $$

where $$ \gamma_2 = \frac{\lambda_3\varsigma_1^2}{1+\mu_3\varsigma_1^2} $$, and $$ \lambda_3 $$ and $$ \mu_3 $$ are integers satisfying the constraints: $$ 0<\mu_3<1 $$, $$ 1+\mu_3<\lambda_3 $$. The parameters $$ l_{21} $$ and $$ l_{22} $$ are positive constants with $$ l_{21}>k_{2m} $$ and $$ l_{22}>0 $$.

Theorem 3For the simplified third-order subsystem (20), a fixed-time disturbance observer is developed in the form of

$$ \begin{equation} \hat{d}_2 = m\Big(l_{21}\operatorname{erf}(\varsigma_2)+l_{22}|\varsigma_2|^{\gamma_2}\operatorname{erf}(\varsigma_2)\Big) \end{equation} $$

then it can estimate $$ d_2 $$ in a fixed time, and the observation error $$ \tilde{d}_2 = d_2-\hat{d}_2 $$ can converge into a small region around the origin within a fixed time $$ t_{d_2} $$.

Proof of Theorem 3 Similar to the proof of Theorem 1.

3.2.2. Fixed-time sliding mode controller

For the third-order subsystem (20), introduce the following auxiliary variable:

$$ \begin{equation} \xi = x_e+\int_{0}^{t}\Big(\lambda_1 \operatorname{erf}(x_e)-\lambda_2 \operatorname{erf}(y_e)+\lambda_3 x_e y_e \operatorname{erf}(y_e)+k_1{|x_e|}^{\gamma_3}\operatorname{erf}\big(\frac{x_e}{\epsilon_3}\big)\Big)\rm{d} \tau \end{equation} $$

where $$ \lambda_1 $$$$ \lambda_2 $$, $$ k_1 $$, and $$ \epsilon_3 $$ are positive constants, $$ \lambda_2\leq\lambda_1\operatorname{erf}(1) $$. Let $$ \gamma_3 = \frac{\lambda_4 {x_e}^{2}}{1+\mu_4{x_e}^{2}} $$, with $$ \lambda_4 $$ and $$ \mu_4 $$ being integers and $$ 0<\mu_4<1 $$, $$ 1+\mu_4<\lambda_4 $$.

Select a fixed-time sliding mode surface as:

$$ \begin{equation} s_2 = \dot{\xi}+\int_{0}^{t}\Big(k_{21}\left(\lceil{\xi}\rfloor^{p_4} + \lceil{\xi}\rfloor^{q_4}\right)+k_{22}\left(\lceil{\dot{\xi}}\rfloor^{p_5} + \lceil{\dot{\xi}}\rfloor^{q_5}\right)\Big)\rm{d} \tau \end{equation} $$

where $$ k_{21} $$ and $$ k_{22} $$ are positive constants, $$ 0<p_i<1 $$ and $$ q_i>1 $$, $$ i = 4, 5 $$ are positive odd integers.

Theorem 4For the third-order subsystem (20), if the fixed-time sliding mode surface is chosen as (25) and the fixed-time controller is designed as (26),

$$ \begin{equation} \begin{aligned} u_2 = &m\Big(\dot{v}_r+\dot{\omega}_r y_e+\omega_r\dot{y}_e+\dot{h}(x_{e}, y_{e})+ k_{21}\left(\lceil{\xi}\rfloor^{p_4} + \lceil{\xi}\rfloor^{q_4}\right)\\ &\mathrel{\phantom{ = }}+k_{22}\left(\lceil{\dot{\xi}}\rfloor^{p_5} + \lceil{\dot{\xi}}\rfloor^{q_5}\right)+\beta_1\lceil{s_2}\rfloor^{p_6} +\beta_2\lceil{s_2}\rfloor^{q_6}+\beta_3\operatorname{erf}(s_2)\Big)+\hat d_2\\ \end{aligned} \end{equation} $$

then the sliding mode surface $$ s_2 $$ is fixed-time stable, which will converge into a small region of origin within settling time $$ t_{s_2}\leq \frac{1}{\alpha_4 2^{\bar{p}_6}\phi_2(1-\bar{p}_6)}+\frac{1}{\alpha_5 2^{\bar{q}_6}\phi_2(\bar{q}_6-1)} $$.

$$ h(x_e, y_e) = \lambda_1 \operatorname{erf}(x_e)-\lambda_2 \operatorname{erf}(y_e)+\lambda_3 x_e y_e \operatorname{erf}(y_e)+k_1{|x_e|}^{\gamma_3}\operatorname{erf}\big(\frac{x_e}{\epsilon_2}\big) $$
, in which $$ \beta_1, \beta_2 $$ are positive constants, and $$ p_6, q_6 $$ are positive odd integers satisfying $$ 0<p_6<1 $$, $$ q_6>1 $$.

Proof of Theorem 4 The proof process will be conducted in 3 steps: (1) After the angular error $$ \theta_{e} $$ converges to zero according to Theorem 2, $$ s_2 $$ and the auxiliary variable $$ \xi $$ can converge into a small region around the origin within a fixed time; (2) The error variables $$ x_e, y_e $$ can converge into a small region around the origin within a fixed time; (3) It should be proved that $$ x_e $$ and $$ y_e $$ do not escape to infinity before the angular error $$ \theta_{e} $$ converges to zero.

Step 1 Select a positive Lyapunov function $$ V_4 = \frac{1}{2} {s}_2^2 $$, differentiating it and substituting (24)-(26) yields to

$$ \begin{equation} \begin{aligned} \dot{V}_4 & = s_2 \dot{s}_2 \\ & = s_2\big(-\beta_1\lceil{s_2}\rfloor^{p_6} -\beta_2\lceil{s_2}\rfloor^{q_6}-\beta_3\operatorname{erf}(s_2)+\frac{1}{m}(\hat{d_2}-d_2)\big)\\ &\leq -\beta_1 |s_2|^{p_6+1}-\beta_2 |s_2|^{q_6+1}-\beta_3s_2\operatorname{erf}(s_2)+\frac{1}{m}s_2\tilde{d}_2\\ &\leq -\beta_1 |s_2|^{p_6+1}-\beta_2 |s_2|^{q_6+1}-\beta_3s_2\operatorname{tanh}(s_2)+\frac{1}{m}s_2\tilde{d}_2\\ &\leq -\beta_1 |s_2|^{p_6+1}-\beta_2 |s_2|^{q_6+1}-\beta_3|s_2|+\frac{1}{m}|s_2||\tilde{d}_2|\\ &\leq -2^{\bar{p}_6}\beta_1 V^{\bar{p}_6}_3 -2^{\bar{q}_6}\beta_2 V^{\bar{q}_6}_3+\bar{\vartheta}_2 \end{aligned} \end{equation} $$

where $$ \bar{p}_6 = \frac{p_6+1}{2}, \bar{q}_6 = \frac{q_6+1}{2} $$, $$ \bar{\vartheta}_2 = \beta_3 \vartheta_2 $$ with $$ \vartheta_2 $$ being a positive constant. Using the Lemma Lemma4, the third-order subsystem (14) is fixed-time stable, and $$ s_2 $$ will converge into a small set $$ \Delta_4 = \Big\{s_2|V(s_2)\leq {\rm{min}}\big\{\big(\frac{c_3}{\alpha_1 2^{\bar{p}_6}}\big)^{\frac{1}{\bar{p}_6}}, \big(\frac{c_3}{\alpha_2 2^{\bar{q}_6}}\big)^{\frac{1}{\bar{p}_6}}\big\}\Big\} $$ around zero in the fixed time $$ t_{s_2} $$, which is determined by

$$ \begin{equation} t_{s_2}\leq \frac{1}{\alpha_4 2^{\bar{p}_6}\phi_2(1-\bar{p}_6)}+\frac{1}{\alpha_5 2^{\bar{q}_6}\phi_2(\bar{q}_6-1)} \end{equation} $$

Then, $$ s_2 $$ will hold in a small region of origin, which guarantees a real sliding mode surface [23]. Therefore, the auxiliary variable $$ \xi $$ and its derivative $$ \dot{\xi} $$ will also converge into the origin along the sliding mode surface [24].

Step 2 According to (25), when the auxiliary $$ \dot{\xi} = 0 $$, one has

$$ \begin{equation} \dot{x}_e = -\lambda_1 \operatorname{erf}(x_e)+\lambda_2 \operatorname{erf}(y_e)-\lambda_3 x_e y_e \operatorname{erf} (y_e)-k_1{|x_e|}^{\gamma_3}\operatorname{erf}\big(\frac{x_e}{\epsilon_3}\big) \end{equation} $$

Choose a Lyapunov function as $$ V_5 = x^{2}_e $$, the time derivative of $$ V_5 $$ is

$$ \begin{equation} \begin{aligned} \dot{V}_5 & = 2 x_e \dot{x}_e\\ & = -2\Big(\lambda_1 x_e \operatorname{erf}(x_e)- \lambda_2 x_e \operatorname{erf}(y_e)+\lambda_3 x^{2}_e y_e \operatorname{erf}(y_e)+k_1{|x_e|}^{\gamma_3}\operatorname{erf}(\frac{x_e}{\epsilon_3})\Big)\\ &\leq -2\Big(\lambda_1 |x_e|-\lambda_2 |x_e|- \lambda_1 \epsilon_4 +k_1 x_e |x_e|^{\gamma_3}\operatorname{erf}(\frac{x_e}{\epsilon_3})\Big)\\ &\leq -2\Big(k_1\epsilon_3 |x_e| |x_e|^{\gamma_3}\operatorname{erf}(|\frac{x_e}{\epsilon_3}|)-\lambda_1 \epsilon_4\Big)\\ \end{aligned} \end{equation} $$

where $$ \epsilon_4 $$ is a positive constant. The rest of the proof is similar to the proof of Theorem 2. There exists a constant $$ 0<\vartheta_3<1 $$ such that the variable $$ x_e $$ will reach and keep in a small region $$ \Delta_5 $$ around the origin within a fixed time $$ T_2 $$:

$$ \begin{equation} T_2\leq \frac{1}{k_1\vartheta_3e^{-\frac{\lambda_4}{2e}} }+\frac{1}{2k_1\epsilon_3(\frac{\lambda_4}{1+\mu_4}-1)}\ \end{equation} $$

Then, it can be obtained that the $$ \dot{x}_e $$ is a uniformly continuous form (29). Employ Barbalat Lemma [25] to prove $$ \dot{x}_e \rightarrow 0 $$ as $$ t \rightarrow \infty $$, then $$ \dot{x}_e $$ is bounded after the variable $$ x_e $$ converges. Hence, there exists a small region $$ \Delta_5 $$ around the origin that $$ y_e $$ can converge into $$ \Delta_5 $$.

Step 3 Before the angular error $$ \theta_{e} $$ converges to zero, $$ \theta_e\neq0 $$, such that subsystem (13) cannot be simplified as (19). It should be proved that system state variables $$ x_e, y_e, $$ and $$ v $$ are bounded before the angular error $$ \theta_{e} $$ converges to zero.

Consider the following bounded function:

$$ \begin{equation} V_6 = \frac{1}{2}x^{2}_e+\frac{1}{2}y^{2}_e+|v| \end{equation} $$

The time derivative of $$ V_6 $$ is

$$ \begin{equation} \begin{aligned} \dot{V}_6 & \leq |x_e||\dot{x}_e|+|y_e||\dot{y}_e|+|\dot{v}| \\ & \leq |x_e||\dot{x}_e|+|y_e||\dot{y}_e|+m\Big(\omega_{r}\dot{y}_e+|\dot{h}(x_e, y_e)|+k_{21}\big(|\xi|^{p_4}+|\xi|^{q_4}\big)+k_{22}(|\dot{\xi}|^{p_5}+|\dot{\xi}|^{q_5}\big) \\ &\mathrel{\phantom{ = }}+\beta_1|s|^{p_6}_2+\beta_2|s|^{q_6}_2\Big)+\frac{|\tilde{d}_2|}{m} \\ & \leq |x_e||\dot{x}_e|+|y_e||\dot{y}_e|+m\Big(|\dot{y}_e|+k_{21}\big(|\xi|^{p_4}+|\xi|^{q_4}\big)+k_{22}\big(|\dot{\xi}|^{p_5}+|\dot{\xi}|^{q_5}\big) +\beta_1|s|^{p_6}_2+\beta_2|s|^{q_6}_2\\ &\mathrel{\phantom{ = }}+\lambda_3\big(|\dot{x}_e||y_e|+|x_e||\dot{y}_e|+|x_e||y_e||\dot{y}_e|\big) +\frac{k_1\lambda_4|x_e||\dot{x}_e|}{1+\mu_0 x^{2}_e} |x_e|^{\gamma_3} \Big(1+\frac{2|x_e|}{1+\mu_4x^2_e}\Big) \\ &\mathrel{\phantom{ = }}+\frac{2k_1}{\epsilon_2\sqrt{\pi}}|x_e|^{\gamma_3}|\dot{x}_e|\Big)+\frac{|\tilde{d}_2|}{m} \\ &\leq |x_e||\dot{x}_e|+|y_e||\dot{y}_e|+m\Big(|\dot{y}_e|+k_{21}\big(|\xi|^{p_4}+|\xi|^{q_4}\big)+k_{22}\big(|\dot{\xi}|^{p_5}+|\dot{\xi}|^{q_5}\big) +\beta_1|s|^{p_6}_2+\beta_2|s|^{q_6}_2\\ &\mathrel{\phantom{ = }}+\lambda_3\big(|\dot{x}_e||y_e|+|x_e||\dot{y}_e|+|x_e||y_e||\dot{y}_e|\big)+\frac{k_1\lambda_4}{\mu_4} (1+2|x_e|)|x_e|^{\frac{\lambda_4}{\mu_4}-1}|x_e|\\ &\mathrel{\phantom{ = }}+\frac{2k_1}{\epsilon_2\sqrt{\pi}}|x_e|^{\frac{\lambda_4}{\mu_4}}|x_e|\Big)+\frac{|\tilde{d}_2|}{m} \\ \end{aligned} \end{equation} $$

Let $$ \eta_1 = \sqrt{x^{2}_e+y^{2}_e+|v|}\geq\eta>1 $$, then one has the following inequalities: $$ |x_e|\leq\eta_1, |y_e|\leq\eta_1, |v|\leq \eta_1 $$. Furthermore, there exist positive constants $$ a_i(i = 3, 4), b_l, c_l, (l = 4, 5, ..., 9) $$, which satisfy $$ |x_e|^{\frac{\lambda_4}{\mu_4}}\leq a_3 \eta_1 $$, $$ |x_e|^{\frac{\lambda_4}{\mu_4}-1}\leq a_4 \eta_1 $$, $$ |s_2|^{p_6}<b_4+c_4\eta_1 $$, $$ |s_2|^{q_6}<b_5+c_5\eta_1 $$, $$ |\xi|^{p_4}<b_6+c_6\eta_1 $$, $$ |\xi|^{q_4}<b_7+c_7\eta_1 $$, $$ |\dot{\xi}|^{p_5}<b_8+c_8\eta_1 $$, $$ |\dot{\xi}|^{q_5}<b_9+c_9\eta_1 $$. According to Theorem 3, the state variable $$ \theta_e, \omega $$ will converge into the origin within a fixed time $$ t_{s_1} $$, then one has $$ \theta_{e}\leq \theta_m $$, $$ |\omega|\leq \omega_m $$, $$ |\tilde{d}_2|<k_d $$. Further $$ |\dot{y}_e|<|\dot{x}_e|\leq v_{1_{\rm{max}}}+\omega_m \eta_1 $$ can be obtained. Then, (33) can be simplified as

$$ \begin{equation} \begin{aligned} \dot{V}_6 &\leq (2\eta_1+2m\lambda_3+\eta^2_1+m)(v_{1_{\rm{max}}}+\omega_{m}\eta_1)+\frac{mk_1\lambda_4}{\mu_4}(v_{1_{\rm{max}}}+\omega_{m}\eta_1)a_4\eta_1 \\ &\mathrel{\phantom{ = }}+\frac{2mk_1}{\epsilon_2\sqrt{\pi}}a_3\eta_1(v_{1_{\rm{max}}}+\omega_{m}\eta_1)+\frac{k_d}{m}\\ &\leq \frac{\eta^2_1}{2}2\Big(2v_{1_{\rm{max}}}+3\omega_{m}+2m\lambda_3+m+\frac{2mk_1a_4\lambda_4}{\mu_4}(v_{1_{\rm{max}}}+\omega_{m})+\frac{2ma_3k_1}{\epsilon_2\sqrt{\pi}}(v_{1_{\rm{max}}}+\omega_{m}) \Big) \\ &\mathrel{\phantom{ = }}+2m\lambda_3v_{1_{\rm{max}}}+\frac{k_d}{m}\\ &\leq KV_6+\varrho_1 \end{aligned} \end{equation} $$

where $$ K $$ and $$ \varrho_1 $$ satisfy the following constraints:

$$ \begin{equation} \begin{aligned} K & = 2\Big(2v_{1_{\rm{max}}}+3\omega_{m}+2m\lambda_3+m+\frac{2mk_1a_4\lambda_4}{\mu_4}(v_{1_{\rm{max}}}+\omega_{m})+\frac{2ma_3k_1}{\epsilon_2\sqrt{\pi}}(v_{1_{\rm{max}}}+\omega_{m}) \Big) \end{aligned} \end{equation} $$

$$ \begin{equation} \begin{aligned} \varrho_1 & = 2m\lambda_3v_{1_{\rm{max}}}+\frac{k_d}{m} \end{aligned} \end{equation} $$

On the contrary, if $$ \eta_1>1 $$, there exists a positive constant $$ \varrho_2 $$, which satisfies $$ \dot{V}_6\leq\varrho_2 $$. One has $$ \dot{V}_6\leq KV_6+\varrho_3 $$ for the state variable $$ x_e, y_e, v $$. Further, before the angular error $$ \theta_{e} $$ converges to zero, one can obtain

$$ \begin{equation} V_6\leq \big(V_6(0) +\frac{\varrho_3}{K}\big)e^{Kt}-\frac{{\varrho_3}}{K} \end{equation} $$

Remark 1 The auxiliary variable $$ \xi $$ in (24) can reduce the order of the third-order subsystem, which simplifies the process of the controller design. In addition, the controller developed in this paper can guarantee that the system state variables converge in a fixed time and the chattering problem is solved by using the error function erf($$ \cdot $$). Furthermore, utilizing the variable exponent coefficient in (24) avoids the common singularity problem.

4. EXPERIMENT RESULTS

To verify the effectiveness of the proposed control scheme, the trajectory tracking experiment is implemented on a Quanser QBot 2e mobile robot platform composed of a QBot 2e mobile robot, an OptiTrack system with 12 infrared cameras, and a computer. The experimental platform is presented in Figure 2. The whole closed-loop experiment structure is as follows: The simulation diagram is compiled on the host computer equipped with MATLAB/Simulink to transform the simulation into an executable file. And the control scheme is written to the Gumstix computer embedded in the QBot 2e through wireless communication protocol. The real-time position information of the QBot 2e is obtained by the OptiTrack positioning system. Then the host computer calculates the information and transmits them to the embedded computer of a WMR for the input of real-time calculation of executable files. So as to complete the trajectory experiment of the mobile robot.

Fixed-time integral sliding mode tracking control of a wheeled mobile robot

Figure 2. The Quanser QBot 2e Mobile Robot Platform.

In the experiment, the physical parameters of the QBot 2e are chosen as follows: $$ m = 4\; \rm{kg} $$, $$ J = 2.5 $$$$ \rm{kg} \cdot \rm{m}^2 $$. The desired reference trajectory is set as $$ x_r = \cos(0.2t)\; {\rm{m}}, y_r = \sin(0.2t)\; {\rm{m}} $$. The initial values of the reference and practical trajectories are $$ [x_r(0), y_r(0), \theta_r(0)]^{\rm{T}} = [1, 0, \pi/2]^{\rm{T}} $$, $$ [x(0), y(0), \theta(0)]^{\rm{T}} = [0.7, -0.02, \pi/6]^{\rm{T}} $$, respectively. The main relevant parameters of the proposed control scheme are as follows: $$ k_1 = 0.001 $$, $$ k_{11} = k_{12} = 0.9 $$, $$ k_{21} = 0.05 $$, $$ k_{22} = 0.06 $$, $$ \epsilon_2 = 0.00001 $$. Choose the parameters $$ \alpha_1 = 2 $$, $$ \alpha_2 = 0.5 $$, $$ \beta_1 = \beta_2 = 1 $$ for the sliding mode surface $$ s_1 $$ in (17) and $$ s_2 $$ in (25), respectively.

It is obvious that the WMR trajectory tracking mission can be achieved by the designed control method as plotted in the red track in Figure 3. The time response curves of the sliding mode surfaces, $$ s_1 $$ and $$ s_2 $$, are shown in Figure 4, which converge very quickly. To illustrate the excellence of the proposed control method, a comparative experiment on the trajectory tracking of WMRs is conducted between this work and reference [26]. The control inputs of the designed control scheme and reference are shown in Figure 5, which are nonsingular and continuous. Figure 6 illustrates the tracking errors in this experiment, which have a big fluctuation due to the influence of external disturbances. In the experiment, the external disturbance is from the experimental environment, such as uneven ground. The observed disturbance values are shown in Figure 7, which indicates the effectiveness of the proposed disturbance observer in this work. From the experimental results, it can be concluded that the designed control scheme has the robustness against the external disturbance and high tracking accuracy.

Fixed-time integral sliding mode tracking control of a wheeled mobile robot

Figure 3. The comparative trajectory tracking experiment results of the WMR. WMR: wheeled mobile robot.

Fixed-time integral sliding mode tracking control of a wheeled mobile robot

Figure 4. Comparative results of sliding mode surfaces in the experiment.

Fixed-time integral sliding mode tracking control of a wheeled mobile robot

Figure 5. Comparative results of control torques in the experiment.

Fixed-time integral sliding mode tracking control of a wheeled mobile robot

Figure 6. Comparative results of tracking errors $$ x_e $$, $$ y_e $$, $$ \theta_{e} $$ in the experiment.

Fixed-time integral sliding mode tracking control of a wheeled mobile robot

Figure 7. Disturbance estimation $$ \hat{d_1} $$ and $$ \hat{d_2} $$ in the experiment.

5. DISCUSSION

In this paper, a universal control scheme for fixed-time trajectory tracking based on ISMC is put forward. The dynamic model of the WMR has been transformed into two error subsystems. Then utilizing the fixed-time technology and ISMC, a new fixed-time disturbance observer has been proposed and applied on the two error subsystems. Furthermore, an observer-based tracking control method has been proposed to achieve a trajectory tracking mission for the WMR, and guarantee the tracking error converges within a fixed time. Finally, the proposed control approach has been verified by a mobile robotic platform, and the experimental results show fine control performances. Our future work will focus on how to realize the formation tracking control of multi-wheeled mobile robots in both theory and experiment.

DECLARATIONS

Acknowledgments

The authors would like to thank the editors and reviewers for their valuable comments dedicated to this article.

Authors' contributions

Made significant contributions in writing, methodology, review: Li B

Executed writing-original draft, experiment verification: Ma L

Methodology validation: Wang C

Data result analysis: Ge C

Supervision and modification: Liu H

Availability of data and materials

Not applicable.

Financial support and sponsorship

This work was supported in part by the National Natural Science Foundation of China (62073212), Natural Science Foundation of Shanghai (23ZR1426600), and Innovation Fund of Chinese Universities Industry-University-Research (2021ZYB05004).

Conflicts of interest

All authors declared that there are no conflicts of interest.

Ethical approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Copyright

© The Author(s) 2023.

REFERENCES

1. Panahandeh P, Alipour K, Tarvirdizadeh B, Hadi A. A kinematic Lyapunov-based controller to posture stabilization of wheeled mobile robots. Mech Syst Signal Pr 2019;134:1-19.

2. Huang H, Li Y, Bai Q. An improved a star algorithm for wheeled robots path planning with jump points search and pruning method. Complex Eng Syst 2022;2:11.

3. Kanayama Y, Kimura Y, Miyazaki F, Noguchi T. A stable tracking control method for an autonomous mobile robot. In: Proceedings., IEEE International Conference on Robotics and Automation; 1990;1. pp. 384-89.

4. Fierro R, Lewis F. Control of a nonholonomic mobile robot: backstepping kinematics into dynamics. In: Proceedings of 1995 34th IEEE Conference on Decision and Control; 1995;4. pp. 3805-10.

5. Bloch A. Nonholonomic mechanics and control. Interdisciplinary Applied Mathematics New York, NY: Springer; 2003.

6. Murray R, Sastry S. Nonholonomic motion planning: steering using sinusoids. IEEE Trans Automat Contr 1993;38:700-16.

7. Tayebi A, Tadjine M, Rachid A. Invariant manifold approach for the stabilization of nonholonomic chained systems: Application to a mobile robot. Nonlinear Dynam 2001;24:167-81.

8. Wang X, Zhang G, Neri F, et al. Design and implementation of membrane controllers for trajectory tracking of nonholonomic wheeled mobile robots. Integr Comput-Aid E 2016;23:15-30.

9. Zhai J, Song Z. Adaptive sliding mode trajectory tracking control for wheeled mobile robots. Int J Control 2019;92:2255-62.

10. Ou M, Sun H, Zhang Z, Li L. Fixed-time trajectory tracking control for multiple nonholonomic mobile robots. T I Meas Control 2021;43:1596-608.

11. Li B, Zhang H, Xiao B, Wang C, Yang Y. Fixed-time integral sliding mode control of a high-order nonlinear system. Nonlinear Dynam 2022;107:909-20.

12. Liu Q, Cai Z, Chen J, Jiang B. Observer-based integral sliding mode control of nonlinear systems with application to single-link flexible joint robotics. Complex Eng Syst 2021;1:8.

13. Zhang Z, Leibold M, Wollherr D. Integral sliding-mode observer-based disturbance estimation for euler–lagrangian systems. IEEE Trans Contr Syst T 2020;28:2377-89.

14. Li B, Hu Q, Yang Y. Continuous finite-time extended state observer based fault tolerant control for attitude stabilization. Aerosp Sci Technol 2019;84:204-13.

15. Zhang H, Li B, Xiao B, Yang Y, Ling J. Nonsingular recursive-structure sliding mode control for high-order nonlinear systems and an application in a wheeled mobile robot. ISA T 2022;130:553-64.

16. Chevillard S. The functions erf and erfc computed with arbitrary precision and explicit error bounds. Inform Comput 2012;216:72-95.

17. Eltayeb A, Rahmat M, Basri MAM, Mahmoud MS. An improved design of integral sliding mode controller for chattering attenuation and trajectory tracking of the quadrotor UAV. Arab J Sci Eng 2020;45:6949-61.

18. Ba D, Li Y, Tong S. Fixed-time adaptive neural tracking control for a class of uncertain nonstrict nonlinear systems. Neurocomputing 2019;363:273-80.

19. Bagul Y, Chesneau C. Sigmoid functions for the smooth approximation to the absolute value function. MJPAA 2021;7:12-19.

20. Polycarpou M, Ioannou P. A robust adaptive nonlinear control design. Automatica 1996;32:423-27.

21. Jiang Z, Nijmeijer H. Tracking control of mobile robots: A case study in backstepping. Automatica 1997;33:1393-99.

22. Wang C, Wen G, Peng Z, Zhang X. Integral sliding-mode fixed-time consensus tracking for second-order non-linear and time delay multi-agent systems. J Franklin I 2019;356:3692-710.

23. Plestan F, Shtessel Y, Brégeault V, Poznyak A. New methodologies for adaptive sliding mode control. Int J Control 2010;83:1907-19.

24. Song T, Fang L, Wang H. Model‐free finite‐time terminal sliding mode control with a novel adaptive sliding mode observer of uncertain robot systems. Asian J of Control 2022;24:1437-51.

25. Li B, Gong W, Yang Y, Xiao B, Ran D. Appointed fixed time observer-based sliding mode control for a quadrotor UAV under external disturbances. IEEE Trans Aero Elec Sys 2022;58:290-303.

26. Tian B, Liu L, Lu H, et al. Multivariable finite time attitude control for quadrotor UAV: Theory and experimentation. IEEE Trans Ind Electron 2018;65:2567-77.

Cite This Article

Export citation file: BibTeX | RIS

OAE Style

Ma L, Wang C, Ge C, Liu H, Li B. Fixed-time integral sliding mode tracking control of a wheeled mobile robot. Complex Eng Syst 2023;3:10. http://dx.doi.org/10.20517/ces.2023.14

AMA Style

Ma L, Wang C, Ge C, Liu H, Li B. Fixed-time integral sliding mode tracking control of a wheeled mobile robot. Complex Engineering Systems. 2023; 3(2): 10. http://dx.doi.org/10.20517/ces.2023.14

Chicago/Turabian Style

Ma, Ling, Chenghu Wang, Cheng Ge, Hui Liu, Bo Li. 2023. "Fixed-time integral sliding mode tracking control of a wheeled mobile robot" Complex Engineering Systems. 3, no.2: 10. http://dx.doi.org/10.20517/ces.2023.14

ACS Style

Ma, L.; Wang C.; Ge C.; Liu H.; Li B. Fixed-time integral sliding mode tracking control of a wheeled mobile robot. Complex. Eng. Syst. 2023, 3, 10. http://dx.doi.org/10.20517/ces.2023.14

About This Article

© The Author(s) 2023. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, sharing, adaptation, distribution and reproduction in any medium or format, for any purpose, even commercially, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Data & Comments

Data

Views
630
Downloads
141
Citations
0
Comments
0
5

Comments

Comments must be written in English. Spam, offensive content, impersonation, and private information will not be permitted. If any comment is reported and identified as inappropriate content by OAE staff, the comment will be removed without notice. If you have any queries or need any help, please contact us at support@oaepublish.com.

0
Download PDF
Cite This Article 14 clicks
Like This Article 5 likes
Share This Article
Scan the QR code for reading!
See Updates
Contents
Figures
Related
Complex Engineering Systems
ISSN 2770-6249 (Online)

Portico

All published articles are preserved here permanently:

https://www.portico.org/publishers/oae/

Portico

All published articles are preserved here permanently:

https://www.portico.org/publishers/oae/