CN114488790B - Omnidirectional mobile robot self-adaptive sliding film control method based on nominal model - Google Patents

Omnidirectional mobile robot self-adaptive sliding film control method based on nominal model Download PDF

Info

Publication number
CN114488790B
CN114488790B CN202111551149.5A CN202111551149A CN114488790B CN 114488790 B CN114488790 B CN 114488790B CN 202111551149 A CN202111551149 A CN 202111551149A CN 114488790 B CN114488790 B CN 114488790B
Authority
CN
China
Prior art keywords
mobile robot
omnidirectional mobile
representing
equation
omnidirectional
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
CN202111551149.5A
Other languages
Chinese (zh)
Other versions
CN114488790A (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.)
Southwest University of Science and Technology
Original Assignee
Southwest University of Science and Technology
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 Southwest University of Science and Technology filed Critical Southwest University of Science and Technology
Priority to CN202111551149.5A priority Critical patent/CN114488790B/en
Publication of CN114488790A publication Critical patent/CN114488790A/en
Application granted granted Critical
Publication of CN114488790B publication Critical patent/CN114488790B/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 an omnidirectional mobile robot self-adaptive synovial membrane control method based on a nominal model, which belongs to the technical field of tracking system control, and comprises the following steps: establishing an omnidirectional mobile robot kinematics equation; establishing an omnidirectional mobile robot dynamics equation according to the omnidirectional mobile robot dynamics equation; calculating to obtain a slip form surface s (t) and a slip film control rate tau based on an omnidirectional mobile robot kinematic equation and a dynamic equation; setting a controller for adjusting the self-adaptive gain of the nominal model according to the sliding mode surface and the sliding film control rate, obtaining a dynamic equation of the controller omnidirectional mobile robot based on the self-adaptive gain adjustment of the nominal model, and completing the self-adaptive sliding film control of the omnidirectional mobile robot based on the nominal model; the problem that the track tracking error is large because the omnidirectional mobile robot can not completely eliminate interference under the condition of interference is solved.

Description

Omnidirectional mobile robot self-adaptive sliding film control method based on nominal model
Technical Field
The invention belongs to the technical field of tracking system control, and particularly relates to an omnidirectional mobile robot self-adaptive synovial membrane control method based on a nominal model.
Background
With the rapid development of industry, research on omni-directional mobile robots is also receiving more and more attention, and omni-directional wheel mobile robots are used in various scenes due to advantages of better operability, zero turning radius in any direction and the like; compared with the traditional omni-directional wheel type mobile robot, the Mecanum wheel type mobile robot can move transversely and even can move along a curved path and rotate in situ; however, with the use of independent motor control, its motion control becomes with uncertainty when there is a disturbance; many scholars have derived kinematic and kinetic modeling about omni-wheel mobile robots; mechanical characteristics of a single wheel are studied in detail from the incision of the single omnidirectional wheel, the pose error of the robot is obtained by establishing an omnidirectional mobile robot kinematics equation, and a dynamics equation is established by analyzing influence factors of friction force on the four wheels; however, the simplest PD control is used in the control process, so that the shake of the robot cannot be eliminated; many scholars deduce a dynamic model of the omnidirectional wheel mobile robot under the condition of friction force, an active disturbance rejection controller with friction compensation is designed by establishing the dynamic model of the friction force, a friction compensation item is firstly carried out, then a nonlinear controller is adopted for track tracking, and the influence of dead zone nonlinearity of the friction force on the control precision of a system is greatly reduced; however, only the limited input condition is considered, and the conditions such as gradual stability, exponential stability and the like are not considered, so that the stability is poor, and meanwhile, the influence of dead zone nonlinearity of friction force on a system cannot be completely eliminated.
The motion of the omni-wheel mobile robot under the disturbance condition has better tracking effect compared with a linear controller (PID) by using a nonlinear controller; the existing slide film control research on track tracking of the wheeled mobile robot aims at the problem of external interference track tracking control of the wheeled mobile robot, a fuzzy rapid double-power approach law slide film control strategy is designed, the robustness of a track tracking control system is improved, the buffeting phenomenon of the robot is eliminated, and the self-adaptive slide film control is realized on the wheeled mobile robot; in order to improve the robustness of path tracking control of the wheeled mobile robot, an inversion self-adaptive sliding mode control strategy is provided, a method of inversion control is adopted by a kinematic control law, linear speed and angular speed are used as auxiliary quantity control input, so that the track approaches to an ideal track as much as possible, a dynamic model adopts a self-adaptive sliding film control law, an input torque controller of the mobile robot is obtained, and the speed approaches to an expected speed as much as possible, but the controller is not well applied to the Mecanum wheel mobile robot.
Disclosure of Invention
Aiming at the defects in the prior art, the self-adaptive sliding film control method of the omnidirectional mobile robot based on the nominal model combines the advantages of a sliding mode control system and a self-adaptive control system, and solves the problem that the omnidirectional mobile robot cannot completely eliminate interference under the condition of interference, so that the track tracking error is larger.
In order to achieve the aim of the invention, the invention adopts the following technical scheme:
the invention provides an omnidirectional mobile robot self-adaptive synovial membrane control method based on a nominal model, which comprises the following steps:
s1, establishing an omnidirectional mobile robot kinematics equation;
s2, establishing an omnidirectional mobile robot dynamics equation according to the omnidirectional mobile robot dynamics equation;
s3, calculating to obtain a sliding mode surface S (t) and a sliding film control rate tau based on an omnidirectional mobile robot kinematics equation and a dynamics equation;
and S4, setting a controller for adjusting the self-adaptive gain of the nominal model according to the sliding mode surface and the sliding film control rate, obtaining a dynamic equation of the controller omnidirectional mobile robot based on the self-adaptive gain adjustment of the nominal model, and completing the self-adaptive sliding film control of the omnidirectional mobile robot based on the nominal model.
The beneficial effects of the invention are as follows: according to the self-adaptive synovial membrane control method of the omnidirectional mobile robot based on the nominal model, which is provided by the invention, the robustness is very strong under the interference condition according to the synovial membrane control, the self-adaptive control can enable the system to work in an optimal or near-optimal motion state, and the accurate control performance is obtained.
Further, the step S1 includes the steps of:
s11, setting the geometric center of the omnidirectional mobile robot to coincide with the mass center, and defining a global coordinate system O f X f Y f Moving coordinate system O m X m Y m And a Mecanum wheel coordinate system O wi X wi Y wi Wherein i is equal to 1,2,3,4;
s12, decomposing the axle center speed of each Mecanum wheel of the omnidirectional mobile robot to obtain the speed of each Mecanum wheel respectively;
s13, obtaining the relation between the axle center speed of each Mecanum wheel and the main body center speed of the omnidirectional mobile robot;
s14, obtaining an inverse kinematics equation of the omnidirectional mobile robot according to the relation between the axle center speed of each Mecanum wheel and the central speed of the main body of the omnidirectional mobile robot and the structural characteristics of the omnidirectional mobile robot;
and S15, obtaining an omnidirectional mobile robot kinematics equation according to the omnidirectional mobile robot inverse kinematics equation, and completing establishment of the omnidirectional mobile robot kinematics equation.
The beneficial effects of adopting the further scheme are as follows: the kinematics equation of the omnidirectional mobile robot is established by analyzing the omnidirectional mobile robot main body and the Mecanum wheel, and a foundation is provided for establishing the dynamics equation of the omnidirectional mobile robot.
Further, the expression of the kinematic equation of the omnidirectional mobile robot in the step S15 is as follows:
Figure BDA0003417234180000041
wherein,,
Figure BDA0003417234180000042
first derivative representing lateral displacement of omni mobile robot,/->
Figure BDA0003417234180000043
First derivative representing longitudinal displacement of omni mobile robot,/->
Figure BDA0003417234180000044
Representing the first derivative of the rotation angle of an omnidirectional mobile robot, R representing the radius of the Mecanum wheel, l 1 And l 2 Respectively represent a Mecanum wheel center distance movement coordinate system x m Axes and y m Distance of axis>
Figure BDA0003417234180000045
And
Figure BDA0003417234180000046
the angular velocities of the first, second, third and fourth mecanum wheels of the omnidirectional mobile robot are represented, respectively.
The beneficial effects of adopting the further scheme are as follows: the omnidirectional mobile robot kinematics equation is provided, and a foundation is provided for establishing the omnidirectional mobile robot kinematics equation.
Further, the step S2 includes the steps of:
s21, constructing an omnidirectional mobile robot kinetic energy equation according to the omnidirectional mobile robot kinetic energy equation;
s22, calculating to obtain an omnidirectional mobile robot kinetic equation according to the Lagrangian kinematic equation and the omnidirectional mobile robot kinetic equation, and completing establishment of the omnidirectional mobile robot kinetic equation.
The beneficial effects of adopting the further scheme are as follows: when the robot dynamics equation is established, the condition of external disturbance is considered, and a better accurate dynamics model can be obtained under the condition of disturbance.
Further, the expression of the kinetic equation of the omnidirectional mobile robot in the step S22 is as follows:
Figure BDA0003417234180000047
wherein M represents an inertial matrix of the omnidirectional mobile robot,
Figure BDA0003417234180000048
represents the angular acceleration of the wheels of the Mecanum robot,/->
Figure BDA0003417234180000049
The friction moment of the Mecanum wheel is represented, tau represents the total moment of the omnidirectional mobile robot, and d represents the external disturbance moment.
The beneficial effects of adopting the further scheme are as follows: and establishing an omnidirectional mobile robot dynamics equation under the condition of considering external disturbance factors, providing an omnidirectional mobile robot dynamics equation, and providing a foundation for the omnidirectional mobile robot to run on a sliding mode surface by using a controller.
Further, the step S3 includes the following steps:
s31, defining a tracking error e (t) of the omnidirectional mobile robot according to a kinematic equation and a dynamic equation of the omnidirectional mobile robot;
s32, calculating to obtain a sliding mode surface and a sliding film control rate according to the tracking error of the omnidirectional mobile robot.
The beneficial effects of adopting the further scheme are as follows: the linear sliding hyperplane is selected, so that the tracking error is gradually converged to zero after the system reaches the sliding mode, and the convergence speed can be adjusted by selecting a control rate matrix of the sliding mode surface.
Further, the expression of the tracking error e (t) of the omnidirectional mobile robot in the step S31 is as follows:
Figure BDA0003417234180000051
wherein,,
Figure BDA0003417234180000052
indicating the desired angular velocity of the Mecanum wheel, < >>
Figure BDA0003417234180000053
The actual angular velocities of the respective mecanum wheels are represented, where i=1, 2,3,4, respectively represent the angular velocities of the first, second, third and fourth mecanum wheels.
The beneficial effects of adopting the further scheme are as follows: the tracking error equation of the omnidirectional mobile robot is provided, and error reference is provided for controlling and adjusting the omnidirectional mobile robot to run on the sliding mode surface.
Further, the expressions of the slip plane S (t) and the slip film control rate τ' in the step S32 are as follows:
Figure BDA0003417234180000054
Figure BDA0003417234180000055
wherein,,
Figure BDA0003417234180000056
derivative, κ, representing tracking error 4×4 Representing a constant positive definite matrix, e (t) representing tracking error,/->
Figure BDA0003417234180000057
Nominal parameters representing the inertial matrix of an omnidirectional mobile robot,/->
Figure BDA0003417234180000058
Represents the angular acceleration of the wheels of the Mecanum robot,/->
Figure BDA0003417234180000059
Representing the derivative of the slide face,
Figure BDA00034172341800000510
Nominal parameter representing moment of inertia,/->
Figure BDA00034172341800000511
Representing the angular velocity of the Mecanum wheel, ρ representing the slip-mode control rate gain and sgn representing the step function.
The beneficial effects of adopting the further scheme are as follows: the sliding film control convergence speed is high, the sliding film control convergence speed can be adjusted by selecting a sliding mode surface parameter matrix, but the stability is poor, so that the controller for nominal model self-adaptive gain adjustment is adopted to control the omnidirectional mobile robot.
Further, the step S4 includes the steps of:
s41, setting a nominal model self-adaptive gain adjustment controller according to a sliding mode surface and a sliding film control rate to obtain a control rate tau' of the nominal model self-adaptive gain adjustment controller;
s42, obtaining a controller omnidirectional mobile robot dynamic model based on the nominal model adaptive gain adjustment according to the nominal parameters of the inertia matrix of the omnidirectional mobile robot, the nominal parameters of the friction coefficient of the omnidirectional mobile robot, the control rate tau' of the controller of the nominal model adaptive gain adjustment and the omnidirectional mobile robot dynamic model, and completing the omnidirectional mobile robot adaptive synovial membrane control based on the nominal model.
The beneficial effects of adopting the further scheme are as follows: because the controlled object is greatly influenced by external disturbance, the system can work in an optimal or near-optimal motion state by adopting self-adaptive control, and the accurate control performance is obtained.
Further, the controller omnidirectional mobile robot dynamics model expression based on the nominal model adaptive gain adjustment in the step S42 is as follows:
Figure BDA0003417234180000061
Figure BDA0003417234180000062
Figure BDA0003417234180000063
u r =-ρsgn(s(t))
wherein M represents an inertial matrix of the omnidirectional mobile robot,
Figure BDA0003417234180000064
represents the angular acceleration of the wheels of the Mecanum robot,/->
Figure BDA0003417234180000065
Representing the friction moment of the Mecanum wheel, +.>
Figure BDA0003417234180000066
Nominal parameters representing the inertial matrix of an omnidirectional mobile robot,/->
Figure BDA0003417234180000067
Nominal parameters representing moment of inertia, d representing external disturbance torque, τ' representing slip film control rate, +.>
Figure BDA0003417234180000068
Indicating jerk, κ of the wheels of a Mecanum robot 4×4 Representing a constant positive definite matrix, ">
Figure BDA0003417234180000069
Representing the derivative of the tracking error, e (t) representing the tracking error,/and (c)>
Figure BDA00034172341800000610
Represents the angular velocity of the Mecanum wheel, +.>
Figure BDA00034172341800000611
Representing the nominal value of the external disturbance moment, u θ Representing the friction compensation term of the omnidirectional mobile robot, u r The slip film approach law is represented, T represents time, T represents period, ρ represents slip mode control rate gain, sgn represents step function.
The beneficial effects of adopting the further scheme are as follows: the controller omnidirectional mobile robot dynamics model based on the nominal model self-adaptive gain adjustment is provided, and is verified according to the Leiprov function derivation.
Drawings
Fig. 1 is a flowchart of steps of a method for controlling an adaptive synovial membrane of an omni-directional mobile robot based on a nominal model according to an embodiment of the present invention.
Fig. 2 is a schematic structural diagram of an omnidirectional mobile robot in an embodiment of the present invention.
Fig. 3 (a) is a position tracking graph of the tracking process of the omni-directional mobile robot in the embodiment of the present invention.
Fig. 3 (b) is a partial enlarged view of a position tracking curve of the tracking process of the omni-directional mobile robot in the embodiment of the present invention.
Fig. 4 (a) is a graph comparing X-axis displacement errors in the tracking process of the omni-directional mobile robot in the embodiment of the present invention.
Fig. 4 (b) is a graph comparing X-axis direction speed errors in the tracking process of the omni-directional mobile robot in the embodiment of the present invention.
Fig. 5 (a) is a graph comparing Y-axis displacement errors in the tracking process of the omni-directional mobile robot in the embodiment of the present invention.
Fig. 5 (b) is a graph comparing Y-axis directional velocity errors in the tracking process of the omni-directional mobile robot in accordance with the embodiment of the present invention.
Fig. 6 is an input torque diagram of a first mecanum wheel in the process of tracking the track of the omnidirectional mobile robot in an embodiment of the present invention.
Detailed Description
The following description of the embodiments of the present invention is provided to facilitate understanding of the present invention by those skilled in the art, but it should be understood that the present invention is not limited to the scope of the embodiments, and all the inventions which make use of the inventive concept are protected by the spirit and scope of the present invention as defined and defined in the appended claims to those skilled in the art.
Example 1
As shown in fig. 1 and fig. 2, the invention provides an omnidirectional mobile robot self-adaptive synovial membrane control method based on a nominal model, which comprises the following steps:
s1, establishing an omnidirectional mobile robot kinematics equation;
the step S1 includes the steps of:
s11, setting the geometric center of the omnidirectional mobile robot to coincide with the mass center, and defining a global coordinate system O f X f Y f Moving coordinate system O m X m Y m And a Mecanum wheel coordinate system O wi X wi Y wi Wherein i is equal to 1,2,3,4;
s12, decomposing the axle center speed of each Mecanum wheel of the omnidirectional mobile robot to obtain the speed of each Mecanum wheel respectively;
the omnidirectional mobile robot Mecanum wheel comprises a first Mecanum wheel, a second Mecanum wheel, a third Mecanum wheel and a fourth Mecanum wheel;
the speed expressions of the first Mecanum wheel and the third Mecanum wheel are as follows:
Figure BDA0003417234180000081
wherein v is ix Representing the speed of the Mecanum wheel of the omnidirectional mobile robot along the x direction, v iy Representing the speed of the Mecanum wheel of the omnidirectional mobile robot along the y direction, R represents the radius of the Mecanum wheel, phi represents the axis of the Mecanum robot body and the coordinate system x of the wheels wi An included angle between the shafts;
Figure BDA0003417234180000082
representing the angular velocity of the wheels of the mic robot; v oi Representing the speed of the mic robot body; i is equal to 1,3, wherein i represents the first Mecanum wheel when i is equal to 1 and the third Mecanum wheel when i is equal to 3;
the speed expressions of the second wheel and the fourth wheel are as follows:
Figure BDA0003417234180000091
wherein i is equal to 2,4, wherein i is equal to 2 for the second Mecanum wheel and i is equal to 4 for the fourth Mecanum wheel;
s13, obtaining the relation between the axle center speed of each Mecanum wheel and the main body center speed of the omnidirectional mobile robot;
in the step S13, the expression of the relationship between the axle center speed of the wheel of the omnidirectional mobile robot and the center speed of the main body of the omnidirectional mobile robot is as follows:
Figure BDA0003417234180000092
q m =[x m y m φ] T
wherein x is m Represents the transverse displacement of the omnidirectional mobile robot, y m Representing the longitudinal displacement of the omni-directional mobile robot,
Figure BDA0003417234180000093
represents the first derivative of the rotation angle phi, J represents a Jacobian matrix,/and>
Figure BDA0003417234180000094
represents a microphone robot body center position vector q m First derivative of l 1 And l 2 Respectively represent a coordinate system X for the movement of the center distance of wheels of the Mecanum robot m Axes and Y m The distance of the axis;
s14, obtaining an inverse kinematics equation of the omnidirectional mobile robot according to the relation between the axle center speed of each Mecanum wheel and the central speed of the main body of the omnidirectional mobile robot and the structural characteristics of the omnidirectional mobile robot;
the inverse kinematics model expression of the mic robot in the step S14 is as follows:
Figure BDA0003417234180000101
wherein J represents a Jacobian matrix,
Figure BDA0003417234180000102
represents a microphone robot body center position vector q m First derivative of l 1 And l 2 Respectively represent a coordinate system X for the movement of the center distance of wheels of the Mecanum robot m Axes and Y m The distance of the axis;
s15, obtaining an omnidirectional mobile robot kinematics equation according to the omnidirectional mobile robot inverse kinematics equation, and completing establishment of the omnidirectional mobile robot kinematics equation;
the expression of the kinematic equation of the omnidirectional mobile robot in the step S15 is as follows:
Figure BDA0003417234180000103
wherein,,
Figure BDA0003417234180000104
first derivative representing lateral displacement of omni mobile robot,/->
Figure BDA0003417234180000105
First derivative representing longitudinal displacement of omni mobile robot,/->
Figure BDA0003417234180000106
Representing the first derivative of the rotation angle of an omnidirectional mobile robot, R representing the radius of the Mecanum wheel, l 1 And l 2 Respectively represent a Mecanum wheel center distance movement coordinate system x m Axes and y m Distance of axis>
Figure BDA0003417234180000107
And
Figure BDA0003417234180000108
respectively representing the angular speeds of a first Mecanum wheel, a second Mecanum wheel, a third Mecanum wheel and a fourth Mecanum wheel of the omnidirectional mobile robot;
s2, establishing an omnidirectional mobile robot dynamics equation according to the omnidirectional mobile robot dynamics equation;
the step S2 includes the steps of:
s21, constructing an omnidirectional mobile robot kinetic energy equation according to the omnidirectional mobile robot kinetic energy equation;
the kinetic energy model expression of the omnidirectional mobile robot in the step S21 is as follows:
Figure BDA0003417234180000111
Figure BDA0003417234180000112
Figure BDA0003417234180000113
wherein E represents the kinetic energy of the omnidirectional mobile robot, m represents the mass of the omnidirectional mobile robot, v x Representing the transverse displacement speed, v, of an omnidirectional mobile robot y Indicating the longitudinal displacement speed of the omnidirectional mobile robot, I z Representing moment of inertia of the omnidirectional mobile robot around the z axis, I o Representing the moment of inertia of the omni-directional mobile robot about a center point,
Figure BDA0003417234180000114
representing the angular velocity of the wheels of the mic robot;
s22, calculating to obtain an omnidirectional mobile robot kinetic equation according to the Lagrangian kinematic equation and the omnidirectional mobile robot kinetic equation, and completing establishment of the omnidirectional mobile robot kinetic equation;
the expression of the dynamic equation of the omnidirectional mobile robot in the step S22 is as follows:
Figure BDA0003417234180000115
Figure BDA0003417234180000116
Figure BDA0003417234180000117
Figure BDA0003417234180000118
wherein M represents an inertial matrix of the omnidirectional mobile robot,
Figure BDA0003417234180000119
represents the angular acceleration of the wheels of the Mecanum robot,/->
Figure BDA00034172341800001110
The method comprises the steps of representing the friction moment of a Mecanum wheel, tau representing the total moment of the omnidirectional mobile robot, d representing the external disturbance moment, alpha representing a first inertia matrix parameter, and beta representing a second inertia matrix parameter.
S3, calculating to obtain a sliding mode surface S (t) and a sliding film control rate tau based on an omnidirectional mobile robot kinematics equation and a dynamics equation;
the step S3 includes the steps of:
s31, defining a tracking error e (t) of the omnidirectional mobile robot according to a kinematic equation and a dynamic equation of the omnidirectional mobile robot;
the expression of the tracking error e (t) of the omnidirectional mobile robot in the step S31 is as follows:
Figure BDA0003417234180000121
wherein,,
Figure BDA0003417234180000122
indicating the desired angular velocity of the Mecanum wheel, < >>
Figure BDA0003417234180000123
The actual angular velocities of the respective mecanum wheels are represented, where i=1, 2,3,4, respectively represent the angular velocities of the first, second, third and fourth mecanum wheels.
S32, calculating to obtain a sliding mode surface and a sliding film control rate according to the tracking error of the omnidirectional mobile robot;
the expressions of the slip-form surface S (t) and slip-film control rate τ' in step S32 are as follows:
Figure BDA0003417234180000124
Figure BDA0003417234180000125
wherein,,
Figure BDA0003417234180000126
derivative, κ, representing tracking error 4×4 Representing a constant positive definite matrix, e (t) representing tracking error,/->
Figure BDA0003417234180000127
Nominal parameters representing the inertial matrix of an omnidirectional mobile robot,/->
Figure BDA0003417234180000128
Indicating the angular acceleration of the wheels of the mecanum robot,
Figure BDA0003417234180000129
representing the derivative of the slide face,
Figure BDA00034172341800001210
Nominal parameter representing moment of inertia,/->
Figure BDA00034172341800001211
Representing the angular velocity of the Mecanum wheel, ρ representing the slip-mode control rate gain, sgn representing the step function;
defining a first Lyapunov function, deriving the function, and verifying the input torque jitter of the omnidirectional mobile robot;
the expression of the first lyapunov function is as follows:
Figure BDA00034172341800001212
wherein V is 1 Representing a first lyapunov function value;
the expression of the first Lyapunov function after derivation is as follows:
Figure BDA00034172341800001213
Figure BDA00034172341800001214
wherein,,
Figure BDA0003417234180000131
representing the function value of the first Lyapunov function value after derivation, <>
Figure BDA0003417234180000132
Representing the difference value between the nominal parameter of the inertia matrix of the omnidirectional mobile robot and the inertia matrix of the omnidirectional mobile robot;
after the first lyapunov function is derived, when the sliding mode control rate gain rho is large enough, the stability of the robot can be ensured, which can cause serious input torque jitter and requires the controller to adjust;
s4, setting a controller for adjusting the self-adaptive gain of the nominal model according to the sliding mode surface and the sliding film control rate, obtaining a dynamic equation of the controller omnidirectional mobile robot based on the self-adaptive gain adjustment of the nominal model, and completing the self-adaptive sliding film control of the omnidirectional mobile robot based on the nominal model;
the step S4 includes the steps of:
s41, setting a nominal model self-adaptive gain adjustment controller according to a sliding mode surface and a sliding film control rate to obtain a control rate tau' of the nominal model self-adaptive gain adjustment controller;
s42, obtaining a controller omnidirectional mobile robot dynamic model based on the nominal model adaptive gain adjustment according to the nominal parameters of the inertia matrix of the omnidirectional mobile robot, the nominal parameters of the friction coefficient of the omnidirectional mobile robot, the control rate tau' of the controller of the nominal model adaptive gain adjustment and the omnidirectional mobile robot dynamic model, and completing the omnidirectional mobile robot adaptive synovial membrane control based on the nominal model;
the dynamic model expression of the controller omnidirectional mobile robot based on the nominal model adaptive gain adjustment in the step S42 is as follows:
Figure BDA0003417234180000133
Figure BDA0003417234180000134
Figure BDA0003417234180000135
u r =-ρsgn(s(t))
wherein M represents an inertial matrix of the omnidirectional mobile robot,
Figure BDA0003417234180000136
represents the angular acceleration of the wheels of the Mecanum robot,/->
Figure BDA0003417234180000137
Representing the friction moment of the Mecanum wheel, +.>
Figure BDA0003417234180000138
Nominal parameters representing the inertial matrix of an omnidirectional mobile robot,/->
Figure BDA0003417234180000139
Nominal parameters representing moment of inertia, d representing external disturbance torque, τ' representing slip film control rate, +.>
Figure BDA0003417234180000141
Indicating jerk, κ of the wheels of a Mecanum robot 4×4 Representing a constant positive definite matrix, ">
Figure BDA0003417234180000142
Representing the derivative of the tracking error, e (t) representing the tracking error,/and (c)>
Figure BDA0003417234180000143
Representing Mecanum wheelsAngular velocity (V/V)>
Figure BDA0003417234180000144
Representing the nominal value of the external disturbance moment, u θ Representing the friction compensation term of the omnidirectional mobile robot, u r Representing a slip film approach law, T representing time, T representing a period, ρ representing a slip mode control rate gain, sgn representing a step function;
designing a second Lyapunov function, deriving the second Lyapunov function, and verifying convergence of a dynamic model of the controller omnidirectional mobile robot based on the nominal model adaptive gain adjustment:
the expression of the second lyapunov function is as follows:
Figure BDA0003417234180000145
wherein V represents the value of the second lyapunov function;
the expression of the second Lyapunov function after derivation is as follows:
Figure BDA0003417234180000146
Figure BDA0003417234180000147
wherein,,
Figure BDA0003417234180000148
a derivative value representing the second lyapunov function,/->
Figure BDA0003417234180000149
Representing the difference between the nominal parameter of the external disturbance torque and the external disturbance torque;
the derivative value of the second Lyapunov function
Figure BDA00034172341800001410
The controller based on the nominal model adaptive gain adjustment has good convergence, obtains high tracking precision, and simultaneously suppresses jitter of control input torque.
Example 2
As shown in fig. 3 (a) and 3 (b), in a practical example of the present invention, the present embodiment adopts a circular model as the trajectory tracking, and the circular model trajectory tracking equation expression is as follows:
Figure BDA0003417234180000151
wherein x represents lateral position tracking and y represents longitudinal position tracking;
setting the omnidirectional mobile robot to run from the origin, wherein the expression of the initial position of the omnidirectional mobile robot is as follows:
Figure BDA0003417234180000152
setting the disturbance d of each Mecanum wheel i Differences between nominal parameters of inertial matrix of omnidirectional mobile robot and inertial matrix of omnidirectional mobile robot
Figure BDA0003417234180000153
Nominal parameter of moment of inertia->
Figure BDA0003417234180000154
Difference between sense parameter and external disturbance moment +.>
Figure BDA0003417234180000155
The expressions are as follows:
d i =2sin(t+1)
Figure BDA0003417234180000156
Figure BDA0003417234180000157
Figure BDA0003417234180000158
setting a nominal model-based adaptive synovial membrane control related parameter: omnidirectional mobile robot mass m=6kg, moment of inertia I about z axis z =6.9210kg·m 2 Moment of inertia I about a point of center o =0.0945kg·m 2 Radius r=0.075 m of the mecanum wheel, center distance of the mecanum wheel moves coordinate system x m Distance of axis l 1 =0.26m, mecanum wheel center distance movement coordinate system y m Distance of axis l 2 =0.135 m, coefficient of friction β of each Mecanum wheel i Constant positive definite matrix κ=0.02 4×4 =diag (5 5 5 5), sliding mode control rate gain ρ=diag (30 30 30 30);
the error comparison of the PID controller, the traditional synovial membrane TSC and the NMASC controller under the conditions of sinusoidal interference, step interference and no interference is obtained through simulation, and the error comparison is shown in the table 1:
TABLE 1
Figure BDA0003417234180000159
Figure BDA0003417234180000161
The NMASC controller adopted by the scheme has the minimum error in track tracking effect.
As shown in fig. 4 (a), 4 (b), 5 (a) and 5 (b), in order to analyze the track tracking effect of the three controllers more clearly, tracking simulation is performed on the condition that sinusoidal interference exists on the omnidirectional mobile robot, according to fig. 4 (a) and 5 (a), the error of the PID controller in the X axis and the Y axis is obviously larger than the error of the adaptive slip film controller NMASC and the conventional slip film controller TSC based on the nominal model, and the PID controller generates larger jitter in the X axis and the Y axis directions, and the PID controller cannot have better convergence in the X axis and the Y axis; according to fig. 4 (b) and fig. 5 (b), the conventional sliding film TSC can achieve a better tracking effect on the desired track, but the conventional sliding film TSC generates obvious jitter on the X and Y axes, and the convergence speed of the NMASC model proposed herein is obviously better than that of the other two controllers.
As shown in fig. 6, the input torque of the wheel 1 at the initial time is greatly changed, and the input torque gradually becomes stable at [ -20 ] but after the adaptive adjustment, and remains within the range of [ -5 5 ].

Claims (5)

1. The self-adaptive synovial membrane control method of the omnidirectional mobile robot based on the nominal model is characterized by comprising the following steps of:
s1, establishing an omnidirectional mobile robot kinematics equation;
s2, establishing an omnidirectional mobile robot dynamics equation according to the omnidirectional mobile robot dynamics equation;
s3, calculating to obtain a sliding mode surface S (t) and a sliding film control rate tau' based on an omnidirectional mobile robot kinematic equation and a dynamic equation;
s4, setting a controller for adjusting the self-adaptive gain of the nominal model according to the sliding mode surface and the sliding film control rate, obtaining a dynamic equation of the controller omnidirectional mobile robot based on the self-adaptive gain adjustment of the nominal model, and completing the self-adaptive sliding film control of the omnidirectional mobile robot based on the nominal model;
the step S3 includes the steps of:
s31, defining a tracking error e (t) of the omnidirectional mobile robot according to a kinematic equation and a dynamic equation of the omnidirectional mobile robot;
s32, calculating to obtain a sliding mode surface and a sliding film control rate according to the tracking error of the omnidirectional mobile robot;
the expression of the tracking error e (t) of the omnidirectional mobile robot in the step S31 is as follows:
Figure FDA0004191701470000011
wherein,,
Figure FDA0004191701470000012
indicating the desired angular velocity of the Mecanum wheel, < >>
Figure FDA0004191701470000013
Representing the actual angular velocity of each mecanum wheel, wherein i=1, 2,3,4, represents the angular velocities of the first, second, third and fourth mecanum wheels, respectively;
the expressions of the slip-form surface S (t) and slip-film control rate τ' in step S32 are as follows:
Figure FDA0004191701470000014
Figure FDA0004191701470000015
wherein,,
Figure FDA0004191701470000016
derivative, κ, representing tracking error 4×4 Representing a constant positive definite matrix, e (t) representing tracking error,/->
Figure FDA0004191701470000021
Nominal parameters representing the inertial matrix of an omnidirectional mobile robot,/->
Figure FDA0004191701470000022
Represents the angular acceleration of the wheels of the Mecanum robot,/->
Figure FDA0004191701470000023
Representation ofDerivative of the sliding surface,
Figure FDA0004191701470000024
Nominal parameter representing moment of inertia,/->
Figure FDA0004191701470000025
Representing the angular velocity of the Mecanum wheel, ρ representing the slip-mode control rate gain, sgn representing the step function;
the step S4 includes the steps of:
s41, setting a nominal model self-adaptive gain adjustment controller according to a sliding mode surface and a sliding film control rate to obtain a control rate tau' of the nominal model self-adaptive gain adjustment controller;
s42, obtaining a controller omnidirectional mobile robot dynamic model based on the nominal model adaptive gain adjustment according to the nominal parameters of the inertia matrix of the omnidirectional mobile robot, the nominal parameters of the friction coefficient of the omnidirectional mobile robot, the control rate tau' of the controller of the nominal model adaptive gain adjustment and the omnidirectional mobile robot dynamic model, and completing the omnidirectional mobile robot adaptive synovial membrane control based on the nominal model;
the dynamic model expression of the controller omnidirectional mobile robot based on the nominal model adaptive gain adjustment in the step S42 is as follows:
Figure FDA0004191701470000026
Figure FDA0004191701470000027
Figure FDA0004191701470000028
u r =-ρsgn(s(t))
wherein M represents an omnidirectional mobile robotThe inertia matrix is used to determine the mass of the material,
Figure FDA0004191701470000029
represents the angular acceleration of the wheels of the Mecanum robot,/->
Figure FDA00041917014700000210
Representing the friction moment of the Mecanum wheel, +.>
Figure FDA00041917014700000211
Nominal parameters representing the inertial matrix of an omnidirectional mobile robot,/->
Figure FDA00041917014700000212
Nominal parameters representing moment of inertia, d representing external disturbance torque, τ' representing slip film control rate, +.>
Figure FDA00041917014700000213
Indicating jerk, κ of the wheels of a Mecanum robot 4×4 Representing a constant positive definite matrix, ">
Figure FDA00041917014700000214
Representing the derivative of the tracking error, e (t) representing the tracking error,/and (c)>
Figure FDA00041917014700000215
Represents the angular velocity of the Mecanum wheel, +.>
Figure FDA00041917014700000216
Representing the nominal value of the external disturbance moment, u θ Representing the friction compensation term of the omnidirectional mobile robot, u r The slip film approach law is represented, T represents time, T represents period, ρ represents slip mode control rate gain, sgn represents step function.
2. The method for controlling the adaptive synovial membrane of the omni-directional mobile robot based on the nominal model according to claim 1, wherein said step S1 comprises the steps of:
s11, setting the geometric center of the omnidirectional mobile robot to coincide with the mass center, and defining a global coordinate system O f X f Y f Moving coordinate system O m X m Y m And a Mecanum wheel coordinate system O wi X wi Y wi Wherein i is equal to 1,2,3,4;
s12, decomposing the axle center speed of each Mecanum wheel of the omnidirectional mobile robot to obtain the speed of each Mecanum wheel respectively;
s13, obtaining the relation between the axle center speed of each Mecanum wheel and the main body center speed of the omnidirectional mobile robot;
s14, obtaining an inverse kinematics equation of the omnidirectional mobile robot according to the relation between the axle center speed of each Mecanum wheel and the central speed of the main body of the omnidirectional mobile robot and the structural characteristics of the omnidirectional mobile robot;
and S15, obtaining an omnidirectional mobile robot kinematics equation according to the omnidirectional mobile robot inverse kinematics equation, and completing establishment of the omnidirectional mobile robot kinematics equation.
3. The method for adaptive synovial membrane control of an omni mobile robot based on a nominal model according to claim 2, wherein the expression of the kinematic equation of the omni mobile robot in step S15 is as follows:
Figure FDA0004191701470000031
wherein,,
Figure FDA0004191701470000032
first derivative representing lateral displacement of omni mobile robot,/->
Figure FDA0004191701470000033
First derivative representing longitudinal displacement of omni mobile robot,/->
Figure FDA0004191701470000034
Representing the first derivative of the rotation angle of an omnidirectional mobile robot, R representing the radius of the Mecanum wheel, l 1 And l 2 Respectively represent a Mecanum wheel center distance movement coordinate system x m Axes and y m Distance of axis>
Figure FDA0004191701470000035
And->
Figure FDA0004191701470000036
The angular velocities of the first, second, third and fourth mecanum wheels of the omnidirectional mobile robot are represented, respectively.
4. The method for controlling the adaptive synovial membrane of the omni-directional mobile robot based on the nominal model according to claim 1, wherein said step S2 comprises the steps of:
s21, constructing an omnidirectional mobile robot kinetic energy equation according to the omnidirectional mobile robot kinetic energy equation;
s22, calculating to obtain an omnidirectional mobile robot kinetic equation according to the Lagrangian kinematic equation and the omnidirectional mobile robot kinetic equation, and completing establishment of the omnidirectional mobile robot kinetic equation.
5. The method for adaptive synovial membrane control of an omni mobile robot based on a nominal model according to claim 4, wherein the expression of the kinetic equation of the omni mobile robot in step S22 is as follows:
Figure FDA0004191701470000041
wherein M represents an inertial matrix of the omnidirectional mobile robot,
Figure FDA0004191701470000042
represents the angular acceleration of the wheels of the Mecanum robot,/->
Figure FDA0004191701470000043
The friction moment of the Mecanum wheel is represented, tau represents the total moment of the omnidirectional mobile robot, and d represents the external disturbance moment.
CN202111551149.5A 2021-12-17 2021-12-17 Omnidirectional mobile robot self-adaptive sliding film control method based on nominal model Active CN114488790B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111551149.5A CN114488790B (en) 2021-12-17 2021-12-17 Omnidirectional mobile robot self-adaptive sliding film control method based on nominal model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111551149.5A CN114488790B (en) 2021-12-17 2021-12-17 Omnidirectional mobile robot self-adaptive sliding film control method based on nominal model

Publications (2)

Publication Number Publication Date
CN114488790A CN114488790A (en) 2022-05-13
CN114488790B true CN114488790B (en) 2023-07-04

Family

ID=81494529

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111551149.5A Active CN114488790B (en) 2021-12-17 2021-12-17 Omnidirectional mobile robot self-adaptive sliding film control method based on nominal model

Country Status (1)

Country Link
CN (1) CN114488790B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115635485B (en) * 2022-11-09 2024-03-15 嘉兴智康机器人有限公司 Real-time human-computer interaction force control method of mobile rehabilitation robot

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112091829A (en) * 2020-08-31 2020-12-18 江苏大学 Sand blasting and rust removing parallel robot friction force mutation compensating fuzzy self-adaptive sliding mode control method

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2866585C (en) * 2012-03-06 2021-02-23 Sirius Xm Radio Inc. Systems and methods for audio attribute mapping
CN105607122B (en) * 2015-12-23 2018-04-10 西南科技大学 A kind of earthquake texture blending and Enhancement Method based on full variation geological data decomposition model
CN106647584A (en) * 2017-01-17 2017-05-10 南京航空航天大学 Fault tolerant control method of four-rotor-wing aircraft based on optimal sliding mode
CN108196607B (en) * 2017-12-29 2019-11-05 西南科技大学 The temprature control method of the electromagnetic induction heating system of uncertain stochastic Time-Delay
CN108614426B (en) * 2018-06-06 2021-03-16 黄山学院 Multi-mobile-robot formation robust control method based on disturbance observer
CN112068428B (en) * 2020-08-31 2022-05-17 五邑大学 Design method and system of double-pendulum PI type Terminal sliding mode controller of bridge crane
CN112180729B (en) * 2020-10-09 2023-03-03 电子科技大学中山学院 Self-adaptive trajectory tracking control method for mechanical arm
CN112904872B (en) * 2021-01-19 2021-12-10 哈尔滨工程大学 Benthic AUV (autonomous Underwater vehicle) fixed-time rapid trajectory tracking control method based on disturbance approximation
CN113697131B (en) * 2021-08-27 2022-04-22 哈尔滨工业大学(深圳) Unwinding-resistant sliding mode attitude tracking control method and system for rigid spacecraft

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112091829A (en) * 2020-08-31 2020-12-18 江苏大学 Sand blasting and rust removing parallel robot friction force mutation compensating fuzzy self-adaptive sliding mode control method

Also Published As

Publication number Publication date
CN114488790A (en) 2022-05-13

Similar Documents

Publication Publication Date Title
CN108919837B (en) Second-order sliding mode control method of automatic driving vehicle based on visual dynamics
Yong et al. Balance control of two-wheeled self-balancing mobile robot based on TS fuzzy model
CN108459605B (en) Trajectory tracking control method based on AGV system
CN109857100B (en) Composite track tracking control algorithm based on inversion method and fast terminal sliding mode
CN110083061A (en) A kind of Control of Wheeled Mobile Robots method
CN111158376B (en) Rocker rocker arm type planet car trajectory tracking coordination control method in soft and rugged terrain
CN111665855A (en) Robot presetting formation control method
CN111208830B (en) Three-closed-loop formation track tracking control method for wheeled mobile robot
CN114488790B (en) Omnidirectional mobile robot self-adaptive sliding film control method based on nominal model
CN112606843A (en) Intelligent vehicle path tracking control method based on Lyapunov-MPC technology
Liu et al. Direction and trajectory tracking control for nonholonomic spherical robot by combining sliding mode controller and model prediction controller
CN115202213A (en) Four-rotor aircraft control method based on active disturbance rejection control
CN112936277A (en) Underwater robot-manipulator system fixed time trajectory tracking method
CN115480583A (en) Visual servo tracking and impedance control method of flying operation robot
Ni et al. Path tracking motion control method of tracked robot based on improved LQR control
CN112008728A (en) Omnidirectional mobile robot track tracking control method based on selective disturbance compensation
CN115268475B (en) Robot fish accurate terrain tracking control method based on finite time disturbance observer
Huo et al. Design and validation of a novel adaptive motion control for a pendulum spherical robot
Ma et al. Finite-time trajectory tracking control of quadrotor UAV via adaptive RBF neural network with lumped uncertainties
CN109648566A (en) The Trajectory Tracking Control method of the unknown all directionally movable robot of the parameter of electric machine
CN113867136A (en) Steel structure damage detection robot attitude control method based on genetic algorithm and fuzzy fractional order PID
Bai et al. Trajectory tracking control of hovercraft based on Shutting-Backstepping Sliding Model
Li et al. Global fast terminal sliding mode control for a quadrotor UAV
CN117850433B (en) Collision avoidance formation control method and system for disturbed communication limited mobile robot
CN115129046B (en) Automatic driving vehicle path tracking method based on sliding mode neural network control

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