Scholarly article on topic 'Fault-Tolerant Region-Based Control of an Underwater Vehicle with Kinematically Redundant Thrusters'

Fault-Tolerant Region-Based Control of an Underwater Vehicle with Kinematically Redundant Thrusters Academic research paper on "Mathematics"

Share paper
Academic journal
Mathematical Problems in Engineering
OECD Field of science

Academic research paper on topic "Fault-Tolerant Region-Based Control of an Underwater Vehicle with Kinematically Redundant Thrusters"

Hindawi Publishing Corporation Mathematical Problems in Engineering Volume 2014, Article ID 527315, 12 pages

Research Article

Fault-Tolerant Region-Based Control of an Underwater Vehicle with Kinematically Redundant Thrusters

Zool H. Ismail,1 Ahmad 'A. Faudzi,1 and Matthew W. Dunnigan2

1 Centre for Artificial Intelligence & Robotics, Universiti Teknologi Malaysia, Jalan Semarak, 54100 Kuala Lumpur, Malaysia

2 School of Engineering and Physical Sciences, Heriot-Watt University, Edinburgh EH14 4AS, UK

Correspondence should be addressed to Zool H. Ismail;

Received 20 January 2014; Revised 2 May 2014; Accepted 20 May 2014; Published 11 June 2014

Academic Editor: Mingcong Deng

Copyright © 2014 Zool H. Ismail et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This paper presents a new control approach for an underwater vehicle with a kinematically redundant thruster system. This control scheme is derived based on a fault-tolerant decomposition for thruster force allocation and a region control scheme for the tracking objective. Given a redundant thruster system, that is, six or more pairs of thrusters are used, the proposed redundancy resolution and region control scheme determine the number of thruster faults, as well as providing the reference thruster forces in order to keep the underwater vehicle within the desired region. The stability of the presented control law is proven in the sense of a Lyapunov function. Numerical simulations are performed with an omnidirectional underwater vehicle and the results of the proposed scheme illustrate the effectiveness in terms of optimizing the thruster forces.

1. Introduction

For several decades, extensive efforts have been dedicated to the development of autonomous underwater vehicles (AUVs) for subsea intervention/inspection (primarily the oil and gas industries), scientific exploration, and military purposes. The nature of an AUV, that is, self-contained and capable of performing multipurpose missions, permits it to operate over a wide area without an umbilical cable and human intervention. University and research laboratories became interested in these vehicles and a survey of the research area was given in [1].

Numerous tracking control schemes for an AUV have been reported over the past few years including adaptive control [2, 3], robust control [4], and learning control [5-7]. The aforementioned research works focused on motion errors to keep the vehicle on its desired trajectory. It is interesting to note that a small disturbance such as a wave or sea current would pull the vehicle off the actual trajectory and hence high control effort is needed to keep the vehicle on its trajectory as well as maintaining the velocity of the vehicle. Recently, an adaptive region tracking scheme for an AUV was proposed

by Li et al. [8] to ensure that the vehicle is kept in a specific region. As opposed to a conventional control method, the presented controller sends command signals for propulsion only when the vehicle is outside the desired region. Therefore, the propellers require less energy to ensure that the vehicle is kept inside the region. Therefore, an arbitrarily small region size implies that a high control effort is required for the AUV.

In the unstructured and hazardous ocean environment, an underwater vehicle is liable to experience faults or failures while executing underwater tasks. Typically, thrusters are known to be one of the most likely sources of faults. In certain cases, the existence of a thruster fault results in the termination of the ongoing mission. The implication of these faults could be very expensive and time consuming. Therefore, research and development to produce a fault-tolerant system for the AUV has gained much attention over the past few years. The study of thruster faults on an underwater vehicle with a kinematically redundant thruster system was reported in [9-13]. This type of AUV has no control surface and as an alternative has more thrusters than the required degrees of freedoms (DOFs) to allow hovering and motion in any direction. Additionally, the thruster redundancy allows

the hovering AUV to complete the given mission and to return to the support vessel, even though one or more thrusters have malfunctioned. In [9, 10], a minimum norm solution with a task space control technique is introduced to solve the problem of optimal thruster force allocation for an overactuated underwater vehicle. The excessive number of thrusters is exploited in order to improve the control performance with the thruster fault. This approach resolves the thruster redundancy in the Cartesian space and allows the AUV to track the task space trajectories with asymptotic reduction of the task space errors. Even though previous works have proposed a range of criterions in the thrust allocation, there is still a significant requirement to maintain high thrust demand for a particular thruster since the vehicle must be kept at various positions on the trajectory. In fact, the presented solution for the control allocation problem and tracking method, that is, computed torque controller [9, 10] and PID control [11], do not necessarily reduce the required forces of the individual thrusts and thus possibly generate thrust demands that may exceed an individual thruster's saturation point. In [13], a fuzzy credit assignment cerebellar model articulation controller (FCA-CMAC) neural network with the weighted pseudoinverse approach was recently presented for fault identification and fault-tolerant control of an actuated underwater vehicle. Nevertheless, the proposed method only aims to maintain the propulsion forces and moments at the desired set point with the unanticipated faults.

In this paper, an adaptive region tracking control scheme is proposed for an underwater vehicle with thruster redundancy. Within the proposed control framework, a minimum norm criterion for thruster force distribution is utilized to take account of the thruster faults during the tracking motion. Due to the redundancy, the AUV can still perform the tracking task even though one or more thrusters have malfunctioned. The proposed dynamic region control scheme, where it is formulated in taskspace, aims to reduce the energy consumed by the redundant thrusters and to simplify the allocation of thruster forces. Unlike the conventional tracking control method, the proposed controller activates and sends commands to the functioning thrusters only when the AUV is outside the desired region and hence it significantly reduces energy consumption. It is impossible to obtain an exact knowledge of the dynamic parameters of an underwater vehicle. Therefore, an adaptive term is introduced in the proposed scheme to overcome the parameter uncertainties. To estimate the unknown parameters of the dynamical system while ensuring the robustness of the system against external disturbances, a combined adaptive sliding mode controller can be utilized [14]. In the robust nonlinear control literature [15,16], this method was proven to overcome the drawbacks of conventional adaptive or sliding mode control although it will lead to the computational burden and increase the complexity of system due to large number of unknown parameters. The rest of the paper is organized as follows. Section 2 describes the kinematics and dynamic properties of an underwater vehicle with kinematically redundant thrusters. Section 3 states the adaptive region control law with fault-tolerant decomposition. The stability analysis in terms of

Earth-fixed framed

Figure 1: An AUV expressed in body-fixed and earth-fixed reference frames.

the Lyapunov technique is also given in this section. Simulation results are given in Section 4. Finally, Section 5 contains concluding remarks.

2. Kinematics and Dynamics

In this section, the structure and properties of an underwater vehicle with kinematically redundant thrusters are briefly discussed. The underwater vehicle state vectors in this paper are conveniently represented according to the Society of Naval Architects and Marine Engineers (SNAME) [17]. The coordinate frames are illustrated in Figure 1.

2.1. Kinematic Model. The relationship between inertial- and body-fixed vehicle velocities can be described using the Jacobian matrix in the following form:

h (%) 03x3

h (thJ

n = l (h) v, (1)

where = [x y z] e R3 and ^ = [00^] e R3 denote the position and the orientation of the vehicle, respectively, expressed in the inertial-fixed frame. /1 and /2 are the transformation matrices expressed in terms of the Euler angles. The linear and angular velocity vectors, v1 = [u v w] e R3 and v2 = [p q r] e R3, respectively, are described in terms of the body-fixed frame.

2.2. Dynamic Model. The dynamic equation of motion for an underwater vehicle has been previously investigated in detail [18]. Due to hydrodynamic effects acting on the system such as added mass, drag, lift, and buoyancy forces, the equation becomes highly nonlinear and coupled. Let the velocity state vector with respect to the body-fixed frame be defined by

v e R6; the underwater vehicle dynamic equation can be expressed in closed form as

Mv + C(v) v + D(v) v + g (tf) = t,

where M is the inertia matrix including the added mass term, C(v) represents the matrix of the Coriolis and centripetal forces including the added mass term, D(v) denotes the hydrodynamic damping and lift force, and g(rj) is the restoring force. For an AUV that has redundant thrusters such that n > 6, the relationship between the vector of generalized force acting on the vehicle and individual thruster forces can be given as follows:

T = EFt

where E is the thruster configuration matrix and Fth is the desired vector of thruster forces. The matrix E represents the geometry of the AUV and its thruster locations to transform the individual thruster force into generalized forces in the body-fixed frame of the AUV. It is interesting to note that the matrix E is a constant matrix whose elements are the coordinates of the location of each thruster with respect to the center of mass of the AUV. In order to allow the AUV to generate forces for three translational motions and moments for rotational motions, the matrix E is kept at rank 6. Therefore, it is assumed that the AUV in this research work is capable of performing the specific motion in 6-dimensional task space provided that the AUV possesses more degrees of freedom than the minimum number required to execute a given task.

The dynamic equation in (2) preserves the following properties [2,18].

Property 1. The inertia matrix M is symmetric and positive definite such that M = MT > 0.

Property 2. C(v) is the skew-symmetric matrix such that C(v) = -Ct(v).

Property 3. The hydrodynamic damping matrix D(v) is positive definite; that is, D(v) = Dt(v) > 0.

Property 4. The dynamic model as described in (2) is linear in a set of dynamic parameters dd e R"p and can be written as

MV + C(v) v + D(v) v + g(n) = Yd (n, v, V)ed, (4)

where Yd(rj, v, V) e R6x"f is a known regression matrix; np is the total number of physical parameters. It is assumed that if the arguments of Yd(-) are bounded, then Yd(-) is bounded.

3. Fault-Tolerant Region Tracking Control Scheme

In this section, a new control law is proposed to achieve the objective of redistribution of the thruster forces between the functioning thrusters in such a way that the AUV is able to track the trajectory of a desired region even if one or

more thrusters are faulty. Within the region-based control framework, the desired moving target is specified by a region [19]. A fault-tolerant region tracking control scheme for an AUV is formulated as follows.

First, a dynamic region with specific shape is defined so that the vehicle converges into the region. This can be viewed as a main objective of the underwater vehicle. The objective function for the desired region of the AUV is defined by the following inequality:


where SrjB = B(t] - t]d) e R6 are the continuous first partial derivatives of the dynamic region, rjd(t) is the time-varying reference point inside the geometric shape, and B(t) is a time-varying and nonsingular scaling factor. It is assumed that rjd(t) and B(t) are bounded functions of time. To achieve the scaling formation, that is, if the scaling factor increases then the size of a desired region also increases, a nonsingular matrix is defined as follows:

B1 0 0 B,

where B1 is the scaling matrix of rj1 and B2 is the scaling matrix of rj2. The scaling of the orientation of the AUV is not required in general, so B2 can be set as an identity matrix. Thus, the scaling matrix B1 is given as

sx 0 0 0 0 0 0

where sx(t), sy(t), and sz(t) are scaling factors. This function is useful when the AUV needs to adapt the moving region, depending on the situation and environment. The corresponding potential energy function for the desired region described in (5) can be specified as

= | [max (0,f(StlB))]2

0, 2 1

where kp is a positive scalar.

f(SnB)<0 f(SnB)>0,

Remark 1. The main contribution of this paper is that an approach called the region function approach is formulated to minimize the required forces amongst functioning thrusters rather than the conventional method of tracking control. The region function approach activates the thrusters only when the AUV is outside the desired region and hence it suggests that less energy is needed to perform the given task. To illustrate this, consider an inequality function f(Sxy) = (x - xd)2 + (y - yd)2 < r2, where sx = sy = 1 and r is a regional bound, expressed in 2D task space. Then, the associated potential energy can be illustrated in Figure 2(a) where P(8qB) e C2 is lower bounded by zero. Next, invoking

Side view 1 ......■ ■ ■

i : '•. : ' ■. 0.8

Figure 2: An illustration of the required energy in 2D Cartesian space when r is set to a positive constant: (a) potential energy and (b) kinetic


the following linear relationship to obtain the required forces of two onboard thrusters:

fth = k\\f(öxy

where Fth is the net thruster forces, k represents a force constant, || • ||2 is the l2-norm, and f(Sxy) is the displacement in a planar plane. If kp is set to 1, then the respective forces can be depicted in Figure 2(b). Note that the functioning thrusters consume energy only when the AUV is outside the desired region. In contrast, Figure 3 shows the potential energy and

net thruster forces when r is set to be arbitrarily small or the region bound becomes a point. In this case, the onboard thrusters will always be active and, therefore, more energy is needed for each thruster to ensure that the AUV tracks the desired point.

Now, differentiating (8) with respect to gives

= kp max (0,f(S^B))


Now, let (10) be represented as the region error 6B in the following form:

eB = max (0, f (8^B))

dfjSfle tyß

From the definition of the Cartesian region function, note that max(0,f(d^B))(df(S^B)/d^B)T = 0 when the AUV is outside the Cartesian region. Once the vehicle reaches the region, this term becomes zero.

Next, a vector vr that is useful is defined as v, = Tl {iid-B-lBS%)-arlB-%

where a is a positive constant. The matrix J-1 represents the inverse of the Jacobian matrix, B is the time derivative of B, and B-1 is the inverse of the scaling matrix. From the arguments of trigonometric functions, these matrices are all bounded. Based on the structure of (11) and (12) and the subsequent stability analysis, a filtered tracking error vector for an underwater vehicle is defined as

r(t) = v - J-1 (yd - B-1BS^b) + aJ-1B-1eB. (13)

In general, the development of the open-loop error system for r(t) can be obtained by premultiplying the inertia matrix with the time derivative of r(t) to yield

Mr + C(V)r + D(V)r + Yd (-)dd = EFt

Yd (■) dd = Mvr + C (v) vr + D (v) vr + g (n), (15) and the derivative of vr in (12) is given as Vr = j-1 {Vd -B-1BStlE)

+ J-1 {fjd - B-lBSnE - {B-1B + B-1 B) SVE) (16) - aj-1 £-1 eE - aj-1 {B-1eE + B-1eE) ,

where ^d(^), tfd(t), and B(t) are all assumed to be bounded functions of time. Next, substituting (13) into (2), the following relationship between the derivative of filtered tracking error vector and thruster forces is obtained:

r = M-1EFth - M-1 [C (v)r + D (v) r + Yd (■) dd]. (17) Replacing r with V - Vr leads to

i]-Jv -JVr = JM 1EFi

-JM-1 [C(v)r + D(v)r + Yd (■)dd],

where q = Jv+jv and Vr is defined in (16). Hence, the mapping between task-space acceleration and the thruster forces can be given as

r\ = uFth - JM-1 [C (v)r + D (v) r + Yd (■) dd ] +}v+ JVr,

where u = JM-1E is a nonsquare thruster control matrix. To find the l2-norm solution for the desired thruster force Fth, the weighted pseudoinverse matrix is utilized as a generalized inverse matrix of u:

Fth = uw* + uw* JM-1 [C (v)r + D (v) r + Yd (■) dd] -uw* {jv + JVr),

where t^* = W-1 uT(uW-1 uT) ; W is a weighted matrix that is used to identify the fault information of each thruster. It is assumed that a faulty thruster(s) has been detected and the detected faulty thruster has completely malfunctioned. Thus, it is also known as the thruster fault matrix [9]. The obtained thruster forces from (20) minimize (l^F^^WF^. In the case of no fault in any of the thrusters, the weighted matrix W is an identity matrix. Hence, uw* = u* is a Moore-Penrose pseudoinverse matrix and (1/2)FthFth will be minimized. Based on the error system development, the l2-norm solution for thruster force in (20), and the subsequent stability analysis, the proposed control law for AUV can be formulated as

Fth = uw* JM-1 [-aJTBeE - Kvr + Yd (■) , (21) where KV is positive constant matrix.

Remark 2. The £m-norm can be used to optimize the individual thrust force instead of minimizing the l2-norm. The £m-norm provides the exact representation of the feasible thrust solution space, while the l2-norm provides an approximation of the feasible solution space. This suggests that the lTO-norm always finds a feasible solution as long as a solution exists within the thrusters' saturation limits [20].

The estimated parameters dd are updated using the following parameter update law:

= -LdYTd (-)r,

where Ld is a symmetric positive definite matrix. Substituting (21) into (19) produces a closed-loop dynamic equation for r(t) as follows:

"wt ij + uwt JM-1 [C (v)r + D (v) r + Yd (■) dd]

(jv + JVr) = uwtJM-1 [-aJTBeE - Kvr + Yd (■) dd] ■

Simplifying (23) leads to

n + JM-1 [C (v)r + D (v) r + Yd (■) dd] - {jv + JVr) = JM-1 [-aJTBeE -Kvr + Yd or

Mr + C(v)r + D(v)r + KVr + aJT (tj) BeE - Yd (■) 6d = 0:

where dd = dd - dd denotes the parameter estimation error. Next, the following nonnegative function is introduced to analyze the stability of the proposed control law:

v = 2rTMr + + aP (8nB). (26)

Differentiating V with respect to time and using the update law (22) yield

V = rTMr - GTyT (■) r + a{BStiE + B8nE)TeE. (27)

Utilizing (17), (20), and (25) and cancelling the common terms lead to

V = -rTD (v) r - rTKVr - aeTEeE < 0, (28)

where Property 3 is used. Now, a new theorem can be stated as follows.

Theorem 3. Given a closed-loop underwater vehicle with kinematically redundant n thrusters in (25), the proposed adaptive control law (21) and parameter update law (22) guarantee the convergence of q into a dynamic region f(8yB) < 0 in the sense that eE ^ 0 and r ^ 0,as t ^ rn.

Proof. In view of the fact that M is uniformly positive definite,

V in (26) is positive definite in r, 9d. Since V > 0 and

V < 0, V is bounded. Hence, r, 9d, and f(SqB) are bounded. The boundedness of f(8qB) ensures the boundedness of (df(8qB)/dqB). Hence, eB in (11) and its derivative are bounded. The boundedness of v guarantees the boundedness of y since y = l(rf)v and J(q) is a trigonometric function. The boundedness of J(q) and J- (q) leads to the boundedness of u(q) and (q) in (19) and (20), respectively. Moreover, Vr in (16) remains a bounded signal by employing standard signal chasing arguments. From the closed-loop equation, it can be concluded that r in (17) is also bounded using the fact that the kinematic and dynamic terms denoted by u(q), (r(), C(v),D(v),and g(q) are bounded for all possible q(t) and v(i). Thus, V is bounded since r = V - Vr. Differentiating (28) with respect to time gives

V = -rTD (v) r - 2rTD (v) r - 2rTKVr - 2oc%eB. (29)

Hence, V is bounded since r, r, eB, and eB are all bounded. Therefore, V is uniformly continuous. Applying Barbalat's lemma [21] gives V 0 and r —> 0. From (11), eB

0 which also indicates that еБ ^ 0 indicates that ек

is approaching zero. Therefore, y converges to the moving desired scalable region. □

Remark 4. To keep the thruster force within the saturation limits, the diagonal elements of the weight matrix can be formulated as [10]

= diag [h1,h2,...,hn]

where n is the number of thrusters and ht is given as follows:

h = 1 +

Note that the second term on the right hand side of (31) is obtained from differentiation of the thruster force variable H with respect to the thrust force of the ith thruster. The same concept was presented in [22], where a similar function based on joint position variables was used to avoid joint limits of a manipulator.

4. Simulation Results

In this section, two simulation studies are carried out to assess the effectiveness of the proposed fault-tolerant region-based control law for an underwater vehicle. The ODIN [9, 23] AUV is known as a near-spherical omnidirectional vehicle mounted with four horizontal thrusters and four vertical thrusters and was chosen for these numerical simulations. The estimated parameters in (22) are 80% to 90% of the original values in order to verify the adaptive performance of the proposed controller. Values for vehicle parameters such as hydrodynamic damping matrices, added masses, and the thruster configuration matrix are given in the appendix. Hydrodynamic damping is assumed to be the summation

-2-, 5

Figure 4: The desired trajectory for the region control scheme where "x" marks the initial position of the AUV.

of linear and quadratic damping whose effects are noncou-pled with respect to other velocities. The performance of conventional tracking control and the proposed technique is observed concerning two cases. The first case is with no thruster fault and the second case is where any one of the four horizontal thrusters is faulty and any one of the four vertical thrusters is faulty.

In these simulation studies, the following inequality function is defined:

f (s1b) = sx (x- xo)2 + sy (у- уо)2

+ sz(z - Zof + s<b((ï - фоу

+ se(e - eoY + sv(f - fo)2 ^ K2

where the elements of , sy, sz, $ф, se, s^j are the components of the time-varying scaling matrix Б and кг is a scalar tolerance. In these simulations, the matrix Б is defined as the identity matrix. Note that (32) can also be represented as the root mean square error for all axes in the following form:

(Sx (x - xo)2 + Sy {у- Уо)2 + Sz(z - zo)2 + s<p (Ф - Фо)2 +s9{e - e0)2 + sY{f - f0)2)m < Kr.

The ODIN vehicle is required to track a predefined trajectory as illustrated in Figure 4 where the green (cross-sectional lines) path is the horizontal basis position initialized at the position [1.5 0 -1.2]T m. Moreover, the vehicle is initialized at the same position ^1(0) = [1.5 0 - 1.2] m while the initial values for attitude are set to be ц2{0) = [0 0 0] degrees. The attitude is kept constant for the entire motion. The control gains are set to the following:

kp = 28; a= 1;

KV = diag ([1.08 1.08 1.08 1.1 1.1 1.1]).

Table 1: Simulation results for Case 2: required thruster forces.

Figure 5: Three-dimensional view for conventional tracking control (xr is set to be arbitrarily small: Kr = 0.001) with no fault.


Figure 6: Three-dimensional view for region tracking control (xr 0.25) with no fault.

The results from Case 1, where the weighted matrix is an identity matrix, are depicted in Figures 5 and 6. The vehicle represented by solid blue lines initially navigates to the desired line and desired region as can be seen in Figures 5 and 6, respectively. The abrupt changes of errors at the beginning are caused by vehicle parameters' errors before it stays at the desired position even though the region is moving or scaling. The position of the AUV is indicated with solid blue lines. Meanwhile, the forces for all the thrusters of Case 1 using region tracking control law are given in Figure 7.

In Case 2, two thrusters on the AUV are faulty; one of the four horizontal thrusters and one of the four vertical thrusters stop functioning after 120 s of motion. The results are plotted in Figures 8 and 9. Even though the fault case happens during the motion, the thruster allocation forces are different from that of Case 1 and the position of the AUV is still kept inside the region. Figure 10 depicts the forces of functioning and faulty thrusters using region tracking control law at various time instants. The total required thruster forces are listed in Table 1, where all values are computed based on the Euclidean

Conventional Proposed control

tracking controller scheme

Horizontal thrusters, (N) 487.42 483.29

Vertical thrusters, (N) 229.53 218.43

Total input, (N) 716.95 701.72

norm for the entire time. The term "Total input" in the table is calculated by numerical summation of the norm of

thruster forces and it can be interpreted as the total energy

requirement. From Table 1, the region function formulation under the proposed control method needs significantly less energy than the conventional controller.

5. Conclusion

In this paper, a new adaptive fault-tolerant region-based tracking control scheme is proposed for an autonomous underwater vehicle with redundant thrusters. Two cases have been considered. The first case is with no thruster fault and the second case is where any one of the four horizontal thrusters is faulty and any one of the four vertical thrusters is faulty. Although some thrusters are faulty during task execution, the AUV is still able to track a desired moving region based on the thruster allocation forces. Moreover, it has been shown that the proposed technique can minimize the required thruster forces unlike the conventional adaptive tracking control. A Lyapunov-like function has been proposed for stability analysis. Simulation results have been presented to demonstrate the performance of the proposed controller.


The omnidirectional intelligent navigator (ODIN) is a near-spherical AUV designed by the University of Hawaii. The dynamic model of ODIN is given by [9]

[MRB + MA] v + [CRB (v) + CA (v)] v + D (v) v + g(n) = r,

where the subscripts RB and A represent the rigid body and added mass terms of the relevant parameters, respectively. The numerical values for the matrices of the vehicle dynamic equation (A.1) are given as

m 0 0 0 mzG 0

0 m 0 -mzG 0 0

0 0 m 0 0 0

0 -mzG 0 ^-xx 0 0

mzG 0 0 0 Iyy 0

0 0 0 0 0 Izz

where Ixx = Iyy = Izz = I = (8/15)npvr{ moments of inertia about the principle axes.

{S -20

'S 0 3

fS -20

fS -20

fS -20

50 100 150 200 250 300

50 100 150 200 250 300

50 100 150 200 250 300

50 100 150 200 250 300

Time (s)

fS -20

fS -50

100 150 200 250

100 150 200 250

50 100 150 200 250 300

50 100 150 200 250 300

Time (s)

Figure 7: Thruster forces for Case 1 using proposed tracking control law. (a) Horizontal thrusters and (b) vertical thrusters.


Figure 8: Three-dimensional view for conventional tracking control (xr is set to be arbitrarily small: Kr = 0.001) with fault.

-2 -L 5'

Figure 9: Three-dimensional view for region tracking control (xr = 0.25) with fault.


Xud 0 0 0 0 0

0 Yvd 0 0 0 0

0 0 Zwd 0 0 0

0 0 0 0 0 0

0 0 0 0 0 0

0 0 0 0 0 0

Crb (v) =


-m (u2 - zGu4) m (ui + zGu5) 0

-mu m (u2 - znu, 0 -Iu( IU5

ca (v) =

D(v) =

dti \ui 1 0 0 0 0

0 -dti \u-21 0 0 0

0 0 -dti \u31 0 0

0 0 0 -dri 1W41 dr2 0

0 0 0 0 -dri \U5

0 0 0 0 0

g(n) =

mg - 3nr pgj sin (d)

mg - jnr pg ) cos (d) sin (ifi)

mg - ^nr pgj cos (0) cos ($)

zGmgcos (d) sin ($) zGmg sin (d)

cGu6 Ui + Znu

0 0 0 0 mu3

0 0 0 -mu3 0 0

0 0 0 mu2 -mui

0 mu3 -mu2 0 0 0

-mu3 0 mui 0 0 0

mu2 -mui 0 0 0

mu, 0 -Iu,

r1 r*6

Provided that rODIN = 0.31 m is the radius of ODIN, m = 125.0 kg is the mass of ODIN, zG = 0.05 m is the distance of the center of gravity from the geometric center, pV = 965 kg/m3 is the average density of the ODIN AUV, and p = 1000 kg/m3 is the density of fresh water, g = 9.81 m/s2. The hydrodynamic derivatives

are given by Xud = Yvd = Zwd = (2/3)npr30Dm, the translational quadratic damping factor dn = -248 N(s/m)2, the angular quadratic damping factor dr1 = -280 Ns2/m, and the angular linear damping factor dr2 = -230 Ns2/m. The thruster configuration matrix is defined as follows:

■s 0

fS -20

£ -20

fS -50

IS 0 a

fS -50

100 150 200 250

100 150 200 250

100 150 200 250

50 100 150 200 250

Time (s)

fS -20

fS -50

100 150 200 250

100 150 200 250

100 150 200 250

50 100 150 200 250 300

Time (s)

Figure 10: Thruster forces for Case 2 using proposed tracking control law. (a) Horizontal thrusters and (b) vertical thrusters.

s -s -s s 0 0 0 0

s s -s -s 0 0 0 0

0 0 0 0 -1 -1 -1 -1

0 0 0 0 to to -to -to

0 0 0 0 to -to -to to

Rz Rz 0 0 0 0

where s = sin((1/4)rc), R = 0.381 m denotes the distance from the center of the vehicle to the vertical thrusters center, and Rz = 0.508 m represents the radial distance from the center of the vehicle to the horizontal thruster's center.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


The authors would like to thank Offshore Geo-Surveys (OGS) Sdn. Bhd. for their invaluable help and advice on the underwater vehicle. This work was supported in part by the Ministry of Higher Education, Malaysia, and Universiti Teknologi Malaysia under Grant no. R.J130000.7823.4F101 and Q.J130000.2509.05H54.


[1] J. Yuh, "Design and control of autonomous underwater robots: a survey," Autonomous Robots, vol. 8, no. 1, pp. 7-24, 2000.

[2] G. Antonelli, Underwater Robots: Motion and Force Control of Vehicle-Manipulator Systems, Springer, Germany, Berlin, 2003.

[3] O.-E. Fjellstad and T. I. Fossen, "Position and attitude tracking of AUV's: a quaternion feedback approach," IEEE Journal of Oceanic Engineering, vol. 19, no. 4, pp. 512-518,1994.

[4] D. R. Yoerger and J.-J. E. Slotine, "Robust trajectory control of underwater vehicles," IEEE Journal of Oceanic Engineering, vol. 10, no. 4, pp. 462-470,1985.

[5] K. P. Venugopal, R. Sudhakar, and A. S. Pandya, "On-line learning control of autonomous underwater vehicles using feedforward neural networks," IEEE Journal of Oceanic Engineering, vol. 17, no. 4, pp. 308-319,1992.

[6] J. Yuh, "Learning control for underwater robotic vehicles," IEEE Control Systems Magazine, vol. 14, no. 2, pp. 39-46,1994.

[7] F. Wang, Y. Xu, L. Wan, and Y. Li, "Modified learning of TS fuzzy neural network control for autonomous underwater Vehicles," in Proceedings of the International Conference on Information Technology and Computer Science (ITCS '09), pp. 361-365, July 2009.

[8] X. Li, S. P. Hou, and C. C. Cheah, "Adaptive region tracking control for autonomous underwater vehicle," in Proceedings of the 11th International Conference on Control, Automation, Robotics and Vision (ICARCV '10), pp. 2129-2134, Singapore, December 2010.

[9] T. K. Podder and N. Sarkar, "Fault-tolerant control of an autonomous underwater vehicle under thruster redundancy," Robotics and Autonomous Systems, vol. 34, no. 1, pp. 39-52,2001.

[10] N. Sarkar, T. K. Podder, and G. Antonelli, "Fault-accommodating thruster force allocation of an AUV considering thruster redundancy and saturation," IEEE Transactions on Robotics and Automation, vol. 18, no. 2, pp. 223-233, 2002.

[11] D. Zhu, Q. Liu, and Y. Yang, "An active fault-tolerant control method ofunmanned underwater vehicles with continuous and uncertain faults," International Journal of Advanced Robotic Systems, vol. 5, no. 4, pp. 411-418, 2008.

[12] J.-K. Choi and H. Kondo, "On fault-tolerant control of a hovering AUV with four horizontal and two vertical thrusters," in Proceedings of the Conference on Oceans (OCEANS '10), Sydney, Australia, May 2010.

[13] Q. Liu, D. Zhu, and S. X. Yang, "Unmanned underwater vehicles fault identification and fault-tolerant control method based on FCA-CMAC neural networks, applied on an actuated vehicle," Journal of Intelligent and Robotic Systems: Theory and Applications, vol. 66, no. 4, pp. 463-475, 2012.

[14] M. Motamedi, M. T. Ahmadian, G. Vossoughi, S. M. Rezaei, and M. Zareinejad, "Adaptive sliding mode control of a piezo-actuated bilateral teleoperated micromanipulation system," Precision Engineering, vol. 35, no. 2, pp. 309-317, 2011.

[15] M. Zeinali and L. Notash, "Adaptive sliding mode control with uncertainty estimator for robot manipulators," Mechanism and Machine Theory, vol. 45, no. 1, pp. 80-90, 2010.

[16] M. Roopaei, B. R. Sahraei, and T.-C. Lin, "Adaptive sliding mode control in a novel class of chaotic systems," Communications in Nonlinear Science and Numerical Simulation, vol. 15, no. 12, pp. 4158-4170, 2010.

[17] The Society of Naval Architects and Marine Engineers (SNAME), "Nomenclature for treating the motion of a submerged body through a fluid," in Technical and Research Bulletin, pp. 1-15,1950.

[18] T. I. Fossen, Guidance and Control of Ocean Vehicles, John Wiley and Sons, New York, NY, USA, 1st edition, 1994.

[19] Z. H. Ismail and M. W. Dunnigan, "Tracking control scheme for an underwater vehicle-manipulator system with single and multiple sub-regions and sub-task objectives," IET Control Theory and Applications, vol. 5, no. 5, pp. 721-735, 2011.

[20] S. Soylu, B. J. Buckham, and R. P. Podhorodeski, "A chattering-free sliding-mode controller for underwater vehicles with fault-tolerant infinity-norm thrust allocation," Ocean Engineering, vol. 35, no. 16, pp. 1647-1659, 2008.

[21] J. J. E. Slotine and W. Li, Applied Nonlinear Control, Prentice Hall, 1991.

[22] T. F. Chan and R. V. Dubey, "Weighted least-norm solution based scheme for avoiding joint limits for redundant joint manipulators," IEEE Transactions on Robotics and Automation, vol. 11, no. 2, pp. 286-292, 1995.

[23] S. K. Choi, J. Yuh, and N. Keevil, "design of omni-directional underwater robotic vehicle," in Proceedings of the Conference on Oceans (OCEANS '10), pp. I192-I197, October 1993.

Copyright of Mathematical Problems in Engineering is the property of Hindawi Publishing Corporation and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.