Next Article in Journal
Multimodal Power Management Based on Decision Tree for Internet of Wearable Things Systems
Next Article in Special Issue
Ship Defense Strategy Using a Planar Grid Formation of Multiple Drones
Previous Article in Journal
Dynamics Modeling and Load-Sharing Performance Optimization of Concentric Face Gear Split-Torque Transmission Systems
Previous Article in Special Issue
Three-Dimensional Rendezvous Controls of Multiple Robots with Amplitude-Only Measurements in Cluttered Underwater Environments
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hierarchical Sliding Mode Control Combined with Nonlinear Disturbance Observer for Wheeled Inverted Pendulum Robot Trajectory Tracking

School of Automation, Beijing Information Science & Technology University, Beijing 100192, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(7), 4350; https://doi.org/10.3390/app13074350
Submission received: 28 February 2023 / Revised: 27 March 2023 / Accepted: 27 March 2023 / Published: 29 March 2023
(This article belongs to the Special Issue Advances in Robot Path Planning, Volume II)

Abstract

:

Featured Application

The research presented herein could be used mostly in warehouse logistics transport activities in smart manufacturing.

Abstract

A proposed optimized model for the trajectory tracking control of a wheeled inverted pendulum robot (WIPR) system is presented in this study, which addresses the problem of poor trajectory tracking performance in the presence of unknown disturbances due to the nonlinear and underactuated characteristics of the system. First, a kinematic controller was used to track a reference trajectory and generate a control law that specifies the desired forward and rotation speeds of the system. Next, a nonlinear disturbance observer (NDO) was designed to enhance the system’s robustness to external disturbances and improve its tracking performance. Then, the coupled system state variables were decoupled into two subsystems: a forward rotation subsystem and a tilt angle velocity subsystem. An improved hierarchical sliding mode controller was designed to control these subsystems separately. Finally, simulation experiments were conducted to compare the proposed method with a common sliding mode control approach. The simulation results demonstrate that the proposed method achieves better tracking performance in the presence of unknown disturbances.

1. Introduction

With the rapid development of technology, human society has simultaneously achieved increased convenience and comfort [1]. In today’s factory warehouses and production lines, a variety of robots add to the possibilities of Industry 4.0. In the warehouses of more advanced companies, transport robots can be found everywhere, replacing traditional manpower and eliminating the need for workers to carry out repetitive lifting and carrying. These “smart” robots can accomplish a task as long as they can follow the requirements of a given transport trajectory. This paper studies the trajectory tracking of a mobile wheeled inverted pendulum on a given reference trajectory to achieve perfect tracking of the ideal motion trajectory of the robots, meeting the requirements of factory transport robots and providing a powerful source of assistance to realizing smart factories [2,3,4,5].
Mobile wheeled inverted pendulum models, such as WIPRs, have attracted much attention because of their special advantages, such as compactness, mobility, and human-like functions. WIPRs are widely used to verify the effectiveness of nonlinear underactuated control methods, and compared with the traditional inverted pendulum, WIPRs have more applications than traditional inverted pendulum vehicles, especially in unknown, dynamic, and nonlinear environments, and are commonly used in logistics transportation, commuting, and navigation, as well as in the aforementioned application in the environment of factory transportation. However, a WIPR is classified as a typical model of nonlinear underactuated systems with two input torques driving two wheels and three degrees of freedom (forward, rotation, and tilt angle of the pendulum), and achieving its high-performance motion control is still a challenging task for the control community [6,7,8,9].
On the one hand, when a WIPR moves, it is always assumed that the ground can provide enough friction to prevent the robot from side-slipping and wheel-sliding (i.e., the robot is guaranteed to move with purely rolling wheels without skidding phenomena), which is a non-complete constraint at this point. On the other hand, consider that the underdriven inverted pendulum body needs to use the input torque of two driving wheels to control the three degrees of freedom of WIPR forward movement, rotation, and the angle of the inverted pendulum. If we want to use MWIPR to track the trajectory, we need to drive a WIPR in real time to control the three form variables with two input variables, which is a typical underactuated problem. Finally, in the real world, WIPRs operate in factories or other similar environments and always encounter various unknown disturbances that interfere with the system. Therefore, the three problems of incomplete constraints and underactuated and unknown perturbations are the main challenges faced by this particular mobile robot for trajectory tracking control [10,11,12,13].
The three issues mentioned above are of importance for the following reasons. First, the incomplete constraint will lead to the WIPR being unable to follow any trajectory movement, especially in the case of high-speed heavy load; if the robot’s incomplete constraints are not considered in motion planning, this is likely to lead to untimely obstacle avoidance and unreachable trajectory. Second, underdriven robots often have excellent dynamic performance or price advantages in terms of drive cost, but their biggest problem is the higher requirements in controller design. Finally, unknown disturbances will affect the control accuracy of the system to a certain extent, and more seriously, will affect the stability of the control system [14,15].
Many researchers and practitioners have proposed several control algorithms to overcome the difficulties faced by the problems associated with WIPR systems. One of the widely used methods is fuzzy control, which is an empirical, rule-based control technique that can effectively control nonlinear systems. By establishing a dynamic model of the WIPR, a fuzzy-logic-based controller can be designed to take the position and angle information of the WIPR as the input and output control signals to control its motion state. For instance, Jian Huang, in [16], proposed an Integral Interval Type 2 Fuzzy Logic (IT2FL) method that can maintain the MTWIP equilibrium while obtaining the desired position and orientation to make it work in an uncertain environment. However, the disadvantages of fuzzy control include low control accuracy, strong dependence on control rules, and difficulty in designing control rules.
The second control algorithm type is neural network control. Chenguang Yang [17] decomposed the underdriven WIPR model into two subsystems. The approximation characteristics of the neural network were used for motion control of the fully driven subsystem, and the sub-fully driven system was used to indirectly control the tilt angle motion of the pendulum. However, the method requires a large number of wavelet coefficient vectors, making the neural network computationally intensive.
Finally, sliding mode control, as the most typical robust control method, shows good tracking performance and strong robustness, which support its wide use in linear and nonlinear systems. For underactuated systems, various sliding mode control methods have been proposed by researchers to achieve different control effects, such as integral sliding mode control, terminal sliding mode control, and hierarchical sliding mode control [18,19,20]. Among them, the application of hierarchical sliding mode control in practical underdriven systems is receiving more and more attention, such as balancing control of a double-inverted pendulum and trajectory tracking control of a wheeled inverted pendulum. Nabanita Adhikary, in [21], proposed an integral inverse-step sliding mode controller for underdriven system control. A feedback control law was designed based on the backpropagation method, and a sliding surface was introduced in the final stage of the algorithm. Jian Huang, in [22], designed two terminal sliding mode controllers to control the speed and braking of a UW-Car based on the dynamic model and the terminal sliding mode control method. He Ping [23] proposed a hierarchical sliding mode controller (HSMC) developed to simultaneously perform speed control and balance control of a two-wheeled self-balancing vehicle (TWSBV).
Hierarchical sliding mode control is a control strategy based on sliding mode control, which divides the sliding surface into two layers. In the first layer, a high-speed sliding surface is introduced, and the control system approaches the desired state quickly. In the second layer, a low-speed sliding surface is introduced, and the control system stabilizes near the desired state. The layered sliding mode control can improve the control accuracy and stability and also has a good effect on the response speed and robustness of the system. Therefore, hierarchical sliding mode control has the same drawback in that it is insensitive to disturbances, which can easily cause the “jitter” phenomenon of the system. To address the shortcomings of sliding mode control, this paper proposes an improved hierarchical sliding mode control method with adaptive exponential convergence law, which can adaptively adjust the control convergence law according to the control state and smooth the sign function, thus effectively improving the problem of the strong jitter of the traditional sliding mode control, and combining the nonlinear disturbance observer (NDO), which is the most powerful method for the control of sliding mode. The NDO can effectively solve the negative impact caused by the unknown disturbance and make the system more robust, and achieve an ideal control effect on the trajectory tracking ability of the WIPR system [24,25,26,27,28,29,30].
Overall, this paper includes the following four aspects: the first part constructs the dynamic model of the WIPR system, decouples the multi-coupled state variables, and facilitates the subsequent controller design; the second part establishes the kinematic trajectory tracking controller of the system and solves to obtain the desired speed of the dynamic control system. In the third part, an optimization model of the WIPR system combining nonlinear disturbance observer and hierarchical sliding mode control is designed, and the convergence of the nonlinear disturbance observer and the stability of the improved hierarchical sliding mode controller is demonstrated. The fourth part constructs the simulation model using the MATLAB/Simulink platform and conducts numerical simulation comparison experiments.
The contributions of this paper are as follows:
(1)
A wheeled inverted pendulum robot with a transport platform is envisioned for use in warehouses or other application scenarios to move goods.
(2)
The convergence law of hierarchical sliding mode control is improved to mitigate the jitter phenomenon of the sliding mode control system, and an adaptive function is introduced to minimize the system jitter.
(3)
By combining a nonlinear disturbance observer and hierarchical sliding mode control to estimate unknown external disturbances as input compensation, the system is made to control more accurately.

2. Materials and Methods

2.1. WIPR Model

A WIPR is a wheeled inverted pendulum transport robot with a placement table, as illustrated in Figure 1. Its left and right wheels are independent drive wheels that control the robot’s movement speed, rotation direction, and tilt angle of the pendulum using the principle of differential drive to manage the position and posture of WIPR. The generalized world coordinate system is denoted Σ O X Y Z while x , y , representing the center coordinate of the robot wheels. The robot’s forward velocity and rotational angular velocity are denoted as v and w , respectively. The angle of the robot’s direction of motion concerning the X -axis is represented by θ , while α is the tilt angle of the pendulum concerning the Z -axis. M refers to the total weight of the transport platform plus the pendulum, whereas m denotes the weight of each drive wheel. The distance between the two wheels is represented by d , while τ r , and τ l are the torque of the right wheel and the left wheel, respectively. The rotational inertia of each driven wheel is denoted by I w and I M represents the rotational inertia of the transport platform and the pendulum together. The length of the pendulum is represented by L . Detailed introduction of robot parameters can be seen in Table 1.
Remark 1.
The forward velocity of MWIPR is  x v , and  x v = x ˙ cos θ + y ˙ sin θ .
Assumption 1.
The tires of the MWIPR do not experience any skidding, and there is no potential for lateral deflection during its motion.
According to Assumption 1, the incomplete constraint equation of WIPR in Equation (1) can be listed as follows:
x ˙ sin θ y ˙ cos θ = 0 ,
The position and posture of the WIPR in the world coordinate system are represented by q = x , y , θ , α T . As the Lagrangian modeling method does not require the inclusion of internal forces within the system, it is a quick and straightforward method of building a model. This property makes it particularly well-suited for constructing multivariable and nonlinear dynamic models for the WIPR, as demonstrated in this paper. By dividing q into q m and α , the position of the robot in the coordinate system is denoted by q m , while the angle of the pendulum is represented by α . Therefore, the Lagrangian method [31] is employed to establish the dynamic model of the WIPR, and the resulting mathematical model is presented below as Equation (2).
M m q M m α q M α m q M α q q ¨ m α ¨ + C m q C m α q C α m q C α q q ˙ m α ˙ + G m G α = B m q m τ m 0 + A T q m λ 0 + d m d α ,
By defining A q m = sin θ , cos θ , 0 , the incompleteness constraint of Equation (1) yields the following result:
A q m q m = 0 ,
The WIPR system’s incomplete constraint force is A T q m λ , where λ is the Lagrange Multiplier. To eliminate the constraint forces in the system, we seek to find a matrix S q m 3 × 2 that satisfies S T q m A T q m = 0 .
By defining ν = v , w T , therefore, Equation (4) can be deduced.
q ˙ m = S q ν ,
To eliminate the incompetent constraint forces, a new vector ζ ˙ = ζ ˙ 1 , ζ ˙ 2 , ζ ˙ 3 T = v , w , α ˙ T is defined and used to transform the equation. The transformation involves multiplying both sides of the equation by a scalar S T q m , resulting in Equation (5):
M q ζ ¨ + C q , q ˙ ζ ˙ + G q = τ + τ d ,
The dynamics of the system can be described using the following equation, in which M q 3 × 3 represents the inertia matrix, C q , q ˙ q ˙ 3 × 3 is the Coriolis force matrix, G q 3 × 1 is the gravity matrix, τ is the control input matrix, and τ d is the total unknown disturbance. The detailed expressions of each vector or matrix are presented below.
M q = S T M m S S T M m α M α m S M α = m 11 0 m 13 0 m 22 0 m 31 0 m 33 ,   C q , q ˙ = S T C m S + S T M m S ˙ C m α C α m S + M α m S ˙ C α 0 0 c 13 0 c 22 c 23 0 c 32 0 ,   G q = S T G m G α = 0 0 g 3 = 0 0 M g L sin α , τ = S T B m q m τ m 0 = τ 1 τ 2 0 τ d = S T d m d α = d 1 d 2 d 3 .
The value of each variable in the expression is indicated as m 11 = 2 m + 2 I M / r 2 + M , m 13 = m 31 = M L cos α , m 22 = d 2 m / 2 + I M d 2 / 2 r 2 + I ω + M L 2 sin 2 α , m 33 = M L 2 + I M , c 22 = 1 / 2 M L 2 α ˙ sin 2 2 α , c 23 = 1 / 2 ω M L 2 sin 2 α , c 13 = M L α ˙ sin α , c 32 = 1 / 2 ω M L 2 sin 2 α .
Due to the coupling of the state variables in the system, Equation (5) is decoupled into Equation (6).
m 11 m 33 m 13 m 31 ζ ¨ 1 m 13 c 32 ζ ˙ 2 + m 33 c 13 ζ ˙ 3 m 13 g 3 = m 33 τ 1 + d 1 m 13 d 3 m 22 ζ ¨ 2 + c 22 ζ ˙ 2 + c 23 ζ ˙ 3 = τ 2 + d 2 m 11 m 33 m 13 m 31 ζ ¨ 3 + m 11 c 32 ζ ˙ 2 m 31 c 13 ζ ˙ 3 + m 11 g 3 = m 11 d 3 m 31 τ 1 + d 1 ,
By defining ξ 1 = ζ 1 ,   ζ 2 ,   ζ 3 T ,   ξ 2 = ζ ˙ 1 , ζ ˙ 2 , ζ ˙ 3 T , Equation (6) is converted to Equation (7).
ξ ˙ 1 = ξ 2 ξ ˙ 2 = f ξ + g ξ τ + D ,
where
f ξ = f 1 f 2 f 3 = 1 Ω m 13 c 32 ζ ˙ 2 m 33 c 13 ζ ˙ 3 m 13 g 3 1 m 22 c 22 ζ ˙ 2 c 23 ζ ˙ 3 1 Ω m 11 c 32 ζ ˙ 2 m 31 c 13 ζ ˙ 3 + m 11 g 3 ,   g ξ = m 33 Ω 1 m 22 m 31 Ω ,   D = D 1 D 2 D 3 = m 33 d 1 m 13 d 3 Ω d 2 m 22 m 11 d 3 m 31 d 1 Ω .
where Ω = m 11 m 33 m 13 m 31 .

2.2. The Design of the Kinematic Control Law

In kinematic trajectory tracking control for a WIPR, the system can be simplified to a general two-wheeled non-complete mobile robot for trajectory tracking. The process involves utilizing a reference trajectory state vector q m r = x r , y r , θ r T , an actual state vector q m = x , y , θ T , and a control objective designed to manage linear and angular velocities. The objective is to ensure that the actual robot travel trajectory aligns with the reference trajectory, even if the trajectory error q m e = e x , e y , e θ T = x r x , y r y , θ r θ T approaches zero. To meet the requirement expressed above, the control laws for v d and w d can be devised as Equation (8).
lim t q m e = 0 ,
The Lyapunov function is selected as Equation (9).
V 1 = 1 2 x e 2 + 1 2 y e 2 + 1 cos θ e ,
While the error persists, the value V 1 remains greater than zero, thereby rendering the function positive definite. Equation (10) describes the derivative of V 1 .
V ˙ 1 = e x e ˙ x + e y e ˙ y + e ˙ θ sin e θ = e x w e y v + v r cos e θ + e y w e x + v r sin e θ + w r w sin e θ = v e x w sin e θ + e x v r cos e θ + e y v r sin e θ + w r sin e θ 0 ,
Lyapunov’s stability theorem [32] establishes that the system can achieve asymptotic stability if the function is negative definite, i.e., if V ˙ 1 0 . Accordingly, the sought control law is as follows (11).
v d = v r cos e θ + k 1 e x w d = w r + v r e y + k 2 sin e θ ,
Both k 1 and k 2 are positive constants.
Substituting the control law into V ˙ 1 gives the following equation.
V ˙ 1 = e x v r cos e θ + k 1 e x w r + v r e y + k 2 sin e θ sin e θ + e x v r cos e θ + e y v r sin e θ + w r sin e θ = e x v r cos e θ + e x v r cos e θ k 1 e x 2 w r sin e θ + w r sin e θ v r e y sin e θ + e y v r sin e θ k 2 sin e θ 2 = k 1 e x 2 k 2 sin e θ 2 0 .
Therefore, V ˙ 1 0 can be proven.
So far, the desired velocity required for the design of the dynamical system is shown in Equation (11), and the velocity tracking problem of the dynamical system and the angle tracking problem of the pendulum will be solved next.

2.3. The Design of NDO

A nonlinear disturbance observer is developed to estimate the actual disturbance in the system for an unknown disturbance D, thereby strengthening the system’s robustness. To address practical considerations, it is assumed that any disturbance is bounded as follows [33,34,35].
Lemma 1.
For initial conditions that are bounded, a Liapunov function is also uniformly bounded  x t  if there exists a continuous positive definite Liapunov function  V x  satisfying the following conditions:
δ 1 x V x δ 2 x ,   V ˙ x κ V x + c ,
where  δ 1 ,   δ 2 : n  is the  V  class function, and  κ ,   c  all are positive constants.
Assumption 2.
Since no disturbance can be infinite in the real world, we assume that the perturbations in the WIPR system studied in this paper are all bounded, and their first-order derivatives and second-order derivatives are assumed to be bounded; thus, the following equations can be obtained.
D ˙ ε 1 ,   D ¨ ε 2 ,   ε 1 > 0 , ε 2 > 0
represents the Euclidean norm of the vector. The NDO is designed as in Equation (14).
D ^ = Z 1 + P 1 ξ 2 Z ˙ 1 = L 1 ξ 2 f ξ + g ξ τ + D ^ + D ˙ ^ D ˙ ^ = Z 2 + P 2 ξ 2 Z ˙ 2 = L 2 ξ 2 f ξ + g ξ τ + D ^ ,
D ^ and D ˙ ^ represent the estimates of the total perturbation and its derivative, respectively, while Z 1 and Z 2 are intermediate variables in the observer. To meet the requirements of the system, the self-designed nonlinear functions P 1 ξ 2 and P 2 ξ 2 are used and must satisfy the conditions L 1 ξ 2 = P 1 ξ 2 ξ 2 and L 2 ξ 2 = P 2 ξ 2 ξ 2 .
Property 1.
The error between the estimated and actual values of the disturbance is represented by  D ˜ = D D ^ , while D ˙ ˜ = D ˙ D ˙ ^ represents the error between the derivative of the actual value of a perturbation and the derivative of the estimated value of the same perturbation.
The derivation of D ˜ and D ˙ ˜ substitution of Equations (7) and (14) into the above equation leads to results D ˜ ˙ and D ˙ ˜ , which are the equations of the NDO, (15) and (16), respectively.
D ˜ ˙ = D ˙ D ^ ˙ = D ˙ Z ˙ 1 P ξ 2 ξ 2 d ξ 2 d t = D ˙ Z ˙ 1 L 1 ξ 2 ξ ˙ 2 = D ˙ + L 1 ξ 2 f ξ + g ξ τ + D ^ D ˙ ^ L 1 ξ 2 ξ ˙ 2 D ˙ + L 1 ξ 2 f ξ + g ξ τ + D ^ D ˙ ^ L 1 ξ 2 f ξ + g ξ τ + D = L 1 ξ 2 D ˜ + D ˙ ˜ ,
D ˙ ˜ ˙ = D ¨ D ˙ ^ ˙ = D ˙ Z ˙ 2 P ξ 2 ξ 2 d ξ 2 d t = D ¨ Z ˙ 2 L 2 ξ 2 ξ ˙ 2 = D ¨ + L 2 ξ 2 f ξ + g ξ τ + D ^ L 2 ξ 2 f ξ + g ξ τ + D = L 2 ξ 2 D ˜ + D ¨ ,
By letting E = D ˜ T ,   D ˙ ˜ T T and substituting the appropriate Equations (15)–(17) can be obtained.
E ˙ = L E + δ D ¨ ,
where L = L 1 ξ 2 I 3 L 2 ξ 2 0 , δ = 0 I 3 . The observer’s stability is examined, and a Lyapunov function is selected to make sure that it can reliably predict the system state despite any nonlinear disturbances. By selecting an appropriate Liapunov function, we rigorously prove the stability of the observer and the precision of its estimation precisely.
Property 2.
L  is a skew-symmetric matrix.
V 2 = 1 2 E T E ,
The proof of the derivative of V 2 can be expressed as follows:
V ˙ 2 = E T E ˙ = E T L E + δ D ¨ E T L + 0.5 δ 2 I 6 E + 0.5 ε 2 2 ,
The above design of an NDO for a WIPR can be summarized in the following theorem.
Theorem 1.
For the existence of an unknown disturbance in a WIPR system, the perturbation estimation error is bounded for the observer designed according to Equation (14).
Proof. 
The design parameters L 1 ξ 2 = P 1 ξ 2 ξ 2 and L 2 ξ 2 = P 2 ξ 2 ξ 2 are such that L + 0.5 δ 2 I 6 is a negative definite matrix, according to Lemma 1, then E is bounded. □

2.4. The Design of Improved Slide Mode Control

The sliding mode control algorithm consists of two key elements: (i) the design of the sliding mode surface; and (ii) the design of the convergence rate. The design of the sliding mode surface is mainly based on the system structure as well as the control objective. As for the design of convergence law, there are four different convergence laws: the isokinetic convergence law, exponential convergence law, power convergence law, and general convergence law. In this paper, based on the optimal control objective of WIPPR to cope with nonlinearity and underdrive, as well as unknown disturbances, the traditional exponential convergence law is improved by introducing an adaptive control function, and an improved sliding mode control based on the adaptive exponential convergence law is proposed, which can weaken the system jitter while speeding up the system response, making the sliding mode control more suitable for tracking the reference trajectory of the WIPR system under the action of unknown disturbances [36]. The design for the traditional exponential convergence law is shown in Equation (20).
s ˙ = η sgn s μ s ,
where: s denotes the slip surface function; the parameters η , μ denote the convergence coefficient; and sgn s denotes the sign function.
In the traditional exponential convergence law, the isokinetic term is denoted η sgn s , and the exponential term is denoted μ s . When the state of the system is far from the slip surface, the exponential term and the isokinetic term in the convergence law act simultaneously to help the system move toward the slip surface, and the magnitude of the isokinetic term and the exponential term are mainly determined by the reference η , μ . The exponential term is small, and the isokinetic term acts mainly when the system is moving close to the surface.
This paper makes a corresponding improvement based on the traditional exponential convergence law and introduces the adaptive o s function to adjust the convergence law in accordance with the control state of the system, as shown in Figure 2, which can accelerate the convergence speed of the sliding mode and weaken the overshoot phenomenon. This allows the sliding mode control to reduce the jitter phenomenon of the system. Following the inclusion of the adaptive function o s , the new exponential convergence law is as follows:
s ˙ = η o s sgn s μ s o s = | s | + 1 log a | s | + a | s | + log a | s | + a + b ,
where a > 0 , b > 0 , η > 0 , μ > 0 .
Through the analysis, relative to not adding the adaptive function (i.e., o(s) = 1), it can be found that when the system motion point is far away from the sliding surface (namely, when s is far away from the origin 0), the adaptive function o s will increase the convergence law, which will speed up the system convergence speed, shorten the system state convergence time to the target state, and reduce the control time; when the system motion point is close to the sliding surface, s will converge to 0 and o s will be less than 1. The role of o s here is to suppress the jitter amplitude and weaken the state variable fluctuation problem after the system is stabilized, and the suppression effect will be more obvious as the parameter b increases. To further weaken the jitter problem of the stabilized system, the smoothing process is carried out for the symbolic function sgn s in this paper, which is known as the traditional symbolic function [37], as shown in the following equation.
sgn s = 1     s > 0 0     s = 0 1    s < 0 ,
The symbolic function after the smoothing process is shown below.
sgn s = s s + β ,   β = 0.01 ,

2.5. The Design of the Forward-Rotation Subsystem

For convenience, the system has been reorganized into the following form (24).
ζ ¨ 1 = 1 Ω m 13 c 32 ζ ˙ 2 m 33 c 13 ζ ˙ 3 m 13 g 3 + m 33 Ω τ 1 + D 1 ζ ¨ 2 = 1 m 22 c 22 ζ ˙ 2 c 23 ζ ˙ 3 + 1 m 22 τ 2 + D 2 ,
The state variables in the system described by Equation (24) are highly coupled. To address this issue and to expand the system’s asymptotic stability domain, a hierarchical sliding mode controller was designed. The controller’s primary objective is to utilize an input control law that can simultaneously control both system variables ζ 1 and ζ 2 , thereby, mitigating the problem of system coupling [38].
Having obtained the expected forward velocity ( v d ) and angular velocity ( w d ) from Equation (11), the error between the actual and expected values can be defined as follows:
e x v = ζ 1 d ζ 1 = x r cos θ r + y r sin θ r x v e v = ζ ˙ 1 d ζ ˙ 1 = v d v e θ = ζ 2 d ζ 2 = θ r θ e w = ζ ˙ 2 d ζ ˙ 2 = w d w ,
To design the sliding mode control error tracking scheme for the v - w subsystem, two mutually independent first-layer sliding mode surfaces were initially constructed. The equations used to create these slide surfaces are as follows:
s 1 = γ 1 e x v + e v ,   γ 1 > 0 ,
s 2 = γ 2 e θ + e w ,   γ 2 > 0 ,
The results of deriving Equations (26) and (27) are presented below.
s ˙ 1 = γ 1 e ˙ x v + e ˙ v = γ 1 e ˙ x v + ζ ¨ 1 d ζ ¨ 1 = 0 s ˙ 2 = γ 2 e ˙ θ + e ˙ w = γ 2 e ˙ θ + ζ ¨ 2 d ζ ¨ 2 = 0 ,
According to Filippov’s equivalent control theory, the equivalent control laws for ζ 1 and ζ 2 are as follows:
τ e q 1 = Δ m 33 γ 1 e ˙ x v + ζ ¨ 1 d + D ^ 1 + 1 m 33 m 13 c 32 ζ ˙ 2 m 33 c 13 ζ ˙ 3 m 13 g 3 τ e q 2 = m 22 γ 2 e ˙ θ + ζ ¨ 2 d + c 22 ζ ˙ 2 + c 23 ζ ˙ 3 D ^ 2
The second sliding surface can be expressed as a linear combination of the first sliding surface.
s 3 = γ 3 s 1 + γ 4 s 2 , γ 3 , γ 4 > 0 ,
To control ζ 1 and ζ 2 , the equivalent control law must be included at the same time to control and enter their designed sliding surface, respectively. Therefore, the total control law is shown in the following equation.
τ u = τ e q 1 + τ e q 2 + τ s w ,
τ s w is the switching law of the converging slide surface phase, and the expressions are as follows.
τ s w = γ 3 m 22 Ω τ e q 2 + γ 4 1 m 22 τ e q 1 + η 1 o s 3 sgn s 3 + μ 1 s 3 γ 3 m 22 Ω + γ 4 1 m 22 ,
To mitigate the jitter phenomenon of the system, the isokinetic and exponential terms of the sliding mode control are improved, where η 1 and μ 1 are the isokinetic and exponential terms of the previous design convergence law, and both are positive constants; o s 3 and sgn s 3 are the adaptive and symbolic functions designed in the previous paper. To prove that the designed controller is stable, the Lyapunov function is chosen as follows.
V 3 = 1 2 s 3 2 > 0
The derivative of V 3 for time is given by the following expression.
V ˙ 3 = s 3 s ˙ 3 = s 3 γ 3 s ˙ 1 + γ 4 s ˙ 2 = s 3 γ 3 γ 1 e ˙ x v + e ˙ v + γ 4 γ 2 e ˙ θ + e ˙ w = s 3 γ 3 γ 1 e ˙ x v + v ˙ d v ˙ + γ 4 γ 2 e ˙ θ + w ˙ d w ˙ = s 3 γ 3 γ 1 e ˙ x v + v ˙ d f 1 m 33 Ω τ e q 1 + τ e q 2 + τ s w D 1 + γ 3 γ 2 e ˙ θ + w ˙ d f 2 1 m 22 τ e q 1 + τ e q 2 + τ s w D 2 = s 3 η 1 o s 3 sgn s 3 μ 1 s 3 2 + γ 1 D ^ 1 D 1 + γ 2 D ^ 2 D 2 η 1 o s 3 s 3 μ 1 s 3 2 + s 3 γ 1 D ˜ 1 + γ 2 D ˜ 2
By design η 1 o s 3 > γ 1 D ˜ 1 + γ 2 D ˜ 2 , the result of the following equation can be obtained.
V ˙ 3 η 1 s 3 2 η 1 2 V 3
From Lemma 1 in [38], the following equation can be obtained.
V 3 t e η 1 2 t t 0 V t 0
It can be seen that the V 3 t index converges to 0, and the rate of convergence depends on η 1 .
As demonstrated by the preceding equation, the error state can attain the slip surface in a finite amount of time. Subsequently, the first layer of slip surfaces s 1 and s 2 can converge asymptotically to zero, leading to the convergence of both the rotational and forward velocities of WIPR to the desired values.

2.6. The Design of the Tilt-Angle Subsystem

The system discussed in the previous section can achieve complete tracking of ζ 1 and ζ 2 within a finite time, which enables us to transform it into the following form:
ζ ¨ 3 = 1 Ω m 11 c 32 ζ ˙ 2 d m 31 c 13 ζ ˙ 3 + m 11 g 3 m 31 Ω τ 3 + D 3 ,
As WIPR aims to maintain a vertical and stable direction of the pendulum during its motion, all relevant parameters ( α d , α ˙ d , α ¨ d ) can be set to zero. As such, the following definitions can be employed:
e α = α d α = α e α ˙ = α ˙ d α ˙ = α ˙ ,
Let the sliding mode surface be defined as Equation (39), with its derivative expressed as Equation (40).
s 4 = γ 5 e α + e α ˙ ,
s ˙ 4 = γ 5 e ˙ α + e ˙ α ˙ = γ 5 α ˙ α ¨ = η 2 o s 4 sgn s 4 μ 2 s 4 ,
After substituting Equation (37), the control law for the tilt angle subsystem can be derived as presented in Equation (41).
τ α = Φ m 31 γ 5 α ˙ + D ^ 3 1 m 31 ( m 11 c 12 ζ ˙ 2 m 31 c 13 ζ ˙ 3 + m 11 g 3 η 2 o s 4 sgn s 4 ) μ 2 s 4
To prove the stability of the designed system, the Lyapunov function is chosen as follows.
V 4 = 1 2 s 4 2
The derivative of V 4 for time is given by the following expression.
V ˙ 4 = s 4 s ˙ 4 = s 4 γ 5 α ˙ α ¨ = s 4 γ 5 α ˙ f 3 1 m 22 τ α D 3 = s 4 η 2 o s 4 sgn s 4 μ 2 s 4 2 + γ 5 D ^ 3 D 3 η 2 o s 4 s 4 μ 2 s 4 2 + γ 5 D ˜ 3
By choosing η 2 o s 4 > γ 5 D ˜ 3 , it enables V ˙ 4 < 0 to hold, indicating that the system achieves asymptotic stability.
Figure 3 displays the schematic block diagram of the control system.

3. Simulation

The focus of this section is to discuss the trajectory-tracking effect of the system in a simulation environment, and to verify the feasibility of the proposed control scheme and what the advantages of the proposed method are compared with other control systems in this paper. Next, the simulation results of different control systems in the face of the same disturbance will be compared to verify the control effectiveness of each system. The parameters in the system are shown in the following Table 2.
The simulation experiments in Matlab/Simulink verified the high-precision trajectory tracking capability of the system and the stability of the pendulum in robot motion. During the simulation study, the initial position was set as q m = 0 ,   0 ,   0 T , and the initial position of the reference trajectory was set as q m r = 0 , 0 , π / 2 T , where the desired tracking velocity = 1 m/s and the angular velocity = 1 rad/s. Therefore, the trajectory of the robot should be a circle with a radius of 1 m, and the center of the circle is 0 , 0 . Therefore, the time function of the reference trajectory was chosen as x = sin t , y = cos t .
An external perturbation was added, as shown in the following equation.
D 1 = 0.4 sin 0.4 t N m D 2 = 0.5 cos 0.5 t N m D 3 = 0.6 sin 0.6 t N m ,
To demonstrate the superiority of the proposed method in this paper, three comparative experiments were conducted under the given disturbance conditions: the first experiment involved the simulation results of the unimproved HSMC method, the second experiment involved the simulation results of the improved IHSMC method with adaptive law but without nonlinear disturbance observer, and the third experiment involved the simulation results of the proposed method in this paper (referred to as PC).
First of all, by observing Figure 4, Figure 5 and Figure 6, it can be concluded that the proposed method in this paper is better than the other two control methods in terms of both the speed of convergence of the error to the steady state and the magnitude of the fluctuation of the error after reaching the steady state when compared with the other two methods. This undoubtedly reflects the effectiveness of the method in this paper, which can track the given reference trajectory very accurately.
Figure 7 and Figure 8 give the tracking of the desired speed of the WIPR system under the three control methods. Compared with the other two methods, firstly, the control method in this paper can track the desired speed more rapidly, reaching the effect of tracking the desired speed at 0.7 s, whereas the other two methods track the desired speed in more than 1 s, which is much slower than the method in this paper, and the fluctuation frequency is high, which may affect the stability of the WIPR. As can be seen from Figure 8, the present method exceeds the other two methods in the tracking effect of rotational velocity relative to the forward velocity, for one. The convergence speed is fast, and more importantly, the proposed method is very stable after the velocity tracking reaches the steady state, which can be regarded as showing no fluctuation compared with the other two methods.
The HSMC method with general exponential convergence law has more frequent angle oscillations, and the system is more unstable, as can be seen from the angle change of the WIPR pendulum shown in Figure 9, whereas the IHSMC improved convergence law method’s pendulum has smoother oscillations after reaching stability, and the control effect is obviously stronger than that of the HSMC with general exponential convergence law. In terms of response time and maximum overshoot, the suggested method outperforms the other two ways, and it can continue to operate smoothly and without oscillations once it has reached the stabilization point.
The variograms of the input torque for the three control methods are presented in Figure 10, Figure 11 and Figure 12. The results indicate that when the convergence law of HSMC follows the general exponential convergence law, the jitter vibration of the input torque for the left and right wheels of WIPR is evident, which adversely affects the output of the actuator (i.e., affects the output of the drive motors of the left and right wheels). In contrast, Figure 11 illustrates that the improved convergence law significantly reduces the jitter phenomenon, resulting in a more beneficial improvement for the actuator.
Figure 12 presents the variation of input torque under the proposed control method. It can be observed that the input torque obtained by this method is smoother than IHSMC, and the jitter suppression effect is more satisfactory. This approach achieves a better torque input graph, making it the most effective method among the three for actuator benefits. Therefore, the proposed control method demonstrates superior performance in terms of reducing jitter and enhancing actuator benefits compared with the other two control methods, making it a better solution.
Figure 13 shows the trajectory tracking diagrams of different control systems. From an intuitive point of view, the proposed scheme is also significantly better than the other two schemes. As shown by the simulation comparison experiment, the method proposed in this paper is feasible, and its effect is excellent.

4. Conclusions

The purpose of this paper was to study the trajectory-tracking problem for WIPRs and propose a hierarchical sliding mode controller with a nonlinear perturbation observer to achieve accurate control of the reference trajectory and maintain the pendulum stability during motion. A nonlinear disturbance observer was designed to make the system more robust to unknown external disturbances. The underdriven coupling of WIPRs was addressed by dividing the system into two subsystems through the decoupling of its control state variables. The hierarchical sliding mode control method with an improved convergence law was then applied to control the system and suppress the “jitter” phenomenon. Finally, the Lyapunov function was chosen to verify the stability of the system mathematically.
The feasibility of the control system was verified using simulation software. However, considering the complexity of the real-world environment and external uncertainty, future work will focus on building a hardware system for the robot to study the real effects of the control method of WIPRs in the real world.

Author Contributions

All of the nominated authors initially made significant contributions to the paper. Conceptualization, M.H. and X.Z.; methodology, X.Z.; software, X.Z.; validation, X.Z., D.C. and M.H.; formal analysis, M.H.; investigation, M.H.; resources, X.Z.; data curation, X.Z., D.C. and Z.X; writing—original draft preparation, M.H.; writing—review and editing, X.Z., D.C. and Z.X.; visualization, D.C.; supervision, M.H.; project administration, M.H.; funding acquisition, M.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (61971048) and the Key Incubation Project of Beijing University of Information Science and Technology (5212110927).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to internal laboratory confidentiality agreements.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kim, S.; Kwon, S. Nonlinear Optimal Control Design for Underactuated Two-Wheeled Inverted Pendulum Mobile Platform. IEEE/ASME Trans. Mechatron. 2017, 22, 2803–2808. [Google Scholar] [CrossRef]
  2. Chen, M. Robust tracking control for self-balancing mobile robots using disturbance observer. IEEE/CAA J. Autom. Sin. 2017, 4, 458–465. [Google Scholar] [CrossRef]
  3. Kim, Y.; Kwon, S. Robust Stabilization of Underactuated Two-Wheeled Balancing Vehicles on Uncertain Terrains with Nonlinear-Model-Based Disturbance Compensation. Actuators 2022, 11, 339. [Google Scholar] [CrossRef]
  4. Liu, J.; Vazquez, S.; Wu, L.; Marquez, A.; Gao, H.; Franquelo, L.G. Extended State Observer-Based Sliding-Mode Control for Three-Phase Power Converters. IEEE Trans. Ind. Electron. 2017, 64, 22–31. [Google Scholar] [CrossRef] [Green Version]
  5. Yang, C.; Li, Z.; Li, J. Trajectory Planning and Optimized Adaptive Control for a Class of Wheeled Inverted Pendulum Vehicle Models. IEEE Trans. Cybern. 2013, 43, 24–36. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Mirzaeinejad, H. Optimization-based nonlinear control laws with increased robustness for trajectory tracking of non-holonomic wheeled mobile robots. Transp. Res. Part C Emerg. Technol. 2019, 101, 1–17. [Google Scholar] [CrossRef]
  7. Wang, F.-C.; Chen, Y.-H.; Wang, Z.-J.; Liu, C.-H.; Lin, P.-C.; Yen, J.-Y. Decoupled Multi-Loop Robust Control for a Walk-Assistance Robot Employing a Two-Wheeled Inverted Pendulum. Machines 2021, 9, 205. [Google Scholar] [CrossRef]
  8. Yue, M.; Ning, Y.; Yu, S.; Zhang, Y. Composite following control for wheeled inverted pendulum vehicles based on hu-man-robot interaction. Sci. China Inf. Sci. 2019, 62, 50206. [Google Scholar]
  9. Watson, M.T.; Gladwin, D.T.; Prescott, T.J.; Conran, S.O. Dual-Mode Model Predictive Control of an Omnidirectional Wheeled Inverted Pendulum. IEEE/ASME Trans. Mechatron. 2019, 24, 2964–2975. [Google Scholar] [CrossRef]
  10. Albert, K.; Phogat, K.S.; Anhalt, F.; Banavar, R.N.; Chatterjee, D.; Lohmann, B. Structure-Preserving Constrained Optimal Trajectory Planning of a Wheeled Inverted Pendulum. IEEE Trans. Robot. 2020, 36, 910–923. [Google Scholar] [CrossRef]
  11. Ginoya, D.; Shendge, P.D.; Phadke, S.B. Sliding Mode Control for Mismatched Uncertain Systems Using an Extended Disturbance Observer. IEEE Trans. Ind. Electron. 2014, 61, 1983–1992. [Google Scholar] [CrossRef]
  12. Pathak, K.; Franch, J.; Agrawal, S. Velocity and position control of a wheeled inverted pendulum by partial feedback linearization. IEEE Trans. Robot. 2005, 21, 505–513. [Google Scholar] [CrossRef]
  13. Takei, T.; Imamura, R.; Yuta, S. Baggage Transportation and Navigation by a Wheeled Inverted Pendulum Mobile Robot. IEEE Trans. Ind. Electron. 2009, 56, 3985–3994. [Google Scholar] [CrossRef]
  14. Gong, S.; Zhang, A.; She, J.; Zhang, X.; Liu, Y. Trajectory Design and Tracking Control for Nonlinear Underactuated Wheeled Inverted Pendulum. Math. Probl. Eng. 2018, 2018, e6134764. [Google Scholar] [CrossRef]
  15. Huang, C.-F.; Yeh, T.-J. Anti Slip Balancing Control for Wheeled Inverted Pendulum Vehicles. IEEE Trans. Control Syst. Technol. 2019, 28, 1042–1049. [Google Scholar] [CrossRef]
  16. Yang, C.; Li, Z.; Cui, R.; Xu, B. Neural Network-Based Motion Control of an Underactuated Wheeled Inverted Pendulum Model. IEEE Trans. Neural Netw. Learn. Syst. 2014, 25, 2004–2016. [Google Scholar] [CrossRef] [PubMed]
  17. Ren, C.; Li, X.; Yang, X.; Ma, S. Extended State Observer-Based Sliding Mode Control of an Omnidirectional Mobile Robot With Friction Compensation. IEEE Trans. Ind. Electron. 2019, 66, 9480–9489. [Google Scholar] [CrossRef]
  18. Chen, L.; Wang, H.; Huang, Y.; Ping, Z.; Yu, M.; Zheng, X.; Ye, M.; Hu, Y. Robust hierarchical sliding mode control of a two-wheeled self-balancing vehicle using perturbation estimation. Mech. Syst. Signal Process. 2020, 139, 106584. [Google Scholar] [CrossRef]
  19. Zhou, Y.; Wang, Z. Robust motion control of a two-wheeled inverted pendulum with an input delay based on optimal integral sliding mode manifold. Nonlinear Dyn. 2016, 85, 2065–2074. [Google Scholar] [CrossRef]
  20. Adhikary, N.; Mahanta, C. Integral backstepping sliding mode control for underactuated systems: Swing-up and stabilization of the Cart–Pendulum System. ISA Trans. 2013, 52, 870–880. [Google Scholar] [CrossRef]
  21. Ri, S.; Huang, J.; Wang, Y.; Kim, M.; An, S. Terminal Sliding Mode Control of Mobile Wheeled Inverted Pendulum System with Nonlinear Disturbance Observer. Math. Probl. Eng. 2014, 2014, e284216. [Google Scholar] [CrossRef] [Green Version]
  22. Ping, H.; Hai, W.; Linfeng, L.; Huifang, K.; Ming, Y.; Canghua, J.; Zhihong, M. A novel hierarchical sliding mode control strategy for a two-wheeled self-balancing vehicle. In Proceedings of the 2017 36th Chinese Control Conference (CCC), Dalian, China, 26–28 July 2017; pp. 3731–3736. [Google Scholar] [CrossRef]
  23. Moness, M.; Mahmoud, D.; Hussein, A. Real-time Mamdani-like fuzzy and fusion-based fuzzy controllers for balancing two-wheeled inverted pendulum. J. Ambient. Intell. Humaniz. Comput. 2020, 13, 3577–3593. [Google Scholar] [CrossRef]
  24. Sun, W.; Su, S.-F.; Xia, J.; Wu, Y. Adaptive Tracking Control of Wheeled Inverted Pendulums with Periodic Disturbances. IEEE Trans. Cybern. 2020, 50, 1867–1876. [Google Scholar] [CrossRef] [PubMed]
  25. Palm, R. Sliding mode fuzzy control. In Proceedings of the 1992 IEEE International Conference on Fuzzy Systems, San Diego, CA, USA, 8–12 March 1992; pp. 519–526. [Google Scholar] [CrossRef]
  26. Jmel, I.; Dimassi, H.; Hadj-Said, S.; M’Sahli, F. An adaptive sliding mode observer for inverted pendulum under mass variation and disturbances with experimental validation. ISA Trans. 2020, 102, 264–279. [Google Scholar] [CrossRef] [PubMed]
  27. Huang, J.; Guan, Z.; Matsuno, T.; Fukuda, T.; Sekiyama, K. Sliding-Mode Velocity Control of Mobile-Wheeled Inverted-Pendulum Systems. IEEE Trans. Robot. 2010, 26, 750–758. [Google Scholar] [CrossRef]
  28. Fukushima, H.; Muro, K.; Matsuno, F. Sliding-Mode Control for Transformation to an Inverted Pendulum Mode of a Mobile Robot with Wheel-Arms. IEEE Trans. Ind. Electron. 2015, 62, 4257–4266. [Google Scholar] [CrossRef] [Green Version]
  29. Chen, X.; Komada, S.; Fukuda, T. Design of a nonlinear disturbance observer. IEEE Trans. Ind. Electron. 2000, 47, 429–437. [Google Scholar] [CrossRef]
  30. Saradagi, A.; Muralidharan, V.; Krishnan, V.; Menta, S.; Mahindrakar, A.D. Formation Control and Trajectory Tracking of Nonholonomic Mobile Robots. IEEE Trans. Control Syst. Technol. 2017, 26, 2250–2258. [Google Scholar] [CrossRef]
  31. Yoshida, K.; Sekikawa, M.; Hosomi, K. Nonlinear analysis on purely mechanical stabilization of a wheeled inverted pendulum on a slope. Nonlinear Dyn. 2016, 83, 905–917. [Google Scholar] [CrossRef] [Green Version]
  32. Chen, W.-H.; Yang, J.; Guo, L.; Li, S. Disturbance-Observer-Based Control and Related Methods—An Overview. IEEE Trans. Ind. Electron. 2016, 63, 1083–1095. [Google Scholar] [CrossRef] [Green Version]
  33. Huang, J.; Ri, S.; Liu, L.; Wang, Y.; Kim, J.; Pak, G. Nonlinear Disturbance Observer-Based Dynamic Surface Control of Mobile Wheeled Inverted Pendulum. IEEE Trans. Control Syst. Technol. 2015, 23, 2400–2407. [Google Scholar] [CrossRef]
  34. Xie, L.; Yu, X. State Observer Based Robust Backstepping Fault-Tolerant Control of the Free-Floating Flexible-Joint Space Manipulator. Appl. Sci. 2023, 13, 2634. [Google Scholar] [CrossRef]
  35. Yue, M.; Wei, X.; Li, Z. Adaptive sliding-mode control for two-wheeled inverted pendulum vehicle based on zero-dynamics theory. Nonlinear Dyn. 2014, 76, 459–471. [Google Scholar] [CrossRef]
  36. Chang, W.-J.; Hsu, F.-L. Sliding mode fuzzy control for Takagi–Sugeno fuzzy systems with bilinear consequent part subject to multiple constraints. Inf. Sci. 2016, 327, 258–271. [Google Scholar] [CrossRef]
  37. Xu, J.-X.; Guo, Z.-Q.; Lee, T.H. Design and Implementation of Integral Sliding-Mode Control on an Underactuated Two-Wheeled Mobile Robot. IEEE Trans. Ind. Electron. 2014, 61, 3671–3681. [Google Scholar] [CrossRef]
  38. Liu, Y.; Jing, Y.W.; Liu, X.P.; Li, X.H. Survey on finite-time control for nonlinear systems. Control Theory Appl. 2020, 37, 1–12. [Google Scholar]
Figure 1. WIPR system.
Figure 1. WIPR system.
Applsci 13 04350 g001
Figure 2. The adaptive function o s .
Figure 2. The adaptive function o s .
Applsci 13 04350 g002
Figure 3. The control system.
Figure 3. The control system.
Applsci 13 04350 g003
Figure 4. The error of x.
Figure 4. The error of x.
Applsci 13 04350 g004
Figure 5. The error of y.
Figure 5. The error of y.
Applsci 13 04350 g005
Figure 6. The error of θ.
Figure 6. The error of θ.
Applsci 13 04350 g006
Figure 7. WIPR forward velocity v .
Figure 7. WIPR forward velocity v .
Applsci 13 04350 g007
Figure 8. WIPR rotation velocity w .
Figure 8. WIPR rotation velocity w .
Applsci 13 04350 g008
Figure 9. The angle of the WIPR pendulum α .
Figure 9. The angle of the WIPR pendulum α .
Applsci 13 04350 g009
Figure 10. Input torque under HSMC.
Figure 10. Input torque under HSMC.
Applsci 13 04350 g010
Figure 11. Input torque under IHSMC.
Figure 11. Input torque under IHSMC.
Applsci 13 04350 g011
Figure 12. Input torque under PC.
Figure 12. Input torque under PC.
Applsci 13 04350 g012
Figure 13. Tracking of circular trajectories.
Figure 13. Tracking of circular trajectories.
Applsci 13 04350 g013
Table 1. Parameter descriptions.
Table 1. Parameter descriptions.
ParameterDescription
m w Mass of each wheel
M The total weight of the transport platform plus the pendulum
I w The rotational inertia of each driven wheel
I M The rotational inertia of the transport platform and the pendulum
d The distance between the two wheels
L The length of the pendulum
τ l The torque of the left wheel
τ r The torque of the right wheel
v WIPR forward velocity
w WIPR rotation velocity
θ WIPR Yaw angle
α The tilt angle of the pendulum
Table 2. The value of each parameter variable in the system.
Table 2. The value of each parameter variable in the system.
Parameter (Unit)Value
M kg 8
m kg 0.5
I M kg m 2 5
I w kg m 2 0.3
d m 0.5
L m 0.5
r m 0.1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hou, M.; Zhang, X.; Chen, D.; Xu, Z. Hierarchical Sliding Mode Control Combined with Nonlinear Disturbance Observer for Wheeled Inverted Pendulum Robot Trajectory Tracking. Appl. Sci. 2023, 13, 4350. https://doi.org/10.3390/app13074350

AMA Style

Hou M, Zhang X, Chen D, Xu Z. Hierarchical Sliding Mode Control Combined with Nonlinear Disturbance Observer for Wheeled Inverted Pendulum Robot Trajectory Tracking. Applied Sciences. 2023; 13(7):4350. https://doi.org/10.3390/app13074350

Chicago/Turabian Style

Hou, Ming, Xuedong Zhang, Du Chen, and Zheng Xu. 2023. "Hierarchical Sliding Mode Control Combined with Nonlinear Disturbance Observer for Wheeled Inverted Pendulum Robot Trajectory Tracking" Applied Sciences 13, no. 7: 4350. https://doi.org/10.3390/app13074350

APA Style

Hou, M., Zhang, X., Chen, D., & Xu, Z. (2023). Hierarchical Sliding Mode Control Combined with Nonlinear Disturbance Observer for Wheeled Inverted Pendulum Robot Trajectory Tracking. Applied Sciences, 13(7), 4350. https://doi.org/10.3390/app13074350

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