Open access peer-reviewed chapter - ONLINE FIRST

Autonomous Underwater Vehicle Guidance, Navigation, and Control

By Timothy Sands and Kevin Bollino

Submitted: May 9th 2018Reviewed: July 13th 2018Published: November 5th 2018

DOI: 10.5772/intechopen.80316

Downloaded: 317

Abstract

A considerable volume of research has recently blossomed in the literature on autonomous underwater vehicles accepting recent developments in mathematical modeling and system identification; pitch control; information filtering and active sensing, including inductive sensors of ELF emissions and also optical sensor arrays for position, velocity, and orientation detection; grid navigation algorithms; and dynamic obstacle avoidance among others. In light of these modern developments, this article develops and compares integrative guidance, navigation, and control methodologies for the Naval Postgraduate School’s Phoenix, a submerged autonomous vehicle. The measure of merit reveals how well each of several methodologies cope with known and unknown disturbance currents that can be constant or harmonic while maintaining safe passage distance from underwater obstacles, in this case submerged mines.

Keywords

  • submersible vehicles
  • ocean research
  • obstacle avoidance
  • guidance
  • navigation
  • and control
  • linear quadratic optimal control
  • approximated optimal control
  • reduced-order observers
  • MIMO
  • SIMO

1. Introduction

The Naval Postgraduate School’s consortium for robotics and unmanned systems education and research (CRUSER) uses three autonomous underwater vehicles, the Remus, Aries [1], and Phoenix [2] vehicles to enhance education and research. The oldest vehicle, Phoenix [3] is used in this study to investigate integrated methodologies [4] for vehicle guidance, navigation, and control through a field of obstacles amidst unknown ocean currents that can be approximated by steady state, fixed disturbance ocean velocities, and can also be represented by harmonically oscillating velocities. This integrated approach is a natural extension of the recent innovations. The Phoenix vehicle’s nominal mathematical modeling was articulated in the 1988 article [5] using surge motion to perform system identification. Recent innovations [6, 7, 8, 9, 10] have extended and improved the nominal system identification resulting in high-confidence mathematically modeling in computer simulations. Such simulations permitted Wu et al. [11] to redesign the L1 adaptive control architecture for pitch-control with anti-windup compensation based on solutions to the Riccati equation to guarantee robust and fast adaption of the underwater vehicle with input saturation and coupling disturbances and the approach was applied to the pitch channel alone. Stability was emphasized in the single-channel approach to emphasize dynamic nonlinearities and measurement errors. The Riccati equation is also utilized in this research and proves effective when applied to all six degrees of freedom per [4], where the approach is applied to instances of disturbances that are constant with simultaneous harmonic disturbances simulating unknown ocean currents and waves. In addition to these recent achievements in control, improvements have also been made to guidance and navigation. In recent years, Bo He et al. [12] demonstrated in simulations and open water experiments, the ability to overcome weak data links and sparse navigation data using a technique called extended information filter (EIF) applied to simultaneous localization and mapping (i.e. “SLAM”) that proved computationally easier to implement than traditional extended Kalman filter (EKF) SLAM. Low computational cost is emphasized here to keep the vehicle size low, but also to exaggerate the laudable goal of achieving optimal or near optimal results with methods that are simple. Such is an overt goal of the new research presented here.

Just last year, Yan et al. [13] integrated the navigation system using a modified fuzzy adaptive Kalman filter (MFAKF) to combine traditional strap-down inertial navigation with OCTANS and Doppler velocity log (DVL) to navigate the challenging polar regions where rapidly converging earth meridians and challenging ocean environments filled with submersed obstacles. This benchmark achievement requires the research here to utilize similar challenging ocean conditions, and provide the motivation for selection of simultaneous steady-state ocean currents together with sinusoidal varying unknown wave conditions amidst an ocean filled with obstacles (where here the non-polar ocean is used, so mines are added to fulfill the role of malignant submersed obstacles). Furthermore, simplified waypoint guidance is derived, based on the onboard-calculated distance from the vehicle to a submerged obstacle. The simplified waypoint guidance is proven effective, and should be considered in situations where onboard operation of a modified fuzzy adaptive Kalman filter proves to be computationally prohibitive. The distance to an underwater obstacle was measured by Wang et al. [14] with a novel method: measuring extremely low frequency (ELF) emissions with onboard inductive sensors. Such emissions are produced by ship hulls with relatively pronounced amplitudes compared to small subsurface obstacles, but the harmonic line spectra and fundamental signal frequency relate directly to the closing speed of approach to the obstacle. Experiments proved that even such small signals were detectable at long range with high sensitivity and low-noise sensors of the current state of the art, thus closing distance to obstacles may now be presumed to be known passively, permitting the simplified waypoint guidance proposed in this manuscript. Particularly after ELF queuing, position, orientation, and velocity of obstacles may be monitored optically as developed by Eren et al. [15], and these states may be used as feedback signals together with the waypoint guidance (desired trajectory) permitting augmentation with linear quadratic Gaussian techniques, as done in this manuscript where full-order state observers are together optimized with attitude controller gains, followed by demonstration that reduced-order observers may also be optimized allowing vehicle operators to compensate for individually failed or degraded sensors, or instances where optimally estimated signals are superior to sensor signals in an individual or multiple channel.

Integrating these latest technological developments was demonstrated last year by Wei et al. [16], who integrated the Doppler velocity methods for obstacle monitoring into a dynamic obstacle avoidance scheme for collision avoidance. Following data fusion, a collision risk assessment model is used to avoid collisions, and claims to be effective in unknown dynamic environments, although the experiments did not go so far as to stipulate near-constant ocean currents in addition to harmonic wave actions. These challenging dynamic environments are addressed in this manuscript as a natural extension of the current state of the art.

Autonomous vehicle angular momentum control of rotational mechanics may be achieved using control moment gyroscopes, one potential momentum exchange actuator with a long, historic legacy actuating space vehicles, where mathematical singularities have just recently been overcome [17, 18, 19, 20, 21, 22, 23], permitting use of the actuator for underwater vehicles as done recently achieved by Thorton et al. [24, 25] including combined attitude and energy storage control. These developments suffice to reveal that attitude control is not controversial, and thus the remainder of this manuscript focuses on guidance and navigation with a residual necessity to implement nominal, effective pitch and yaw control.

2. Materials and methods

Submersible vehicles require control systems to guide the vehicle around obstacles that can present dangers to vehicle health and safety in the presence of ocean currents. The challenge addressed here is to navigate the Naval Postgraduate School’s Phoenix submersible vehicle (Figure 1) through a minefield whose dimensions are 200 m × 5100 m in the presence of 0.5 m/s ocean currents. The field will contain at least 30 mines placed at locations using a random number generator. The resulting controller structure has an inner-outer loop structure, and several technologies will be described including pole-placement designs, linear-optimal (quadratic) Gaussian techniques, full and partial order observers for online disturbance identification for ocean currents (both constant lateral underwater ocean currents and also sinusoidal varying currents), tracking systems and feedforward control designed to counter open ocean currents, in addition to integral control. The outer loop controller uses Line-of-Site (LOS) guidance to provide a heading command to the inner loop. The inner loop controller uses output heading feedback to track heading commands. The vehicle is simulated to traverse the minefield and successfully travels no closer than 5 m from any mine and arrive within one half meter from the commanded destination autonomously.

Figure 1.

Submersible vehicle sample and notional minefield [1]: (a) field of randomly placed submersed mines to be avoided by autonomous vehicle and (b) Aries submersible in open ocean.

2.1. System dynamics

The equations of motion used to simulate the dynamic behavior of the autonomous submersible vehicle in a horizontal plane are listed in Eqs. (1)(4). All variables in these equations are assumed to be in nondimensional form with respect to the vehicle length (7.3′) and constant forward speed (∼3 ft./s). The vehicle weighs 435 lbs. and is neutrally buoyant. Time is nondimensionalized such that 1 s represents the time it takes to travel one vehicle length (Figure 2).

mYν´ν´Yr´mxGr´=Yr´ν+Yrmr+Yδsδs+YδbδbE1
mxGNν´ν´Nr´Izr´=Nvν+NrmxGr+Nδsδs+NδbδbE2
ψ´=rE3
y´=sinψ+νcosψE4
In addition to the following dependent equationx´=cosψνsinψE5

Figure 2.

Vehicle geometry and reference axes: (a) Phoenix in open ocean [1] and (b) vehicle geometry and reference axis.

where

The variable areThe constants areYδs = 0.01241.
νlateral (sway) velocitym = 0.0358Yδb = 0.01241
rturning rate (yaw)Iz = 0.0022Nr´ = − 0.00047
ψheading angle (degrees)xG = 0.0014Nν´ = − 0.00178
ylateral deviation (cross-track error)Yr´ = − 0.00178Nr = − 0.00390
δsstern rudder deflectionYv´ = − 0.03430Nv = − 0.00769
δbbow rudder deflectionYr = 0.01187Nδs = − 0.0047
Yv = − 0.10700Nδb = 0.0035

The constant definitions in the mass m, mass moment of inertia with respect to a vertical axis that passes through the vehicle’s geometric center (amidships) Iz, position of the vehicle’s center of gravity (measured positive forward of amidships) xG, with the remaining terms referred to as the hydrodynamic coefficients. These constants are all presented in non-dimensional form.

Defining the state vector xνrψyTand the control uδsδbTand assuming small angles, the dynamics expressed in Eqs. (1)(4) may be expressed in state space form as x´=Ax+Buwhere

A=1.47760.3083001.86731.26820001100010B=0.22710.14541.91591.21120000E6

The system may also be expressed in a transfer function ratio of outputs divided by inputs in Laplace form using Eq. (7) where observer matrix [C] is merely a proper identity matrix to this point of the manuscript. Eq. (7) yields two transfer function relationships between each of the two possible rudder inputs as seen in Eqs. (8) and (9). Notice that both transfer functions have poles and zeros at the origin, while pole-zero cancelation is possible in the case of the stern rudder. On the other hand, even after pole-zero cancelation in the bow rudder Eq. (9), there remains an open loop pole at the origin that must be dealt with during control design, since it represents a potentially unstable element (at the very least, in the instance where the estimated constants are exactly correct, and these equations of motion exactly describe the system, an oscillatory element exists that will not decay). Nonetheless, the dynamics accord to nature. Consider trying to steer a row-boat using the rear rudder. It is much more stable than trying to steer the rowboat using a rudder in the front. This analogy applies to the submersible vehicle and is verified in these results.

GS=CsIA1BE7
GSδsYsδss=0.2271s3+0.875s2s4+2.746s3+1.298s2=s20.2271s+0.875s2s2+2.746s+1.298E8
GSδbYsδbs=1.211s2+1.518ss4+2.746s3+1.298s2=s1.211s+1.518s2s2+2.746s+1.298E9

In Figure 3, the uncontrolled system is analyzed by merely performing a circular turn with each (and then both) rudders. The bow and stern rudders alone are each compared to the combined use of both bow and stern rudders. The bow rudder was deflected +15° for about 21 s, while the stern rudder was deflected for −15° for about 11 s. When both rudders were deflected the maneuver was completed in roughly 8 seconds. Two initial conditions for the sway velocity were investigated (ν0=0and then ν0=8). In all cases, the bow rudder alone performed the poorest, with the stern rudder alone performing the turn in a smaller radius and shorter time. Furthermore, the combined use of both rudders resulting in tightest maneuver.

Figure 3.

Analysis of uncontrolled system: comparison of rudder performance: (a) counter-clockwise turn, ν 0 = 0 and (b) initial sway velocity ν 0 = 8 .

Two simulation methodologies were used to investigate sensitivities to integration method. MATLAB was used with Euler integration, while SIMULINK was used with Runge-Kutta integration with identical timesteps, Δt = 0.1 s. The results were nearly negligible and are displayed in Table 1, from which insensitivity to integration approach is established.

Rudder deflectedEuler: x-distanceRunge-Kutta: x-distanceEuler: y-distanceRunge-Kutta: y-distance
Bow6.54716.54696.86476.8646
Stern3.16653.16653.57683.5768
Both2.45462.45462.65672.6567

Table 1.

Comparison of simulation integration methodology.

Distances calculated to traverse one circular path.


2.2. Control law design

In the system analysis, the optimal rudder implementation scheme was determined to be the application of both rudders, where the rudders were slaved to the same maneuver angle magnitude with the opposite sign, i.e. a “scissored-pair” per Eq. (10). In the case where only variable y is to be measured, the new state space formulation of the system equation components are in Eq. (11). Under the assumption of rudders constrained to behave as a scissored-pair the transfer function from rudder input to output y is given by Eq. (12) whose poles and zeros are listed in Eq. (13), with Eq. (14) revealing the system’s eigenvalues, noting the values are identical to the location of the poles in accordance with theory. The controllability and observability matrices ([CO] and [OB] respectively) are listed in Eq. (15) (whose matrix product [OC] is in Eq. (16)) verifying these system equations are both controllable and observable, since these matrices are full rank, while the determinant of the controllability matrix is 63.1778, a large value with a small value of the matrix condition number, 13.4513. The nonzero determinant of the controllability matrix proves controllability, but to see how close the system is to being uncontrollable, the matrix condition number proves more useful. These two figures of merit indicate the system equations are highly controllable, and accordingly this manuscript will investigate and compare several options for navigation control: pole placement, linear quadratic optimal control, linear quadratic Gaussian, and time optimal control. The same holds true for observability, and thus linear quadratic Gaussian. The matrix product [OC] is the same for every definition of state variables for the given system.

δb=δsE10
A=1.47760.3083001.86731.26820001100010,B=0.08163.127100,C=0001,D=0E11
GSδYsδs=0.08164s22.06s4.773s4+2.746s3+1.298s2E12
polesat:s=0,0,0.6070,2.1388;zerosat:s=6.1279e13,near0,near0E13
eigA=λ=0,0,0.6070,2.1388E14
CO=0.08160.84332.42165.55443.12713.81326.410512.6514003.12710.08163.81326.41052.28381.3916OB=000110101.47760.89170.69170.42170000E15
OC=00.08162.8381.39160.08162.28381.39160.85612.28381.39161.39160.85610.85610.54410.54410.3825E16

Diagonalizing the original system [A] matrix, the spectral decomposition TΛ=ATΛ=T1ATin Eq. (17) may be used to verify a diagonal matrix of eigenvalues [Λ], and then write the system of equations in normal-coordinate form x´=Ax+Bu;y=Cxusing the following transformation: A=Λ=T1AT, B=T1B, and C=CTwhose results are in Eq. (18).

Λ=0.46630.10740010.3033000.46760.00060.4996100111T11.47760.3083001.86731.26820001100010A0.46630.10740010.3033000.46760.00060.499610011TE17
A=2.138800000.60700000000000,B=1.25020286.18888068.4625e78.4625e7,C=0.0006111E18
ubaseline=δ=KυυKrrKψψKyyE19

For the pole placement proportional-derivative (PD) controller articulated in Eq. (19), the poles are set to have roughly the same time constant, while avoiding exactly coincident poles. Gains are iterated for various time constants as displayed in Figure 4, but the following rule of thumb is asserted as well to quickly achieve performance that closely mimics the performance of linear-quadratic optimal (LQR) gains where the control effort and tracking error are equally weighted in the cost function of the optimization.

Figure 4.

Gain values for each state iterated for various time constants.

RULE OF THUMB: Select unity time-constant tcto roughly locate closed-loop poles per Eq. (20). Then place other poles at slightly different locations (e.g. sp=s1±0.01p)

Pole:s1=1tcE20

The gains achieved using the rule of thumb KR.O.T. = {0.5070 –0.3687 -0.7157 -0.1972} have quite different values compared to the gains calculated through the matrix Riccati equation in the linear-quadratic optimization KLQR = {−0.0939 –1.2043 -2.2138 -1}, but nonetheless the resulting behaviors are indeed very similar.

Next, the initial feedback control design was evaluated in simulations where the ship is initially located off the desired track by one ship’s length port side with zero heading, and rudder deflection was limited to 0.4 radians (∼23°). Next, another simulation was performed to test an initial heading angle of 30° starboard where the initial y(0) = 0. The results are displayed in Figure 5(a) and (b) respectively. All state variations were plotted in Figure 4, highlighting the fact that y converges to zero along with the other states. Furthermore, the results of rudder-limited simulations are displayed in Figure 6 and Figure 7 for both scenarios (Table 2).

Figure 5.

Simulations testing the initial baseline feedback controller in two scenarios: (a) initially one ship’s length port side and (b) initial heading 30o starboard.

Figure 6.

State variations for both scenarios simulated using pole-placement gains via rule of thumb: (a) initially one ship’s length port side and (b) initial heading 30o starboard.

Time constantKνKrKψKy
0.5−1.5135−1.7005−5.1508−3.22524
10.5070−0.3687−0.7157−0.1972
21.12480.2870−0.0906−0.0116
LQR−0.0939−1.2043−2.2138−1

Table 2.

Gains for various time constants and also solution to linear quadratic optimization.

Reminder: state definition xνrψyT.


2.3. Observer design

To design a state observer, the system must be observable [4], verifiable through examination of the observability matrix [OB] per Eq. (21), where [C] = [ν r ψ y] = [0 1 1 1]. The condition of the observability matrix reveals the degree of observability, and it is defined by the ratio of maximum to minimal singular values.

OB=CCACA2CAn1E21

Figure 7.

Rudder-limited trajectory track using pole-placement gains via rule of thumb and LQR: (a) initially one ship’s length port side and (b) initial heading 30o starboard.

2.3.1. Full-order observer design

x̂´=Ax̂+Bu+LyCx̂E22
x´x̂´=AxAx̂LCxCx̂E23
e´x´x̂´=ALCxx̂E24
e´=ALCeE25

Assuming that only νmeasurements are available, a mathematical model of the estimated system is in Eq. (22) with a full order observer design using the observer error Eq. (23) leading to the error vector in Eq. (24) allowing the re-expression of Eq. (22) as Eq. (25), where the dynamic behavior of the error vector is determined by the eigenvalues of matrixALC, where Lgains of the observer may be chosen as desired for systems that prove observable, such that the error vector will converge to zero for any stable AKeC. In the following paragraphs, Lis designed by solving the matrix Ricatti equation leading to linear quadratic optimal gains, and also by solving the rule of thumb relationship between gains and time constant as done for the controller gains (Table 3).

Multiple of controller time constant used for observerObserver gain matrix
12tc0.74641.80778.82705.1942T
10tc1031.59093.41211.59530.0020T

Table 3.

Full-order observer gains designed by rule of thumb for various time constants as multiple of controller time constant, tc.

Reminder: 12tcis used in subsequent simulations.


Figure 8 displays the results of simulations revealing the accuracy of state estimation when Lis calculated by the rule of thumb, where the time constant is chosen to be half (tc = 1/2) the time constant of the controller (tc = 1) and the simulation is initialized with the heading angle 30° off, while Figure 9 displays the simulation initialized one boat-length starboard position.

Figure 8.

Simulations starting 30° off heading gains via rule of thumb state observer gains, (a) true and estimated sway velocity, ν(t), (b) true and estimated turning rate, r(t), (c) true and estimated heading angle, ψ(t), (d) true and estimated cross track, y(t).

Figure 9.

Simulations starting 1 boat-length starboard with gains via rule of thumb, (a) true and estimated sway velocity, ν(t), (b) true and estimated turning rate, r(t), (c) true and estimated heading angle, ψ(t), (d) true and estimated cross track, y(t).

2.3.2. Reduced-order observer design

Assuming that some measurements are available from sensors, this paragraph describes the possible iterations and reveals states that are relatively more important to measure with sensors. Four possible output matrices are used to investigate observability. Four options for output matrices Cifor i = 1…4 result in four reduced-order observers OBifor i = 1…4 are detailed in Eqs. (26)(29). Output matrix C1produces an observability matrix OB1with rank = 4 (observable) and determinant not nearly equal to zero. Output matrix C2produces an observability matrix OB2with rank = 4 (observable) and determinant not nearly equal to zero. Output matrix C3produces an observability matrix OB3with rank = 4 (observable) and determinant nearly equal to zero. The matrix condition number is very high indicating the system is barely observable. Output matrix C4produces an observability matrix OB4with rank = 3 (not observable) and determinant equal to zero with a matrix condition number equal to infinity. This means if all other states are measured by sensors, it is not possible to use an observer (even an optimal observer) to determine lateral deviation (cross-track error), y. It is a key state to measure with sensors. The sensor combinations that include y are observable. Using every other sensor, (except y) results in a system that is not observable. Furthermore, measuring y alone results in a barely observable system.

C1=νrψy=0111OB1=CCACA2CAn1=01110.86730.2682101.78235.63521.60742.58790000E26
C2=νrψy=0101OB2=CCACA2CAn1=01010.86731.2682103.649610.76242.87564.77170000E27
C2=νrψy=0001OB3=CCACA2CAn1=000110101.47760.89170.69170.42170000E28
C2=νrψy=1110OB4=CCACA2CAn1=11103.34500.5764106.9019012.18431.76214.09000000E29

Assuming yis to be measured by a sensor, Table 4 reveals that measuring νin addition to y produces the most observable system, and is recommended for designing reduced-order observers. The drawback is measuring ν requires a Doppler sonar, which may not always be available. If all states are measureable except ν the resulting reduced-order observer merely estimates νusing gains on the measureable states displayed in Table 5. Figure 10 reveals very good estimation of νwhen all other states are sensed, and this estimated value of νwas fed to the motion controller in addition to the measured states (the poorly estimated states were neglected instead favoring the more-accurate measurements). State convergence to zero is achieved in the instance of state initialization 30° off-heading. Figure 11 displays similar results for the instance of state initialization one boat-length starboard.

Sensors used to measure statesObservability matrix condition number
yand ν8.8456
yand r21.1306
yand ψ31.2919

Table 4.

Observability matrix condition number for options to supplement y measurement.

Reminder: high condition number is less observable system.


Multiple of controller time constant used for observerObserver gain matrix
110tc0.217400.1164T
2tc0.406900.2179T
10tc0.594100.3182T

Table 5.

Reduced-order observer gains designed by rule of thumb for various time constants as multiple of controller time constant, tc.

Relatively faster 110tcis used in subsequent simulations.


Figure 10.

Simulations starting 30° off heading gains via rule of thumb reduced-order state observer gains: (a) true and estimated sway velocity, ν(t) versus time (seconds), (b) true and estimated turning rate, r(t) versus time (seconds), (c) true and estimated heading angle, ψ(t) versus time (seconds), (d) true and estimated cross track, y(t) versus time (seconds).

Figure 11.

Simulations starting one boat-length starboard with gains via rule of thumb reduced-order state observer gains: (a) true and estimated sway velocity, ν(t) versus time (seconds), (b) true and estimated turning rate, r(t) versus time (seconds), (c) true and estimated heading angle, ψ(t) versus time (seconds), (d) true and estimated cross track, y(t) versus time (seconds).

2.3.3. Gain margin and phase margin

Figure 12 compares the loop gains of the system with and without a compensator via gain margin and phase margin with full-state feedback, while Figure 13 displays the loop gains when output-feedback via observers is used. Each has relative strengths. Full state (theoretical) feedback yields infinite gain margin, yet relatively lower phase margin (usually consider more important of the two), while output feedback (real-world) yields good (but lesser) gain margin with increased phase margin.

Figure 12.

Infinite gain margin and 80.2o phase margin using full state feedback via full-ordered observer with rule of thumb controller gains: (a) root locus real Axis, (b) bode plot frequency (rad/sec).

Figure 13.

61.4° gain margin and 145° phase margin using reduced-order observer (both rule of thumb gains for half-controller tc = 0.5, and compensator with rule of thumb gains (tc = 1), (a) root locus, real axis, (b) bode plot, frequency (rad/sec).

2.4. Tracking systems and feedforward control in the presence of constant disturbance currents

This section evolves the earlier developed system equations and performance analysis by adding non-quiescent conditions, in particular introduction of a lateral underwater ocean current with an absolute velocity, υ0, requiring a modification of the system equations to add the lateral current to Eq. (4) resulting in Eq. (30).

y´=sinψ+νcosψ+υ0E30

2.4.1. Analysis of disturbed system in ocean currents via state equations and simulations

Using the controller (Eq. (19)) and the modified system equations where Eq. (4) is replaced by Eq. (30), and applying the final value theorem: fttsFss0, a steady state value 1/ω + 1 has some variable quantity added to unity for various υ0. Thus, steady-state errors exist in all cases with such disturbances, which are verified by simulations depicted in Figure 14 using gain values from the rule of thumb (ROT) for unity time constant. The steady-state errors are directly proportional to the disturbance magnitude. Figure 15 displays max rudder deflection for the maximal lateral ocean current in the study (to verify the control design continues to remain less than 0.4 radians) where we learned any current greater than 0.4 cannot be eliminated; therefore we next investigate feedforward control and integral control.

Figure 14.

Steady-state position error for various lateral underwater ocean currents.

Figure 15.

Feedback alone unable to counter constant lateral underwater ocean currents, (a) rudder deflection, υ 0 = 0.5 , (b) steady state error vs. υ 0 .

2.4.2. Elimination of steady-state error using feedforward control

Modify the control law to ufeedforward=δ=K1υK2rK3ψK4yK0in order to eliminate the steady-state error, where K0is chosen to insure zero steady-state error, where the feedback gains are chosen by the rule of thumb (Figures 16 and 17).

Figure 16.

Feedforward element included to counter constant lateral underwater ocean currents, (a) rudder deflection, υ 0 = 0.5 , (b) all states when υ 0 = 0.5 .

Figure 17.

Comparison: feedback control with and without feedforward ( υ 0 = 0.5 ).

2.5. Disturbance estimation with reduced-order observer and integral control

Section 2.4 demonstrated feedforward control effectively countered the disturbance currents, but the current was presumed to be known. In to truly be effective, the reduced order observer is next augmented to include estimation of the unknown disturbance current velocity ν̂c, where the observer now estimates the disturbance current velocity, the lateral sway velocity, ν, the lateral deviation (cross-track error), y,and the heading angle ψ. Figure 18a and b display the estimates of the unknown current for two current velocity conditions: ν̂c1=νest1=0.1and ν̂c2=νest2=0.5respectively, while Figure 18c and d display the yand ψstates for each current velocity conditions. Notice how large rudder deflections modify the heading angle to the command-tracking value which counters the disturbance current (sometimes referred to as “crabbing”), and after establishing the crab heading angle, the rudder deflection goes towards zero, illustrating the effectiveness of command tracking.

Figure 18.

Reduced-order observer state estimates versus time (seconds) for two disturbance currents υ c 0 = 0.1 0.5 , where Δ is the rudder deflection using these estimates when the worst-case disturbance current is applied. (a) Sway velocity, (b) disturbance current, (c) lateral deviation (cross-track error), (d) heading angle.

Figure 19 displays all the states versus time in seconds and also the trajectory when a worst-case unknown disturbance current νc=0.5is applied and estimated by the reduced-order observer where the observer gains are solutions to the linear quadratic Gaussian optimization. Meanwhile Figure 20 displays the results in cases utilizing command tracking with reduced order observer and with command: ψ = −0.5 and sinusoidal disturbance current υc0= Asin(0.1 t) but no disturbance estimation or feedforward, while Figure 21 uses disturbance estimation and feedforward and rule of thumb gains. Lastly, Figure 22 displays the performance of reduced-order observers, which is especially useful in instances of limited at-sea computational capabilities.

Figure 19.

Performance with disturbance estimation and command tracking using LQR and rule of thumb gains in reduced order observer, and command tracking to ψ = − 0.5 amidst constant disturbance current υ c = 0.5 ., (a) states, (b) trajectory.

Figure 20.

Utilization of command tracking with reduced order observer, with command: ψ = − 0.5 and sinusoidal disturbance current ν c 0 = Asin 0.1 t but no disturbance estimation or feedforward, (a) all states vs. time (seconds), (b) trajectory.

Figure 21.

Utilization of command tracking with reduced order observer, with command: ψ = −0.5 and sinusoidal disturbance current ν c 0 = Asin 0.1 t and disturbance estimation and feedforward and rule of thumb gains, (a), (b).

Figure 22.

Utilization of command tracking with reduced order observer, with command: ψ = −0.5 and sinusoidal disturbance current ν c 0 = Asin 0.1 t , (a) with disturbance estimation (and feedforward), reduced order observer, (b) with integral control but no disturbance estimation or feedforward.

2.6. Waypoint guidance

A simple line-of-sight guidance routine was employed based on fixing waypoints through a minefield in order to navigate to a specified point and safely return home. The coordinates are fed to a logic determining when to turn per Eq. (31), where d is the distance to the waypoint, and the heading command was autonomously calculated per Eq. (32).

Turnif:xcx2+ycy2dE31
ψcommand=Ktan1ycyxcxE32

Particular attention is brought to the inverse tangent calculation, since quadrant must be preserved in the calculation, since the vehicle will navigate in 360°.

3. Results

The following paragraphs mirror Section 2. Above to provide a concise and precise description of the experimental results, their interpretation as well as the experimental conclusions that can be drawn in each sub-topic introduced and developed so far. Some new development naturally follows in the paragraphs of results, in response to the lessons learned.

3.1. System dynamics

Some basics lessons come from a brief analysis of the uncontrolled system dynamics. The open loop plant equations are potentially unstable (at least persistently oscillatory) with respect to only the bow rudder, while the relationship can be stable with respect to the stern rudder alone. Can be stable is exaggerated to emphasize the presence of pole-zero cancelation, which is an unwise practice (especially in this instance with both poles and zeros at the origin on the stability boundary) unless the estimates for the constants in the system equations are very well known. The analysis of the dynamics also revealed the bow rudder was least relatively-effective at maneuvering alone when compared to the stern rudder, however the bow rudder does enhance vehicle maneuverability when used together with the stern rudder as a “scissored-pair” where the sign of the maneuver angle is opposite for each rudder. This “scissored-pair” constraint simplified the MIMO control design, allowing the design engineer to treat the system as a SISO design, since one rudder’s deflection become a dependent variable constrained to the other rudder’s deflection.

3.2. Control law design

Baseline proportional-derivative control designs effectively stabilized the dynamics, but were ineffective in the presence of a constant lateral open ocean current. Gains selected by rule of thumb performed similar to the linear-quadratic optimal control designs, so this underwater vehicle control could be designed at sea with rudimentary math in instances when higher level computational abilities are not available. Augmentation of the control in include gains tuned to reject the constant current proved effective, but required the current to be measure to permit the control component to be properly tuned. Furthermore, when the lateral disturbance current had sinusoidal variation, the controller was rendered ineffective rejecting the disturbance.

3.3. Observer design

The submersible vehicle’s system equations were verified observable by calculation of a full-ranked observability matrix in Section 2.3. A full sate observer was designed first to permit vehicle control with “full state feedback”, yet without directly measuring velocity. Observer gains may be tuned using classical methods in the general spirit of duality between controller and observers. Their dual nature also permits the matrix Riccati equation to produce optimal gains for a linear-quadratic cost function that exclusively emphasizes state estimation error, unlike the controller optimization where the cost function balanced control effort with state error. State observers permit the vehicle operator to have smooth calculated estimates of all states at all times, which proves useful in the event of sensor interruptions or failures, and reduced-ordered observers may be used in instances where computations on-board the vehicle must be limited, for example to minimize computer size, weight, and/or power.

Especially in light of naturally occurring (roughly) sinusoidal variations in ocean current, the system equations were augmented to include the presumed-unknown disturbance as a state.

3.4. Tracking systems and feedforward control in the presence of disturbance currents

Simple feedforward control elements proved effective against known or estimated constant lateral disturbance currents by allowing the vehicle to autonomously perform “set-and-drift” principles where a highly trained helmsman would turn the bow of a ship into a current, but the simple feedforward elements were ineffective at countering currents with sinusoidal variation. In the set and drift principle the heading is de facto non-zero, so the vehicle cannot simultaneously maintain center-pointing while countering the disturbance. If such a requirement were added, designers must decouple the scissored-pair rudder constraint and design the rudder commands separately to simultaneously counter the disturbance while maintaining centerline pointing.

3.5. Disturbance estimation and integral control

Full-ordered observers effectively estimated constant and sinusoidal disturbance currents and proved useful in the control designs for feedforward control, but furthermore reduced-ordered observer were applied in cases where disturbances were forces and moments and feedforward control was not used. Integral control was used instead to drive steady-state error to zero where sufficiently large time-constants were used for the integrator, i.e. the fifth pole in the pole placement control must be less negative than the other poles.

3.6. Fully assembled system demonstration

In light of all these results, a fully assembled control system was used to navigate the proper mathematical models of the Phoenix autonomous submersible vehicle through a simulated 200 m × 500 m minefield in the presence of unknown ocean currents. The field was populated randomly with 30+ mines, and vehicle successfully traversed the minefield in the presence of an unknown 0.5 m/s current with a miss distance from the nearest mine not less than 5 m, navigating from the starting point to pass within 0.5 m of a commanded en route point at sea, and then return to the start point. The outer loop controller used line-of-sight guidance to provide heading commands to the inner loop, and the inner loop controller was an output-feedback heading controller. Two control strategies both proved effective: Linear-quadratic Gaussian, and approximate optimal pole-placement by rule of thumb. In the linear-quadratic Gaussian case, both the controller gains and observer gains were selected by optimization of the respective matrix Riccati equation. Figure 23 displays the completed maneuver where each dot displays the location of a randomly placed mine. Full state feedback was achieved with state observers via the certainly equivalence principle and the states were utilized in a proportional-derivative-integral feedback control architecture. Detailed outputs and figures of merit are plotted in Figures 24-28 including performance of a second transit of the minefield for validation purposes.

Figure 23.

Navigation through simulated field of 30 randomly placed mines in −0.5 m/s current with linear quadratic Gaussian PID controller and full-state observer.

Figure 24.

Continuous distance (m) to closest mine with linear quadratic Gaussian optimized PID controller and full-state observer versus time (s).

Figure 25.

Linear quadratic (Gaussian) optimal observer convergence with actual value in light-pink near zero, while estimates are depicted oscillating in blue, (a) state ν , (b) state r .

Figure 26.

Linear quadratic (Gaussian) optimal observer convergence, (a) state ψ , (b) command tracking (radians) versus time (s).

Figure 27.

Linear quadratic (Gaussian) optimal observer convergence of y, (c) state y , (d) state.

Figure 28.

Validation trajectory through simulated field of 30 randomly placed mines in −0.5 m/s current with linear quadratic Gaussian optimized PI controller and reduced-ordered observer, (a) second trajectory (results validation), (b) heading command tracking.

4. Discussion

The results of this study establish both classical and modern control paradigms to guide autonomous submersible vehicles through obstacles in unknown ocean currents. Both elegant and simplified autonomous controls proved effective, making this technology immediately assessable to low-end technology implementations. The results are consistent with the significant body of literature on motion mechanics in the presence of unknown disturbances with the added complication of restricted path planning due to randomly placed obstacles, where mines were used in this study driving an additional requirement of minimum safe distance for obstacle passage. This consistency with the current literature leads to a natural direction for future research, since recent innovations in nonlinear idealized (and sometimes also adaptive) methods have recently proven to be natural extensions of technology in these fields.

A natural sequel to this manuscript would utilize the aforementioned methods ([26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39] in particular), which comprise nonlinear mathematical amplifications of the linear methods utilized here. The sequel should include an investigation of idealized nonlinear and adaptive methods with a direct comparison to the current state-of-the art including time-optimal control methods.

Acknowledgments

Conceptualization by Distinguished Professor Healey foremost, and then Professor Kevin Bollino and Sands afterwards; methodology by Healey; Software by Bollino and Sands; validation by Professor Isaac Kaminer, Sands, and Bollino; formal analysis by Sands and Bollino; writing-original draft preparation by Sands; writing-review & editing by Sands; research supervision by Kaminer. Authorship has been limited to those who have contributed substantially to the work reported. Professor Sands is not a tenured faculty member of Stanford University; rather he is the Associate Dean of the Naval Postgraduate School’s Graduate School of Engineering and Applied Sciences. In order to avoid legal jeopardy, Dr. Sands publishes government-funded research under his association with the Naval Postgraduate School, while publishing non-government funded research under his continuing associations with Stanford and Columbia Universities.

Conflict of interest

The authors declare no conflict of interest.

Download

chapter PDF

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Timothy Sands and Kevin Bollino (November 5th 2018). Autonomous Underwater Vehicle Guidance, Navigation, and Control [Online First], IntechOpen, DOI: 10.5772/intechopen.80316. Available from:

chapter statistics

317total chapter downloads

1Crossref citations

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us