Next Article in Journal
Rapid Synthesis of Silver Nanowires in the Polyol Process with Conventional and Microwave Heating
Next Article in Special Issue
Causality Analysis with Different Probabilistic Distributions Using Transfer Entropy
Previous Article in Journal
Lyapunov-Guided Energy Scheduling and Computation Offloading for Solar-Powered WSN
Previous Article in Special Issue
Vibration Fatigue Life Reliability Cable Trough Assessment by Using Weibull Distribution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Controller for an Asymmetric Underactuated Hovercraft in Terms of Quasi-Velocities

by
Przemyslaw Herman
Institute of Automatic Control and Robotics, Poznan University of Technology, ul. Piotrowo 3a, 60-965 Poznan, Poland
Appl. Sci. 2023, 13(8), 4965; https://doi.org/10.3390/app13084965
Submission received: 13 March 2023 / Revised: 5 April 2023 / Accepted: 13 April 2023 / Published: 14 April 2023
(This article belongs to the Special Issue Recent Advances in Nonlinear Vibration and Control)

Abstract

:
In this paper, a nonlinear controller for tracking a desired trajectory for an underactuated hovercraft is considered. It is a modification of a method known from the literature. However, the control algorithm considered here has two important features that differ from the mentioned control strategy. First, it concerns the case when the center of mass does not coincide with the geometric center, which results in additional forces and moments of force. The lack of symmetry causes the original trajectory tracking method not to take this fact into account, while the proposed approach is a generalization of the known concept. Here, a diagonalization of the inertia matrix has been applied, by means of a velocity transformation, which made it possible to reduce the symmetric matrix to a diagonal form. Secondly, the transformed quasi-velocity equations of motion allow some insight into the dynamics of the vehicle as it moves, which was not shown in the source work. The offered approach was verified by numerical tests for a hovercraft model with three DOF and for two desired trajectories. The method can be useful in preliminary simulation studies at the controller selection stage without experimental validation.

1. Introduction

Underactuated vehicles, i.e., those that do not have all the control signals, have enjoyed growing interest for many years. This is due to the fact that their cost can be lower than fully actuated vehicles and their design can be simplified. To this group of vehicles belong hovercrafts. The control methods used can be divided from various points of view. They can be distinguished, for example, on the basis of the mathematical description, design or control task being carried out. Mostly considered is a model with three DOF for a vehicle moving horizontally. One of the criteria adopted here is the implementation of the controller into either the vehicle earth-fixed (inertial) frame or body-fixed body frame. Of the works concerning the first group, the following can be mentioned: [1] (nonlinear receding horizon control with real-time optimization), [2] (stabilization, fuzzy control), [3,4] (tracking control with two loops), [5] (LQR trajectory tracking), [6,7] (sliding mode control) and [8] (flatness-based tracking control). Trajectory tracking and path-following control schemes are discussed in the following: [9,10] (Lyapunov-based approach with backstepping) and [11] (Lyapunov-based approach with adaptive switching supervisory control). Despite the universality of the schemes, they have been applied to the hovercraft described in the inertial frame.
However, there are also many works belonging to the second group. An article on the stabilization of the hovercraft in [12] can be pointed out here. A modified genetic algorithm together with a fuzzy proportional integral derivative controller was presented in [13]. Trajectory tracking schemes based on a differential flatness property were proposed in [14,15,16]. Other methods used a flatness property with generalized proportional integral (GPI) control [17], backstepping [18], active disturbance rejection control, differential flatness or reinforcement learning [19]. Modified by disturbances, the model was used for tracking control in [20] (higher-order sliding mode control). A Lyapunov-based backstepping technique tracking algorithm was proposed in [9]. A strategy composed of the line-of-sight guidance law with a neural network-based adaptive dynamic inversion for the underactuated hovercraft following a prescribed path was found in [21]. The solution to the tracking control problem with simultaneous parameter identification was given in [22]. In [23], trajectory tracking and path following algorithms were developed for a hovercraft vehicle in the presence of uncertain parameters and unknown external disturbances, using the backstepping technique. An article [24] discussed and verified the design and trajectory tracking controller for an underactuated hovercraft using an observer, with unmeasured linear velocity and subject to time-varying disturbances. Another approach to trajectory tracking, namely based on the state transformation technique and cascade method, combined with a transverse function was proposed in [25]. In [26], it was shown that the quantitative feedback theory (QFT) technique provided an approach as well as a feasible control design methodology to solve the tracking problem of a nonlinear multivariable underactuated hovercraft, by using a manageable approach to nonlinear QFT synthesis based on the local linearization of the nonlinear plant about appropriate closed-loop outputs. It is also possible to adapt known control methods related to marine vehicles in horizontal motion. However, this requires careful verification because under some operating conditions such algorithms are unreliable, as was examined in [27].
Control methods can also be divided according to the design of the vehicle, that is, when only motors or motors and a rudder are used. If the vehicle has only thrusters, then schemes, e.g., from [4,5,6,9,17,19,20,25], can be proposed. On the other hand, when the hovercraft is equipped with a thruster and rudder, the relevant methods are as presented in [8,22,23,24]. The design of the vehicle is not a barrier to the applicability of methods with thrusters when equipped with a rudder because there are formulas for converting forces and moments of force as given in [23], recalling work from [28].
In this work, the trajectory tracking control scheme in the body-fixed frame for underactuated hovercrafts is considered. The proposed algorithm is a modification of the known controller developed in [22]. This modification extends the original scheme to non-symmetric vehicles. In the classical dynamics model, the vehicle is assumed to be symmetric which means that any inaccuracies in the model must be reduced using a controller. Furthermore, mutual interference between variables is also neglected in such a model. With this approach, we lose insight into the dynamics of the system, but we can effectively accomplish the task of trajectory tracking. The proposed algorithm addresses the trajectory tracking desired for a non-symmetric vehicle model. Such a model is closer to the real vehicle because it takes into account the fact of inaccuracy of the position of the center of mass as well as the external interactions on the individual variables. This is the case when the vehicle is not fully balanced or is carrying additional mass, which shifts the center of symmetry. The dynamics model is based on a velocity transformation which leads to a diagonalized inertia matrix. As a consequence, the new inertia matrix has a diagonal form and the unaccounted quantities are shifted to the other components of the dynamics equation. In addition, the control method by such a model can be used for simulation tests, which will help in later experimental studies. It can also be useful at the design phase when the model parameters can still be changed.
The contributions of this paper can be given as follows:
(1)
The proposed algorithm, unlike the scheme in [22] and many other algorithms, is suitable for asymmetric vehicles in the horizontal plane. Thus, it is an extension of the known control concept to a more general class of vehicles.
(2)
The theoretical development concerns the inclusion in the control algorithm of components related to the fact that the geometric center does not coincide with the center of mass. These components are then introduced into the stability proof of the system. In addition, the dynamic model of the vehicle is transformed to velocity space and, as a result of this transformation, new variables include the parameters of this model. Consequently, the combination consisting of a velocity transformation and a control algorithm yields a useful tool for estimating the behavior of a vehicle during trajectory tracking.
(3)
The simulation verification of the proposed trajectory tracking algorithm for estimating the effect of mutual coupling between velocities and evaluating the reduction in kinetic energy by each quasi-velocity (which includes coupling with other variables).
(4)
Pointing out some performance differences between the IQV-based control scheme and two other selected algorithms that use a diagonal inertia matrix model.
The remainder of this paper is organized as follows. In Section 2, the hovercraft model in terms of quasi-velocities is shown. The tracking control algorithm is presented in Section 3. Section 4 contains the simulation results of the confirmed effectiveness of the proposed controller. The conclusions are summarized in Section 5. At the end of the work, there is Appendix A.

2. Hovercraft Model

The model of a hovercraft moving horizontally that is considered in this work is shown in Figure 1.
Using the North-East-Down (NED) frame (the inertial frame), the position and orientation of the vehicle is η = [ x , y , ϕ ] T . Based on [22,29], the equations of motion of the hovercraft moving horizontally can be written as follows (symbols are explained in Appendix A):
η ˙ = R ( ϕ ) ν ,
M ν ˙ = C ( ν ) D ( ν ) + τ ,
where ν = [ u , v , r ] T denotes the vector of velocities in the body frame (containing the surge velocity, the sway velocity and the yaw velocity). If the hovercraft is equipped with thrusters and a rudder, then the τ vector can be written as follows [22]:
τ = [ b T T cos ψ , b T T sin ψ , a b T T sin ψ ] T . The remaining matrices and vectors have the following forms:
M = m 0 m y g 0 m m x g m y g m x g J , R ( ϕ ) = cos ϕ sin ϕ 0 sin ϕ cos ϕ 0 0 0 1 , C ( ν ) = m v r m x g r 2 m u r + m y g r 2 m x g u r m y g v r , D ( ν ) = d u 0 sign u + d u u + d u r r d v 0 sign v + d v v + d v r r d r 0 sign r + d u r u + d v r v + d r r ,
in which m, J are the vehicle mass and inertia, respectively; y g , x g are the mass displacements; d 0 u , d 0 v , d 0 r mean the friction coefficients; and coefficients d u , d v , d r , d u r , d v r represent linear drag.
Comment. In Equations (1)–(3), there are mass components, coefficients of friction linear drag and displacements. In the proposed method, it is not necessary to know the exact displacements (after transformation their influence is small but rather the mass components m x g , m y g are important). Although these are difficult to determine for a real vehicle, there are methods by which such information can be obtained. For example, one can use the parameter identification methods described in [22,30] for a hovercraft (a more detailed description of the general method is also known from [31]).
After velocity transformation, e.g., [32], the new equations instead of (2) are of the form
N ζ ˙ = Υ T C ( ν ) Υ T D ( ν ) + Υ T τ ,
ν = Υ ζ , Υ = 1 0 Υ 13 0 1 Υ 23 0 0 1 , N = d i a g { N 1 , N 2 , N 3 } ,
where ζ = [ ζ 1 , ζ 2 , ζ 3 ] T is the vector of inertial quasi-velocities. Other symbols mean N 1 = m , N 2 = m , N 3 = J m ( x g 2 + y g 2 ) , Υ 13 = y g , Υ 23 = x g .
Equations of motion instead of (4) and (5) can be written as follows:
ζ 1 = u Υ 13 r ,
ζ 2 = v Υ 23 r ,
ζ 3 = r ,
N 1 ζ ˙ 1 = m v r + m x g r 2 d u 0 sign u d u u d u r r + b T T cos ψ ,
N 2 ζ ˙ 2 = m u r m y g r 2 d v 0 sign v d v v d v r r + b T T sin ψ ,
N 3 ζ ˙ 3 = m x g u r + m y g v r d u r u d v r v d r 0 sign r d r r a b T T sin ψ + Υ 13 b T T cos ψ + Υ 23 b T T sin ψ .
Finally, after dividing into two sub-systems, the set of Equations (6)–(11) can be given in the compact vector matrix form:
p ˙ = R 1 ν t r ,
R ˙ 1 = R 1 S ζ 3 ,
m ζ ˙ t r = m S ν t r + f ;
N 3 ζ ˙ 3 = τ r ,
where v t r = [ u v ] T is the linear velocity vector, ζ ˙ t r = [ ζ 1 ζ 2 ] is the corresponding quasi-velocity vector, R 1 = cos ϕ sin ϕ sin ϕ cos ϕ , S = 0 1 1 0 and
f = m x g r 2 d u 0 sign u d u u d u r r + b T T cos ψ m y g r 2 d v 0 sign v d v v d v r r + b T T sin ψ ,
τ r = m x g u r + m y g v r d u r u d v r v d r 0 sign r d r r a b T T sin ψ + Υ 13 b T T cos ψ + Υ 23 b T T sin ψ .
Assumption 1. 
The presence of a discontinuous function s i g n ( a r g ) ( a r g denotes the argument of the function) in the system model makes the whole system a discontinuous system. Thus, in order to ensure the continuity condition of the system model, the condition s i g n ( a r g ) t a n h ( c · a r g ) (where c is a constant value that provides a good enough approximation) is assumed to be acceptable.
Assumption 2. 
Considering Assumption 1, it is also assumed that the functions used for the controller design can also be approximated by t a n h ( c · a r g ) functions. This approach to the approximation of signals of the controller is known from the literature, for example, [33,34]. In the t a n h ( c · a r g ) functions, the constant c has the same value for both the model and the controller.
Comment. The notation of the s i g n function in the following work has been retained because of the symbolism used in the original paper [22]. Here, however, such notation should be understood as the t a n h function.

3. Control Scheme Using Quasi-Velocities

The proposed controller is a modification of the one coming from [22] and is an extension to vehicles with a model that takes into account the asymmetry caused by the shift of the center of mass. Thus, the goal is to also track the desired trajectory p d , which has a fixed point defined in the body frame δ R 2 . The disturbances are not taken into account. This means that the algorithm only considers the friction model as presented in Section 2, while it does not handle other types of external disturbances (deterministic or stochastic).
The tracking error z 1 is introduced as a measure of accuracy. Calculating its time derivative, one has
z 1 = p p d + R 1 δ , z 2 = z ˙ 1 = R 1 v t r p ˙ d + R 1 S δ ζ 3 , δ = [ δ x 0 ] T .
Recalling now (6) and (7), denoting Υ 3 = [ Υ 13 Υ 23 ] T , one obtains ν t r = ζ t r + Υ 3 ζ 3 and its time derivative ν ˙ t r = ζ ˙ t r + Υ 3 ζ ˙ 3 . Thus, the time derivative of z 2 has the following form:
z ˙ 2 = R 1 S v t r ζ 3 + m 1 R 1 ( m S v t r ζ 3 + f + m Υ 3 N 3 1 τ r ) p ¨ d + R 1 S 2 δ ζ 3 2 + R 1 S δ N 3 1 τ r .
The following LFC is proposed:
V = 1 2 k 1 z 1 T z 1 + 1 2 z 2 T z 2 + μ z 1 T z 2 ,
where k 1 > 0 , μ > 0 are some coefficients. Determining analogously as in [22] the time derivative of V, one obtains
V ˙ = A ( z 1 , z 2 ) + ( μ z 1 + z 2 ) T ( z ˙ 2 + k 1 z 1 + k 2 z 2 ) ,
with
A ( z 1 , z 2 ) = k 1 μ z 1 T z 1 + k 2 μ z 1 T z 2 + ( k 2 μ ) z 2 T z 2 ,
and k 2 > 0 . Moreover, in order to guarantee that A ( z 1 , z 2 ) and V ( z 1 , z 2 ) are positively definite, it is assumed that [22]
μ < k 1 , μ < 4 k 1 k 2 4 k 1 + k 2 2 .
Using z ˙ 2 and grouping the terms, Equation (21) can be written in the form:
V ˙ = A ( z 1 , z 2 ) + ( μ z 1 + z 2 ) T R 1 · ( B T cos ψ sin ψ + m 1 f 0 R 1 T p ¨ d δ x 0 ζ 3 2 + Υ 13 Υ 23 + δ x N 3 1 τ r 0 + R 1 T ( k 1 z 1 + k 2 z 2 ) ) ,
where
B = b T β 11 β 12 β 21 β 22 = b T m 1 + N 3 1 Υ 13 2 N 3 1 Υ 13 ( a Υ 23 ) ( Υ 23 + δ x ) N 3 1 Υ 13 m 1 ( Υ 23 + δ x ) N 3 1 ( a Υ 23 ) ,
f 0 = m x g r 2 d u 0 sign ( u ) d u u d u r r m y g r 2 d v 0 sign ( v ) d v v d v r r ,
τ r 0 = m x g u r + m y g v r d u r u d v r v d r 0 sign ( r ) d r r .
Under such conditions, a stabilizing controller should be designed in order to track the system error.
Theorem 1. 
It is assumed that the equations of motion of a hovercraft expressed in terms of quasi-velocities can be written as (6)–(11), with the external force (16) and (17) taking into consideration the states (18) and the bounded desired trajectory p d and bounded its time derivative. If an invertible matrix B (25) is selected, this means that
β 11 β 22 β 12 β 21 0 ,
and the control thrust T as well as the rudder angle ψ are used in the form
T cos ψ sin ψ = B 1 ( m 1 f 0 R 1 T p ¨ d δ x 0 ζ 3 2 + Υ 13 Υ 23 + δ x N 3 1 τ r 0 + R 1 T ( k 1 z 1 + k 2 z 2 ) ) ,
with k 1 > 0 , k 2 > 0 , and then origin of the dynamical system of the system is globally exponentially stable.
Proof. 
Taking into account (23), the LFC (20) is radially unbounded and positive definite function. Applying the controller (29) into (24), the closed-loop time derivative as well as the LFC lead to
V ˙ = A ( z 1 , z 2 ) .
Thus, from Lyapunov’s direct method (a simple convergence lemma [35] (p. 91)), one obtains
V ˙ ( t ) α V ( t ) ,
V ( t ) V ( 0 ) e α t ,
where α is a real positive number. Hence, it can be concluded that the origin of the error system is globally exponentially stable. □
Zero dynamics question. Because ζ 3 = r , the analysis of the zero dynamics is thus analogous as in [22] but not direct. So, let us consider this issue.
The angular kinetic energy of the hovercraft is given in the form
V r = 1 2 r .
Using (8), (11), the controller (29) and calculating the time derivative of V ˙ r , one obtains
V ˙ r = r N 3 1 ( m x g u r + m y g v r d u r u d v r v d r 0 sign r d r r a b T T sin ψ + Υ 13 b T T cos ψ + Υ 23 b T T sin ψ ) .
Taking into account z 2 (18), δ , and denoting e 2 = [ 0 1 ] T , it can be written as follows:
v = e 2 T R 1 T ( z 2 + p ˙ d ) δ x r .
In addition, the inverse matrix to matrix B is of the form
B 1 = b T 1 β ¯ 11 β ¯ 12 β ¯ 21 β ¯ 22 ,
where β ¯ 11 = β 22 / P , β ¯ 12 = β 12 / P , β ¯ 21 = β 21 / P , β ¯ 22 = β 11 / P ( P = β 11 β 22 β 12 β 21 ).
The sway velocity v (35) is applied in the time derivative V ˙ r . Thus, (34) can be rewritten as follows:
V ˙ r = N 3 1 r d r 0 sign r + ( N 3 1 d r r 2 + N 3 1 r ( m x g u r 2 + m y g v r 2 d u r u d v r v ) N 3 1 r [ 1 1 ] Υ 13 0 0 ( Υ 23 a ) β ¯ 11 β ¯ 12 β ¯ 21 β ¯ 22 ( m 1 m x g r 2 d u 0 sign ( u ) d u u d u r r m y g r 2 d v 0 sign ( v ) d v v d v r r + R 1 T ( k 1 z 1 + k 2 z 2 p ¨ d ) δ x 0 r 2 + Υ 13 Υ 23 + δ x N 3 1 ( m x g u r + m y g v r d u r u d v r v d r 0 sign ( r ) d r r ) ) ) .
Denoting now
Ξ = N 3 1 d r r 2 + N 3 1 r ( m x g u r 2 + m y g v r 2 d u r u d v r v ) N 3 1 r [ 1 1 ] Υ 13 0 0 ( Υ 23 a ) β ¯ 11 β ¯ 12 β ¯ 21 β ¯ 22 ( m 1 m x g r 2 d u 0 sign ( u ) d u u d u r r m y g r 2 d v 0 sign ( v ) d v v d v r r + R 1 T ( k 1 z 1 + k 2 z 2 p ¨ d ) δ x 0 r 2 + Υ 13 Υ 23 + δ x N 3 1 ( m x g u r + m y g v r d u r u d v r v d r 0 sign ( r ) d r r ) ) ,
one obtains
V ˙ r = N 3 1 r d r 0 sign r + Ξ .
The following theorem can now be proposed.
Theorem 2. 
Consider the closed loop consisting of the system described by (12)–(15) together with the trajectory tracking controller (29). For any length of the arm a, there is a constant ρ > 0 such that if the control point meets the condition | | δ x | | < ρ , inequality (28) is true, and additionally Ξ < 0 , then the linear velocity v t r = [ u v ] T and the vehicle angular velocity r are bounded.
Proof. 
Based on Theorem 1, it can be concluded that tracking errors z 1 and z 2 , vehicle parameters in (2) and (3) (the mass, the inertia moment, the displacements and the drag coefficients) and the reference trajectory and its derivatives are bounded. If condition Ξ < 0 is satisfied, then V ˙ r 0 and the boundedness of r is ensured. However, because component (38) has a complex form, it is necessary to pay attention to the design parameters and model parameters to guarantee this condition. Design parameters a , δ x can be chosen small enough to reduce the impact of unwanted components. In addition, the values of the components for vehicle asymmetry m x g , m y g must not be too large (because it affects Υ 13 and Υ 23 ). The values of the drag due to asymmetry d u r , d v r and components m 1 d u and m 1 d v should also be small. To sum up this consideration, it can be said that the influence of components due to asymmetry should be limited. When such conditions are met, it is possible to guarantee that Ξ < 0 . In this case, r will be bounded. Given the definition of tracking error z 2 and the boundedness of z 2 and r, one can deduce that the linear velocities u, v are also bounded. □

4. Simulation Results

Proposed verification method. The simulation study was based on the following assumptions:
(1)
A hovercraft model with parameters similar to those of a vehicle known from the literature was selected;
(2)
Two curvilinear trajectories were selected for tracking;
(3)
Comparative tests were carried out for the proposed control algorithm based on quasi-velocities, the algorithm using the TSMC method and the original algorithm from [22];
(4)
The test results were compared on the basis of the assumed set of criteria and indexes.
Comment. The thrust saturation effect is sometimes taken into account in the control algorithm, such as [36,37]. However, when the thruster saturation effect is not significant, then this problem is ignored [38,39]. Such a situation occurs when the desired trajectories have an uncomplicated mathematical form, and the thrust saturation effect is only a temporary or sporadic phenomenon.

4.1. Hovercraft Model and Test Conditions

The purpose of the numerical simulations was to verify the effectiveness of the control algorithm under near-real conditions for two vehicles with different dynamics. This meant that the operating conditions had to be chosen so as not to exceed the maximum forces and torques that could be obtained from the thrusters. In addition, the velocity values could not exceed those that are achievable by each vehicle.
The parameters of the hovercraft adopted for testing are similar to those contained in [21]. They are collected in Table 1. The simulations were performed using the software from [40] (Matlab/Simulink source codes) but modified for an algorithm containing an asymmetric inertia matrix and for quasi-velocities. In this program, in order to avoid discontinuities in the system, the s i g n function was approximated by t a n h functions as used in the literature, e.g., [34,41].
The vehicle models under consideration have vastly different dynamics, which makes it possible to verify the suitability of the control algorithm for different vehicles.
Moreover, the center of mass is shifted, i.e., x g = 0.1 m, y g = 0.1 m, and the desired trajectory position profiles p d = [ x d , y d ] T were selected as follows:
p d 1 = 0.2 t , 7 sin ( 0.05 t ) T ,
p d 2 = 0.2 t , 0.5 t sin ( 0.02 t ) + 0.3 t cos ( 0.12 t ) + 0.22 t T .
Other parameters applied are as follows: a = 0.5 m, b T = 0.5 m and δ x = 0.005 m. The first is the sine trajectory, whereas the second is a complex trajectory. The initial points for the hovercraft were p 0 = [ 0 , 1.0 ] T . The saturation of the control signal was | T cos ψ | 30 N, | T sin ψ | 30 N. A comparison of the simulation results for two different set trajectories can be useful at the stage of trajectory design to perform the set task. The simulations were performed using Matlab/Simulink, the method ODE 3 (Bogacki–Shampine), the time of motion t = 30 s and the time step Δ t = 0.05 s.

4.2. Performance Indexes and Evaluation Criteria

A performance test and a test of changes in motion dynamics were used to verify the proposed control scheme. The following quantities were applied to estimate the performance:
(1)
The mean integrated absolute error, i.e., M I A = 1 t f t 0 t 0 t f | χ e ( t ) | d t , where χ e = Δ p x , Δ p y , Δ ζ 1 , Δ ζ 2 ;
(2)
The mean integrated absolute control, i.e., M I A C = 1 t f t 0 t 0 t f | P ( t ) | d t , where P = T cos ψ , T , sin ψ ;
(3)
The root mean square of the tracking error, i.e., R M S = 1 t f t 0 t 0 t f e ( t ) 2 d t , with e ( t ) = Δ p x 2 + Δ p y 2 ( Δ p x , Δ p y mean the position error in the reference frame).
The estimation of the effects of the dynamics is carried out by means of the following quantities [42,43]:
(1)
The kinetic energy for the total vehicle and corresponding to each quasi-velocity, i.e., K T = 1 2 ν T M ν = 1 2 ζ T N ζ = 1 2 i = 1 3 N i ζ i 2 = i = 1 3 K i .
(2)
The mean kinetic energy relating to each quasi-velocity and the total kinetic energy: K m = m e a n ( K ) = i = 1 3 m e a n ( K i ) .
(3)
The deformation of each velocity resulting from dynamic coupling, i.e., Δ ζ i = ζ i ν i , where i = 1 , 2 , 3 .

4.3. Simulation Tests Conducted

The parameters of the controllers were determined based on the method given in [27]. According to this procedure, the parameters are grouped according to the function played in the control algorithm. Their tuning and analysis of the results consists of visualizing the position errors of the values of the parameters in the same group (the other parameters have unit values of the values, and if they are too large, they should be reduced so that they are close to zero). If the values of several parameters (or one when it is enough) lead partially to the goal, e.g., the convergence of the position errors is acceptable, then the effect of other enhancements on the effectiveness of the controller can be studied. If the quality of the control improves, tuning is performed. With promising results, we change the values of the next group of parameters and check whether the position errors have improved.

4.3.1. Test of Quasi-Velocity Based Controller

The set of control parameters to ensure acceptable tracking results was (for both kind of desired trajectories)
k 1 = 5.0 , k 2 = 5.0 .
These gains were obtained by tuning the control algorithm, i.e., increasing the values from 1 to such values that guaranteed acceptable trajectory tracking accuracy (taking into account the other operating conditions).
Sine trajectory. As is evident in Figure 2a, the set trajectory is tracking correctly and the position errors (presented in Figure 2b) in the x and y directions are tending to near zero after about 5 s. As can be observed in Figure 2c, the non-steady state for the control signals disappears after about 15 s of movement. Figure 2d demonstrates that at a similar time, this condition disappears for hovercraft velocities. Additional information about the effect of dynamic parameters on movement is obtained from changes in the value of Δ ζ 1 , Δ ζ 2 as seen in Figure 2e. Because the displacement values of the center of mass position are the same in the x and y directions, Δ ζ 1 , Δ ζ 2 are therefore symmetrical. Figure 2f presents the evolution of kinetic energy when the vehicle is in motion. Using this graph, it is possible to determine which variable causes more or less kinetic energy consumption, assuming dynamic parameters and a desired trajectory. This can be used to design the vehicle and for planning its movement.
Complex trajectory. Figure 3 shows the results of the analogous tests conducted for the complex trajectory. It can refer, for example, to a situation in which the vehicle should avoid a fixed obstacle. As can be seen in Figure 3a, the algorithm works properly and the desired trajectory is tracked efficiently. This observation is also confirmed in Figure 3b. However, from Figure 3c, it is apparent that the values of the applied force components do not reach a steady state but are constantly changing and, moreover, are larger than for a sine trajectory. A similar phenomenon can be seen in Figure 3d, which shows the velocities of the vehicle. The changes in the Δ ζ 1 , Δ ζ 2 values are larger than before but also symmetrical (for the same reason as previously), as observed in Figure 3e. The kinetic energy changes, as shown in Figure 3f, have larger values than for the sinus trajectory, and their waveforms are considerably different from the previous ones. So, the selection of the trajectory significantly affects the consumption of this energy.
In order to obtain more objective information, other indexes were used to evaluate the achieved results. The corresponding values of the controller’s performance are summarized in Table 2.
From the values of the MIA indexes, it can be seen that the tracking accuracy is comparable for both selected trajectories. The Δ ζ 1 , Δ ζ 2 indexes have larger values for the complex trajectory, which means that the velocities are more deformed for such a trajectory. The observations regarding the larger values of the applied force components when using a complex trajectory instead of a sinus trajectory were confirmed. However, the RMS value was lower for this trajectory; therefore, it is somewhat more accurately tracked. On the other hand, the kinetic energy consumption was indeed higher for the complex trajectory (the observations from the graphs were confirmed). However, it is also interesting to know which quasi-variable is the reason for this energy consumption. For the sine trajectory, it was calculated that K 1 = 0.4321 , K 2 = 0.1265 , K 3 = 0.1469 J, and for the complex trajectory K 1 = 0.6854 , K 2 = 0.2661 , K 3 = 0.2278 J. As can be seen from the calculations, the most kinetic energy is consumed by the variable Δ ζ 1 which corresponds to the movement in the longitudinal direction.
Robustness issue. In addition, the robustness of the proposed algorithm to changes in the hovercraft parameters was examined. The inaccuracy of the model parameters was assumed as W = 0.5 (50%).
It turned out that the controller was equally effective in this case (only small index values were obtained), except for kinetic energy, which had higher values (Figure 4) due to the fact that the vehicle had a greater mass and moment of inertia, or a greater displacement of the center of mass. The increase in the value of the kinetic energy is proportional to the increase in the value of the dynamic parameters. For the sine trajectory, it was calculated that K m = 1.0511 , K 1 = 0.6599 , K 2 = 0.11784 , K 3 = 0.2128 J, and for the complex trajectory K m = 1.7586 , K 1 = 1.0355 , K 2 = 0.3873 , K 3 = 0.3357 J.

4.3.2. Comparison with Terminal Sliding Mode Controller

In order to compare the performance and capabilities of the proposed control scheme, the simulation tests were conducted for the same hovercraft model and the same desired trajectories but using a different control method. An algorithm based on the terminal sliding mode control (TSMC) method proposed in [44] was chosen. The following considerations were followed:
(1)
According to the cited reference, the TSMC method provides robustness against unmodeled dynamics, model uncertainty and external interference from ocean currents and waves. It is also known to be superior to the conventional sliding mode control technique in terms of finite-time convergence and high-precision steady-state tracking. It has been successfully applied in many works.
(2)
Despite the fact that the approach based on TSMC is designed for underactuated underwater vehicles, it has been successfully implemented for the hovercraft in [27].
(3)
In the cited article, several algorithms based on different methods were investigated through simulation, and the mentioned TSMC method proved successful in tracking different trajectories while two others proved ineffective.
The comparative of the TSMC algorithm tests was performed taking into account three issues:
(1)
The model with diagonal inertia matrix M d i a g and without friction.
(2)
The model with diagonal inertia matrix M d i a g with friction.
(3)
The model with non-diagonal inertia matrix M s y m .
During the testing, it was noted that the TSMC algorithm is very sensitive to changes in control parameters, and numerical problems also arise (e.g., the integration step of the integration method). For this reason, the values of a set of these parameters were assumed similar to those given in [27] (the model of the tested vehicle has dynamics very similar to the model used in that paper). On the other hand, changing the gain of the controller led to a loss of trajectory tracking efficiency.
The control scheme of the TSMC algorithm was introduced and studied in [44]. The simulation investigation used software from [45] (Matlab/Simulink source codes). The following set of control parameters was assumed for the trials, which provided satisfactory trajectory tracking results:
l x = l y = 1.1 , k x = k y = 1.0 , W 1 = W 2 = 1.0 ,
β 1 = β 2 = 0.6 , q 1 / p 1 = 1 / 2 , q 2 / p 2 = 3 / 5 ,
and e ˜ u ( 0 ) = 0.05 , e ˜ v ( 0 ) = 0 . The saturation of the control signals was limited to | τ u | 30 N and | τ r | 30 Nm to obtain comparable constraints with the previous simulations.
Sine trajectory without friction and using M d i a g . Due to numerical problems, acceptable results could not be obtained as in the case of the algorithm using IQV for the same step time. Therefore, this time was reduced to Δ t = 0.003 s (ODE 3 integration method). As can be seen in Figure 5a, the trajectory tracking is correct and the position errors are close to zero after about 10 s, as shown in Figure 5b. However, the time to obtain the error convergence is longer than with an IQV-type algorithm. The values of the force and torque decrease significantly after a short time (Figure 5c), but the torque suddenly increases and decreases after about 8 s, which is not good for the performance of the task. The velocity values are acceptable, as can be seen in Figure 5d.
Sine trajectory with friction and using M d i a g . The same step time Δ t = 0.003 s and ODE 3 integration method were applied. From Figure 6a,b, it can be observed that the sine trajectory tracking is not as accurate when frictional forces appear in the dynamics model. However, in this case, the force and torque have smaller values and the suddenly jumping value of the torque is clearly smaller, as shown in Figure 6c. The velocity values are now also slightly lower, as can be seen in Figure 6d.
Complex trajectory without friction and using M d i a g . For this desired trajectory, the numerical problems were even greater than for the previous one. The ODE 3 method was adopted, but the time step was reduced 10 times, i.e., Δ t = 0.0003 s. As can be noted in Figure 7a, and which is also confirmed by the results in Figure 7b, the algorithm works properly, ensuring trajectory tracking. The time to reach steady state is about 10 s, so the controller works quickly. Figure 7c shows that the force and moment values are acceptable. However, after about 15 s, there is a sharp change in the value of the moment due to a change in the shape of the desired trajectory. This also indicates the sensitivity of the algorithm when this situation occurs.
Complex trajectory without friction and using M d i a g . The algorithm did not provide any acceptable results even when reducing the integration step.
Sine and complex trajectory without friction and using M s y m . The algorithm proved to be completely useless and the task of trajectory tracking remained unsolved.
The performance parameters of the TSMC controller are shown in Table 3.

4.3.3. Comparison with Original Controller

The IQV-based control algorithm was also compared with the original algorithm in terms of performing the trajectory tracking task. The same trajectories, operating conditions and control coefficients (42) were assumed.
The performance parameters of the original controller are shown in Table 4. Comparing the parameter values in Table 2 and Table 4, it can be observed that the original algorithm gave better MIA results on the x direction than the IQV algorithm. However, this is an error difference of only 0.2 mm (sine) and 2.5 mm (complex). In the y direction, the error difference was 2.3 mm (sine). For the complex trajectory, the IQV controller delivered a lower value of 2.0 mm. The RMS index in both test cases had lower values after the IQV controller (2.7 mm—sine and 0.9 mm—complex). The smaller values of the MIAC and KE indexes are due to the fact that the original algorithm does not take into account the effect of vehicle asymmetry (additional mass components).
It is also worth noting that the favorable results from the original algorithm cannot be directly compared with those from the IQV algorithm, because the latter works under more difficult conditions resulting from vehicle asymmetry.

4.3.4. Discussion of Results

Tests were conducted on one hovercraft model and two desired trajectories selected for tracking. The goal was to compare the performance in both cases based on runs of different quantities, and to indicate the information possible with the criteria and indicators introduced. Based on the results obtained, the following was determined regarding the control algorithm:
  • It provides information on the asymmetric movement of the hovercraft (with a displacement in the center of mass caused, for example, by an additional load);
  • It makes it possible by simulating the effect of a set trajectory on the vehicle motion, velocity deformation and kinetic energy consumption for the vehicle and for the individual quasi-velocity corresponding to motion in each direction;
  • A small effect of the inaccuracy of the model parameters (with a 50% increase in their values) on the tracking accuracy was noticed, except for the kinetic energy values;
  • If possible for the realization of the set task of trajectory tracking, it is better to choose trajectories described by less complicated mathematical functions.
The fact that no differences were observed in the accuracy of tracking the set trajectory, but the dynamics indexes changed, can be explained by the reason that the proposed control scheme is based on the control of kinematic quantities, but at the same time it allows the detection of changes caused by modifications in vehicle dynamics.
Comment. Only an approximate friction model is included in the proposed control algorithm, which limits the controller’s applicability under conditions where disturbances may be negligible. Such a situation is the case, for example, in laboratory conditions or under favorable atmospheric conditions that will not disturb the work of the controller. In addition, a t a n h continuous function approximation was used in the simulations to avoid model discontinuities.
Based on the comparative tests conducted, some observations can be made about the performance of the proposed IQV algorithm and the TSMC method. The results relating to trajectory tracking errors are given in Table 5. The percentage increase UP was calculated according to the following formula:
U P = 1 W I Q V W T S M C · 100 % .
As can be seen from this table, the IQV controller yielded lower values than the TSMC controller for all cases considered.
Table 6 shows a comparison of the percentages of the indexes used in the tests. The minus symbol for the IQV controller indicates that the values were worse. Although in 3 cases the IQV controller gave worse index values, these are position errors of 0.2 mm to 2.5 mm. It is also significant that the initial position in the y direction was 1 m.
Benefits of IQV controller compared to TSMC algorithm.
  • The TSMC method can give satisfactory results but for a diagonal inertia matrix. For a symmetric matrix, the problem of tracking the desired trajectory, as shown by the tests performed, was not solved at all. This means that if the model errors are significant, the algorithm will be ineffective.
  • For the proposed IQV algorithm, no such significant numerical problems arise as in the TSMC method.
  • The applicability of the proposed IQV algorithm is greater than that of the TSMC scheme, because the former is suitable for both a model with a diagonal inertia matrix and a symmetric matrix, and it also takes into account functions representing friction.
  • The selection of the control parameters in the case of the offered algorithm is intuitive, while the controller parameters in the TSMC method must be selected with a lot of time unless another selection method is added. Even then, when a set of parameters is selected correctly, a small change in them can lead to a loss of performance of the algorithm. Such an inconvenience is not present in the control scheme proposed in this paper.
Benefits of IQV controller compared to original algorithm.
  • An IQV controller is suitable for asymmetric vehicles in the longitudinal and lateral directions, so the geometric center can be different from the center of mass.
  • During the realization of the tracking task of the desired trajectory, information on changes in the behavior of the vehicle due to the displacement of the center of mass is obtained simultaneously. Thus, it is a combination of a control algorithm and a scheme for analyzing vehicle dynamics.
  • It is possible to estimate the effect of dynamic parameters on changes in vehicle velocity and kinetic energy reduced by each quasi-velocity (taking into account couplings). In the classical algorithm, the effect of couplings is ignored.
Advantages, drawbacks and limitations of the IQV controller
  • Because only two controller parameters and three design parameters are changed, it may therefore be difficult to achieve very accurate tracking (however, this may be because according to the mathematical proof only uniform stability is guaranteed).
  • The control scheme is suitable for studying hovercraft with rudders, while so far there is no way to directly extrapolate to marine vehicle models with rudders.
The simulation verification involved testing a control algorithm using a model based on a velocity transformation containing components due to the fact that the geometric center does not coincide with the center of mass. New variables including the parameters of this model were incorporated into the control algorithm resulting in a useful tool for estimating vehicle behavior during trajectory tracking. Based on the comparison with the TSMC-type algorithm, better values of the adopted indicators were obtained, while in comparison with the original algorithm the results did not always yield a better performance. Nevertheless, a feature of the IQV controller is that it is a useful tool for estimating the impact of dynamic parameters during trajectory tracking.

5. Conclusions

This paper proposes a modified algorithm to track the desired trajectory. The original control scheme is designed for a hovercraft model described by a diagonal inertia matrix, whereas the design developed in this paper is for controlling a vehicle with an asymmetric dynamic model (caused, for example, by uneven load distribution). It is therefore more general than the classical version but does not allow the identification of model parameters because it is used to study the effect of dynamic parameters on the movement of the hovercraft. To achieve this, a velocity transformation using decomposition of the inertia matrix is utilized which leads to new equations but with a diagonal inertia matrix, and then only the control scheme is designed. It is also shown what mathematical conditions such a scheme must meet, and the performance criteria are adapted to evaluate its performance. The simulations performed for two different set trajectories allow comparing the performance of the algorithm under the assumed operating conditions. The presented approach can be helpful at the stage of motion planning when it is necessary to decide on the route of the vehicle, taking into account the kinetic energy consumption. Analogous tests as presented in this paper can be performed for different vehicle models.

Funding

This work was supported by the Poznan University of Technology Grant No. 0211/SBAD/0122.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Author would like to thank Dariusz Pazderski and Jakub Bernat for their comments on conducting stability proofs.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A

Table A1. Explanation of symbols and abbreviations.
Table A1. Explanation of symbols and abbreviations.
SymbolExplanation
DOFdegree of freedom
IQVinertial quasi-velocities
LFCLyapunov function candidate
R ( ϕ ) transformation matrix between the body frame and the inertial frame
Mvehicle inertia matrix
C ( ν ) Coriolis and centripetal forces vector
D ( ν ) drag vector
τ vector of control inputs (the total force and torque)
b T input scaling coefficient
alength of the arm from the measured center of mass to the surface of the rudder
Ndiagonal vehicle inertia matrix in terms of IQV
Υ velocity transformation matrix
fvector of total forces
τ r total torque
ζ vector of IQV including the vehicle velocities and the inertial parameters
z 1 tracking error in the inertial frame
z 2 = z ˙ 1 time derivative of the tracking error in the inertial frame
z ˙ 2 = z ¨ 1 second time derivative of the tracking error in the inertial frame
δ arbitrary vector for the point to be controlled in the body frame
Bcontrol signal transformation matrix

References

  1. Seguchi, H.; Ohtsuka, T. Nonlinear receding horizon control of an underactuated hovercraft. Int. J. Robust Nonlinear Control 2003, 13, 381–398. [Google Scholar] [CrossRef]
  2. Tanaka, K.; Iwasaki, M.; Wang, H.O. Switching Control of an R/C Hovercraft: Stabilization and Smooth Switching. IEEE Trans. Syst. Man Cybern. Part Cybern. 2001, 31, 853–863. [Google Scholar] [CrossRef] [PubMed]
  3. Chaos, D.; Moreno-Salinas, D.; Aranda, J. Control no lineal de un aerodeslizador no holonomico con acciones de control limitadas [Nonlinear control of an underactuated nonholonomic hovercraft with limited control inputs]. RIAI Rev. Iberoam. Autom. Inform. Ind. 2013, 10, 402–412. [Google Scholar] [CrossRef] [Green Version]
  4. Chaos, D.; Moreno-Salinas, D.; Munoz-Mansilla, R.; Aranda, J. Nonlinear Control for Trajectory Tracking of a Nonholonomic RC-Hovercraft with Discrete Inputs, Hindawi Publishing Corporation. Math. Probl. Eng. 2013, 2013, 589267. [Google Scholar] [CrossRef] [Green Version]
  5. Cremean, L.; William, B.; Dunbar, W.B.; van Gogh, D.; Hickey, J.; Klavins, E.; Meltzer, J.; Murray, R.M. The Caltech Multi-Vehicle Wireless Testbed. In Proceedings of the 41st IEEE Conference on Decision and Control, Las Vegas, NV, USA, 10–13 December 2002; pp. 86–88. [Google Scholar]
  6. Jeong, S.; Chwa, D. Coupled Multiple Sliding-Mode Control for Robust Trajectory Tracking of Hovercraft with External Disturbances. IEEE Trans. Ind. Electron. 2018, 65, 4103–4113. [Google Scholar] [CrossRef]
  7. Jeong, S.; Eom, M.; Chwa, D. Disturbance-Estimation-Based Hierarchical Sliding Mode Control of Hovercraft with Wind Disturbance. In Proceedings of the 2015 15th International Conference on Control, Automation and Systems (ICCAS 2015), in BEXCO, Busan, Republic of Korea, 13–16 October 2015; pp. 532–537. [Google Scholar]
  8. Pröhl, L.; Aschemann, H. Nonlinear Observer-Based Control of an Under-Actuated Hovercraft Vehicle. IFAC PapersOnLine 2020, 53-2, 8973–8978. [Google Scholar] [CrossRef]
  9. Aguiar, A.P.; Cremean, L.; Hespanha, J.P. Position tracking for a nonlinear underactuated hovercraft: Controller design and experimental results. In Proceedings of the 42nd IEEE Conference on Decision and Control, Maui, HI, USA, 9–12 December 2003; pp. 3858–3863. [Google Scholar]
  10. Aguiar, A.P.; Hespanha, J.P. Trajectory-Tracking and Path-Following of Underactuated Autonomous Vehicles with Parametric Modeling Uncertainty. IEEE Trans. Autom. Control 2007, 52, 1362–1379. [Google Scholar] [CrossRef] [Green Version]
  11. Aguiar, A.P.; Hespanha, J.P. Logic-Based Switching Control for Trajectory-Tracking and Path-Following of Underactuated Autonomous Vehicles with Parametric Modeling Uncertainty. In Proceedings of the 2004 American Control Conference, Boston, MA, USA, 30 June–2 July 2004; pp. 3004–3010. [Google Scholar]
  12. Fantoni, I.; Lozano, R.; Mazenc, F.; Pettersen, K.Y. Stabilization of a nonlinear underactuated hovercraft. Int. J. Robust Nonlinear Control 2000, 10, 645–654. [Google Scholar] [CrossRef]
  13. Tran, H.K.; Son, H.H.; Duc, P.V.; Trang, T.T.; Nguyen, H.N. Improved Genetic Algorithm Tuning Controller Design for Autonomous Hovercraft. Processes 2020, 8, 66. [Google Scholar] [CrossRef] [Green Version]
  14. Sira-Ramirez, H. Dynamic Second-Order Sliding Mode Control of the Hovercraft Vessel. IEEE Trans. Control Syst. Technol. 2002, 10, 860–865. [Google Scholar] [CrossRef] [Green Version]
  15. Sira-Ramirez, H.; Ibanez, C.A. The Control of the Hovercraft System: A Flatness Based Approach. In Proceedings of the 2000 IEEE International Conference on Control Applications, Anchorage, AK, USA, 25–27 September 2000; pp. 692–697. [Google Scholar]
  16. Sira-Ramirez, H.; Ibanez, C.A. On the Control of the Hovercraft. Syst. Dyn. Control 2000, 10, 151–163. [Google Scholar] [CrossRef]
  17. Morales, R.; Sira-Ramirez, H.; Somolinos, J.A. Linear active disturbance rejection control of the hovercraft vessel model. Ocean Eng. 2015, 96, 100–108. [Google Scholar] [CrossRef] [Green Version]
  18. de Souza, W.F.; Rafikova, E.; Meza, M.E.M.; Gafurov, S. Backstepping Trajectory Tracking of Underactuated Hovercraft. In Proceedings of the 2018 Global Fluid Power Society PhD Symposium (GFPS), Samara, Russia, 18–20 July 2018. [Google Scholar] [CrossRef]
  19. Kong, X.; Xia, Y.; Hu, R.; Lin, M.; Sun, Z.; Dai, L. Trajectory Tracking Control for Under-Actuated Hovercraft Using Differential Flatness and Reinforcement Learning-Based Active Disturbance Rejection Control. J. Syst. Sci. Complex. 2022, 35, 502–521. [Google Scholar] [CrossRef]
  20. Defoort, M.; Floquet, T.; Kokosy, A.; Perruquetti, W. A novel higher order sliding mode control scheme. Syst. Control Lett. 2009, 58, 102–108. [Google Scholar] [CrossRef] [Green Version]
  21. Kim, K.; Lee, Y.K.; Oh, S.; Moroniti, D.; Mavris, D.; Vachtsevanos, G.J.; Papamarkos, N.; Georgoulas, G. Guidance, Navigation, and Control of an Unmanned Hovercraft. In Proceedings of the 2013 21st Mediterranean Conference on Control & Automation (MED), Platanias-Chania, Crete, Greece, 25–28 June 2013; pp. 380–387. [Google Scholar]
  22. Cabecinhas, D.; Batista, P.; Oliveira, P.; Silvestre, C. Hovercraft Control With Dynamic Parameters Identification. IEEE Trans. Control Syst. Technol. 2018, 26, 785–796. [Google Scholar] [CrossRef]
  23. Xie, W.; Cabecinhas, D.; Cunha, R.; Silvestre, C. Robust Motion Control of an Underactuated Hovercraft. IEEE Trans. Control Syst. Technol. 2019, 27, 2195–2208. [Google Scholar] [CrossRef]
  24. Xie, W.; Cabecinhas, D.; Cunha, R.; Silvestre, C. Global Practical Tracking for a Hovercraft with Unmeasured Linear Velocity and Disturbances. IFAC PapersOnLine 2020, 53, 8959–8964. [Google Scholar] [CrossRef]
  25. Yan, L.; Ma, B.; Xie, W. Robust practical tracking control of an underactuated hovercraft. Asian J. Control 2021, 23, 2201–2213. [Google Scholar] [CrossRef]
  26. Munoz-Mansilla, R.; Chaos, D.; Aranda, J.; Diaz, J.M. Application of quantitative feedback theory techniques for the control of a non-holonomic underactuated hovercraft. IET Control Theory Appl. 2012, 6, 2188–2197. [Google Scholar] [CrossRef]
  27. Herman, P. Assessing the Feasibility of Underwater Vehicle Controllers in Underactuated Hovercrafts via Simulations. Int. J. Model. Identif. Control 2021, 37, 45–56. [Google Scholar] [CrossRef]
  28. Sanders, R.M.W. Control of a Model Sized Hovercraft; DCT Rapporten; Technische Universiteit Eindhoven; University of New South Wales: Kensington, Australia, 2003; Volume 27. [Google Scholar]
  29. Do, K.D.; Pan, J. Control of Ships and Underwater Vehicles; Springer: London, UK, 2009. [Google Scholar]
  30. Patelski, R.; Pazderski, D. Extended State Observer Based Parameter Identification of the Hovercraft System. In Proceedings of the 2022 26th International Conference on Methods and Models in Automation and Robotics (MMAR), Międzyzdroje, Poland, 22–25 August 2022; pp. 330–335. [Google Scholar] [CrossRef]
  31. Patelski, R.; Pazderski, D. Novel Adaptive Extended State Observer for Dynamic Parameter Identification with Asymptotic Convergence. Energies 2022, 15, 3602. [Google Scholar] [CrossRef]
  32. Loduha, T.A.; Ravani, B. On First-Order Decoupling of Equations of Motion for Constrained Dynamical Systems. Trans. ASME J. Appl. Mech. 1995, 62, 216–222. [Google Scholar] [CrossRef]
  33. Do, K.D. Global tracking control of underactuated ODINs in three-dimensional space. Int. J. Control 2013, 86, 183–196. [Google Scholar] [CrossRef]
  34. Elmokadem, T.; Zribi, M.; Youcef-Toumi, K. Trajectory tracking sliding mode control of underactuated AUVs. Nonlinear Dyn. 2016, 84, 1079–1091. [Google Scholar] [CrossRef] [Green Version]
  35. Slotine, J.J.; Li, W. Applied Nonlinear Control; Prentice Hall: Englewood Cliffs, NJ, USA, 1991. [Google Scholar]
  36. Huang, J.; Wen, C.; Wang, W.; Song, Y.D. Global stable tracking control of underactuated ships with input saturation. Syst. Control Lett. 2015, 85, 1–7. [Google Scholar] [CrossRef]
  37. Zhang, M.; Liu, X.; Yin, B.; Liu, W. Adaptive terminal sliding mode based thruster fault tolerant control for underwater vehicle in time-varying ocean currents. J. Frankl. Inst. 2015, 352, 4935–4961. [Google Scholar] [CrossRef]
  38. Qiao, L.; Zhang, W. Adaptive non-singular integral terminal sliding mode tracking control for autonomous underwater vehicles. IET Control Theory Appl. 2017, 11, 1293–1306. [Google Scholar] [CrossRef]
  39. Qiao, L.; Zhang, W. Double-Loop Integral Terminal Sliding Mode Tracking Control for UUVs with Adaptive Dynamic Compensation of Uncertainties and Disturbances. IEEE J. Ocean. Eng. 2019, 44, 29–53. [Google Scholar] [CrossRef]
  40. Kicki, P.; Kozlowski, J. Software: Hovercraft Control With Dynamic Parameters Identification; Unpublished Project; Poznan University of Technology: Poznan, Poland, 2018. [Google Scholar]
  41. Mitchell, A.J.; McGookin, E.W.; Murray-Smith, D.J. Implementation of Sliding Mode Observer Based Reconfiguration in an Autonomous Underwater Vehicle. Acta Polytech. 2005, 45, 18–24. [Google Scholar] [CrossRef] [PubMed]
  42. Herman, P. Application of nonlinear controller for dynamics evaluation of underwater vehicles. Ocean Eng. 2019, 179, 59–66. [Google Scholar] [CrossRef]
  43. Herman, P. Vehicle Dynamics Study Based on Nonlinear Controllers Conclusions and Perspectives. In Inertial Quasi-Velocity Based Controllers for a Class of Vehicles; Springer Tracts in Mechanical Engineering; Springer: Cham, Switzerland, 2022. [Google Scholar] [CrossRef]
  44. Elmokadem, T.; Zribi, M.; Youcef-Toumi, K. Terminal sliding mode control for the trajectory tracking of underactuated Autonomous Underwater Vehicles. Ocean Eng. 2017, 129, 613–625. [Google Scholar] [CrossRef]
  45. Gryczka, N.; Rychlewicz, K. Software: Terminal Sliding Mode Control for the Trajectory Tracking of Underactuated Autonomous Underwater Vehicles; Unpublished Project; Poznan University of Technology: Poznan, Poland, 2018. [Google Scholar]
Figure 1. Underactuated hovercraft model sketch.
Figure 1. Underactuated hovercraft model sketch.
Applsci 13 04965 g001
Figure 2. Simulation results for sine trajectory: (a) desired ( p d ) and realized (p) trajectory; (b) position errors Δ p x , Δ p y (in inertial frame); (c) input signals components (applied force) T cos ψ , T sin ψ ; (d) actual velocities u , v , r ; (e) changes in quasi-velocity ( Δ ζ 1 , Δ ζ 2 ); (f) changes in kinetic energy ( K 1 , K 2 , K 3 , K T ) during motion.
Figure 2. Simulation results for sine trajectory: (a) desired ( p d ) and realized (p) trajectory; (b) position errors Δ p x , Δ p y (in inertial frame); (c) input signals components (applied force) T cos ψ , T sin ψ ; (d) actual velocities u , v , r ; (e) changes in quasi-velocity ( Δ ζ 1 , Δ ζ 2 ); (f) changes in kinetic energy ( K 1 , K 2 , K 3 , K T ) during motion.
Applsci 13 04965 g002
Figure 3. Simulation results for complex trajectory: (a) desired ( p d ) and realized (p) trajectory; (b) position errors Δ p x , Δ p y (in inertial frame); (c) input signals components (applied force) T cos ψ , T sin ψ ; (d) actual velocities u , v , r ; (e) changes in quasi-velocity ( Δ ζ 1 , Δ ζ 2 ); (f) changes in kinetic energy ( K 1 , K 2 , K 3 , K T ) during motion.
Figure 3. Simulation results for complex trajectory: (a) desired ( p d ) and realized (p) trajectory; (b) position errors Δ p x , Δ p y (in inertial frame); (c) input signals components (applied force) T cos ψ , T sin ψ ; (d) actual velocities u , v , r ; (e) changes in quasi-velocity ( Δ ζ 1 , Δ ζ 2 ); (f) changes in kinetic energy ( K 1 , K 2 , K 3 , K T ) during motion.
Applsci 13 04965 g003
Figure 4. Changes in kinetic energy ( K 1 , K 2 , K 3 , K T ) during motion—results with 50% parameter increase in parameter values: (a) sine trajectory; (b) complex trajectory.
Figure 4. Changes in kinetic energy ( K 1 , K 2 , K 3 , K T ) during motion—results with 50% parameter increase in parameter values: (a) sine trajectory; (b) complex trajectory.
Applsci 13 04965 g004
Figure 5. Simulation results for TSMC and sine trajectory (without friction): (a) desired ( p d ) and realized (p) trajectory; (b) position errors Δ p x , Δ p y (in inertial frame); (c) input signals components (applied force and torque) τ u , τ r ; (d) actual velocities u , v , r .
Figure 5. Simulation results for TSMC and sine trajectory (without friction): (a) desired ( p d ) and realized (p) trajectory; (b) position errors Δ p x , Δ p y (in inertial frame); (c) input signals components (applied force and torque) τ u , τ r ; (d) actual velocities u , v , r .
Applsci 13 04965 g005
Figure 6. Simulation results for TSMC and sine trajectory (with friction): (a) desired ( p d ) and realized (p) trajectory; (b) position errors Δ p x , Δ p y (in inertial frame); (c) input signals components (applied force and torque) τ u , τ r ; (d) actual velocities u , v , r .
Figure 6. Simulation results for TSMC and sine trajectory (with friction): (a) desired ( p d ) and realized (p) trajectory; (b) position errors Δ p x , Δ p y (in inertial frame); (c) input signals components (applied force and torque) τ u , τ r ; (d) actual velocities u , v , r .
Applsci 13 04965 g006
Figure 7. Simulation results for TSMC and complex trajectory (without friction): (a) desired ( p d ) and realized (p) trajectory; (b) position errors Δ p x , Δ p y (in inertial frame); (c) input signals components (applied force and torque) τ u , τ r ; (d) actual velocities u , v , r .
Figure 7. Simulation results for TSMC and complex trajectory (without friction): (a) desired ( p d ) and realized (p) trajectory; (b) position errors Δ p x , Δ p y (in inertial frame); (c) input signals components (applied force and torque) τ u , τ r ; (d) actual velocities u , v , r .
Applsci 13 04965 g007
Table 1. Parameters of hovercraft.
Table 1. Parameters of hovercraft.
SymbolValueUnit
m12kg
J1.65kgm 2
d u 0 1.0kg/s
d u 4.5kg/s
d v 0 1.0kg/s
d v 4.5kg/s
d r 0 0.5kg/s
d r 0.8kg/s
d u r 1.0kg/s
d v r 0.5kg/s
Table 2. Performance for IQV controller (with friction).
Table 2. Performance for IQV controller (with friction).
SineComplex
Index TrajectoryTrajectory
MIA Δ p x 0.00450.0072
Δ p y 0.03980.0371
Δ ζ 1 0.02440.0327
Δ ζ 2 0.02440.0327
MIAC T cos ψ 5.14265.9542
T sin ψ 2.49812.9009
RMS | | e | | 0.13650.1307
KE K m 0.70551.1793
Table 3. Performance for TSMC controller.
Table 3. Performance for TSMC controller.
Sine TrajectorySine TrajectoryComplex Trajectory
Index without Frictionwith Frictionwithout Friction
MIA Δ p x 0.07130.12090.0260
Δ p y 0.09550.10710.0487
MIAC τ u 1.57902.45332.0831
τ r 0.35390.41370.3547
RMS | | e | | 0.25670.26570.1593
KE K m 0.62930.58660.9914
Table 4. Performance for classical controller (with friction).
Table 4. Performance for classical controller (with friction).
SineComplex
Index TrajectoryTrajectory
MIA Δ p x 0.00430.0047
Δ p y 0.03750.0391
MIAC T cos ψ 4.76065.5258
T sin ψ 1.82472.7876
RMS | | e | | 0.13920.1316
KE K m 0.68131.0801
Table 5. Percentage increase in IQV tracking performance over TSMC.
Table 5. Percentage increase in IQV tracking performance over TSMC.
Sine TrajectorySine TrajectoryComplex Trajectory
Index without Frictionwith Frictionwithout Friction
MIA Δ p x 93.7%96.3%72.3%
Δ p y 58.3%62.8%23.8%
RMS | | e | | 46.8%48.6%18.0%
Table 6. IQV tracking performance vs. original algorithm [22].
Table 6. IQV tracking performance vs. original algorithm [22].
Sine TrajectoryComplex Trajectory
Index with Frictionwith Friction
MIA Δ p x −4.7%−53.2%
Δ p y −6.1%5.1%
RMS | | e | | 1.9%0.7%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Herman, P. Controller for an Asymmetric Underactuated Hovercraft in Terms of Quasi-Velocities. Appl. Sci. 2023, 13, 4965. https://doi.org/10.3390/app13084965

AMA Style

Herman P. Controller for an Asymmetric Underactuated Hovercraft in Terms of Quasi-Velocities. Applied Sciences. 2023; 13(8):4965. https://doi.org/10.3390/app13084965

Chicago/Turabian Style

Herman, Przemyslaw. 2023. "Controller for an Asymmetric Underactuated Hovercraft in Terms of Quasi-Velocities" Applied Sciences 13, no. 8: 4965. https://doi.org/10.3390/app13084965

APA Style

Herman, P. (2023). Controller for an Asymmetric Underactuated Hovercraft in Terms of Quasi-Velocities. Applied Sciences, 13(8), 4965. https://doi.org/10.3390/app13084965

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop