CN113377016B - Slip-considered multi-mobile-robot cooperative formation buffeting-free sliding mode control method - Google Patents

Slip-considered multi-mobile-robot cooperative formation buffeting-free sliding mode control method Download PDF

Info

Publication number
CN113377016B
CN113377016B CN202110805055.XA CN202110805055A CN113377016B CN 113377016 B CN113377016 B CN 113377016B CN 202110805055 A CN202110805055 A CN 202110805055A CN 113377016 B CN113377016 B CN 113377016B
Authority
CN
China
Prior art keywords
formation
mobile robot
formula
mobile
robot
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110805055.XA
Other languages
Chinese (zh)
Other versions
CN113377016A (en
Inventor
郭一军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huangshan University
Original Assignee
Huangshan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Huangshan University filed Critical Huangshan University
Priority to CN202110805055.XA priority Critical patent/CN113377016B/en
Publication of CN113377016A publication Critical patent/CN113377016A/en
Application granted granted Critical
Publication of CN113377016B publication Critical patent/CN113377016B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

The invention discloses a slip-form control method for multi-mobile-robot cooperative formation without buffeting in consideration of slipping, and solves the formation control problem of multi-mobile-robots under the condition of slipping. The method comprises the following specific steps: establishing a mathematical model of a multi-mobile-robot formation control system considering a slipping condition; designing a formation buffeting-free sliding mode controller based on the established formation system mathematical model, and realizing that each mobile robot in the formation can form a formation with a pilot robot according to a preset formation; and (5) stability analysis of the formation controller. The multi-mobile-robot formation control method can eliminate the adverse effect of wheel slip on the formation control effect, and can ensure the rapid and stable control of formation even under the condition of slip.

Description

Slip-considered multi-mobile-robot cooperative formation buffeting-free sliding mode control method
Technical Field
The invention relates to the field of multi-mobile robot control, in particular to a slip-form control method for multi-mobile robot cooperative formation without buffeting, which considers slip.
Background
In order to complete related tasks, the efficiency and the success rate of task completion can be greatly improved through the cooperative cooperation of the multiple mobile robots. For example, when complex tasks are executed, the tasks can be decomposed and then completed by different mobile robots, so that the tasks can work in parallel, and the time for completing the tasks can be greatly shortened. In addition, when a certain mobile robot or some mobile robots are in fault in the task execution process, other mobile robots can still complete the given task through mutual cooperation. Based on this, the problem of cooperative control of multiple mobile robots has become a hot research problem in the field of robot application. As one of the important research directions of cooperative control, formation control of multiple mobile robots also has an important research value.
With the wide application of mobile robots, the working environment is more complicated, and the multi-mobile robot formation system is interfered by more uncertainties. The problem of slipping of the wheels of the mobile robot is a problem to be considered, especially when the convoy system is applied to icy or slippery road surfaces.
The sliding mode control has strong robustness on the uncertainty of system matching. However, the conventional sliding mode control has a serious buffeting problem, which not only is easy to accelerate the abrasion of the system actuator, but also has a serious negative effect on the control precision of the system.
Disclosure of Invention
The invention aims to solve the problem that the wheel slip of a mobile robot affects the formation control effect of multiple mobile robots, provides a slip-mode control method for the cooperative formation without buffeting of the multiple mobile robots considering the slip, and realizes the rapid and stable control of the formation of the multiple mobile robots even under the condition that the wheel slip exists.
The invention adopts a continuous approach law to design a multi-mobile robot formation controller, and solves the buffeting problem of sliding mode control signals. In addition, the interference observer can be designed to effectively estimate the uncertainty of the system, so that the control effect of the formation controller can be further improved.
The invention is realized by the following technical scheme:
and (1) establishing a mathematical model of the multi-mobile-robot formation control system considering the slipping condition.
(1.1) for under-actuated wheeled mobile robots, the under-actuated wheeled mobile robots are generally composed of two driving wheels and one following wheel, and the kinematic model of the under-actuated wheeled mobile robots can be expressed as follows without considering wheel slip
Figure BDA0003166184290000021
In the formula, q i =[x i ,y ii ] T ∈R 3 A pose vector representing the i (i ═ 1,2,3, …, N) th mobile robot; [ v ] of ii ] T ∈R 2 The control input velocity vector indicating the i-th mobile robot is composed of a control input linear velocity and a control input angular velocity of the mobile robot.
(1.2) when the mobile robot works on wet and slippery ground and other complex application occasions, the wheels of the mobile robot skid with the ground, and the actual control input speed of the mobile robot is not equal to the control input speed of the mobile robot, so that a certain deviation exists between the actual control input speed and the control input speed. When considering the case where the wheels of the mobile robot slip, the kinematic model of the ith mobile robot can be described as the following equation:
Figure BDA0003166184290000022
in the formula, ρ 1i ,ρ 2i The wheel slip parameter directly acting on the control input linear speed and the control input angular speed can be understood as the actual control input speed and the control input speed of the mobile robotThe ratio of degrees.
(1.3) further, the control input linear velocity and the control input angular velocity of the mobile robot and the angular velocity omega output by the two driving wheels li 、ω ri The relationship between can be expressed as:
Figure BDA0003166184290000023
Figure BDA0003166184290000024
in the formula, r i Radius of a driving wheel of the mobile robot; b is a mixture of i The distance between two driving wheels of the mobile robot; omega li The angular velocity of the left driving wheel of the mobile robot; omega ri Is the angular velocity of the right drive wheel of the mobile robot.
(1.4) the piloting mobile robot plays a role in piloting in the formation control, and the motion track of the piloting mobile robot is described by the following equation:
Figure BDA0003166184290000031
in the formula, q 0 =[x 0 ,y 00 ] T ∈R 3 The pose coordinates of the piloted mobile robot are obtained; v. of 0 Is the linear velocity of the piloted mobile robot; omega 0 Is the angular velocity of the piloted mobile robot.
(1.5) defining the expected relative distance and direction angle of formation as l id ,θ id And further, an expected track which is formed into a certain formation and is required to be tracked by following a reference point of the mobile robot i can be obtained through coordinate transformation:
Figure BDA0003166184290000032
wherein [ x ] id ,y id ] T ∈R 2 Representing the coordinate representation in the global coordinate system XOY following the desired trajectory of the mobile robot i.
(1.6) selecting a local coordinate system x of the following mobile robot bi A point p in the forward direction of the axis i As a reference point to follow the mobile robot, its coordinates may be expressed in the global coordinate system XOY as
Figure BDA0003166184290000033
L i Representing the distance between the reference point and the origin of the local coordinate system of the following mobile robot.
(1.7) to p xi ,p yi Derivation can be obtained:
Figure BDA0003166184290000034
in the formula (I), the compound is shown in the specification,
Figure BDA0003166184290000035
for systematic uncertainty, assume f i And its first derivative satisfy the relationship: | f i ||≤d 1i
Figure BDA0003166184290000036
Wherein d is 1i And d 2i As unknown normal numbers.
(1.8) the error equation for the formation control can be defined as
Figure BDA0003166184290000041
In the formula, E i =[e xi ,e yi ] T
(1.9) the error dynamic equation for the formation control of the multiple mobile robots can be obtained from the equations (5), (6) and (7):
Figure BDA0003166184290000042
in the formula (I), the compound is shown in the specification,
Figure BDA0003166184290000043
(1.10) designing a buffeting-free sliding-mode formation controller under the condition that the wheels of the formation mobile robot slip with the ground, so that the formation error of the multi-mobile robot system can overcome the influence of wheel slip, a more expected formation control effect is obtained, and meanwhile, the control signal of the formation mobile robot is ensured to be smooth.
And (2) designing a formation buffeting-free sliding mode controller in order to realize a formation control target.
(2.1) defining an integral slip-form surface as follows:
Figure BDA0003166184290000044
in the formula, S i =[s 1i ,s 2i ] T ;k 1i And more than 0 is an integral sliding mode surface design parameter.
(2.2) the derivation of equation (10) can be:
Figure BDA0003166184290000045
(2.3) the approximation rule is defined as
Figure BDA0003166184290000046
In the formula, the design parameter k of the approach law 2i ,k 3i ,k 4i ,α 1i And alpha 2i The relationship needs to be satisfied: k is a radical of 2i >0,k 3i >0,k 4i >0,0<α 1i <1,α 2i >1;
Figure BDA0003166184290000051
Figure BDA0003166184290000052
sign (·) is a sign function.
(2.4) is obtained from the formulae (11) and (12):
Figure BDA0003166184290000053
(2.5) the first part of the controller may be designed to
Figure BDA0003166184290000054
In the formula (I), the compound is shown in the specification,
Figure BDA0003166184290000055
is f i An estimate of (d). Uncertainty term f i The estimated value of (c) can be obtained by designing the disturbance observers equation (15) and equation (16):
Figure BDA0003166184290000056
Figure BDA0003166184290000057
in the formula, K oi =diag{k 01i ,k 02i The value of 0 is a gain matrix of the interference observer to be designed; beta is a i Is the intermediate vector of the disturbance observer.
Defining a system uncertainty estimation error of
Figure BDA0003166184290000058
From equations (7), (15) and (16), the dynamic equation of the estimation error can be obtained:
Figure BDA0003166184290000059
thus, the uncertainty estimation error can be determined by simply selecting the appropriate observer gain matrixConverge to a very small area near the origin, i.e. satisfy the relation
Figure BDA00031661842900000510
ε i An upper bound for the error is estimated for the disturbance observer.
(2.6) the second part of the controller may be designed to
Figure BDA0003166184290000061
(2.7) the complete formation controller expression is obtained from equations (14) and (18):
Figure BDA0003166184290000062
and (3) analyzing the stability of the formation controller.
(3.1) designing a Lyapunov function:
Figure BDA0003166184290000063
(3.2) is obtained from the formula (11), the formula (20):
Figure BDA0003166184290000064
(3.3) substitution of formula (19) for formula (21) gives:
Figure BDA0003166184290000065
in the formula, beta i =min{k 2i ,k 3i ,k 4i }。
(3.4) rewriting the formula (22) as
Figure BDA0003166184290000071
(3.5) As is clear from formula (23), provided that s ji One of the following conditions is satisfied:
Figure BDA0003166184290000072
the following equation can be made:
Figure BDA0003166184290000073
(3.6) therefore, the system is stable as seen by Lyapunov's stability theory and can be controlled by selecting the appropriate controller parameter k 2i ,k 3i ,k 4i Can make the sliding mode variable s ji And the method converges to a tiny area near the zero point, so that a good formation error control effect is obtained.
The invention has the advantages that: the method can eliminate the adverse effect of wheel slip on the formation control effect, can ensure the rapid and stable control of formation even under the condition of wheel slip, and simultaneously, the control signal output by the formation controller is smooth.
Drawings
Fig. 1 is a diagram illustrating a formation structure of a multi-mobile robot considering a slip situation according to the present invention.
Fig. 2 shows the formation error of the following mobile Robot Follower Robot1 according to the present invention.
Fig. 3 shows the formation error of the following mobile Robot Follower Robot2 according to the present invention.
FIG. 4a shows the uncertainty f of the following mobile Robot Follower Robot1 disturbance observer to the system 11 Is estimated.
FIG. 4b shows uncertainty f of following mobile Robot Follower Robot1 disturbance observer to system 21 Is estimated.
FIG. 5a shows the uncertainty f of the following mobile Robot Follower Robot2 disturbance observer to the system 12 Is estimated.
FIG. 5b shows uncertainty f of following mobile Robot Follower Robot2 disturbance observer to system 22 Is estimated.
Fig. 6 shows control signals of a following mobile Robot, namely a Follower Robot 1.
Fig. 7 shows control signals of a following mobile Robot, namely a Follower Robot 2.
Fig. 8 is a schematic structural diagram of the controller according to the present invention.
Detailed Description
In order to more intuitively illustrate the technical solution and technical advantages of the present invention, the following description of the technical solution of the present invention with reference to specific embodiments is made clearly and completely with reference to fig. 2 to 7.
A slip-considered multi-mobile-robot collaborative formation buffeting-free sliding mode control method comprises the following specific technical steps:
and (1) establishing a mathematical model of the multi-mobile-robot formation control system considering the slipping condition.
(1.1) for the wheel-type mobile robot under-actuated, the mobile robot is generally composed of two driving wheels and one following wheel, and the kinematic model of the mobile robot can be expressed as the following wheel without considering the wheel slip
Figure BDA0003166184290000081
In the formula, q i =[x i ,y ii ] T ∈R 3 A pose vector representing the i (i ═ 1,2,3, …, N) th mobile robot; [ v ] of ii ] T ∈R 2 The control input velocity vector representing the ith mobile robot is composed of a control input linear velocity and a control input angular velocity of the mobile robot.
(1.2) when the mobile robot works on wet and slippery ground and other complex application occasions, the wheels of the mobile robot skid with the ground, and the actual control input speed of the mobile robot is not equal to the control input speed of the mobile robot, so that a certain deviation exists between the actual control input speed and the control input speed. When considering the case where the wheels of the mobile robot slip, the kinematic model of the ith mobile robot can be described as the following equation:
Figure BDA0003166184290000082
in the formula, ρ 1i ,ρ 2i The wheel slip parameter, which is indicative of the direct action on the control input linear and control input angular velocities, is understood to be the ratio of the actual control input velocity of the mobile robot to its control input velocity.
(1.3) further, the control input linear velocity and the control input angular velocity of the mobile robot and the angular velocity omega output by the two driving wheels li 、ω ri The relationship between can be expressed as:
Figure BDA0003166184290000091
Figure BDA0003166184290000092
in the formula, r i Radius of a driving wheel of the mobile robot; b i The distance between two driving wheels of the mobile robot; omega li The angular velocity of the left driving wheel of the mobile robot; omega ri Is the angular velocity of the right drive wheel of the mobile robot.
(1.4) the piloting mobile robot plays a role in piloting in the formation control, and the motion track of the piloting mobile robot is described by the following equation:
Figure BDA0003166184290000093
in the formula, q 0 =[x 0 ,y 00 ] T ∈R 3 The pose coordinates of the piloted mobile robot are obtained; v. of 0 Is the linear velocity of the piloted mobile robot; omega 0 To pilot the angular velocity of a mobile robot。
(1.5) defining the expected relative distance and direction angle of formation as l id ,θ id And further, expected tracks which form a certain formation and need to be tracked along a reference point of the mobile robot i can be obtained through coordinate transformation:
Figure BDA0003166184290000094
wherein [ x ] id ,y id ] T ∈R 2 Representing the coordinate representation in the global coordinate system XOY following the desired trajectory of the mobile robot i.
(1.6) selecting a local coordinate system x of the following mobile robot bi A point p in the forward direction of the axis i As a reference point to follow the mobile robot, its coordinates may be expressed in the global coordinate system XOY as
Figure BDA0003166184290000101
L i Representing the distance between the reference point and the origin of the local coordinate system of the following mobile robot.
(1.7) to p xi ,p yi The derivation can be:
Figure BDA0003166184290000102
in the formula (I), the compound is shown in the specification,
Figure BDA0003166184290000103
for systematic uncertainty, assume f i And its first derivative satisfy the relationship: i f i ||≤d 1i
Figure BDA0003166184290000104
Wherein d is 1i And d 2i As unknown normal numbers.
(1.8) the error equation for the formation control can be defined as
Figure BDA0003166184290000105
In the formula, E i =[e xi ,e yi ] T
(1.9) the error dynamic equation for the formation control of the multiple mobile robots can be obtained from the equations (5), (6) and (7):
Figure BDA0003166184290000106
in the formula (I), the compound is shown in the specification,
Figure BDA0003166184290000107
(1.10) designing a buffeting-free sliding-mode formation controller under the condition that the wheels of the formation mobile robot slip with the ground, so that the formation error of the multi-mobile robot system can overcome the influence of wheel slip, a more expected formation control effect is obtained, and meanwhile, the control signal of the formation mobile robot is ensured to be smooth.
And (2) designing a formation buffeting-free sliding mode controller for realizing a formation control target, as shown in fig. 8.
(2.1) defining an integral slip-form surface as follows:
Figure BDA0003166184290000111
in the formula, S i =[s 1i ,s 2i ] T ;k 1i And more than 0 is an integral sliding mode surface design parameter.
(2.2) derivation of equation (10) can be found:
Figure BDA0003166184290000112
(2.3) the approximation rule is defined as
Figure BDA0003166184290000113
In the formula, the design parameter k of the approach law 2i ,k 3i ,k 4i ,α 1i And alpha 2i The relationship needs to be satisfied: k is a radical of 2i >0,k 3i >0, k 4i >0,0<α 1i <1,α 2i >1;
Figure BDA0003166184290000114
Figure BDA0003166184290000115
sign (·) is a sign function.
(2.4) is obtained from the formulae (11) and (12):
Figure BDA0003166184290000116
(2.5) the first part of the controller may be designed to be
Figure BDA0003166184290000117
In the formula (I), the compound is shown in the specification,
Figure BDA0003166184290000118
is f i An estimate of (d). Uncertainty term f i The estimated value of (c) can be obtained by designing the disturbance observers equation (15) and equation (16):
Figure BDA0003166184290000119
Figure BDA0003166184290000121
in the formula, K oi =diag{k 01i ,k 02i The value of 0 is a gain matrix of the interference observer to be designed; beta is a beta i Is the intermediate vector of the disturbance observer;
defining a system uncertainty estimation error as
Figure BDA0003166184290000122
From equations (7), (15) and (16), the dynamic equation of the estimation error can be obtained:
Figure BDA0003166184290000123
therefore, the uncertainty estimation error can be converged to a very small area near the origin, i.e. the relation is satisfied, by selecting a suitable observer gain matrix
Figure BDA0003166184290000124
ε i An upper bound for the error is estimated for the disturbance observer.
(2.6) the second part of the controller may be designed to
Figure BDA0003166184290000125
(2.7) the complete formation controller expression is obtained from equations (14) and (18):
Figure BDA0003166184290000126
and (3) analyzing the stability of the formation controller.
(3.1) designing a Lyapunov function:
Figure BDA0003166184290000127
(3.2) is obtained from the formula (11), the formula (20):
Figure BDA0003166184290000131
(3.3) the formula (19) may be substituted for the formula (21):
Figure BDA0003166184290000132
in the formula, beta i =min{k 2i ,k 3i ,k 4i }。
(3.4) rewriting formula (22) as
Figure BDA0003166184290000133
(3.5) As is clear from formula (23), provided that s ji One of the following conditions is satisfied:
Figure BDA0003166184290000134
the following equation can be made:
Figure BDA0003166184290000135
(3.6) therefore, the system is stable as known from Lyapunov's stability theory and can be controlled by selecting the appropriate controller parameter k 2i ,k 3i ,k 4i Namely, the sliding mode variable s ji And the method converges to a tiny area near the zero point, so that a good formation error control effect is obtained.
In the embodiment, three mobile robots are provided, namely a Leader Robot, a Follower Robot1 and a Follower Robot 2.
Linear velocity v of piloted mobile robot 0 Is set as v 0 0.3m/s, angular velocity ω 0 Set to ω 0 0.1rad/s, with initial pose coordinate set to q 0 =[0,0,0] T . The physical parameter of the following mobile Robot Robot1 is set as r 1 =0.08m,b 1 =0.3m,L 1 1m, the reference point initial coordinate is set to
Figure BDA0003166184290000141
The physical parameter of the following mobile Robot2 is set as r 2 =0.08m, b 2 =0.3m,L 2 1m, the reference point initial coordinate is set to
Figure BDA0003166184290000142
The controller parameters for the Follower Robot1 are set to: k is a radical of 11 =2,k 21 =0.4,k 31 =0.1, k 41 =0.1,α 11 =0.8,α 21 1.8; observer gain matrix parameter set to K o1 Diag { -6, -6 }. The controller parameters for the Follower Robot1 are set to: k is a radical of 12 =2,k 22 =0.4,k 32 =0.1,k 42 =0.1,α 12 =0.8,α 22 1.8; observer gain matrix parameter set to K o2 =diag{-6,-6}。
Simulation experiment results of the control method are shown in fig. 2-7, and fig. 2 and 3 are schematic diagrams of formation errors of following mobile robots, namely, a Follower Robot1 and a Follower Robot2 respectively; FIGS. 4 and 5 are schematic diagrams of the following mobile robots Follower Robot1 and Follower Robot2 for estimating the system uncertainty by the disturbance observer, respectively; fig. 6 and 7 are control signals of following mobile robots, Follower Robot1 and Follower Robot2, respectively. As can be seen from fig. 2 and 3, the controller can well realize the control of the formation system of the multiple mobile robots with the slipping condition between the wheels and the ground, when the slipping condition occurs between the wheels and the ground within 15 seconds, the formation error fluctuates to a certain extent, but the formation error can still obtain a good control effect due to the timely intervention of the designed interference observer. As can be seen from fig. 6 and 7, the control signal of the controller is smooth, which is advantageous for practical application of the controller.

Claims (1)

1. A slip-considered multi-mobile-robot collaborative formation buffeting-free sliding mode control method is characterized by comprising the following steps: the method specifically comprises the following steps:
(1) establishing a mathematical model of a multi-mobile-robot formation control system considering a slipping condition;
(2) designing a formation buffeting-free sliding mode controller based on the established mathematical model, and realizing that each mobile robot in the formation can form a formation with a pilot robot according to a preset formation;
(3) analyzing the stability of the designed formation buffeting-free sliding mode controller;
establishing a mathematical model of the multi-mobile-robot formation control system considering the slipping condition in the step (1), which comprises the following specific steps:
(1.1) the under-actuated wheeled mobile robot is composed of two driving wheels and one following wheel, and its kinematic model is expressed as the following wheel without considering the wheel slip
Figure FDA0003679326750000011
In the formula, q i =[x i ,y ii ] T ∈R 3 A pose vector representing the ith mobile robot, i ═ 1,2,3, …, N; [ v ] of ii ] T ∈R 2 The control input velocity vector of the ith mobile robot is represented and is composed of a control input linear velocity and a control input angular velocity of the mobile robot;
(1.2) when the mobile robot works on a slippery ground, slipping occurs between wheels of the mobile robot and the ground, the actual control input speed of the mobile robot is not equal to the control input speed of the mobile robot, and a certain deviation exists between the actual control input speed and the control input speed, and when the situation that the wheels of the mobile robot slip is considered, the kinematic model of the ith mobile robot is described as the following equation:
Figure FDA0003679326750000012
in the formula, ρ 1i ,ρ 2i Indicating wheel slip parameters acting directly on control input linear and control input angular velocities, i.e. mobile machinesThe ratio of the actual human control input speed to its control input speed;
(1.3) control input linear velocity and control input angular velocity of mobile robot and angular velocity omega output by two driving wheels li 、ω ri The relationship between them is expressed as:
Figure FDA0003679326750000021
Figure FDA0003679326750000022
in the formula, r i Radius of a driving wheel of the mobile robot; b i The distance between two driving wheels of the mobile robot; omega li Is the angular velocity of the left drive wheel of the mobile robot; omega ri The angular velocity of the right driving wheel of the mobile robot;
(1.4) the motion trajectory of the piloted mobile robot is described by the following equation:
Figure FDA0003679326750000023
in the formula, q 0 =[x 0 ,y 00 ] T ∈R 3 The pose coordinates of the piloted mobile robot are obtained; v. of 0 The linear velocity of the piloted mobile robot; omega 0 Is the angular velocity of the piloted mobile robot;
(1.5) defining the expected relative distance and direction angle of formation as l id ,θ id And obtaining an expected track which forms a certain formation and needs to be tracked along with a reference point of the mobile robot by coordinate transformation:
Figure FDA0003679326750000024
wherein [ x ] id ,y id ] T ∈R 2 A coordinate representation in the global coordinate system XOY representing a desired trajectory for following the mobile robot i;
(1.6) selecting a local coordinate system x of the following mobile robot bi A point p in the forward direction of the axis i As a reference point to follow the mobile robot, its coordinates are represented in the global coordinate system XOY as
Figure FDA0003679326750000025
L i Representing the distance between the reference point and the origin of the local coordinate system of the following mobile robot;
(1.7) to p xi ,p yi And (5) obtaining a derivative:
Figure FDA0003679326750000031
in the formula (I), the compound is shown in the specification,
Figure FDA0003679326750000032
for systematic uncertainty, assume f i And its first derivative satisfy the relationship: | f i ||≤d 1i
Figure FDA0003679326750000033
Wherein d is 1i And d 2i Is an unknown normal number;
(1.8) formation control error equation is defined as
Figure FDA0003679326750000034
In the formula, E i =[e xi ,e yi ] T
(1.9) error dynamic equation for formation control of multiple mobile robots by equation (5), equation (6) and equation (7):
Figure FDA0003679326750000035
in the formula (I), the compound is shown in the specification,
Figure FDA0003679326750000036
designing a formation buffeting-free sliding mode controller based on the established mathematical model in the step (2), and realizing that each mobile robot in the formation can form a formation with a pilot robot according to a preset formation, wherein the formation comprises the following specific steps:
(2.1) defining an integral slip-form surface as follows:
Figure FDA0003679326750000037
in the formula, S i =[s 1i ,s 2i ] T ;k 1i More than 0 is an integral sliding mode surface design parameter;
(2.2) deriving equation (10):
Figure FDA0003679326750000038
(2.3) the approximation rule is defined as
Figure FDA0003679326750000041
In the formula, the design parameter k of the approach law 2i ,k 3i ,k 4i ,α 1i And alpha 2i The relationship needs to be satisfied: k is a radical of formula 2i >0,k 3i >0,k 4i >0,0<α 1i <1,α 2i >1;
Figure FDA0003679326750000042
Figure FDA0003679326750000043
sign (·) is a sign function;
(2.4) is obtained from the following formulae (11) and (12):
Figure FDA0003679326750000044
(2.5) the first part of the formation buffeting-free sliding mode controller is designed to be
Figure FDA0003679326750000045
In the formula (I), the compound is shown in the specification,
Figure FDA0003679326750000046
is f i Estimated value of (a), uncertainty term f i The estimated value of (c) is obtained by the designed disturbance observers equation (15) and equation (16):
Figure FDA0003679326750000047
Figure FDA0003679326750000048
in the formula, K oi =diag{k 01i ,k 02i The value of 0 is a gain matrix of the interference observer to be designed; beta is a beta i Is the intermediate vector of the disturbance observer;
defining a system uncertainty estimation error as
Figure FDA0003679326750000049
The dynamic equation of the estimation error is obtained from equations (7), (15) and (16):
Figure FDA00036793267500000410
satisfy the relationship
Figure FDA00036793267500000411
ε i Estimating an upper bound of error for the disturbance observer;
(2.6) second part of the formation buffeting-free sliding mode controller is designed to be
Figure FDA0003679326750000051
(2.7) the complete formation controller expression is obtained from equations (14) and (18):
Figure FDA0003679326750000052
analyzing the stability of the designed formation buffeting-free sliding mode controller in the step (3), and concretely comprising the following steps:
(3.1) designing a Lyapunov function:
Figure FDA0003679326750000053
(3.2) is obtained from the following formulae (11) and (20):
Figure FDA0003679326750000054
(3.3) substitution of formula (19) for formula (21):
Figure FDA0003679326750000055
in the formula, beta i =min{k 2i ,k 3i ,k 4i };
(3.4) rewriting formula (22) as
Figure FDA0003679326750000061
(3.5) As long as | s is given by the formula (23) ji L satisfies one of the following conditions:
Figure FDA0003679326750000062
that is, the following equation holds:
Figure FDA0003679326750000063
(3.6) the system is stable from Lyapunov's stability theory and by choosing the appropriate controller parameter k 2i ,k 3i ,k 4i So that the sliding mode variable s ji And the method converges to a tiny area near the zero point, so that a good formation error control effect is obtained.
CN202110805055.XA 2021-07-16 2021-07-16 Slip-considered multi-mobile-robot cooperative formation buffeting-free sliding mode control method Active CN113377016B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110805055.XA CN113377016B (en) 2021-07-16 2021-07-16 Slip-considered multi-mobile-robot cooperative formation buffeting-free sliding mode control method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110805055.XA CN113377016B (en) 2021-07-16 2021-07-16 Slip-considered multi-mobile-robot cooperative formation buffeting-free sliding mode control method

Publications (2)

Publication Number Publication Date
CN113377016A CN113377016A (en) 2021-09-10
CN113377016B true CN113377016B (en) 2022-08-30

Family

ID=77582371

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110805055.XA Active CN113377016B (en) 2021-07-16 2021-07-16 Slip-considered multi-mobile-robot cooperative formation buffeting-free sliding mode control method

Country Status (1)

Country Link
CN (1) CN113377016B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115202349B (en) * 2022-07-14 2024-07-05 广州大学 Multi-mobile robot collaborative formation control method, device, equipment and storage medium based on communication interference

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111665865A (en) * 2020-06-16 2020-09-15 天津大学 Multi-unmanned aerial vehicle formation robust control method based on virtual structure method

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2005177977A (en) * 2003-11-25 2005-07-07 Matsushita Electric Works Ltd Robot hand control device
CN105607636B (en) * 2016-01-21 2018-03-02 浙江工业大学 A kind of wheeled mobile robot Leader-Follower Formation control method based on Integral Sliding Mode algorithm
CN107577145B (en) * 2017-08-25 2020-06-09 湘潭大学 Backstepping sliding mode control method for formation flying spacecraft
CN108614425B (en) * 2018-06-06 2021-03-16 黄山学院 Self-adaptive integral sliding mode control method for mobile robot
CN109032137B (en) * 2018-07-24 2020-11-06 西北工业大学 Distributed tracking control method for multi-Euler-Lagrange system
CN109739094A (en) * 2019-02-28 2019-05-10 哈尔滨理工大学 A kind of mobile robot trace tracking new method based on adaptive sliding-mode observer
CN110083061B (en) * 2019-04-12 2022-04-08 国网宁夏电力有限公司电力科学研究院 Wheel type mobile robot control method
CN110109363B (en) * 2019-05-28 2022-10-21 重庆理工大学 Neural network self-adaptive control method for wheeled mobile robot formation
CN111897358B (en) * 2020-07-30 2022-04-15 南京航空航天大学 Unmanned aerial vehicle formation fault-tolerant control method based on self-adaptive sliding mode

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111665865A (en) * 2020-06-16 2020-09-15 天津大学 Multi-unmanned aerial vehicle formation robust control method based on virtual structure method

Also Published As

Publication number Publication date
CN113377016A (en) 2021-09-10

Similar Documents

Publication Publication Date Title
Li et al. Development of a new integrated local trajectory planning and tracking control framework for autonomous ground vehicles
Rijalusalam et al. Implementation kinematics modeling and odometry of four omni wheel mobile robot on the trajectory planning and motion control based microcontroller
CN108614426B (en) Multi-mobile-robot formation robust control method based on disturbance observer
CN108008720B (en) Fuzzy sliding mode trajectory tracking control and method for wheeled mobile robot
Ji et al. Path planning and tracking for vehicle collision avoidance based on model predictive control with multiconstraints
Fang et al. Trajectory tracking control of farm vehicles in presence of sliding
CN108227491B (en) Intelligent vehicle track tracking control method based on sliding mode neural network
CN109799814B (en) Design method of track unmanned vehicle track tracking controller based on kinematics model
CN108381555B (en) Design method of lower controller of redundancy mechanical arm of flying operation robot
CN113268064B (en) Multi-mobile-robot cooperative formation control method considering communication time delay
CN110045598B (en) Method for tracking and controlling underground advancing path of cantilever type heading machine
CN109857100B (en) Composite track tracking control algorithm based on inversion method and fast terminal sliding mode
Peicheng et al. Intelligent vehicle path tracking control based on improved MPC and hybrid PID
CN109739094A (en) A kind of mobile robot trace tracking new method based on adaptive sliding-mode observer
CN113377016B (en) Slip-considered multi-mobile-robot cooperative formation buffeting-free sliding mode control method
Low A flexible leader-follower formation tracking control design for nonholonomic tracked mobile robots with low-level velocities control systems
CN110109363B (en) Neural network self-adaptive control method for wheeled mobile robot formation
CN107085432B (en) Target track tracking method of mobile robot
Yue et al. Path tracking control of skid-steered mobile robot on the slope based on fuzzy system and model predictive control
CN115202349A (en) Multi-mobile-robot cooperative formation control method, device, equipment and storage medium based on communication interference
CN113848905B (en) Mobile robot track tracking method based on neural network and self-adaptive control
Li et al. Fixed-time sliding mode-based adaptive path tracking control of maize plant protection robot via extreme learning machine
CN113467475A (en) Track tracking sliding mode control method for Mecanum wheel omnidirectional mobile robot
Bacha et al. Autonomous vehicle path tracking using nonlinear steering control and input-output state feedback linearization
Xie et al. Parameter self‐learning feedforward compensation‐based active disturbance rejection for path‐following control of self‐driving forklift trucks

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant