CN115686034A - Unmanned underwater vehicle track tracking control method considering failure of speed sensor - Google Patents

Unmanned underwater vehicle track tracking control method considering failure of speed sensor Download PDF

Info

Publication number
CN115686034A
CN115686034A CN202211453151.3A CN202211453151A CN115686034A CN 115686034 A CN115686034 A CN 115686034A CN 202211453151 A CN202211453151 A CN 202211453151A CN 115686034 A CN115686034 A CN 115686034A
Authority
CN
China
Prior art keywords
formula
control law
observer
speed
longitudinal
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.)
Granted
Application number
CN202211453151.3A
Other languages
Chinese (zh)
Other versions
CN115686034B (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.)
Hefei Jinglong Environmental Protection Technology Co ltd
Original Assignee
Harbin 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 Harbin University of Science and Technology filed Critical Harbin University of Science and Technology
Priority to CN202211453151.3A priority Critical patent/CN115686034B/en
Publication of CN115686034A publication Critical patent/CN115686034A/en
Application granted granted Critical
Publication of CN115686034B publication Critical patent/CN115686034B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

The invention provides a track tracking control method of an unmanned underwater vehicle considering failure of a speed sensor, which comprises the following steps: establishing a navigation mathematical model, and setting an expected track mathematical model; designing a speed observer based on a position signal, and designing a longitudinal moment and transverse moment disturbance observer on the basis of the speed observer; designing virtual control laws of longitudinal speed and transverse speed, introducing an instruction filter to constrain the amplitude and speed of the virtual control laws, and solving the problem of increased computational complexity caused by high-order derivation of the virtual control laws; designing a trajectory tracking sliding mode controller based on the virtual control law and the output of the instruction filter, and compensating the disturbance of longitudinal torque and transverse torque based on the estimation value of the disturbance observer; the technical scheme of the invention solves the problems that the performance failure of the speed sensor of the unmanned underwater vehicle after underwater operation is not considered in the prior art, and the tracking failure of the controller is caused because the torque output of the controller exceeds the maximum output of the propeller of the unmanned underwater vehicle.

Description

Unmanned underwater vehicle track tracking control method considering failure of speed sensor
Technical Field
The invention relates to the technical field of track tracking control of unmanned underwater vehicles, in particular to a track tracking control method of an unmanned underwater vehicle considering failure of a speed sensor.
Background
In recent years, the motion control of an under-actuated unmanned underwater vehicle has become a research hotspot of ocean engineering, and the unmanned underwater vehicle plays an important role in oceanographic observation, submarine resource exploration, hydrological measurement and mapping, submarine rescue and the like, and has wide application prospect. With the improvement of the automation degree, the unmanned submersible vehicle is deployed and recovered by using the unmanned boat without depending on manpower, so that a new research hotspot is formed. In order to complete automatic deployment and recovery of the unmanned underwater vehicle, the unmanned underwater vehicle must realize track tracking, particularly high-precision water surface track tracking. Therefore, research on trajectory tracking during deployment and recovery of unmanned underwater vehicles becomes more and more important. However, currently, the research on trajectory tracking of the unmanned underwater vehicle on the water surface is relatively lacking, and in the prior art, the trajectory tracking of the unmanned underwater vehicle in the automatic deployment and recovery process has the defects that the influence caused by wind, wave and flow in the navigation process or the interference caused by unmodeled parts outside the system cannot be accurately estimated, the interference is not compensated, the accuracy of a speed sensor of the unmanned underwater vehicle after completing underwater operation is greatly influenced or even fails, the expected output of a controller is not limited according to the dynamic performance of the unmanned underwater vehicle, and the like. These disadvantages are more evident especially in environments where sea surface wave disturbances are complex and uncertain.
The existing research on the water surface track tracking method in the laying and recovery process of the unmanned underwater vehicle mainly has the following defects:
1) The accuracy of the speed sensor of the unmanned underwater vehicle is greatly influenced or even fails after the unmanned underwater vehicle finishes underwater operation. The track Tracking Control methods designed by the backstepping Sliding Mode, such as the papers of track Tracking Control of an under-utilized UUV Using a Novel Nonlinear Integral Tracking Mode Surface, and dynamic Tracking Mode Control of the track Tracking of under-utilized underwater vehicles, do not consider that the influence of the sensor precision after the underwater operation of the unmanned underwater vehicle can cause the Control precision to be reduced or even fail. The traditional track tracking control method designed for the unmanned ship on the water surface does not consider the influence of the precision of the speed sensor on track tracking, for example, the limitation on the virtual control law is considered in the article of 'under-actuated ship track tracking instruction filtering sliding mode robust control', but the influence caused by the failure of the speed sensor is not considered.
2) The traditional track tracking control method designed for the unmanned ship on the water surface does not consider the dynamic performance of the unmanned underwater vehicle, does not limit the expected output of the controller, and does not solve the problem of increased computational complexity caused by calculating a high-order derivative by a virtual control law. For example, the control method proposed in article "under-actuated ship adaptive dynamic surface output feedback trajectory tracking control with nonlinear observer" does not impose limitation on the desired output of the controller, and does not avoid the problem of increased computational complexity caused by high-order derivation of the virtual control law in the process of calculating the control law. Similarly, the control method proposed in patent No. CN113835338A, CN110134012a and the like is also not applicable to the water surface trajectory tracking control of the unmanned underwater vehicle because the desired output is not limited.
3) The influence caused by wind, wave and flow or the interference generated by unmodeled parts outside the system is not observed in the sailing process, and the interference is not compensated. For example, the control method proposed in the article "underwater robot position tracking control based on command filtering" does not observe and compensate for various disturbances of the system.
Disclosure of Invention
In order to solve the problems of track tracking control in the process of laying and recovering the unmanned underwater vehicle in the prior art, the invention considers the track tracking problem of the unmanned underwater vehicle under the condition that the unmanned underwater vehicle cannot measure a speed vector, simultaneously considers the interference generated by the internal part, the external part and the unmodeled part of the system, limits the output of the controller by restricting the virtual control law in order to prevent the expected output of the controller from exceeding the upper limit of the thrust of the unmanned underwater vehicle driver, and provides a corresponding control method.
In order to achieve the purpose, the invention is realized by the following technical scheme:
the track tracking control method of the unmanned underwater vehicle considering the failure of the speed sensor is characterized by comprising the following steps of:
s1, establishing a water surface navigation mathematical model of the unmanned underwater vehicle, and setting an expected track mathematical model.
S2, designing a speed vector observer based on the position signal, designing a disturbance observer in the longitudinal moment and transverse moment directions on the basis, and verifying the validity of the disturbance observer through Lyapunov theorem.
And S3, designing a virtual control law of longitudinal speed and transverse speed based on the real-time position error, introducing an instruction filter to constrain the amplitude and the speed of the virtual control law on the basis of the virtual control law, solving the problem that the calculation complexity is increased due to the fact that high-order derivation needs to be carried out on the virtual control law when the actual control law is solved, and verifying the effectiveness of the virtual control law through the Lyapunov theorem.
And S4, designing a track tracking sliding mode controller of the unmanned underwater vehicle based on the virtual control law and the output value of the instruction filter, compensating the disturbance of the longitudinal moment and the transverse moment based on the system disturbance estimation value observed by the disturbance observer, and verifying the effectiveness of the track tracking sliding mode controller through Lyapunov theorem.
And S5, verifying the effectiveness of the proposed method through simulation.
Further, the step S1 specifically includes:
defining a northeast coordinate system OXYZ and an attached coordinate system O respectively B X B Y B Z B And two coordinate systems, wherein a northeast coordinate system XOY is set as an inertial coordinate system, an arbitrary point O on the earth is taken as an origin of the coordinate system, OX points to the true north, OY points to the true east, and OZ is vertical to the plane where OX and OY are located and faces downwards. Will attach body coordinate system O B X B Y B Z B As an inertial coordinate system, assuming that the unmanned underwater vehicle body is symmetrical left and right, the central point is taken as the origin O of the inertial coordinate system B ,O B X B The axis points along the hullDirection of boat bow, O B Y B Directed vertically to starboard, O B Z B Perpendicular to the hull and down. And performing dynamic modeling on the unmanned underwater vehicle in the water surface navigation process to obtain a six-degree-of-freedom model comprising surging (surge), swaying (sway), heaving (heave), rolling (roll), pitching (pitch) and yawing (yaw). Neglecting the hydrodynamic resistance term higher than the second order and the motions of heave, pitch and roll, and establishing a kinematic and dynamic model as follows:
Figure BDA0003952362400000031
in the formula, eta = [ x, y, psi =] T X, y represents the unmanned vehicle motion position vector, psi ∈ [0,2 pi]Representing a heading angle, J (ψ) is a conversion matrix of a northeast coordinate system and an attached body coordinate system, v = [ u, v, r] T The velocity vector of the motion of the unmanned underwater vehicle in the attached coordinate system is represented, and u, v and r respectively represent surging, swaying and yawing velocities.
Figure BDA0003952362400000032
Representing an inertia matrix containing additional mass, C (v) representing a diagonally symmetric matrix,
Figure BDA0003952362400000033
denotes a damping matrix, τ = [ ] uvr ] T Representing moment control input vector, τ u 、τ v 、τ r Respectively representing the control forces of surge, sway and yaw, tau d =[τ dudvdr ] T System disturbance torque vector, τ du 、τ dv 、τ dr Representing the disturbance moments in the surge, sway, and yaw directions, respectively.
Further, the step S2 specifically includes:
firstly, designing a speed vector observer based on a position signal, eliminating a C (v) v term with a speed vector v in a model because the speed vector v cannot be measured, and introducing a variable matrix P (eta, v):
P(η,v)=H(η)v (2)
in the formula, H (eta) 3×3 Is a mapping matrix of velocity vectors v to P (η, v).
The derivation of equation (2) yields:
Figure BDA0003952362400000041
to eliminate the term with velocity vector v in equation (3), the variable Y is introduced:
Figure BDA0003952362400000042
to make Y =0 and thereby eliminate the term containing the velocity vector v in equation (3), let:
Figure BDA0003952362400000043
thus, there are:
P=Y-H(η)M -1 DH -1 (η)P+H(η)M -1 (τ+τ d ) (6)
and then designing a speed observer as follows:
Figure BDA0003952362400000044
in the formula (I), the compound is shown in the specification,
Figure BDA0003952362400000045
is the observed value of eta, P of the observer, K 1 、K 2 A matrix is designed for the gain of the signal,
Figure BDA0003952362400000046
is the observation error of the observer.
To prove the effectiveness of the designed speed vector observer, the observation error of the observer is defined as:
Figure BDA0003952362400000047
in the formula (I), the compound is shown in the specification,
Figure BDA0003952362400000048
in order to observe the error of the observer,
Figure BDA0003952362400000049
is the position error of the observer observations.
Defining the Lyapunov function:
Figure BDA0003952362400000051
in the formula, x 1 、x 2 The matrix is designed for positive definite. The derivation is carried out on the formula (9) and is obtained from the formula (8):
Figure BDA0003952362400000052
and order:
x 1 J(ψ)H -1 (η)=x 2 K 2 (11)
then there are:
Figure BDA0003952362400000053
as defined by Lyapunov's theorem, the observer shown in equation (7) converges asymptotically.
On the basis of designing a speed observer, designing a disturbance observer for an unknown disturbance part, and carrying out online estimation on disturbance in the surge and heading directions according to the formula (6):
Figure BDA0003952362400000054
in the formula (I), the compound is shown in the specification,
Figure BDA0003952362400000055
for a disturbance observer, carrying out unknown disturbance tau on unmodeled parts and inside and outside of the system d Observed value of, K 0 ∈R 3×3 Is a design matrix composed of positive definite parameters.
In the actual situation where the temperature is too high,
Figure BDA0003952362400000056
failure to pass the observed value
Figure BDA0003952362400000057
Directly obtaining, and thus introducing, auxiliary variables
Figure BDA0003952362400000058
And then designing a disturbance observer as follows:
Figure BDA0003952362400000059
in the formula, unknown disturbance tau is generated for the disturbance observer to the unmodeled part and the inside and outside of the system d Observed value of, K 0 ∈R 3×3 Is a design matrix composed of positive definite parameters.
To prove the effectiveness of the designed moment disturbance observer, a Lyapunov function is defined:
Figure BDA00039523624000000510
in the formula (I), the compound is shown in the specification,
Figure BDA00039523624000000511
τ d is the true value of the disturbance of the system,
Figure BDA00039523624000000512
an error is estimated for the parameter data.
The derivation of equation (15) can be:
Figure BDA0003952362400000061
from this, the observation error of the disturbance observer expressed by the equation (14) can be known
Figure BDA0003952362400000062
Converging progressively to zero in a finite time.
Further, the step S3 specifically includes:
and designing a virtual control law of longitudinal speed and transverse speed based on the real-time position error. First, the error between the actual position and the desired trajectory position is defined as:
Figure BDA0003952362400000063
in the formula, x d As the abscissa, y, of the desired trajectory d Is the ordinate of the desired trajectory, x is the abscissa of the actual tracking position, and y is the ordinate of the actual tracking position.
And (3) solving a first derivative of the formula (17), and combining the unmanned submersible vehicle kinematics and the dynamic model formula (1) to obtain:
Figure BDA0003952362400000064
in order to make the position tracking error converge to zero, a vertical virtual control law a is designed ud And a horizontal virtual rate a vd Comprises the following steps:
Figure BDA0003952362400000065
in the formula (I), the compound is shown in the specification,
Figure BDA0003952362400000066
wherein k > 0,c > 0 is a design constant.
To prevent the output of the virtual control from exceeding the maximum value of the torque which can be provided by the propeller of the unmanned underwater vehicle and solve the problem that the output of the virtual control exceeds the maximum value of the torque which can be provided by the propeller of the unmanned underwater vehicleIn the process of solving the actual control law, the high-order derivation of the virtual control law is required to generate the problem of increased calculation complexity, the speed and the amplitude of the virtual control law are limited by using an instruction filter, and the longitudinal virtual control law a is required to be solved ud And a horizontal virtual rate a vd As input to the command filter, the output of the command filter
Figure BDA0003952362400000067
And
Figure BDA0003952362400000068
and the method is used for solving the actual control law.
Definition of instruction filter:
in order to solve the problem of complexity expansion of differentiation which is possibly generated by derivation of the virtual control law when the actual control law is obtained, reduce the workload of reverse calculation and prevent the control law from being invalid due to input saturation, the amplitude and the speed of the virtual control law are limited; the instruction filter can also handle constraints of intermediate state variables, as compared to a first order filter. Therefore, the command filter is introduced to process the virtual control law.
The state space expression for the instruction filter is defined as:
Figure BDA0003952362400000071
in the formula, x id As input to the filter, z 1 For the output of this filter, ζ and ω are the damping and bandwidth, respectively, of the instruction filter. When the effect of amplitude and velocity is not taken into account, x id And z 1 The linear relationship between them is:
Figure BDA0003952362400000072
when inputting x id When bounded, then z is output 1
Figure BDA0003952362400000073
Is also bounded, and in the case where ω is large, | z 1 -x id I tends to zero.
To prove the effectiveness of the virtual control law, a Lyapunov function is defined:
Figure BDA0003952362400000074
derived from formula (19):
Figure BDA0003952362400000075
substituting equation (22) into equation (18) has:
Figure BDA0003952362400000076
when (u-a) ud ) And (v-a) vd ) Towards 0:
Figure BDA0003952362400000077
the derivation of equation (22) in combination with equation (25) yields:
Figure BDA0003952362400000078
from this, it can be seen that the virtual control law shown in equation (19) can control the position error x e 、y e Asymptotically converging to zero. The virtual control law may also be considered a speed expectation for the unmanned vehicle to reach the target position.
Further, the step S4 specifically includes:
designing a longitudinal thrust control law and a transverse thrust control law; firstly, designing a longitudinal thrust control law, and setting an error u between a longitudinal speed and an expected longitudinal speed e Comprises the following steps:
Figure BDA0003952362400000081
in the formula (I), the compound is shown in the specification,
Figure BDA0003952362400000082
for the observation of the longitudinal velocity u by the velocity observer, a ud Is a virtual control law of longitudinal speed.
Designing integral slip-form surface s for longitudinal speed error 1
Figure BDA0003952362400000083
In the formula, mu 1 > 0 is a constant.
In order to further obtain the longitudinal thrust control law, the derivation is carried out on the formula (28) and the formula (1) is combined to obtain:
Figure BDA0003952362400000084
in the formula (I), the compound is shown in the specification,
Figure BDA0003952362400000085
is an estimate of the observer's velocity vector u, v, r, τ u Is the actual control law of the longitudinal thrust,
Figure BDA0003952362400000086
for virtual control law a ud Output value, m, through instruction filter 11 、m 22 、d 11 The coefficients in the inertia matrix M and the damping matrix D of the additional mass in equation (1) are shown, respectively.
Further, the longitudinal thrust control law τ is obtained from the equations (29) and (14) u Comprises the following steps:
Figure BDA0003952362400000087
in the formula, λ 1 >0;
Figure BDA0003952362400000088
Is the observed compensation value of the interference observer to the longitudinal disturbance in the equation (7).
Law of longitudinal thrust control designed to demonstrate u Defining the Lyapunov function:
Figure BDA0003952362400000089
the derivation is obtained by the derivation of equation (31) and the substitution of equations (29) and (30):
Figure BDA00039523624000000810
since when t → ∞ is reached,
Figure BDA00039523624000000811
therefore, the following are provided:
Figure BDA0003952362400000091
thus signal s 1 The agreement is eventually bounded.
Secondly, designing a transverse thrust control law, and setting the error between the transverse speed and the expected transverse speed as follows:
Figure BDA0003952362400000092
because the power system of the unmanned underwater vehicle is an under-actuated system and has no thrust output in the transverse direction, the steering torque tau is required r The following linear sliding mode surfaces are designed:
Figure BDA0003952362400000093
in the formula, mu 2 > 0 is a constant.
To further determine the lateral thrust control law, equation (35) is derived with respect to time:
Figure BDA0003952362400000094
in the formula (I), the compound is shown in the specification,
Figure BDA0003952362400000095
for virtual control law a vd A first derivative output value that passes through the instruction filter,
Figure BDA0003952362400000096
for virtual control law a vd The second derivative output value of the instruction filter. Also, in combination with formula (1):
Figure BDA0003952362400000097
in the formula, τ r Is the actual control law of the transverse direction moment. Combining formula (37) with formula (1), formula (36) can be expressed as:
Figure BDA0003952362400000098
from formula (38) in combination with formula (14):
Figure BDA0003952362400000099
in the formula, λ 2 >0;
Figure BDA00039523624000000910
And (4) compensating the observed lateral disturbance value for the disturbance observer.
To demonstrate the designed transverse direction moment control law tau r Defining the Lyapunov function:
Figure BDA0003952362400000101
the derivation of equation (40) and the substitution of equations (38) and (39) yields:
Figure BDA0003952362400000102
since when t → ∞ is reached,
Figure BDA0003952362400000103
therefore, there are:
Figure BDA0003952362400000104
thus signal s 2 The agreement is eventually bounded.
Compared with the prior art, the invention has the beneficial effects that:
1. the invention fully considers the problems that in the recovery process of the unmanned underwater vehicle after the underwater task is finished, the speed sensor is possibly influenced by pressure and a humid environment to cause the reduction of observation precision, and the reduction of control precision of a control law is caused, even the control law is invalid, and specially designs the speed observer based on the position vector. The control effect is not influenced by pressure intensity and a moist environment.
2. The influence of internal and external disturbance of the system and part of disturbance which is not modeled on the system on the laying and recovery process of the unmanned underwater vehicle is fully considered, a disturbance observer of longitudinal moment and transverse moment is designed based on a speed observer model, and the control law is compensated based on the observation value of the disturbance observer in actual control.
3. In consideration of the performance of the unmanned underwater vehicle driver, the expected output of the control law is limited by the command filter, compared with the prior art, the control law designed by the invention takes longer time for tracking the expected track, but the control law is advantageous in that the control law is prevented from failing because the output torque exceeds the maximum output performance of the driver.
Drawings
In order to more clearly illustrate the technical solutions of the prior art and the embodiments of the present invention, the drawings used in the description of the prior art and the embodiments will be briefly introduced below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
FIG. 1 is a system block diagram of the present invention;
FIG. 2 is a schematic diagram of kinematics and dynamics modeling of the unmanned underwater vehicle in the invention;
FIG. 3 is a schematic diagram of the trajectory tracking effect of the present invention;
FIG. 4 is a graph comparing the expected longitudinal position response to the actual longitudinal position response of the present invention;
FIG. 5 is a graph comparing the expected longitudinal position response to the actual longitudinal position response of the present invention;
FIG. 6 is a comparison of the observed value of longitudinal velocity and the actual value of longitudinal velocity by the observer according to the present invention;
FIG. 7 is a comparison of the observed value of the lateral velocity and the actual value of the lateral velocity by the observer of the present invention;
FIG. 8 is a comparison graph of the observed value of the heading angular velocity and the actual value of the heading angular velocity by the observer according to the present invention;
FIG. 9 is a graph of the variation of the vertical virtual control law according to the present invention;
FIG. 10 is a graph of the variation of the horizontal virtual control law of the present invention;
FIG. 11 is a graph of the change in longitudinal speed control law input according to the present invention;
FIG. 12 is a lateral velocity control law input profile according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. The following description of at least one exemplary embodiment is merely illustrative in nature and is in no way intended to limit the invention, its application, or uses. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As shown in fig. 1, the present invention provides a track following control method for an unmanned underwater vehicle considering the failure of a speed sensor, comprising the following steps:
s1, establishing a water surface navigation mathematical model of the unmanned underwater vehicle, and setting an expected track mathematical model.
S2, designing a speed vector observer based on the position signal, designing a disturbance observer in the longitudinal moment and transverse moment directions on the basis, and verifying the validity of the disturbance observer through Lyapunov theorem.
And S3, designing a virtual control law of longitudinal speed and transverse speed based on the real-time position error, introducing an instruction filter to constrain the amplitude and the speed of the virtual control law on the basis of the virtual control law, solving the problem that the calculation complexity is increased due to the fact that high-order derivation needs to be carried out on the virtual control law when the actual control law is solved, and verifying the effectiveness of the virtual control law through the Lyapunov theorem.
And S4, designing a track tracking sliding mode controller of the unmanned underwater vehicle based on the virtual control law and the output value of the instruction filter, compensating the disturbance of the longitudinal moment and the transverse moment based on the system disturbance estimation value observed by the disturbance observer, and verifying the effectiveness of the track tracking sliding mode controller through Lyapunov theorem.
And S5, verifying the effectiveness of the proposed method through simulation.
The step S1 specifically includes:
respectively defining a northeast coordinate system OXYZ and an attached coordinate system O B X B Y B Z B And two coordinate systems, wherein a northeast coordinate system XOY is set as an inertial coordinate system, an arbitrary point O on the earth is taken as an origin of the coordinate system, OX points to the true north, OY points to the true east, and OZ is vertical to the plane where OX and OY are located and faces downwards. Will enclose body coordinate system O B X B Y B Z B As an inertial coordinate system, assuming that the unmanned underwater vehicle body is symmetrical left and right, the central point is taken as the origin O of the inertial coordinate system B ,O B X B The shaft points to the boat bow direction along the boat body, O B Y B Directed vertically to starboard, O B Z B Perpendicular to the hull and down. The dynamics modeling is carried out on the unmanned underwater vehicle in the water surface navigation process, and a six-degree-of-freedom model comprising surging (surge), swaying (sway), heaving (heave), rolling (roll), pitching (pitch) and yawing (yaw) is obtained and is shown in fig. 2. Neglecting hydrodynamic resistance terms higher than the second order and heave (yaw), pitch (pitch) and roll (roll) motions, and establishing a kinematic and dynamic model as follows:
Figure BDA0003952362400000121
in the formula, eta = [ x, y, psi =] T X, y represents the unmanned vehicle motion position vector, psi ∈ [0,2 pi]Representing the heading angle, J (ψ) is a conversion matrix of the northeast coordinate system and the satellite coordinate system, v = [ u, v, r] T The velocity vector of the motion of the unmanned underwater vehicle in the attached coordinate system is represented, and u, v and r respectively represent surging, swaying and yawing velocities.
Figure BDA0003952362400000122
Representing an inertial matrix containing additional masses, C (v) representing a diagonally symmetric matrix,
Figure BDA0003952362400000123
represents a damping matrix, τ = [ τ ] uvr ] T Representing moment control input vector, τ u 、τ v 、τ r Respectively representing the control forces of surge, sway and yaw, tau d =[τ dudvdr ] T System disturbance torque vector, τ du 、τ dv 、τ dr Representing the disturbance moments in the surging, swaying, and yawing directions, respectively.
The step S2 specifically includes:
firstly, designing a speed vector observer based on a position signal, eliminating a C (v) v term with a speed vector v in a model because the speed vector v cannot be measured, and introducing a variable matrix P (eta, v):
P(η,v)=H(η)v (2)
in the formula, H (eta) 3×3 Is a mapping matrix of velocity vectors v to P (η, v).
The derivation of equation (2) yields:
Figure BDA0003952362400000131
to eliminate the term with velocity vector v in equation (3), the variable Y is introduced:
Figure BDA0003952362400000132
to make Y =0 and thereby eliminate the term containing the velocity vector v in equation (3), let:
Figure BDA0003952362400000133
thus, there are:
P=Y-H(η)M -1 DH -1 (η)P+H(η)M -1 (τ+τ d ) (6)
and then designing a speed observer as follows:
Figure BDA0003952362400000134
in the formula (I), the compound is shown in the specification,
Figure BDA0003952362400000135
is the observed value of eta, P of the observer, K 1 、K 2 A matrix is designed for the gain of the signal,
Figure BDA0003952362400000136
is the observation error of the observer.
To prove the effectiveness of the designed speed vector observer, the observation error of the observer is defined as:
Figure BDA0003952362400000137
in the formula (I), the compound is shown in the specification,
Figure BDA0003952362400000138
in order to observe the error of the observer,
Figure BDA0003952362400000139
is the position error of the observer observations.
Defining the Lyapunov function:
Figure BDA0003952362400000141
in the formula, x 1 、x 2 The matrix is designed for positive definite. The derivation is carried out on the formula (9) and is obtained from the formula (8):
Figure BDA0003952362400000142
and let:
x 1 J(ψ)H -1 (η)=x 2 K 2 (11)
then there are:
Figure BDA0003952362400000143
as defined by Lyapunov's theorem, the observer shown in equation (7) converges asymptotically.
On the basis of designing a speed observer, designing a disturbance observer for an unknown disturbance part, and carrying out online estimation on disturbance in the surge and heading directions according to the formula (6):
Figure BDA0003952362400000144
in the formula (I), the compound is shown in the specification,
Figure BDA0003952362400000145
for a disturbance observer, carrying out unknown disturbance tau on unmodeled parts and inside and outside of the system d Observed value of, K 0 ∈R 3×3 Is a design matrix composed of positive definite parameters.
In the actual situation where the temperature is too high,
Figure BDA0003952362400000146
failure to pass the observed value
Figure BDA0003952362400000147
Directly obtained, thus introducing auxiliary variables
Figure BDA0003952362400000148
And then designing a disturbance observer as follows:
Figure BDA0003952362400000149
in the formula, unknown disturbance tau is generated for the disturbance observer to the unmodeled part and the inside and outside of the system d Observed value of, K 0 ∈R 3×3 Is a design matrix composed of positive definite parameters.
To prove the effectiveness of the designed moment disturbance observer, a Lyapunov function is defined:
Figure BDA00039523624000001410
in the formula (I), the compound is shown in the specification,
Figure BDA00039523624000001411
τ d is the true value of the disturbance of the system,
Figure BDA00039523624000001412
an error is estimated for the parameter data.
The derivation of equation (15) can be:
Figure BDA0003952362400000151
from this, the observation error of the disturbance observer expressed by the equation (14) can be known
Figure BDA0003952362400000152
Converging progressively to zero in a finite time.
The step S3 specifically includes:
and designing a virtual control law of longitudinal speed and transverse speed based on the real-time position error. First, the error between the actual position and the desired trajectory position is defined as:
Figure BDA0003952362400000153
in the formula, x d As the abscissa, y, of the desired trajectory d X is the abscissa of the actual tracking position, and y is the ordinate of the actual tracking position.
And (3) solving a first derivative of the formula (17), and combining the unmanned submersible vehicle kinematics and the dynamic model formula (1) to obtain:
Figure BDA0003952362400000154
in order to make the position tracking error converge to zero, a vertical virtual control law a is designed ud And a horizontal virtual rate a vd Comprises the following steps:
Figure BDA0003952362400000155
in the formula (I), the compound is shown in the specification,
Figure BDA0003952362400000156
whereinK > 0,c > 0 is a design constant.
In order to prevent the output of virtual control from exceeding the maximum value of the moment provided by the propeller of the unmanned underwater vehicle and solve the problem that the high-order derivation of the virtual control law is required in the process of solving the actual control law so as to increase the calculation complexity, an instruction filter is used for limiting the speed and the amplitude of the virtual control law, and a longitudinal virtual control law a is required to be limited ud And a horizontal virtual rate a vd As input to the command filter, the output of the command filter
Figure BDA0003952362400000157
And
Figure BDA0003952362400000158
and the method is used for solving the actual control law.
Definition of instruction filter:
in order to solve the problem of complexity expansion of differentiation which is possibly generated by derivation of the virtual control law when the actual control law is obtained, reduce the workload of reverse calculation and prevent the problem of control law failure caused by input saturation, the amplitude and the speed of the virtual control law are limited; the instruction filter can also handle constraints of intermediate state variables, as compared to a first order filter. Therefore, the command filter is introduced to process the virtual control law.
The state space expression for the instruction filter is defined as:
Figure BDA0003952362400000161
in the formula, x id Is the input of the filter, z 1 For the output of this filter, ζ and ω are the damping and bandwidth, respectively, of the instruction filter. When the effect of amplitude and velocity is not taken into account, x id And z 1 The linear relationship between them is:
Figure BDA0003952362400000162
when inputting x id When bounded, then z is output 1
Figure BDA0003952362400000163
Is also bounded, and in the case where ω is large, | z 1 -x id I tends to zero.
To prove the effectiveness of the virtual control law, a Lyapunov function is defined:
Figure BDA0003952362400000164
obtained by the formula (19):
Figure BDA0003952362400000165
substituting equation (22) into equation (18) has:
Figure BDA0003952362400000166
when (u-a) ud ) And (v-a) vd ) When the ratio is close to 0:
Figure BDA0003952362400000167
the derivation of equation (22) in combination with equation (25) yields:
Figure BDA0003952362400000168
from this, it can be seen that the virtual control law shown in equation (19) can control the position error x e 、y e Asymptotically converging to zero. The virtual control law may also be considered a speed expectation for the unmanned vehicle to reach the target position.
The step S4 specifically includes:
design longitudinal thrust controlLaw and lateral thrust control law; firstly, designing a longitudinal thrust control law, and setting an error u between a longitudinal speed and an expected longitudinal speed e Comprises the following steps:
Figure BDA0003952362400000171
in the formula (I), the compound is shown in the specification,
Figure BDA0003952362400000172
for the observation of the longitudinal velocity u by the velocity observer, a ud Is a virtual control law of longitudinal speed.
Designing integral slip-form surface s for longitudinal speed error 1
Figure BDA0003952362400000173
In the formula, mu 1 > 0 is a constant.
In order to further obtain the longitudinal thrust control law, the derivation is carried out on the formula (28) and the formula (1) is combined to obtain:
Figure BDA0003952362400000174
in the formula (I), the compound is shown in the specification,
Figure BDA0003952362400000175
is an estimate of the observer's velocity vector u, v, r, τ u Is the actual control law of the longitudinal thrust,
Figure BDA0003952362400000176
for the virtual control law a ud Output value, m, through instruction filter 11 、m 22 、d 11 The coefficients in the inertia matrix M and the damping matrix D of the additional mass in equation (1) are shown, respectively.
Further, the longitudinal thrust control law τ is obtained from the equations (29) and (14) u Comprises the following steps:
Figure BDA0003952362400000177
in the formula, λ 1 >0;
Figure BDA0003952362400000178
Is the observed compensation value of the interference observer to the longitudinal disturbance in the formula (7).
Designed to demonstrate longitudinal thrust control law τ u Defining the Lyapunov function:
Figure BDA0003952362400000179
the derivation of equation (31) and the substitution of equations (29) and (30) yields:
Figure BDA00039523624000001710
since when t → ∞ is reached,
Figure BDA00039523624000001711
therefore, the following are provided:
Figure BDA0003952362400000181
thus signal s 1 The agreement is eventually bounded.
Secondly, designing a transverse thrust control law, and setting the error between the transverse speed and the expected transverse speed as follows:
Figure BDA0003952362400000182
because the power system of the unmanned underwater vehicle is an under-actuated system and has no thrust output in the transverse direction, the steering torque tau is required r The following linear sliding mode surfaces are designed:
Figure BDA0003952362400000183
in the formula, mu 2 0 is a constant.
To further determine the lateral thrust control law, equation (35) is derived over time:
Figure BDA0003952362400000184
in the formula (I), the compound is shown in the specification,
Figure BDA0003952362400000185
for virtual control law a vd A first derivative output value that passes through the instruction filter,
Figure BDA0003952362400000186
for virtual control law a vd The second derivative output value of the instruction filter. Also, in combination with formula (1):
Figure BDA0003952362400000187
in the formula, τ r Is the actual control law of the transverse direction moment. Combining formula (37) with formula (1), formula (36) can be expressed as:
Figure BDA0003952362400000188
from formula (38) in combination with formula (14):
Figure BDA0003952362400000189
in the formula of lambda 2 >0;
Figure BDA00039523624000001810
Compensating values for the observation of lateral disturbances by a disturbance observer。
To demonstrate the designed transverse direction moment control law tau r Defining the Lyapunov function:
Figure BDA0003952362400000191
the derivation of equation (40) and the substitution of equations (38) and (39) yields:
Figure BDA0003952362400000192
since when t → ∞ is satisfied,
Figure BDA0003952362400000193
therefore, there are:
Figure BDA0003952362400000194
thus signal s 2 The agreement is eventually bounded.
In the step S5, simulation verification is performed on the trajectory tracking control in the unmanned underwater vehicle deployment and recovery process, and the validity of the verification method specifically includes:
the reference trajectory is defined as follows:
Figure BDA0003952362400000195
setting longitudinal and transverse unmodeled parts and external disturbance tau in the system du And τ dr Comprises the following steps:
Figure BDA0003952362400000196
setting the kinetic equation M, C and D matrix of the unmanned underwater vehicle as follows:
Figure BDA0003952362400000197
simulation related parameter mu 1 =1,μ 2 =0.02,λ 1 =1,λ 2 =1,ζ 1 =0.9,ζ 2 =0.9,ω 1 =20rad/s,ω 2 =20rad/s, design matrix K 0 、K 1 、K 2 Setting as follows:
K 0 =[2,0,2] T ,
Figure BDA0003952362400000198
unmanned underwater vehicle initial position [ x (0), y (0), psi (0)] T =[0m,0m,0.5rad] T Initial state [ u (0), v (0), r (0)] T =[3m/s,0m/s,0rad/s] T
The result of the expected track tracking path and the actual track tracking path pair shown in fig. 3 shows that the control method provided by the invention can well enable the unmanned underwater vehicle to run on the expected path, and keep a stable state all the time, and has good performance.
The unmanned underwater vehicle position tracking variation curves are shown in fig. 4-5, and at 83s, the longitudinal position tracks the upper expected path, and at 189s, the transverse position tracks the upper expected path.
The observed values of the observer on the states of the unmanned underwater vehicle are shown in fig. 6-8, which shows that the velocity vector can be observed quickly and accurately.
The outputs of the horizontal and vertical virtual control laws are shown in fig. 9-10, and the outputs of the vertical and horizontal actual control laws are shown in fig. 11-12, so that it can be seen that the control law curve is smooth and has good robustness.
As can be seen from the variation curves of the vertical and horizontal actual control laws shown in fig. 9-10, there is a strong chattering phenomenon in the interval of 80s-120s, and it can be seen from the combination of equations (30) and (39) and fig. 9-10 that the chattering phenomenon is caused by the strong chattering in the interval of 80s-120s in the virtual control law. The derivation process of the virtual control law combining the formula (17) and the formula (19) can be seen that the essence of the virtual control law is that the unmanned underwater vehicle essentially represents the slave unmanned underwater vehicleThe expected values of the longitudinal speed and the lateral speed in the process from the current actual position to the desired position, and therefore the virtual control law is highly correlated with the current longitudinal and lateral actual positions and the desired longitudinal and lateral positions, in other words, the essential cause of the jitter of the virtual control law is the error value between the actual position and the desired track position shown in equation (17). It can be seen from the relationship between the longitudinal position and the transverse position shown in fig. 4-5 that at the 80 th s, the actual longitudinal position coincides with the desired longitudinal position, but the system is not completely stable at this time, because the actual transverse position does not coincide with the desired transverse position, during which the transverse actual control law is always active, and the magnitude of the transverse velocity is continuously adjusted, but since the unmanned vehicle is an under-actuated system, the longitudinal velocity is affected by the transverse velocity component, and then the longitudinal thrust control law must be affected by the component from the transverse velocity, and the longitudinal thrust control law is designed based on the integral surface of the velocity difference, so that the longitudinal thrust control law switches back and forth near the sliding surface to keep stable, resulting in the buffeting of the output of the longitudinal control law shown in fig. 11. As can be seen in connection with fig. 4 and 6, the longitudinal thrust control law, by switching back and forth near the sliding mode face, ensures that the longitudinal actual position tracks the upper longitudinal desired position exactly from 83 seconds and is not deflected by the influence of the transverse thrust. Similarly, for the lateral thrust control law, before the 80s moment, because the longitudinal speed is not constrained, the lateral thrust control law can freely control and adjust the lateral speed of the unmanned underwater vehicle, so that the lateral actual position of the unmanned underwater vehicle tracks the expected actual position, but after the 80s moment, because the longitudinal speed is constrained by the longitudinal thrust control law and the unmanned underwater vehicle is an under-actuated system, the change of the lateral speed is necessarily influenced by the longitudinal speed component, which partially influences the interference in the switching process of the longitudinal control law near the sliding mode surface, because after the 80s moment, the longitudinal thrust control law ensures the longitudinal tracking effect in the form of buffeting, which partially influences the lateral speed in the form of buffeting, and in the 80s-120s process, the lateral thrust control law is the same as in the process of confirming that the actual lateral position reaches the expected lateral positionWhen the actual transverse position is tracked in the process of the expected transverse position, the output of the transverse thrust control law is gradually reduced, the influence on the longitudinal speed is also reduced, and therefore the buffeting of the transverse thrust control law and the buffeting of the longitudinal thrust control law are synchronously gradually reduced until y in the formula (17) e =0, v in formula (34) e =0, and further s in formula (35) 2 And =0, the buffeting in the transverse control law when the whole system reaches a stable state completely disappears.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.

Claims (1)

1. An unmanned underwater vehicle track tracking control method considering speed sensor failure is characterized by comprising the following steps:
s1, establishing a water surface navigation mathematical model of the unmanned underwater vehicle, and setting an expected track mathematical model;
s2, designing a speed vector observer based on the position signal, designing a disturbance observer in the longitudinal moment and transverse moment directions on the basis, and verifying the validity of the disturbance observer through Lyapunov theorem;
s3, designing a virtual control law of longitudinal speed and transverse speed based on the real-time position error, introducing an instruction filter to constrain the amplitude and speed of the virtual control law on the basis of the virtual control law, solving the problem that the high-order derivation needs to be carried out on the virtual control law when the actual control law is solved to increase the computational complexity, and verifying the effectiveness of the virtual control law through Lyapunov theorem;
s4, designing a track tracking sliding mode controller of the unmanned underwater vehicle based on the virtual control law and the output value of the instruction filter, compensating disturbance of longitudinal moment and transverse moment based on a system disturbance estimation value observed by a disturbance observer, and verifying effectiveness of the track tracking sliding mode controller through Lyapunov theorem;
s5, verifying the effectiveness of the proposed method through simulation;
in S1, establishing a kinematics and dynamics model comprises the following steps:
Figure FDA0003952362390000011
in the formula, eta = [ x, y, psi =] T X, y represents the unmanned vehicle motion position vector, psi ∈ [0,2 pi]Representing the heading angle, J (ψ) is a conversion matrix of the northeast coordinate system and the satellite coordinate system, v = [ u, v, r] T Representing the motion velocity vector of the unmanned underwater vehicle in an attached coordinate system, wherein u, v and r respectively represent surging, swaying and yawing velocities;
Figure FDA0003952362390000012
representing an inertial matrix containing additional masses, C (v) representing a diagonally symmetric matrix,
Figure FDA0003952362390000013
denotes a damping matrix, τ = [ ] uvr ] T Representing moment control input vector, τ u 、τ v 、τ r Respectively representing the control forces of surge, sway and yaw, tau d =[τ dudvdr ] T System disturbance torque vector, τ du 、τ dv 、τ dr Respectively representing disturbance moments in the surge, sway and bow directions;
in S2, a velocity vector observer based on a position signal is designed first, and since the velocity vector v is not measurable, the C (v) v term with the velocity vector v in the model is eliminated, and a variable matrix P (η, v) is introduced:
P(η,v)=H(η)v (2)
in the formula, H (eta) 3×3 Is a mapping matrix of velocity vectors v to P (eta, v);
the derivation of equation (2) yields:
Figure FDA0003952362390000021
to eliminate the term with velocity vector v in equation (3), the variable Y is introduced:
Figure FDA0003952362390000022
thus, there are:
P=Y-H(η)M -1 DH -1 (η)P+H(η)M -1 (τ+τ d ) (5)
and then designing a speed observer as follows:
Figure FDA0003952362390000023
in the formula (I), the compound is shown in the specification,
Figure FDA0003952362390000024
is the observed value of eta, P of the observer, K 1 、K 2 A matrix is designed for the gain of the signal,
Figure FDA0003952362390000025
is the observation error of the observer;
on the basis of designing a speed observer, designing a disturbance observer for an unknown disturbance part, and carrying out online estimation on disturbance in the surge and heading directions according to the formula (5):
Figure FDA0003952362390000026
in the formula (I), the compound is shown in the specification,
Figure FDA0003952362390000027
for a disturbance observer, carrying out unknown disturbance tau on unmodeled parts and inside and outside of the system d Observed value of, K 0 ∈R 3×3 A design matrix consisting of positive definite parameters;
in the actual situation where the temperature is too high,
Figure FDA0003952362390000028
failure to pass the observed value
Figure FDA0003952362390000029
Directly obtained, thus introducing auxiliary variables
Figure FDA00039523623900000210
And then designing a disturbance observer as follows:
Figure FDA00039523623900000211
in the formula (I), the compound is shown in the specification,
Figure FDA00039523623900000212
for a disturbance observer, carrying out unknown disturbance tau on unmodeled parts and inside and outside of the system d Observed value of, K 0 ∈R 3×3 A design matrix consisting of positive definite parameters;
in the S3, a virtual control law of longitudinal speed and transverse speed is designed based on the real-time position error; first, the error between the actual position and the desired trajectory position is defined as:
Figure FDA0003952362390000031
in the formula, x d As the abscissa, y, of the desired trajectory d Is the ordinate of the expected track, x is the abscissa of the actual tracking position, and y is the ordinate of the actual tracking position;
and (3) obtaining a derivative of the formula (9) by combining the unmanned submersible vehicle kinematics and the dynamic model formula (1):
Figure FDA0003952362390000032
in order to make the position tracking error converge to zero, a vertical virtual control law a is designed ud And a horizontal virtual rate a vd Comprises the following steps:
Figure FDA0003952362390000033
in the formula (I), the compound is shown in the specification,
Figure FDA0003952362390000034
wherein k is more than 0,c is more than 0 and is a design constant;
in order to prevent the output of virtual control from exceeding the maximum value of the moment provided by the propeller of the unmanned underwater vehicle and solve the problem that the high-order derivation of the virtual control law is required in the process of solving the actual control law so as to increase the calculation complexity, an instruction filter is used for limiting the speed and the amplitude of the virtual control law, and a longitudinal virtual control law a is required to be limited ud And a horizontal virtual rate a vd As input to the command filter, the output of the command filter
Figure FDA0003952362390000035
And
Figure FDA0003952362390000036
the method is used for solving the actual control law;
in S4, a longitudinal thrust control law and a transverse thrust control law are designed; firstly, designing a longitudinal thrust control law, and setting an error u between a longitudinal speed and an expected longitudinal speed e Comprises the following steps:
Figure FDA0003952362390000037
in the formula (I), the compound is shown in the specification,
Figure FDA0003952362390000038
for the observation of the longitudinal velocity u by the velocity observer, a ud A virtual control law for longitudinal velocity;
designing integral slip-form surface s for longitudinal speed error 1
Figure FDA0003952362390000039
In the formula, mu 1 0 is a constant;
in order to further obtain the longitudinal thrust control law, the derivation is carried out on the formula (13) and the formula (1) is combined to obtain:
Figure FDA0003952362390000041
in the formula (I), the compound is shown in the specification,
Figure FDA0003952362390000042
is an estimate of the observer's velocity vector u, v, r, τ u Is the actual control law of the longitudinal thrust,
Figure FDA0003952362390000043
for virtual control law a ud Output value, m, through instruction filter 11 、m 22 、d 11 Coefficients in an inertia matrix M and a damping matrix D of the additional mass in the formula (1) respectively;
further, the longitudinal thrust control law τ is obtained from the equations (14) and (8) u Comprises the following steps:
Figure FDA0003952362390000044
in the formula, λ 1 >0;
Figure FDA0003952362390000045
Is an observed compensation value of the interference observer in the formula (8) to longitudinal disturbance;
secondly, designing a transverse thrust control law, and setting an error v between the transverse speed and the expected transverse speed e Comprises the following steps:
Figure FDA0003952362390000046
in the formula (I), the compound is shown in the specification,
Figure FDA0003952362390000047
for the observation of the lateral velocity v by the velocity observer, a vd A virtual control law for lateral velocity;
because the power system of the unmanned underwater vehicle is an under-actuated system and has no thrust output in the transverse direction, the steering torque tau is required r Design the following linear sliding form surface s 2
Figure FDA0003952362390000048
In the formula, mu 2 > 0 is a constant;
to further determine the lateral thrust control law, equation (17) is derived with respect to time:
Figure FDA0003952362390000049
in the formula (I), the compound is shown in the specification,
Figure FDA00039523623900000410
for virtual control law a vd A first derivative output value that passes through the instruction filter,
Figure FDA00039523623900000411
for virtual control law a vd Passing through a second derivative output value of the instruction filter; also, in combination with formula (1):
Figure FDA0003952362390000051
in the formula (I), the compound is shown in the specification,
Figure FDA0003952362390000052
is an estimate of the observer's velocity u, v, r, τ r Is the actual control law of the moment in the transverse direction, m 11 、m 22 、m 33 、d 22 、d 33 Coefficients in an inertia matrix M and a damping matrix D of the additional mass in the formula (1) respectively; combining formula (19) with formula (1), formula (18) can be represented as:
Figure FDA0003952362390000053
law of transverse thrust control τ obtained by combining equation (20) with equation (8) r Comprises the following steps:
Figure FDA0003952362390000054
in the formula, λ 2 >0;
Figure FDA0003952362390000055
Is the observed compensation value of the interference observer to the lateral disturbance in the formula (8).
CN202211453151.3A 2022-11-21 2022-11-21 Unmanned submarine trajectory tracking control method considering failure of speed sensor Active CN115686034B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211453151.3A CN115686034B (en) 2022-11-21 2022-11-21 Unmanned submarine trajectory tracking control method considering failure of speed sensor

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211453151.3A CN115686034B (en) 2022-11-21 2022-11-21 Unmanned submarine trajectory tracking control method considering failure of speed sensor

Publications (2)

Publication Number Publication Date
CN115686034A true CN115686034A (en) 2023-02-03
CN115686034B CN115686034B (en) 2023-07-25

Family

ID=85054675

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211453151.3A Active CN115686034B (en) 2022-11-21 2022-11-21 Unmanned submarine trajectory tracking control method considering failure of speed sensor

Country Status (1)

Country Link
CN (1) CN115686034B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117234092A (en) * 2023-11-16 2023-12-15 广州大学 Underwater submerged robot and control method and control device thereof

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109933074A (en) * 2019-04-10 2019-06-25 大连海事大学 A kind of more unmanned boat flocking control device structures and design method having leader
CN110377036A (en) * 2019-07-09 2019-10-25 哈尔滨工程大学 A kind of unmanned water surface ship Track In Track set time control method constrained based on instruction
CN111007854A (en) * 2019-12-18 2020-04-14 哈尔滨工程大学 Under-actuated ship trajectory tracking control system
CN111736617A (en) * 2020-06-09 2020-10-02 哈尔滨工程大学 Speed observer-based benthonic underwater robot preset performance track tracking control method
CN112327883A (en) * 2020-11-27 2021-02-05 哈尔滨工程大学 Under-actuated surface vessel path tracking control method based on extended state observer

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109933074A (en) * 2019-04-10 2019-06-25 大连海事大学 A kind of more unmanned boat flocking control device structures and design method having leader
CN110377036A (en) * 2019-07-09 2019-10-25 哈尔滨工程大学 A kind of unmanned water surface ship Track In Track set time control method constrained based on instruction
CN111007854A (en) * 2019-12-18 2020-04-14 哈尔滨工程大学 Under-actuated ship trajectory tracking control system
CN111736617A (en) * 2020-06-09 2020-10-02 哈尔滨工程大学 Speed observer-based benthonic underwater robot preset performance track tracking control method
CN112327883A (en) * 2020-11-27 2021-02-05 哈尔滨工程大学 Under-actuated surface vessel path tracking control method based on extended state observer

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
DONGHAO XU 等: "Trajectory Tracking of Underactuated Unmanned Surface Vessels: Non-Singular Terminal Sliding Control with Nonlinear Disturbance Observer", vol. 12, no. 6, pages 1 - 17 *
LANPING ZOU 等: "Robust Neural Network Trajectory-Tracking Control of Underactuated Surface Vehicles Considering Uncertainties and Unmeasurable Velocities", vol. 9, pages 29 - 38 *
李宗宣 等: "结合速度和干扰观测的船舶路径跟踪模型预测控制", 上海海事大学学报, vol. 42, no. 1, pages 19 - 25 *
汪慢: "高速欠驱动UUV的三维空间轨迹跟踪控制方法", pages 5 *
高振宇: "自主水下航行器的轨迹跟踪及编队控制", pages 4 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117234092A (en) * 2023-11-16 2023-12-15 广州大学 Underwater submerged robot and control method and control device thereof
CN117234092B (en) * 2023-11-16 2024-01-26 广州大学 Underwater submerged robot and control method and control device thereof

Also Published As

Publication number Publication date
CN115686034B (en) 2023-07-25

Similar Documents

Publication Publication Date Title
CN108803632B (en) Water surface unmanned ship all-state constraint trajectory tracking control method based on saturation compensation technology
CN107168312B (en) Space trajectory tracking control method for compensating UUV kinematic and dynamic interference
CN107024863B (en) UUV trajectory tracking control method for avoiding differential explosion
CN109814392B (en) Fault-tolerant control method for fault robustness of under-actuated underwater robot actuator
CN111736600B (en) Track tracking control method for unmanned surface vehicle under time-lag asymmetric time-varying all-state constraint
Wang et al. Command filter based adaptive neural trajectory tracking control of an underactuated underwater vehicle in three-dimensional space
Do Global robust and adaptive output feedback dynamic positioning of surface ships
CN111580523A (en) Unmanned ship path tracking active disturbance rejection control method based on sideslip angle compensation
CN110618611A (en) Unmanned ship trajectory tracking safety control method based on slew rate constraint
CN104898688A (en) UUV four degree-of-freedom dynamic positioning adaptive anti-interference sliding mode control system and control method
Wang et al. Filter-backstepping based neural adaptive formation control of leader-following multiple AUVs in three dimensional space
CN109739248B (en) Ship-borne three-degree-of-freedom parallel stable platform stability control method based on ADRC
CN108519736A (en) A kind of USV track followings contragradience sliding-mode control
CN113821030B (en) Fixed time track tracking control method for under-actuated unmanned ship
CN113467231A (en) Unmanned ship path tracking method based on sideslip compensation ILOS guidance law
CN115686034A (en) Unmanned underwater vehicle track tracking control method considering failure of speed sensor
CN114967714A (en) Anti-interference motion control method and system for autonomous underwater robot
Dai et al. Dual closed loop AUV trajectory tracking control based on finite time and state observer
CN115079698A (en) Unmanned surface vessel path tracking time trigger control method under deception attack task
KR101690154B1 (en) Waypoint tracking apparatus of autonomous underwater system and method thereof
Shikai et al. Trajectory tracking for underactuated UUV using terminal sliding mode control
Caharija et al. Path following control of underactuated surface vessels in the presence of multiple disturbances
Liu et al. Output feedback control design for station keeping of AUVs under shallow water wave disturbances
CN116820081A (en) Under-actuated unmanned ship path tracking control method based on model-free self-adaptive sliding mode
CN114527772A (en) Design method and system of AUV trajectory tracking controller

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
TR01 Transfer of patent right

Effective date of registration: 20240517

Address after: 230000 Woye Garden Commercial Building B-1017, 81 Ganquan Road, Shushan District, Hefei City, Anhui Province

Patentee after: HEFEI JINGLONG ENVIRONMENTAL PROTECTION TECHNOLOGY Co.,Ltd.

Country or region after: China

Address before: 150040 No. 52, Xuefu Road, Nangang District, Heilongjiang, Harbin

Patentee before: HARBIN University OF SCIENCE AND TECHNOLOGY

Country or region before: China