Next Article in Journal
Deep Learning Techniques for Grape Plant Species Identification in Natural Images
Next Article in Special Issue
A Sensing and Tracking Algorithm for Multiple Frequency Line Components in Underwater Acoustic Signals
Previous Article in Journal
A Subset-Reduced Method for FDE ARAIM of Tightly-Coupled GNSS/INS
Previous Article in Special Issue
Blind Fractionally Spaced Channel Equalization for Shallow Water PPM Digital Communications Links
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Advanced Autonomous Underwater Vehicles Attitude Control with L 1 Backstepping Adaptive Control Strategy

Department of Mechanical Engineering, University of Connecticut, Storrs CT 06269, USA
*
Authors to whom correspondence should be addressed.
Sensors 2019, 19(22), 4848; https://doi.org/10.3390/s19224848
Submission received: 30 September 2019 / Revised: 28 October 2019 / Accepted: 4 November 2019 / Published: 7 November 2019

Abstract

:
This paper presents a novel attitude control design, which combines L 1 adaptive control and backstepping control together, for Autonomous Underwater Vehicles (AUVs) in a highly dynamic and uncertain environment. The Euler angle representation is adopted in this paper to represent the attitude propagation. Kinematics and dynamics of the attitude are in the strict feedback form, which leads the backstepping control strategy serving as the baseline controller. Moreover, by bringing fast and robust adaptation into the backstepping control architecture, our controller is capable of dealing with time-varying uncertainties from modeling and external disturbances in dynamics. This attitude controller is proposed for coupled pitch-yaw channels. For inevitable roll excursions, a Lyapunov function-based optimum linearization method is presented to analyze the stability of the roll angle in the operation region. Theoretical analysis and simulation results are given to demonstrate the feasibility of the developed control strategy.

1. Introduction

With increasing demands for undersea exploration and exploitation, researches in related fields have been given a boost, especially the study of autonomous underwater vehicles (AUVs), which is an undersea system containing its own power and controlling itself while accomplishing a pre-defined task [1]. A fully autonomous underwater robotics network for various missions with high performance is in great demand [2]. Thus, the design and control strategies of AUVs have received considerable attention by researchers.
The attitude control problem is one of the fundamental problems to facilitate the advancement of autonomous underwater networks. However, the attitude dynamics of AUVs have highly-coupled nonlinearity with uncertainties from modeling errors and time-varying external disturbances. Various methods have been proposed for the attitude control. Some work uses feedback controllers separately for each channel neglecting the coupling among angles, while some handle the problem directly by nonlinear control approaches, such as sliding-mode-based control [3,4,5,6], H tools [7,8] and others.
In this paper, based on good study of the nonlinear AUV model presented in [9], which is driven by four water pumps, a precise attitude controller with the combination of backstepping control and L 1 adaptive control is proposed for the yaw and pitch channels. The backstepping approach provides a recursive method for stabilizing the origin of a system in strict-feedback form [10]. The kinematics and dynamics of the attitude problem for mobile system with 6-DOF are in strict feedback form, which is a good application of backstepping control strategy [11]. In order to handle large time-varying uncertainties in dynamics, which are from modeling, complicated hydrodynamics, and external disturbances, we bring the L 1 adaptive control theory into the backstepping control architecture. The key feature of L 1 adaptive control is guaranteed robustness in the presence of fast adaptation [12]. It is a modification of model reference adaptive control (MRAC) and was initially motivated by aerospace applications. With the combination of high gain and low-pass filter, it guarantees fast adaptation and satisfactory transient response. Theoretical results, analysis details and extensions can be found in [12,13,14,15,16,17,18]. Thus, we chose the L 1 backstepping adaptive control architecture to achieve robust attitude control of AUVs.
Moreover, due to physical design, the AUV could only provide sufficient moments on the pitch and yaw channels. The transient response analysis for pitch and yaw control is well studied with inspiration form the works in [11,19]. Unwanted roll excursions are inevitable and dynamically coupled into yaw and pitch motion [8]. The roll motion can be ambiguous without active roll stabilization, especially for smaller AUVs which have a relatively small stabilization moment due to limited vertical distance from the center of gravity to the center of buoyancy. Thus, we introduce a Lyapunov function-based optimum method to analyze the stability of roll angle in the operation region. Simulation results are presented to study and improve the design.
The rest of paper is organized as follows: Section 2 briefly introduces the advanced AUV design presented in [9]. Section 3 presents the nonlinear model and the attitude control. Section 4 introduces attitude controller design based on L 1 backstepping control only considering uncertainties in dynamics. In Section 5, we provide the proof of guaranteed transient response analysis for the coupled pitch-yaw channels. The stability of the roll angle is also discussed here. The simulation results are presented in Section 6. Section 7 gives the conclusion and future work.

2. AUV Mechanical Design

The attitude controller design is based on a full-scale underwater platform depicted in Figure 1, which utilizes a thrust propulsion system powered by four low-cost submersible pumps forming four water-jet thrusters with reducing end nozzles [9]. This provides an optimal thrust to flow rate ratio. By differentiating the power combinations of four motors, the force and moments needed for the propulsion and maneuvering are generated. Figure 1 shows the working mode of AUV, in which the force vector points to the back and two degree of freedoms (DOFs) of the moment to manipulate the pitch and yaw angles. The manufacturing cost is low compared to the deep-water design which would require more control surfaces.
This AUV is aligned to be neutrally buoyant, which means the gravitational force and buoyant force are equal to each other. It is also aligned to be naturally stable, which means two of the three Euler angles, ϕ and θ , are close to zero when it is at rest. The body coordinate frame O is attached to the center of mass of the AUV. The x axis points to the head and y axis goes to the right of the AUV. More modeling details are presented in Appendix A.

3. Problem Formulation and Control Objective

3.1. Kinematics and Dynamics of Attitude

The Euler angle representation is used here to present the attitude of the AUV. Table 1 gives some of the definitions and notation in this paper.
A standard notation for attitude propagation equations and attitude dynamics is adopted here,
Ω ˙ J w ˙ = Ψ ( Ω ) w η b w × J w .
The first line of Equation (1) describes the kinematics of attitude angles, where Ω = ϕ θ ψ T and w = w x w y w z T are state vectors of attitude angles and angle rates expressed in the body frame. ϕ , θ and ψ denote roll, pitch and yaw angle, while w x , w y and w z are angle rates with respect to x, y and z axes of body frame, and
Ψ ( Ω ) = 1 s ϕ t θ c ϕ t θ 0 c ϕ s ϕ 0 s ϕ / c θ c ϕ / c θ .
The second line of Equation (1) describes the dynamics of attitude angles, where the diagonal matrix J R 3 is moment of inertia of the AUV, and η b R 3 is the overall moment applied on the AUV expressed in body frame, which consists of buoyancy–gravity stabilization moment τ G , control input η and fluid dynamic moment τ D , shown as follows,
η b = τ G + η + τ D .
As mentioned before, this could offer two DOFs of moment to manipulate the pitch and yaw angle. Thus, the control signal is η = [ 0 η y η z ] T .
Equation (1) summarizes the overall attitude model of the AUV. The modeling details could be found in [9].
This paper only considers the uncertainties in dynamics. Thus, compared to Equation (1), the attitude dynamics can be represented as follows,
w ˙ = H ¯ ( w ) + H ¯ δ ( w ) + ( K + K δ ) η b + K σ o u t ,
where H ¯ ( w ) = J 1 ( w × J w ) and K = J 1 . The uncertainties, H ¯ δ ( w ) and K δ , are due to the imprecise knowledge of J. σ o u t represents disturbance from environment. This dynamic equation could be
w ˙ = H ¯ ( w ) + K K 1 H ¯ δ ( w ) + K 1 K δ η b + η b + σ o u t .
Expanding the second η b in (5), let H ( w ) = H ¯ ( w ) + τ G + τ D , and defining σ = K 1 K δ η b + K 1 H ¯ δ ( w ) + σ o u t to represent the overall uncertainties, we obtain
w ˙ = H ( w ) + K ( η + σ ) .

3.2. The Trimmed Model for Pitch and Yaw Dynamics

In this section, the original model in Section 3.1 is trimmed for controller design and performance analysis.
Define
I s = 0 1 0 0 0 1 ;
the kinematics Euler angle equations of the pitch and yaw channels are
θ ˙ ψ ˙ = c ϕ s ϕ s ϕ / c θ c ϕ / c θ w y w z , w ˙ y w ˙ z = I s w ˙ = I s H ( w ) + I s K ( η + σ ) .
Define η s = τ y τ z as the control signal. To further simply Equation (8), we introduce the following definitions.
w s = I s w = w y w z , Ω s = I s Ω = θ ψ , σ s = I s σ , K s = I s K I s T , J s = I s J I s T , Ψ s ( Ω ) = I s Ψ ( Ω ) I s T = c ϕ s ϕ s ϕ / c θ c ϕ / c θ , H s ( w ) = I s H ( w ) = J 22 1 w x w z ( J 11 J 33 ) J 33 1 w x w y ( J 22 J 11 ) ,
where J 11 , J 22 and J 33 are the diagonal elements of J.
Thus, Equation (8) is
Ω ˙ s = Ψ s ( Ω ) w s w ˙ s = H s ( w ) + K s ( η s + σ s ) .
All the uncertainties have been added together into σ s R , which is subject to the following assumptions:
Assumption 1.
There exist constants L > 0 and L 0 > 0 such that the following inequalities hold uniformly in t 0 , w 1 , w 2 R 3 and w s 1 , w s 2 R 2 :
σ ( t , w 1 ) σ ( t , w 2 ) L | | w 1 w 2 | | + L 0
Assumption 2.
There exist constants L 1 > 0 , L 2 > 0 and L 3 > 0 such that the following inequalities hold uniformly in t 0 :
σ ˙ ( t ) L 1 w ˙ ( t ) + L 2 w ( t ) + L 3
Assumption 3.
The desired command Ω d , Ω ˙ d and Ω ¨ d R 2 are bounded. Ω and ω exist in Ω ¯ and w ¯ , which are compact subsets of R 3 , where Ψ ( Ω ) and H ( w ) are continuous and bounded, and so are Ω s and ω s .

3.3. Control Objective

The control objective is to design an adaptive attitude controller, which could let attitude angles to track the desired commands Ω d . In this paper, the goal is to let the pitch and yaw angles, Ω s , always track the desired values, Ω d . Meanwhile, the roll angle is self-stable. The self-stability of the roll angle will be proved in the analysis part.

4. Controller Design

4.1. State Predictor

The state predictor is defined as follows,
w ^ ˙ s = Δ H s ( w ) + K s ( η s + σ ^ s ) + A m e ,
where e = Δ w ^ s w s is the prediction error, and A m is a Hurwitz matrix, which defines the desired convergence of e. Using σ ˜ s = σ ^ s σ s , the prediction error is,
e ˙ = A m e + K ( σ ^ s σ s ) .
The Laplace transform of it is,
e ( s ) = ( s I A m ) 1 K σ ˜ s .

4.2. Adaptive Law

Setting the sampling time of the adaptation law by T s , and the prediction error by e, a piecewise constant adaptation law is given by
σ ^ ( i T s ) = J s ϕ ( T s ) u ( i T s ) , t [ i T s , ( i + 1 ) T s ] ,
where ϕ ( T s ) = A m ( I e x p ( A m T s ) ) 1 and u ( i T s ) = e x p ( A m T s ) e ( i T s ) for all i = 1 , 2 , 3 , . . . .

4.3. L 1 Backstepping Euler Angle Controller

Define η s = [ η y η z ] T , and use η b R 2 and η a R 2 to represent control laws coming from the backstepping loop and L 1 adaptive loop respectively. Thus,
η s = η b + η a ,
where
η b = J s [ A m 2 ( w s w d ) H s ( w ) Ψ s T ( Ω ) ( Ω s Ω d ) ]
η a = C ( s ) ( σ ^ s η w ˙ d )
η w ˙ d = J s w d
w s d = Ψ 1 ( A m 1 ( Ω s Ω d ) + Ω ˙ d )
A m 1 , A m 2 are diagonal Hurwitz matrices. K 1 is a positive gain and D 1 ( s ) is a strictly proper transfer function, the value of which ensure C ( s ) = K 1 D 1 1 + K 1 D 1 has unit DC gain.
For the entire system in (1), the overall control law would be
η = 0 η s .

5. Analysis

The performance analysis of this paper has two parts, namely, self-stability analysis in the roll angle channel and response performance analysis in the pitch-yaw angle channels. The roll angle channel relies on the self stabilization mechanism of itself with assumptions of bounded states in the pitch-yaw channels. Then, the pitch-yaw angle channels stability based on the bounded states in the roll channel is analyzed. The overall stability of the system will be discussed at the end.

5.1. Roll Angle Channel Self-Stability Analysis

Expanding Equation (6), the dynamics of the roll angle is described by the following equations,
ϕ ˙ = w x + s θ t θ c ϕ t θ w y w z
w ˙ x = J 11 1 G W d sin ϕ 1 2 C x s i g n ( w x ) w x 2 J 11 1 w y w z ( J 33 J 22 ) + σ 1 ,
where σ 1 = 1 0 0 σ , representing the overall uncertainty in the roll channel. The first two terms in w ˙ x come from τ G and τ D , where d is the distance from the center of gravity and the center of buoyancy, and C x is the damping coefficient.
Define
d 1 = s θ t θ c ϕ t θ w y w z
d 2 = J 11 1 w y w z ( J 33 J 22 ) + σ 1 ,
Equations (22) and (23) can be written as follows,
ϕ ˙ = w x + d 1
w ˙ x = J 11 1 G W d sin ϕ 1 2 C x s i g n ( w x ) w x 2 + d 2 .
Consider the linearization model,
x ˙ = A x + d r o l l ,
where x = ϕ w x is the state vector, and A = 0 1 g 1 g 2 1 2 C x g 3 is a Hurwitz matrix. g 1 = J 11 1 G W d is a positive constant. g 2 and g 3 are positive linearization coefficients. d r o l l R 2 is the uncertainty that
d r o l l = d + d m ,
where d = d 1 d 2 T is the uncertainty from the other two channels, and d m = 0 d 3 T represents the difference between the linearization model and the original model of x. d 1 , d 2 are defined in (24) and (25), while d 3 is defined as follows.
d 3 = g 1 ( g 2 ϕ sin ϕ ) + 1 2 C x ( g 3 w x s i g n ( w x ) w x 2 ) .
Since A is a Hurwitz matrix, there exist matrices P and Q, which are positive definite. Consider the Lyapunov function candidate:
V r o l l ( x ( t ) ) = x T ( t ) P x ( t ) .
The derivative of V r o l l is
V ˙ r o l l ( x ( t ) ) = x T ( t ) Q x ( t ) + 2 x T ( t ) P d r o l l ( t ) x T ( t ) Q x ( t ) + 2 x T ( t ) P d m ( t ) + 2 λ max ( P ) x ( t ) d .
d 1 and d 2 are continuous and bounded in the compact sets Ω ¯ and ω ¯ , such that
d d 1 + d 2 = b d 1 + b d 2 .
Thus, (31) could be written as
V ˙ ( t ) r o l l x T ( t ) Q x ( t ) + 2 x T ( t ) P d m ( t ) + 2 λ max ( P ) x ( t ) ( b d 1 + b d 2 ) ,
which shows the boundary of V ˙ r o l l is effected by the value of the linearization mismatch error d 3 in d m . With the optimum linearization, proper coefficients, g 2 and g 3 , could give the minimum d 3 .
Lemma 1.
For the V r o l l ( x ( t ) ) and V ˙ r o l l ( x ( t ) ) defined in (30) and (31) if there exist b d 1 , b d 2 , b x such that
i d 1 t 1 b d 1 d 2 t 1 b d 2 .
i i V ˙ r o l l ( x ) 0 , x { x | V r o l l ( x ) = b x 2 λ min ( P ) } .
Then the roll channel is bounded and
x t 1 L b x .
where t 1 > 0 is a dummy variable.
Proof of Lemma 1.
Recalling the Lyapunov function in (30), first we prove that
V r o l l t 1 b x 2 λ min ( P )
by contradiction method as follows. Assume the opposite of Equation (37) is true, then there exists some time t [ 0 , t 1 ] that
V r o l l ( t ) > b x 2 λ min ( P ) .
Since V r o l l ( 0 ) < b x λ min ( P ) , and V ( t ) is continuous, there exists a time t [ 0 , t 1 ] , such that
V r o l l ( t ) = b x 2 λ min ( P ) ,
while
V r o l l ( t ) < b x 2 λ min ( P ) .
Equations (39) and (40) imply that V ˙ ( t ) > 0 , which is clearly in contradiction with Condition ii in (35). Thus, Equation (37) is true, holding for all t [ 0 , t 1 ] ,
V r o l l ( t ) b x 2 λ min ( P ) .
Since λ min ( P ) x 2 x T ( t ) P x ( t ) = V ( t ) , then t [ 0 , t 1 ]
x 2 b x 2 λ min ( P ) λ min ( P ) , b x 2 ,
which leads to (36). □
Remark 1.
For a given pair of g 2 and g 3 , with an appropriate choice of P and Q, the Lyapunov function and the derivative of it are specified, which are shown in Figure 2. The blank region in the center indicts where V ˙ ( x ) > 0 , while the colorful region shows where V ˙ ( x ) < 0 . The contours represent V ( x ) = V c , where V c R . Figure 2 shows a bunch of Lyapunov contours with different values of V c . Any contours within the colorful region satisfy Condition ii in Lemma 1.
Optimum Linearization. In what follows, we introduce an optimum procedure to find g 2 and g 3 , such that V ˙ > 0 is a closed region, which means contours V ( x ) = V c exist, and find the minimum value of V c , denoted as V c m i n .
  • For a given pair of ( g 2 , g 3 ) , with the specified P and Q, we have the objective function L ( x ) = V ˙ ( x ) . The optimization problem under such set up would be
    m i n i m i z e L ( x ) = V ˙ ( x ) s u b j e c t t o V ( x ) = V ¯ i ,
    where V ¯ i R . This optimization problem would search along each contour V ( x ) = V ¯ i to find the maximum value of V ˙ in this contour.
  • Base on step 1, set L ( x ) = = 0 to find the V ¯ 0 , which means in this contour V ( x ) = V ¯ 0 the maximum value of V ˙ is 0.
  • Repeat steps 1 and 2 in the compact sets of g 2 and g 3 ; get a set of V ¯ 0 ( g 2 , g 3 ) .
  • Define the boundary function of states, g ( x ) = B x ( V ¯ 0 ) , where x = g ( x ) .
  • Find the m i n ( B x ( V 0 ) ) in the set of V ¯ 0 ( g 2 , g 3 ) , which give the states’ minimum bounds. Then, the optimal linearization coefficients, g 2 and g 3 , are picked up.

5.2. Pitch-Yaw Angle Channel Stability Analysis

Let:
w ˜ s = Δ w s w d
Ω ˜ s = Δ Ω s Ω d
σ ˜ s = Δ σ ^ s σ s
γ 0 ( T s , ρ w ) = Δ n ( I e x p ( A m T s ) ) A m 1 K s ( L ρ w + L 0 )
γ 1 ( T s , ρ w , ρ η ) = Δ λ max [ ( J s A m ) ( J s A m ) ] γ 0 ( T s , ρ w ) + 2 b d σ T s n
γ 2 ( T s , ρ w , ρ η ) = Δ n | | ( s I A m ) 1 K s | | L 1 γ 1 ( T s , ρ w , ρ η )
b d σ ( ρ w ) = Δ L 1 ( H s ( w ) L 1 + K s L 1 ( ρ η + L ρ w + L 0 ) ) + L 2 ρ w + L 3
ρ f = Δ ( C ( s ) 1 ) J s L 1 w s d + 1 C ( s ) L 1 ( L ρ w + L 0 ) + C ( s ) L 1 γ 1 ( T s , ρ w , ρ η )
Lemma 2.
Considering the system described in (9) with the state predictor (12), adaptive law (15) and control law (16), if the truncated norm w s t 1 L ρ w , η t 1 L ρ η , x t 1 L b x for any time t 1 0 , we have
σ ˜ s t 1 γ 1 ( T s , ρ w , ρ η ) , e t 1 γ 2 ( T s , ρ w , ρ η ) .
Proof of Lemma 2.
The solution of system (13) in the time interval [ ( i 1 ) T s , t + ( i 1 ) T s ] , t [ 0 , T s ] is
e ( ( i 1 ) T s + t ) = e x p ( A m t ) e ( i 1 ) T s + ( i 1 ) T s ( i 1 ) T s + t e x p ( A m ( i 1 ) T s + t τ ) K s σ ^ s ( ( i 1 ) T s ) d τ ( i 1 ) T s ( i 1 ) T s + t e x p ( A ( ( i 1 ) T s + t τ ) ) K s σ s ( τ ) d τ .
When t = T s , it follows from Equation (52) that
e ( i T s ) = e x p ( A m T s ) e ( i 1 ) T s + ( i 1 ) T s i T s e x p ( A m ( i 1 ) T s + T s τ ) K s σ ^ s ( i 1 ) T s d τ ( i 1 ) T s i T s e x p ( A m ( i 1 ) T s + T s τ ) K s σ s ( τ ) d τ .
According to the choice of adaptive law in (15), we have
e x p ( A m T s ) e ( i 1 ) T s + ( i 1 ) T s i T s e x p ( A m ( i T s τ ) ) K s σ ^ s ( ( i 1 ) T s ) d τ = 0 .
It follows from Equation (53) that
e ( i T s ) = ( i 1 ) T s i T s e x p ( A ( i T s τ ) ) K s σ s ( τ ) d τ .
The norm of it is
e ( i T S ) = ( i 1 ) T s i T s e x p ( A m ( i T s τ ) ) K s σ s ( τ ) d τ ( i 1 ) T s i T s e x p ( A m ( i T s τ ) ) K s σ s ( τ ) d τ ( i 1 ) T s i T s e x p ( A m ( i T s τ ) ) K s d τ n σ s ( τ ) = e x p ( A m i T s ) e x p ( A m τ ) | ( i 1 ) T s i T s A m 1 K s n σ s ( i T s ) L = ( I e x p ( A m T s ) ) A m 1 K s n σ s ( i T s ) L .
Since σ s = I s σ , σ s is a subset of σ . Using Assumption 1 in (56), we could have
e ( i T s ) ( I e x p ( A m T s ) ) A m 1 K s ( L w s ( i T s ) L + L 0 ) n .
In the condition of this Lemma | | w s t 1 | | L ρ w , for all i while i T s < t 1 , Equation (57) could be
e ( i T s ) ( I e x p ( A m T s ) ) A m 1 K s ( L ρ x + L 0 ) n .
Using the definition of γ 0 ( T s , ρ x ) in (46), finally we could get
e ( i T s ) γ 0 ( T s , ρ x ) .
For all i T s < t 1 , according to Assumption 2, the derivative of uncertainties could be written in
σ ˙ s t 1 L L 1 w ˙ s t 1 L + L 2 w s t 1 L + L 3 L 1 ( H s ( w ) + K s ( η s + σ s ) ) t 1 L + L 2 w s t 1 L + L 3 L 1 ( H s ( w ) L 1 + K s L 1 ( ρ η + L ρ w + L 0 ) ) + L 2 ρ w + L 3 .
Using the definition of b d σ in (49), then
σ ˙ s t 1 L b d σ .
It follows from (54) that
e ( i T s ) = ( I e x p ( A m T s ) ) e ( i T s ) i T s ( i + 1 ) T s e x p ( A m ( ( i + 1 ) T s τ ) ) K s σ ^ s ( i T s ) d τ = i T s ( i + 1 ) T s e x p ( A m ( ( i + 1 ) T s τ ) ) A m K s e ( i T s ) d τ i T s ( i + 1 ) T s e x p ( A ( ( i + 1 ) T s τ ) ) K s σ ^ s ( i T s ) d τ = i T s ( i + 1 ) T s e x p ( A m ( ( i + 1 ) T s τ ) ) K s ( σ ^ s ( i T s ) + A m e ( i T ) ) d τ .
Hence, (55) and (61) imply that
( i 1 ) T s i T s e x p ( A m ( i T s τ ) ) K s σ s ( τ ) d τ = i T s ( i + 1 ) T s e x p ( A m ( ( i + 1 ) T s τ ) ) K s ( σ ^ s ( i T s ) + A m e ( i T s ) ) d τ ,
and there exists t p [ ( i 1 ) T s , i T s ] such that
K s σ ^ s ( i T s ) + A m e ( i T s ) = K s σ s ( t p ) .
For any t < t 1 , there exists t p [ ( i 1 ) T s , i T s ] such that | t t p | 2 T s which satisfies (63),
σ ^ s ( t ) σ s ( t ) σ ^ s ( t ) σ s ( t p ) + σ s ( t ) σ s ( t p ) σ ^ s ( i T s ) σ s ( t p ) + σ s ( t ) σ s ( t p ) J s A m e ( i T s ) + t p t B σ ˙ s ( τ ) d τ .
The bound of σ ˙ s ( t ) is derived in Equation (60). Then we have
σ ^ s ( t ) σ s ( t ) λ max ( ( J s A m ) ( J s A m ) ) γ 0 ( T s , ρ x ) + 2 b d σ T s n .
It follows from the definition of σ ˜ s ( t ) and γ 1 ( T s , ρ w , ρ η ) in (45), (47) that
| | σ ˜ s t 1 | | γ 1 ( T s , ρ x , ρ η ) .
Using the dynamics in Equation (13), we have
e ( s ) = ( s I A m ) 1 K s ( σ ^ s ( s ) σ s ( s ) ) .
Hence,
| | e t 1 | | L | | ( s I A m ) 1 K s | | L 1 | | ( σ ^ s σ s ) t | | L , | | ( s I A m ) 1 K s | | L 1 γ 1 ( T s , ρ w , ρ η ) .
Using the norm property | | e t 1 | | n | | e t 1 | | L , we have
| | e t 1 | | n | | ( s I A m ) 1 K s | | L 1 γ 1 ( T s , ρ w , ρ η )
where n is the dimension of w s ( t ) . The value of n is equal to 2 in this situation. It follows from the definition of γ 2 ( T s , ρ w , ρ η ) in (48),
| | e t 1 | | γ 2 ( T s , ρ w , ρ η )
which completes the proof. □
Lemma 3.
For any given bounded ρ w , ρ η
lim T 0 γ 0 ( T s , ρ w ) 0
lim T 0 γ 1 ( T s , ρ w , ρ η ) 0
lim T 0 γ 2 ( T s , ρ w , ρ η ) 0 .
Proof of Lemma 3.
Recall the definition of γ 0 ( T s , ρ w ) in (46), the insides of the integration are functions with bounded values. As T s 0 we can get γ 0 ( T s , ρ w ) 0 .
Similarly the limits of γ 1 ( T s , ρ w , ρ η ) and γ 2 ( T s , ρ w , ρ η ) going to zero could be proved. The proof is completed. □
Lemma 4.
Considering the system described in (9) with the state predictor (12), adaptive law (15) and control law (16), if the truncated norm w s t 1 L ρ w , η s t 1 L ρ η , x t 1 L b x , Ω s t 1 ρ Ω for any time t 1 0 , we have
Ω ˜ s t 1 C , w ˜ s t 1 C
where C = Δ max ρ f 2 λ min ( A m 2 T A m 2 ) , ρ f 2 λ min ( A m 1 T A m 1 ) λ min ( A m 2 T A m 2 ) .
Proof of Lemma 4.
Considering this Lyapunov function candidate,
V = 1 2 ( Ω ˜ s T Ω ˜ s + w ˜ s T w ˜ s ) ,
first, we prove that
V ( t ) C .
Based on the control laws in (16) and the definition in (20) (43) and (44), the derivative of (75) is:
V ˙ = Ω ˜ ˙ s T Ω ˜ s + w ˜ ˙ s T w ˜ s = ( Ψ s w s Ω ˙ d ) T Ω ˜ s + ( H s ( w ) + K s ( η + σ s ) w ˙ d ) T w ˜ s = ( Ψ s w d Ω ˙ d ) T Ω ˜ s + ( Ψ s w ˜ s ) T Ω ˜ s + ( H s ( w ) + K s ( η s + σ s ) w ˙ d ) T w ˜ = Ω ˜ s T A m 1 Ω ˜ s + Ω ˜ s T ( Ψ s w ˜ s ) + ( H s ( w ) + K s ( η s + σ s ) w ˙ d ) T w ˜ s = Ω ˜ s T A m 1 Ω ˜ s + ( Ψ s T Ω ˜ s + H s ( w ) + K s ( η s + σ s ) w ˙ d ) T w ˜ s = Ω ˜ s T A m 1 Ω ˜ s + w ˜ s T A m 2 w ˜ s + [ ( C ( s ) 1 ) J s w ˙ d + ( 1 C ( s ) ) σ s C ( s ) σ ˜ s ] T w ˜ s = Ω ˜ s T A m 1 Ω ˜ s + w ˜ s T A m 2 w ˜ s + f ( t , w s ) w ˜ s
where A m 1 and A m 2 are Hurwitz matrices, which ensure that the total of the first two quadratic terms would be negative.
For f ( t , w s ( t ) ) = [ ( C ( s ) 1 ) J s w ˙ d + ( 1 C ( s ) ) σ s C ( s ) σ ˜ s ] T , according to Assumption 1 and Lemma 1,
f ( t , w s ( t ) ) ( C ( s ) 1 ) J s s L 1 w d + 1 C ( s ) σ ( t ) s + C ( s ) σ ˜ s ( C ( s ) 1 ) J s s L 1 w d + 1 C ( s ) L 1 ( L ρ w + L 0 ) + C ( s ) L 1 γ 1 ( T s , ρ w , ρ η )
By the definition in (50), thus
f ( t , w s ( t ) ) ρ f .
Next, we prove Equation (76) by contradiction method. Assuming the opposite of Equation (76) is true. Then if at any time t 1 > 0 , one has V ( t 1 ) > C . According to (75), we could have Ω ˜ s > C or w ˜ s > C .
Based on this, two-case discussions under this are presented as follows,
Case 1. For w ˜ s > C ,thus in (77), we could get
V ˙ = Ω ˜ s T A m 1 Ω ˜ s + w ˜ s T A m 2 w ˜ s + K s f ( t , w s ) w ˜ s Ω ˜ s T A m 1 Ω ˜ s + w ˜ s T A m 2 w ˜ s + ρ f w ˜ s Ω ˜ s T A m 1 Ω ˜ s λ min ( A m 2 ) L 2 w ˜ s L 2 2 + ρ f w ˜ s
which could easily find V ˙ ( t 1 ) < 0 .
Case 2. For Ω ˜ s > C
If w ˜ s C , from (77),
V ˙ = Ω ˜ s T A m 1 Ω ˜ s + w ˜ s T A m 2 w ˜ s + K s f ( t , w s ) w ˜ Ω ˜ s T A m 1 Ω ˜ s + w ˜ s T A m 2 w ˜ s + ρ f w ˜ s ( λ min ( A m 1 ) L 2 Ω ˜ s L 2 2 + ρ f w ˜ s ) + w ˜ s T A m 2 w ˜ s ,
based on the conditions settings for this case, we could get V ˙ ( t ) < 0 .
If w ˜ s > C , according to Case 1, we could easily verify V ˙ < 0 .
Hence, if V ( t 1 ) > C , then from the above Cases 1 and 2 we have
V ˙ ( t 1 ) < 0 .
For t = 0 , it means that if V ( 0 ) > C , V ( t ) will keep decreasing until V ( t ) C . If V ( 0 ) C , then V ( t ) C .
Thus, for all t 1 0 , V ( t ) < C . Since w ˜ s t 1 2 = w ˜ s T w ˜ s V ( t ) , Ω ˜ s t 1 2 = Ω ˜ s T Ω ˜ s V ( t ) , we have
Ω ˜ s t 1 C , w ˜ s t 1 C , t 1 0 .
 □
Lemma 5.
For the system in (2), (6) with the L 1 backstepping adaptive controller in (16), if the truncated L norm Ω s t 1 L ρ Ω , w s t 1 L ρ w , x t 1 L b x and η s t 1 L ρ η for any time t 1 0 , then
η s t 1 L C J s A m 2 + J s H s ( w ) L 1 + C J s Ψ s T ( Ω ) L 1 + C ( s ) L 1 ( L ρ w + L 0 ) + C ( s ) J s s L 1 w d + C ( s ) L 1 γ 1 ( T s , ρ w , ρ η )
and
w s t 1 L w d + C
Ω s t 1 L Ω d + C .
Proof of Lemma 5.
According to (16),
η s ( s ) = η b ( s ) + η a ( s ) = J s ( A m 2 w ˜ H s ( w ) Ψ s T ( Ω ) Ω ˜ s ) + C ( s ) ( σ s J s s w d + σ ˜ s ) ,
taking the norm of this equation, we could have
η s t 1 L J s A m 2 w ˜ s L + J s H s ( w ) L 1 + J s Ψ s T ( Ω ) L 1 Ω ˜ s L + C ( s ) L 1 σ s + C ( s ) J s s L 1 w d + C ( s ) L 1 σ ˜ s .
Use Assumption 1, Lemmas 1 and 4, we could get
η t 1 L C J s A m 2 + J s H s ( w ) L 1 + C J s Ψ s T ( Ω ) L 1 + C ( s ) L 1 ( L ρ w + L 0 ) + C ( s ) J s s L 1 w d + C ( s ) L 1 γ 1 ( T s , ρ w , ρ η ) .
Based on the definition of w ˜ s in (43)
w s = w d + w ˜ s .
Thus, we could get
w s L w d + w ˜ s L w d + C .
Same for Ω s ,
Ω s L Ω d + C .
 □
Lemma 6.
There exist ρ w > 0 , ρ Ω > 0 , ρ η > 0 and T s > 0 such that
C J s A m 2 + J s H s ( w ) L 1 + C J s Ψ s T ( Ω ) L 1 + C ( s ) L 1 ( L ρ w + L 0 ) + C ( s ) J s s L 1 w d + C ( s ) L 1 γ 1 ( T s , ρ w , ρ η ) < ρ η
and
w d + C < ρ w
Ω d + C < ρ Ω .
Proof of Lemma 6.
Let us choose ρ η such that
ρ η = C J s A m 2 + J s H s ( w ) L 1 + C J s Ψ s T ( Ω ) L 1 + C ( s ) L 1 ( L ρ w + L 0 ) + C ( s ) J s s L 1 w d + Δ 1
where Δ 1 > 0 is any positive constant. From Lemma 5.2, there exists some T s to make
C ( s ) L 1 γ 1 ( T s , ρ w , ρ η ) < Δ 1 .
For (94) and (95), we could get the value of the left sides. There must exist a set of ρ w and ρ Ω , which could satisfy the inequalities. □
Theorem 1.
For the system in (2), (6) with the L 1 backstepping adaptive controller in (16), choosing T S to make Lemma 6 hold, if x t 1 L b x then
w s L < ρ w , Ω s L < ρ Ω , η s L < ρ η .
Proof of Theorem 1.
For t = 0 , it satisfy that
w s ( 0 ) < ρ w , Ω s ( 0 ) < ρ Ω , η s ( 0 ) < ρ η .
Using proof by contradiction, assume Theorem 1 is not true, since w s ( t ) , Ω s ( t ) and η s ( t ) are continuous. There exists some t 0 where
w s = ρ w o r Ω s = ρ Ω o r η s = ρ η .
Thus,
w s L ρ w , Ω s L ρ Ω , η s L ρ η .
Letting t 1 = t , following from Lemmas 5 and 6 that
w s t L w d + C < ρ w
Ω s t L Ω d + C < ρ Ω
and
η s t L C J s A m 2 + J s H s ( w ) L 1 + C J s Ψ s T ( Ω ) L 1 + C ( s ) L 1 ( L ρ w + L 0 ) + C ( s ) J s s L 1 w d + C ( s ) L 1 γ 1 ( T s , ρ w , ρ η ) < ρ η .
These contradict what we get in (101) with the assumptions in (100). Therefore t does not exist. Thus the statement in Theorem 1 holds for all t > 0 . □

5.3. Overall Stability Analysis

Theorem 2.
If we choose Ω d and design parameters A m 1 , A m 2 to make ρ w and ρ Ω satisfy the following equation,
i d 1 t 1 b d 1 d 2 t 1 b d 2 .
i i V ˙ r o l l ( x ) 0 , x { x | V r o l l ( x ) = b x 2 λ min ( P ) }
i i i x ( 0 ) b x , w s ( 0 ) < ρ w , Ω s ( 0 ) < ρ Ω
where d 1 and d 2 are defined in (24) and (25), then the entire system is stable and
x L b x , w s L < ρ w , Ω s L < ρ Ω .
Proof of Theorem 2.
We prove this theorem by a contradiction argument, for t = 0 , it satisfies that
x ( 0 ) b x , w s ( 0 ) < ρ w , Ω s ( 0 ) < ρ Ω .
Assume Theorem 2 is not true, since x ( t ) , w s ( t ) and Ω s ( t ) are continuous. There exists some t 0 where
w s ( t ) = ρ w o r Ω s ( t ) = ρ Ω w h i l e x ( t ) b x .
Thus,
x t L b x , w s t L ρ w , Ω s t L ρ Ω .
It follows from Theorem 1 that Equation (98) holds and contradicts with Equation (110). Therefore t does not exist. Thus, the statement in Theorem 2 holds for all t > 0 . □
Remark 2.
The conditions (105), (106) and (107) of Theorem 2 are feasible. In the definition of d 1 and d 2 we can readily obtain that d 1 s θ t θ + c ϕ t θ ρ ω and d 1 J 11 1 ( J 33 J 22 ) ρ ω 2 + 2 L ρ ω + L 0 noticing that σ 1 is subject to Assumption 1. By reducing the pitch angle operating range θ and the desired angular velocity ρ ω , we can always satisfy (107). However, this also limits the operation envelop of the vehicle. The control parameters and reference input should not violate the conditions. The simulation results in next section also validating the feasibility condition of Theorem 2.

6. Simulation Results

6.1. Simulation for Roll Angle Stability Analysis

Following the optimum linearization procedure in Remark 1, the numerical solution for optimum linearization is g 2 = 0 . 89 and g 3 = 1 . 85 . The parameter settings for the simulation are
g 1 = 2 Q = 1 0 . 8 0 . 5 1 . 8 d 1 = 0 . 1 , d 2 = 0 . 1 .
The simulation results are shown in Figure 3. The region of V ˙ > 0 has been reduced by picking up the suitable coefficients g 2 and g 3 . The blue contour is the one who has m a x ( V ˙ ) = 0 . The purple line give the bounds of ϕ and w x , which are ϕ [ 0.7392 , 0.7392 ] and w x [ 0.8987 , 0.8987 ] .

6.2. Pith-Yaw Angle Attitude Control

6.2.1. Simulation Structure

The control law generated by the controller is a desired moment which can not be directly taken by the AUV. A dynamic control allocation module is brought into consideration to distribute the moment among these four thrusters. Similar to the study of quadcopters in [20], the total trust and moment provided by this thruster configuration with 4 thrusters can be expressed as the vector sum of the force and moment from each individual thruster. Using w i to represent each rotor’s speed, the trust provided by the i t h thruster is
T i = b ω i 2 ,
where b is the truster coefficient. Thus, the total trust is given by
T b = i = 1 4 T i = b ( ω 1 2 + ω 2 2 + ω 3 2 + ω 4 2 ) .
Pairwise differences in rotors’ speed drive the vehicle to rotate. The torques about the AUV’s y-axis and z-axis are generated by the moments
τ y = l b ( ω 1 2 ω 2 2 + ω 3 2 + ω 4 2 )
τ z = l b ( ω 1 2 ω 2 2 ω 3 2 + ω 4 2 ) ,
where l is the distance from the thruster axis to the center of mass. Thus, the total trust and moments in the body frame are
T b τ y τ z = b b b b l b l b l b l b l b l b l b l b ω 1 2 ω 2 2 ω 3 2 ω 4 2 = A ω 1 2 ω 2 2 ω 3 2 ω 4 2 .
The pseudoinverse is used to calculate the allocation matrix A 1 . Then, by using
ω 1 2 ω 2 2 ω 3 2 ω 4 2 = A 1 T b τ y τ z ,
the desired moments command could be distributed into RPM commands for each thruster. The simulation structure is shown in Figure 4.

6.2.2. Closed-Loop Response

In this section, we will present the performance of two control methods. One is the proposed L 1 backstepping control, while the other is the PID control. Case 1 presents the simulation results with the designed reference inputs shown in Figure 5. Case 2 provides the simulation results with step functions as reference signals. The closed-loop responses of the L 1 backstepping control in Case 2 are shown in Figure 6a,b, while Figure 7a,b are for the PID method.
The control laws, η y and η z , generated by these two controllers are shown in Figure 6c and Figure 7c, respectively, which are the desired moments with respect to axis y and axis z, namely, moment M and N. The RPM commands sent to the four water pumps are shown in Figure 6c and Figure 7c. Comparing Figure 6a,b, we could see the pitch angle uses less time than yaw angle to achieve the goal. When the desired angle of the pitch channel is reached, the commands are switched to put more effort into the yaw angle channel function.

7. Conclusions

In this paper, a robust L 1 backstepping attitude control has been proposed for AUVs in a dynamic environment. Moreover, a Lyapunov function-based optimum linearization method is presented to analyze the stability of the roll angle in the operation region without active stabilization. Simulation results have been provided to show the effectiveness of the proposed approach. Further research will focus on the performance improvement and precise trajectory-following algorithm design, which eventually will be extended to a fully autonomous underwater robotic network.

Author Contributions

Y.L. and J.C. conceived the idea of the study, performed the research, and analyzed the results under the supervision of C.C.; Y.L. wrote the initial draft of the paper; Y.L. and J.C. revised the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. AUV Modeling

The design of the prototype used in this paper is presented in Section 2. In this section, we would give more modeling details.

Appendix A.1. Proof of Concept Testing

The proof of concept system is designed to be passively stable. In initial testing, this criterion has been sufficiently met with a self-righting time of approximately 1 second from a complete roll over. Initial testing of the propulsion system in a swimming pool environment has displayed an acceptable cruising speed and turn radius.
The fluid level and position of the buoyancy control unit is adjustable as shown in Figure A1. With all motors’ power set to zero, we adjust the fluid level to make the buoyancy force equal. Additionally, the position of the buoyancy is also adjusted to make the attitude of the AUV straight in the water.
This static balancing and alignment can make the vehicle naturally closed to the desired system. Although, the control law from the closed-loop controller can compensate these unbalanced and misaligned factors, the initial alignment can still reduce the control efforts since control signals are subject to physical constraints with limited amplitude, like actuator saturation. It also helps to ensure that the vehicle parameters fall into the fully controllable range.
Figure A1. Prototype side view [9].
Figure A1. Prototype side view [9].
Sensors 19 04848 g0a1

Appendix A.2. Dynamic Model of the AUV

The details of dynamic model are presented in three subsections. The general approaches to building the model of a remotely operated vehicle (ROV), an AUV or other underwater vehicles are discussed in [21,22,23]. Some results are adopted, while further simplification of the modeling are derived to help design a more efficient control system.
• Actuator Dynamics
The dynamic relation of the setting value from the speed controller and real-time speed of motor is described in Equation (A1). A single lag model is used to model the real-time speed ω i of the motor i and the setting input value u i .
ω i ( s ) = 1 T m s + 1 u i ( s ) .
• Force and Moment Generation Process
This part analyses how the force and moment are generated and applied to the AUV. There are three terms applied on the vehicle, which are weight and buoyancy force and moment, water pump propulsion force and moment and AUV fluid dynamic force and moment.
(a) Weight and Buoyancy
The gravitational force and buoyant force in terms of body coordinate systems are,
F w b = m g sin θ cos θ sin ϕ cos θ cos ϕ , F B b = ρ g V sin θ cos θ sin ϕ cos θ cos ϕ ,
where g is the gravitational acceleration, ρ is the fluid density and V is the volume of the fluid displaced by the vehicle. The moments generated by these forces can be expressed in terms of the positions of the center of mass C and the center of the buoyancy B [21].
G W = R C × F W , G B = R B × F B ,
where R C and R B are the respective positions of the center of mass and the center of buoyancy in the local coordinate system.
The AUV is aligned to be neutrally buoyant, which means F w = F B . It is also aligned to be naturally stable, which means the two of three Euler angles ϕ , θ are close to zero. The moment generated by the buoyancy force and gravity force can be simplified as follows,
τ x s = F W d sin ( ϕ ) , τ y s = F W d sin ( θ ) ,
where d is the distance from the gravity center and buoyancy center. The more stabilization moments τ x s , τ y s is always trying to maintain stability, which means the larger the d, the more stable the vehicle.
(b) Water Pump Propulsion Force and Moment
As discussed in Section 6.2.1, the resultant force T b and moments, τ y and τ z , of the proposed propulsion system can be expressed as the vector sum of the force and moment form each individual pump in Equation (116).
(c) AUV Fluid Dynamic Force and Moment
The shape of the AUV is complex and therefore the modeling of its behavior is almost impractical. From the perspective of control, the simplified model at the operation points is desired.
The drag equation of fluid is
f D = 1 2 ρ v 2 C d A ,
where ρ is the density of the fluid, v is the speed of the object relative to the fluid, C d is the drag coefficient, A is the reference area. In the vector form,
F D = 1 2 ρ C d x A x 0 0 0 C d y A y 0 0 0 C d z A z s i g n ( v x ) v x 2 s i g n ( v y ) v y 2 s i g n ( v z ) v z 2 .
Similarly, for the rotational moment, it is
τ D = 1 2 C ω x 0 0 0 C ω y 0 0 0 C ω z s i g n ( ω x ) ω x 2 s i g n ( ω y ) ω y 2 s i g n ( ω z ) ω z 2 .
To summarize, the overall force f b and moment η b applied to AUV are
f b = T b + F D , η b = τ G + η + τ D = τ x s τ y s 0 + 0 τ y τ z + τ D ,
where T b , F D , τ x s , τ y s , τ y , τ z and τ D are defined in Equations (116), (A4) and (A6).
• Rigid Dynamics of the AUV Body
Revisiting the notations in [24], the equations of motion for a rigid body subject to body force f b R 3 and torque η b R 3 applied at the center of mass and specified with respect to the body coordinate frame is given by the Newton–Euler equation in the body coordinate which can be written as,
m I 0 0 J v ˙ b ω ˙ + ω × m v b ω × J ω = f b η b ,
where v b R 3 is the body velocity, ω R 3 is the body angular velocity, m R is the mass, I R 3 × 3 is an identity matrix and J R 3 × 3 is an inertial matrix. The effects of added mass [21] will influence the total mass m and total moment of inertial matrix J. There will be large uncertainties in those parameters as discussed in Equations (4) and (5).
The position and velocity of the AUV center of gravity are given by p = x y z T and v = p ˙ R 3 , respectively, expressed to the spatial frame in North-East-Down orientation. Let R S O be the rotation matrix of the body axes relative to the spatial axes and vector. R can be parameterized by the ZYX Euler angles with ϕ , θ and ψ about the x, y and z axes, respectively.
R ( Ω ) = exp ( z ^ ψ ) exp ( y ^ θ ) exp ( x ^ ϕ ) = c θ c ψ s ϕ s θ c ψ c ϕ s ψ c ϕ s θ c ψ + s θ s ψ c θ s ψ s ϕ s θ c ψ + c ϕ c ψ c ϕ s θ s ψ s θ c ψ s θ s ϕ c θ c ϕ c θ ,
where x ^ = 1 0 0 T , y ^ = 0 1 0 T , z ^ = 0 0 1 T . By differentiating R ( Ω ) respect to time, we have the state equations of the Euler angles Ω ,
Ω ˙ = Ψ ( Ω ) ω ,
where Ψ ( Ω ) is give in Equation (2). In the ZYX Euler angle parameterization of rotation matrix, there are singularities at θ = ± π / 2 . For the following discussion, we assume that the trajectory of AUV does not pass through the singularities. The motion equations of a rigid body are
p ˙ m v ˙ Ω ˙ J ω ˙ = v R T ( Ω ) f b Ψ ( Ω ) ω η b ω × J ω .
Equation (A12) summarizes the overall dynamic system of the AUV, where f b and η b are defined in Equation (A8).

References

  1. Blidberg, D.R. The development of autonomous underwater vehicles (auvs): A brief summary. In Proceedings of the IEEE ICRA, Seoul, Korea, 21–26 May 2001; Volume 4. [Google Scholar]
  2. Bovio, E.; Cecchi, D.; Baralli, F. Autonomous underwater vehicles for scientific and naval operations. Annu. Rev. Control 2006, 30, 117–130. [Google Scholar] [CrossRef]
  3. Logan, C.L. A comparison between H-infinity/mu-synthesis control and sliding-mode control for robust control of a small autonomous underwater vehicle. In Proceedings of the IEEE Symposium on Autonomous Underwater Vehicle Technology (AUV’94), Cambridge, MA, USA, 19–20 July 1994; IEEE: Piscataway, NJ, USA, 1995; pp. 399–416. [Google Scholar]
  4. Cui, R.; Zhang, X.; Cui, D. Adaptive sliding-mode attitude control for autonomous underwater vehicles with input nonlinearities. Ocean Eng. 2016, 123, 45–54. [Google Scholar] [CrossRef]
  5. 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. 2018, 44, 29–53. [Google Scholar] [CrossRef]
  6. Qiao, L.; Zhang, W. Adaptive second-order fast nonsingular terminal sliding mode tracking control for fully actuated autonomous underwater vehicles. IEEE J. Ocean. Eng. 2018, 44, 363–385. [Google Scholar] [CrossRef]
  7. Feng, Z.; Allen, R. Reduced order H-infinity control of an autonomous underwater vehicle. Control Eng. Practice 2004, 12, 1511–1520. [Google Scholar] [CrossRef]
  8. Petrich, J.; Stilwell, D.J. Robust control for an autonomous underwater vehicle that suppresses pitch and yaw coupling. Ocean Eng. 2011, 38, 197–204. [Google Scholar] [CrossRef]
  9. Che, J.; Cernio, J.; Prainito, J.; Zuba, M.; Cao, C.; Cui, J.H.; Kazerounian, K. An advanced Autonomous Underwater Vehicle design and control strategy. In Proceedings of the IEEE 2012 Oceans, Hampton Roads, VA, USA, 14–19 October 2012. [Google Scholar]
  10. Khalil, H.K. Nonlinear Systems; Prentice-Hall: Upper Saddle River, NJ, USA, 2002. [Google Scholar]
  11. Mallikarjunan, S.; Nesbitt, B.; Kharisov, E.; Xargay, E.; Hovakimyan, N.; Cao, C. L1 adaptive controller for attitude control of multirotors. In Proceedings of the AIAA Guidance, Navigation, and Control Conference, Minneapolis, MN, USA, 13–16 August 2012; p. 4831. [Google Scholar]
  12. Hovakimyan, N.; Cao, C. L1 Adaptive Control Theory: Guaranteed Robustness with Fast Adaptation; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2010. [Google Scholar]
  13. Cao, C.; Hovakimyan, N. Design and Analysis of a Novel L1 Adaptive Control Architecture With Guaranteed Transient Performance. IEEE Trans. Autom. Control 2008, 53, 586–591. [Google Scholar] [CrossRef]
  14. Cao, C.; Hovakimyan, N. Stability Margins of L1 Adaptive Control Architecture. IEEE Trans. Autom. Control 2010, 55, 480–487. [Google Scholar]
  15. Hovakimyan, N.; Cao, C.; Kharisov, E.; Xargay, E.; Gregory, I.M. L 1 adaptive control for safety-critical systems. IEEE Control Syst. Mag. 2011, 31, 54–104. [Google Scholar]
  16. Luo, J.; Cao, C.; Yang, Q. L1 adaptive controller for a class of non-affine multi-input multi-output nonlinear systems. Int. J. Control 2013, 86, 348–359. [Google Scholar] [CrossRef]
  17. Che, J.; Santone, M.; Cao, C. Adaptive control for systems with output constraints using an online optimization method. J. Optim. Theory Appl. 2015, 165, 480–506. [Google Scholar] [CrossRef]
  18. Ma, T.; Cao, C. L 1 adaptive output-feedback descriptor for multivariable nonlinear systems with measurement noises. Int. J. Robust Nonlinear Control 2019, 29, 4097–4115. [Google Scholar] [CrossRef]
  19. Che, J.; Cao, C. L1 adaptive control of system with unmatched disturbance by using eigenvalue assignment method. In Proceedings of the 2012 IEEE 51st Annual Conference on Decision and Control (CDC), Maui, HI, USA, 10–13 December 2012; pp. 4823–4828. [Google Scholar]
  20. Corke, P. Robotics, Vision and Control: Fundamental Algorithms in MATLAB® Second, Completely Revised; Springer: Berlin, Germany, 2017; Volume 118. [Google Scholar]
  21. Yuh, J. Modeling and control of underwater robotic vehicles. IEEE Trans. Syst. Man Cybern. 1990, 20, 1475–1483. [Google Scholar] [CrossRef]
  22. Abkowitz, M. Stability and Motion Control of Ocean Vehicles; MIT Press: Cambridge, MA, USA, 1969. [Google Scholar]
  23. Clayton, B.; Bishop, R. Mechanics of Marine Vehicles; E. & F.N. Spon: London, UK, 1982. [Google Scholar]
  24. Koo, T.; Ma, Y.; Sastry, S. Nonlinear control of a helicopter based unmanned aerial vehicle model. IEEE Trans. Control Syst. Technol. 2001. [Google Scholar]
Figure 1. Mechanical design and propulsion system.
Figure 1. Mechanical design and propulsion system.
Sensors 19 04848 g001
Figure 2. V ( x ) and V ˙ ( x ) in the roll channel.
Figure 2. V ( x ) and V ˙ ( x ) in the roll channel.
Sensors 19 04848 g002
Figure 3. Roll channel stability.
Figure 3. Roll channel stability.
Sensors 19 04848 g003
Figure 4. Simulation structure.
Figure 4. Simulation structure.
Sensors 19 04848 g004
Figure 5. Case 1: Simulation results of the L 1 backstepping control.
Figure 5. Case 1: Simulation results of the L 1 backstepping control.
Sensors 19 04848 g005
Figure 6. Case 2A: Simulation results of the L 1 backstepping control.
Figure 6. Case 2A: Simulation results of the L 1 backstepping control.
Sensors 19 04848 g006
Figure 7. Case 2B: Simulation results of the PID control.
Figure 7. Case 2B: Simulation results of the PID control.
Sensors 19 04848 g007
Table 1. Definitions and Notation.
Table 1. Definitions and Notation.
SymbolDefinitionSymbolDefinition
ϕ Roll angle σ s Uncertainty in pitch and yaw channel
θ Pitch angle σ ^ Estimation of σ s
ψ Yaw angle η b Moment command from baseline controller
w x Roll rate η a Moment command from adaptive controller
w y Pitch rate Ω s θ ψ T
w z Yaw rate Ω d θ d ψ d T
s * s i n ( * ) w s x y w z T
c * c o s ( * ) JMoment of inertia of the AUV
t * t a n ( * ) K J 1

Share and Cite

MDPI and ACS Style

Liu, Y.; Che, J.; Cao, C. Advanced Autonomous Underwater Vehicles Attitude Control with L 1 Backstepping Adaptive Control Strategy. Sensors 2019, 19, 4848. https://doi.org/10.3390/s19224848

AMA Style

Liu Y, Che J, Cao C. Advanced Autonomous Underwater Vehicles Attitude Control with L 1 Backstepping Adaptive Control Strategy. Sensors. 2019; 19(22):4848. https://doi.org/10.3390/s19224848

Chicago/Turabian Style

Liu, Yuqian, Jiaxing Che, and Chengyu Cao. 2019. "Advanced Autonomous Underwater Vehicles Attitude Control with L 1 Backstepping Adaptive Control Strategy" Sensors 19, no. 22: 4848. https://doi.org/10.3390/s19224848

APA Style

Liu, Y., Che, J., & Cao, C. (2019). Advanced Autonomous Underwater Vehicles Attitude Control with L 1 Backstepping Adaptive Control Strategy. Sensors, 19(22), 4848. https://doi.org/10.3390/s19224848

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