Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 527315, 12 pages

http://dx.doi.org/10.1155/2014/527315

Research Article

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

Zool H. Ismail,¹ Ahmad ’A. Faudzi,¹ and Matthew W. Dunnigan²

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

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

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.

Abstract.

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 task space, 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 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.

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

2.1. Kinematic Model

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

where and denote the position and the orientation of the vehicle, respectively, expressed in the inertial-fixed frame. and are the transformation matrices expressed in terms of the Euler angles. The linear and angular velocity vectors, and , 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 ; the underwater vehicle dynamic equation can be expressed in closed form as

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

where is the thruster configuration matrix and is the desired vector of thruster forces. The matrix 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 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 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 is symmetric and positive definite such that .

Property 2. is the skew-symmetric matrix such that .

Property 3. The hydrodynamic damping matrix is positive definite; that is, .

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

where is a known regression matrix; is the total number of physical parameters. It is assumed that if the arguments of are bounded, then 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 are the continuous first partial derivatives of the dynamic region, is the time-varying reference point inside the geometric shape, and is a time-varying and nonsingular scaling factor. It is assumed that and 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:

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

where , and 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

where is a positive scalar.

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 , where and is a regional bound, expressed in 2D task space. Then, the associated potential energy can be illustrated in Figure 2(a) where is lower bounded by zero. Next, invoking the following linear relationship to obtain the required forces of two onboard thrusters:

where is the net thruster forces, represents a force constant, is the -norm, and is the displacement in a planar plane. If 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 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.

(a)

(b)

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

(a)

(b)

Figure 3: An illustration of the required energy in 2D Cartesian space when is arbitrarily small: (a) potential energy and (b) kinetic energy.

Now, differentiating (8) with respect to gives

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

From the definition of the Cartesian region function, note that when the AUV is outside the Cartesian region. Once the vehicle reaches the region, this term becomes zero.

Next, a vector that is useful is defined as

where is a positive constant. The matrix represents the inverse of the Jacobian matrix, is the time derivative of , and 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

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

where

and the derivative of in (12) is given as

where , , and 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:

Replacing with leads to

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

where is a nonsquare thruster control matrix. To find the -norm solution for the desired thruster force , the weighted pseudoinverse matrix is utilized as a generalized inverse matrix of :

where ; 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 . In the case of no fault in any of the thrusters, the weighted matrix is an identity matrix. Hence, is a Moore-Penrose pseudoinverse matrix and will be minimized. Based on the error system development, the -norm solution for thruster force in (20), and the subsequent stability analysis, the proposed control law for AUV can be formulated as

where is positive constant matrix.

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

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

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

Simplifying (23) leads to

or

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

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

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

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 thrusters in (25), the proposed adaptive control law (21) and parameter update law (22) guarantee the convergence of into a dynamic region in the sense that and , as .

Proof. In view of the fact that is uniformly positive definite, in (26) is positive definite in , . Since and , is bounded. Hence, , , and are bounded. The boundedness of ensures the boundedness of . Hence, in (11) and its derivative are bounded. The boundedness of guarantees the boundedness of since and is a trigonometric function. The boundedness of and leads to the boundedness of and in (19) and (20), respectively. Moreover, in (16) remains a bounded signal by employing standard signal chasing arguments. From the closed-loop equation, it can be concluded that in (17) is also bounded using the fact that the kinematic and dynamic terms denoted by , , , and are bounded for all possible and . Thus, is bounded since . Differentiating (28) with respect to time gives

Hence, is bounded since , , , and are all bounded. Therefore, is uniformly continuous. Applying Barbalat’s lemma [21] gives which also indicates that and . From (11), indicates that is approaching zero. Therefore, 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]

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

Note that the second term on the right hand side of (31) is obtained from differentiation of the thruster force variable with respect to the thrust force of the th 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 of linear and quadratic damping whose effects are noncoupled 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:

where the elements of 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:

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 m. Moreover, the vehicle is initialized at the same position m while the initial values for attitude are set to be degrees. The attitude is kept constant for the entire motion. The control gains are set to the following:

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.

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

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

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

(a)

(b)

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

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 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.

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

Conventional tracking controller
Proposed control scheme

Horizontal thrusters, (N)
487.42
483.29

Vertical thrusters, (N)
229.53
218.43

Total input, (N)
716.95
701.72

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

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

(a)

(b)

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

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.

Appendix

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]

where the subscripts and 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

where are the moments of inertia about the principle axes.

Consider

Provided that is the radius of ODIN, is the mass of ODIN, is the distance of the center of gravity from the geometric center, is the average density of the ODIN AUV, and is the density of fresh water, . The hydrodynamic derivatives are given by , the translational quadratic damping factor , the angular quadratic damping factor , and the angular linear damping factor . The thruster configuration matrix is defined as follows:

where m denotes the distance from the center of the vehicle to the vertical thruster’s center, and 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.

Acknowledgments

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.

References

J. Yuh, “Design and control of autonomous underwater robots: a survey,” Autonomous Robots, vol. 8, no. 1, pp. 7–24, 2000.

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

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.

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.

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.

J. Yuh, “Learning control for underwater robotic vehicles,” IEEE Control Systems Magazine, vol. 14, no. 2, pp. 39–46, 1994.

F. Wang, Y. Xu, L. Wan, and Y. Li, “Modified learning of T-S 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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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

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.

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.

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

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.

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.