# Analysis of Parameter Sensitivity Using Robust Design Techniques for a Flatfish Type Autonomous Underwater Vehicle

Authors: M. Santhakumar; T. Asokan; T. R. Sreeram
Journal: International Journal of Quality, Statistics, and Reliability
 (2009)

---

## Abstract

Hydrodynamic parameters play a major role in the dynamics and control of Autonomous Underwater Vehicles (AUVs). The performance of an AUV is dependent on the parameter variations and a proper understanding of these parametric influences is essential for the design, modeling, and control of high-performance AUVs. In this paper, the sensitivity of hydrodynamic parameters on the control of a flatfish type AUV is analyzed using robust design techniques such as Taguchi's design method and statistical analysis tools such as Pareto-ANOVA. Since the pitch angle of an AUV is one of the crucial variables in the control applications, the sensitivity analysis of pitch angle variation is studied here. Eight prominent hydrodynamic coefficients are considered in the analysis. The results show that there are two critical hydrodynamic parameters, that is, hydrodynamic force and hydrodynamic pitching moment in the heave direction that influence the performance of a flatfish type AUV. A near-optimal combination of the parameters was identified and the simulation results have shown the effectiveness of the method in reducing the pitch error. These findings are significant for the design modifications as well as controller design of AUVs.


---

## Body

## 1. Introduction

Autonomous Underwater Vehicles provide new alternatives for undersea exploration, relieving human divers from the risks and fatigue of working under constrained underwater environments.  These are intelligent robots deployed to carry out predefined underwater tasks without human intervention. Since they are autonomous in nature, dynamics, control, and navigation of these robots are of major concern to the ocean engineering researchers. Developing a highly stable and maneuverable AUV with good payload capability is of major focus in the underwater robotics research today. The dynamic performance and autonomous control of an AUV depend on many factors such as robot shape, weight, buoyancy, propulsion, sensor systems, and control systems. The geometrical shape of the robot determines the hydrodynamic and added mass coefficients and has a very large influence on the robot dynamics [1, 2].There are two major types of AUVs presently in use. The most common type is the torpedo-shaped AUV with a cylindrical hull [3], which is widely used in the design of AUVs. The flatfish type AUV has a geometric shape as shown in Figure 1, which is highly maneuverable and at the same time can carry higher payloads.  This exceptionally stable shape produces low drag to the body and is wide enough to accept a variety of sensors [4].Figure 1
Body-fixed frame and earth-fixed reference frame for AUV.The essential control requirements for any AUV are maneuvering in dive plane, depth control, and station keeping. While all these are critical requirements for an AUV, this paper focuses on the maneuvering in dive plane, especially the pitch angle variation during the forward motion of AUV. The design parameters that influence control requirements are the vehicle buoyancy, vertical distance between centre of gravity (CG) to centre of buoyancy (CB), and the hydrodynamic parameters. Due to safety reasons, most of the AUVs are designed to have some minimum positive buoyancy. Similarly, the vertical distance between CG and CB is always kept constant in order to ensure stability of the vehicle.  Therefore, for a given configuration of AUV, the hydrodynamic parameter variation is the most critical and hence the sensitivity of these parameters in the pitch motion of an AUV is analyzed using Taguchi’s design principles in this paper. The paper is organized in the following manner.  A short introduction to the modeling and simulation of AUV with open-loop and closed-loop controls is presented in Section2.  Parametric model of an AUV is presented in Section 3, followed by a short summary of robust design formulation and an introduction to Taguchi techniques (Section 4). The sensitivity analysis and numerical simulations are presented in Section 5 for a flatfish type AUV. Finally the results are summarized.


## 2. Modeling and Simulation

Underwater robots experience a range of forces while moving in the fluid medium and they are generally referred to as hydrodynamic forces. A detailed discussion on hydrodynamic forces on underwater robots can be found in [5]. Added mass, drag, and lift are the most common reaction forces acting on the AUV and these are related to linear and angular accelerations and velocities. The hydrodynamic parameters relate these forces to the corresponding robot parameters that are highly nonlinear. Many of these parameters are dependent on the shape of AUV. In the case of flatfish type AUV, these hydrodynamic parameters are not easily derived, owing to the complex nature of the geometry. Since the robot performance depends on these parameters, it is necessary to identify the critical parameters and control them to improve vehicle performance. To investigate the dynamic behavior of the AUV, a mathematical model of the robot has been developed. The dynamic motion of AUV can be described in a common way using six degrees of freedom (DOF) nonlinear equations (refer to Figure 1). The equations of motion of the AUV can be written as in [5, 6] and these equations show the various parameters which are involved in the dynamic motion of the robot and are used in the numerical simulation. The simulation method and its details are given in the proceeding section.Surge Motion:(1)(m-Xu̇)u̇-(Xu+Xu|u||u|)u-(Xv+Xv|v||v|)v-(Xw+Xw|w||w|)w-(mxG+Zq̇)q2-(mxG-Yṙ)r2+(Yv̇-m)vr+(m-Zẇ)wq+(W-B)sinθ=(Xδs+Xδs|δs||δs|)δs+(Xδr+Xδr|δr||δr|)δr+τprop.
Sway Motion:(2)(m-Yv̇)v̇-(Yv+Yv|v||v|)v-Yrr+(mxG-Yṙ)ṙ+mur+(mxG+Zq̇)pq+(m-Zẇ)wp-(W-B)cosθsinϕ=(Yδr+Yδr|δr||δr|)δr.
Heave Motion:(3)(m-Zẇ)ẇ-(Zw+Zw|w||w|)w-Zqq+(mxG-Zq̇)q̇-muq+(mxG-Yṙ)rp+(m-Yv̇)vp-(W-B)cosθcosϕ=(Zδs+Zδs|δs||δs|)δs.
Roll Motion:(4)(Ix-Kṗ)ṗ+(Iz-Iy)qr+(Yv̇-Zẇ)vw-Zq̇vq+(Mq̇-Nṙ)rq+Yṙrw-Kpp-zBBcosθsinϕ=0.
Pitch Motion:(5)(Iy-Mq̇)q̇-mxGẇ+mxGuq+(Ix-Iz+Nṙ)rp-mxGvp-(Mw+Mw|w||w|)w-Mqq-zBBsinθ+(xGW-xBB)cosθcosϕ=(Mδs+Mδs|δs||δs|)δs.
Yaw Motion:(6)(Iz-Nṙ)ṙ-(Nv+Nv|v||v|)v-Nrr+mxGur+(Iy-Ix-Mq̇)pq-mxGwp-(xGW-xBB)cosθsinϕ=(Nδr+Nδr|δr||δr|)δr,
where [xyzϕθψ] are positions and orientations with respect to inertial (fixed) frame, x is surge position, y is sway position, z is heave position, ϕ is  roll angle, θ is  pitch angle, and ψ is  yaw angle, [uvwpqr] are linear and angular velocities with respect to body (moving) frame, u is surge velocity, v is sway velocity, w is heave velocity, p is roll rate, q is pitch rate, r is yaw rate. u̇,v̇ and ẇ are linear accelerations in surge, sway, and heave directions, respectively. ṗ,q̇, and ṙ are angular accelerations in roll, pitch, and yaw directions, respectively.  [XYZKMN] are forces and moments (which include hydrodynamic damping and added mass) acting on AUV (such as Xu|u|,Xu,Xu̇), X is resultant force in surge direction, Y is resultant force in sway direction, Z is resultant force in heave direction, K is resultant moment in roll axis, M is  resultant moment in pitch axis, N is  resultant moment in yaw axis, m is mass of the AUV, W is weight of the AUV, B is buoyancy of the AUV, I is inertia terms of the AUV, δs is stern control plane deflection, δr is rudder plane deflection, CG=[xG,yG,zG]= Centre of Gravity, and CB=[xB,yB,zB]= Centre of Buoyancy.


### 2.1. Initial Study

In order to identify the sensitivity of these parameters, an experimental AUV that is being developed at the Indian Institute of Technology, Madras, India, is considered a test case. The dynamic model of the robot is developed using the Newton-Euler formulation [5] and is simulated using the experimental robot (positive buoyant) parameters and the experimentally determined hydrodynamic parameters. Simulation experiments were carried out with the help of MATLAB–Simulink software (step size is 0.01 s) and Runge–Kutta’s fourth-order model (RK4 solver is used) for the above dynamic equations. Repeated sequences of random signals (white noises) are considered the noise effects, which are introduced in the input stage and their amplitudes are limited by ±2 units. The results are shown in Figure 2. It shows a large steady state pitch angle during the forward motion, which is not acceptable from the point of view of control and navigation. This large steady state pitch angle results in the heave and sway motions also (refer to Figure 2), which is, again, not desirable. One option at this stage is to go for closed loop control of the AUV, which will reduce the pitch error. However, the control forces needed to achieve this will be very high and eventually lead to poor overall performance. Hence, it is necessary to identify the causes for pitch variation and rectify the problem through redesign of the AUV. Sensitivity analysis is one of the methods to identify the critical hydrodynamic parameters that contribute to the pitching motion. Robust design methods, developed for general design optimization, can be effectively used here to analyze the sensitivity of the AUV to various hydrodynamic parameters.AUV positions and orientations for constant forward speed (atU = 1.98 m/s).
(a)(b)(c)(d)(e)(f)


## 2.1. Initial Study

In order to identify the sensitivity of these parameters, an experimental AUV that is being developed at the Indian Institute of Technology, Madras, India, is considered a test case. The dynamic model of the robot is developed using the Newton-Euler formulation [5] and is simulated using the experimental robot (positive buoyant) parameters and the experimentally determined hydrodynamic parameters. Simulation experiments were carried out with the help of MATLAB–Simulink software (step size is 0.01 s) and Runge–Kutta’s fourth-order model (RK4 solver is used) for the above dynamic equations. Repeated sequences of random signals (white noises) are considered the noise effects, which are introduced in the input stage and their amplitudes are limited by ±2 units. The results are shown in Figure 2. It shows a large steady state pitch angle during the forward motion, which is not acceptable from the point of view of control and navigation. This large steady state pitch angle results in the heave and sway motions also (refer to Figure 2), which is, again, not desirable. One option at this stage is to go for closed loop control of the AUV, which will reduce the pitch error. However, the control forces needed to achieve this will be very high and eventually lead to poor overall performance. Hence, it is necessary to identify the causes for pitch variation and rectify the problem through redesign of the AUV. Sensitivity analysis is one of the methods to identify the critical hydrodynamic parameters that contribute to the pitching motion. Robust design methods, developed for general design optimization, can be effectively used here to analyze the sensitivity of the AUV to various hydrodynamic parameters.AUV positions and orientations for constant forward speed (atU = 1.98 m/s).
(a)(b)(c)(d)(e)(f)


## 3. AUV Parametric Model

A parametric model of the AUV for the pitch motion is developed to study the pitch angle variation and the parameter sensitivity.


### 3.1. AUV Model for Pitch

With reference to Figure1, the pitch motion is with respect to the Y axis, in the vertical plane. Assuming that (xB,yB,zB) and (xG,yG,zG) be the coordinates of centre of buoyancy and centre of gravity of the vehicle with respect to the body centre, the AUV model equations for pitch are derived as follows [5, 6]:(7)(ẇq̇żθ̇)=M-1[ZwZq-(mzGq+mu)-Xu̇u00Mw+Xu̇uMq+m(xGu+zGw)00cosθ0000100](wqzθ)+M-1((W-B)cosθ(xBB-xGW)cosθ+(zBB-zGW)sinθ-usinθ0)+M-1(ZδMδ00),
where inertia matrix is
(8)M=[m-Zẇ-mxG-Zq̇00-mxGIyy-Mq̇0000100001].
The terms Xu̇,Zẇ,Zq̇,Mq̇,Zw,Zq,Mw, and Mq in (7) are the hydrodynamic parameters (hydrodynamic damping and added mass forces and moments) affecting the pitch motion. Here, Xu̇,Zẇ, and Zq̇ are the added mass forces due to accelerations u̇,ẇ, and q̇, respectively, and Mq̇ is the hydrodynamic added mass moment due to acceleration q̇ along the pitch axis. Similarly Zw and Zq are the hydrodynamic damping forces along the Z axis due to velocities w and q, and Mw and Mq are the hydrodynamic damping moments in the pitch axis due to velocities w and q.


## 3.1. AUV Model for Pitch

With reference to Figure1, the pitch motion is with respect to the Y axis, in the vertical plane. Assuming that (xB,yB,zB) and (xG,yG,zG) be the coordinates of centre of buoyancy and centre of gravity of the vehicle with respect to the body centre, the AUV model equations for pitch are derived as follows [5, 6]:(7)(ẇq̇żθ̇)=M-1[ZwZq-(mzGq+mu)-Xu̇u00Mw+Xu̇uMq+m(xGu+zGw)00cosθ0000100](wqzθ)+M-1((W-B)cosθ(xBB-xGW)cosθ+(zBB-zGW)sinθ-usinθ0)+M-1(ZδMδ00),
where inertia matrix is
(8)M=[m-Zẇ-mxG-Zq̇00-mxGIyy-Mq̇0000100001].
The terms Xu̇,Zẇ,Zq̇,Mq̇,Zw,Zq,Mw, and Mq in (7) are the hydrodynamic parameters (hydrodynamic damping and added mass forces and moments) affecting the pitch motion. Here, Xu̇,Zẇ, and Zq̇ are the added mass forces due to accelerations u̇,ẇ, and q̇, respectively, and Mq̇ is the hydrodynamic added mass moment due to acceleration q̇ along the pitch axis. Similarly Zw and Zq are the hydrodynamic damping forces along the Z axis due to velocities w and q, and Mw and Mq are the hydrodynamic damping moments in the pitch axis due to velocities w and q.


## 4. Robust Design Formulation

Engineering design is increasingly becoming model based, in that its complexity calls for a mathematical model involving multiple quantities, some of which are to be decided by the designer with the purpose of meeting performance specifications, for example, the thrust that an underwater propeller must deliver at a given rpm, under given environment conditions such as ambient temperature, pressure, and so forth. The aim of robust design is to develop products whose performance remains within specifications in the presence of large variations in environment conditions.


### 4.1. Sensitivity Analysis

It is a procedure to determine the sensitivity of the outcomes of an alternative to changes in its parameters (as opposed to changes in the environment). If a small change in a parameter results in relatively large change in the outcomes, the outcomes are said to be sensitive to that parameter. This may mean that the parameter has to be determined very accurately or the alternative is to redesign such that the sensitivity of the parameter is low [6, 7]. Sensitivity Analysis can be used to determine(1)
the model resemblance with the system under study,(2)
the quality of model definition,(3)
factors that mostly contribute to the output variability,(4)
the region in the space of input factors for which the model variation is maximum,(5)
interactions between factors.Parameter sensitivity is usually performed as a series of tests in which the modeler sets different parameter values to see how a change in the parameter causes a change in the dynamic behavior of the system. By showing how the model behavior responds to change in parameter values, sensitivity analysis is a useful tool in model building as well as in model evaluation [8, 9].For sensitivity analysis, Taguchi’s robust design is more appropriate because it is a design and data analysis method and is more engineering oriented than science oriented. The distinct idea of Taguchi’s robust design that differs from the conventional experimental design is that of designing for the simultaneous modeling of both mean and variability [10].The concept of robust design has many aspects such as(i)
finding a set of conditions for design variables which are robust to noise,(ii)
achieving the smallest variation in a product’s function relative to a desired target value,(iii)
minimizing the number of experiments using orthogonal arrays and testing for confirmation.


### 4.2. Robust Design

Robust Design is a method, also called the Taguchi Method, pioneered by Byrne and Taguchi  [11] that greatly improves engineering productivity. Taguchi introduced the concept of parameter design as an effective means of improving the efficiency of a system which involves significant variability or “noise.” A system can perform its intended function at many settings/values of design parameters. The method requires the consideration of the effect of noise factors (environmental variation during the system’s usage, manufacturing variation of the system, and component deterioration) and the cost of failure in the field to ensure the performance requirements of the system. Indeed, it is the most powerful method available to reduce the variation, improve efficiency, and simultaneously reduce development interval [10]. Robust design technique can be applied to many aspects such as optimization, experimental design, sensitivity analysis, parameter estimation, model prediction, and so forth [7, 10].


### 4.3. Taguchi’s Approach

In setting up a framework for robust design, the classification of the quantities is at play in the design task as given below.(i)
Design variables (DVs) are those quantities to be decided by the designer with the purpose of meeting performance specifications under given conditions.(ii)
Design-environment parameters (DEPs) are the quantities over which the designer has no control and define the conditions of the environment under which the designed object will operate.(iii)
Performance functions (PFs) are quantities used to represent the performance of the design in terms of design variables and design-environment parameters.The responses at each setting of parameters are treated as a measure that would be indicative of not only the mean of some quality characteristic but also the variance of the same characteristic. The mean and the variance are combined into a single performance measure known as the signal-to-noise (S/N) ratio [10]. Taguchi classifies robust parameter design problems into different categories depending on the goal of the problem and for each category as follows.Smaller the better: the target value of y, that is, quality variable, is zero. In this situation, S/N ratio is defined as follows: (9)S/NRatio=-10log(1n∑i=1nyi2).Larger the better: the target value of y, that is, quality variable is infinite and S/N ratio is defined as follows:
(10)S/NRatio=-10log(1n∑i=1n1yi2).Nominal the best: the certain target value is given for y value. In this situation, S/N ratio is defined as follows:
(11)S/NRatio=-10log(∑i=1nyi2s2).
In this paper, smaller the better characteristic is used due to the requirement to keep the pitch angle variation to a minimum.Taguchi’s method uses an orthogonal array (OA) and analysis of mean to study the effects of parameters based on statistical analysis of experiments. An OA is a fractional factorial matrix which assures a balanced comparison of levels of any factor or interaction of factors. It is a matrix of numbers arranged in rows and columns where each row represents the level of the factors in each run, and each column represents a specific factor that can be changed from each run. The array is called orthogonal because all columns can be evaluated independently of one another. To compare performances of parameters, the statistical test known as the analysis of variance (ANOVA) is used. Further details and technical merits about robust parameter design can be found in the references of [9, 11–14]. The standard array representation of OA is depicted in Figure 3.  Standard OA is basically derived from linear graph of number of variables and levels; the ones most often used are L4 (23), L8 (27), L9 (34), L25 (56), and L49 (78) [7].Figure 3
Representation of standard orthogonal array.Once the suitable orthogonal array for experimentation is selected, then the experiments are conducted for the identified conditions and the results are analyzed using statistical methods, such as Pareto-ANOVA, ANOVA, or response curves to identify the optimal parameters for performance enhancement [7, 11, 15, 16].


## 4.1. Sensitivity Analysis

It is a procedure to determine the sensitivity of the outcomes of an alternative to changes in its parameters (as opposed to changes in the environment). If a small change in a parameter results in relatively large change in the outcomes, the outcomes are said to be sensitive to that parameter. This may mean that the parameter has to be determined very accurately or the alternative is to redesign such that the sensitivity of the parameter is low [6, 7]. Sensitivity Analysis can be used to determine(1)
the model resemblance with the system under study,(2)
the quality of model definition,(3)
factors that mostly contribute to the output variability,(4)
the region in the space of input factors for which the model variation is maximum,(5)
interactions between factors.Parameter sensitivity is usually performed as a series of tests in which the modeler sets different parameter values to see how a change in the parameter causes a change in the dynamic behavior of the system. By showing how the model behavior responds to change in parameter values, sensitivity analysis is a useful tool in model building as well as in model evaluation [8, 9].For sensitivity analysis, Taguchi’s robust design is more appropriate because it is a design and data analysis method and is more engineering oriented than science oriented. The distinct idea of Taguchi’s robust design that differs from the conventional experimental design is that of designing for the simultaneous modeling of both mean and variability [10].The concept of robust design has many aspects such as(i)
finding a set of conditions for design variables which are robust to noise,(ii)
achieving the smallest variation in a product’s function relative to a desired target value,(iii)
minimizing the number of experiments using orthogonal arrays and testing for confirmation.


## 4.2. Robust Design

Robust Design is a method, also called the Taguchi Method, pioneered by Byrne and Taguchi  [11] that greatly improves engineering productivity. Taguchi introduced the concept of parameter design as an effective means of improving the efficiency of a system which involves significant variability or “noise.” A system can perform its intended function at many settings/values of design parameters. The method requires the consideration of the effect of noise factors (environmental variation during the system’s usage, manufacturing variation of the system, and component deterioration) and the cost of failure in the field to ensure the performance requirements of the system. Indeed, it is the most powerful method available to reduce the variation, improve efficiency, and simultaneously reduce development interval [10]. Robust design technique can be applied to many aspects such as optimization, experimental design, sensitivity analysis, parameter estimation, model prediction, and so forth [7, 10].


## 4.3. Taguchi’s Approach

In setting up a framework for robust design, the classification of the quantities is at play in the design task as given below.(i)
Design variables (DVs) are those quantities to be decided by the designer with the purpose of meeting performance specifications under given conditions.(ii)
Design-environment parameters (DEPs) are the quantities over which the designer has no control and define the conditions of the environment under which the designed object will operate.(iii)
Performance functions (PFs) are quantities used to represent the performance of the design in terms of design variables and design-environment parameters.The responses at each setting of parameters are treated as a measure that would be indicative of not only the mean of some quality characteristic but also the variance of the same characteristic. The mean and the variance are combined into a single performance measure known as the signal-to-noise (S/N) ratio [10]. Taguchi classifies robust parameter design problems into different categories depending on the goal of the problem and for each category as follows.Smaller the better: the target value of y, that is, quality variable, is zero. In this situation, S/N ratio is defined as follows: (9)S/NRatio=-10log(1n∑i=1nyi2).Larger the better: the target value of y, that is, quality variable is infinite and S/N ratio is defined as follows:
(10)S/NRatio=-10log(1n∑i=1n1yi2).Nominal the best: the certain target value is given for y value. In this situation, S/N ratio is defined as follows:
(11)S/NRatio=-10log(∑i=1nyi2s2).
In this paper, smaller the better characteristic is used due to the requirement to keep the pitch angle variation to a minimum.Taguchi’s method uses an orthogonal array (OA) and analysis of mean to study the effects of parameters based on statistical analysis of experiments. An OA is a fractional factorial matrix which assures a balanced comparison of levels of any factor or interaction of factors. It is a matrix of numbers arranged in rows and columns where each row represents the level of the factors in each run, and each column represents a specific factor that can be changed from each run. The array is called orthogonal because all columns can be evaluated independently of one another. To compare performances of parameters, the statistical test known as the analysis of variance (ANOVA) is used. Further details and technical merits about robust parameter design can be found in the references of [9, 11–14]. The standard array representation of OA is depicted in Figure 3.  Standard OA is basically derived from linear graph of number of variables and levels; the ones most often used are L4 (23), L8 (27), L9 (34), L25 (56), and L49 (78) [7].Figure 3
Representation of standard orthogonal array.Once the suitable orthogonal array for experimentation is selected, then the experiments are conducted for the identified conditions and the results are analyzed using statistical methods, such as Pareto-ANOVA, ANOVA, or response curves to identify the optimal parameters for performance enhancement [7, 11, 15, 16].


## 5. Numerical Simulation

As explained in Sections2 and 3, the objective of the work is to identify the critical hydrodynamic parameters influencing the pitch angle (performance function) variation of the AUV. As explained in Section 2.1, there are eight hydrodynamic parameters (design variables) that influence the pitching moment of the robot (AUV). As these parameters are functions of forward speed of the vehicle, this factor is used as the noise factor (design environment parameter) in the analysis. This noise factor (forward velocity of the vehicle) is changed by providing different thrust values in the simulation such as 25%, 50%, and 75% of total thrust available from the actuators. The details of the analysis and the results are presented in the following sections.


### 5.1. Experimental Layout

The first step in the analysis is to choose the orthogonal array for trials. For selecting the orthogonal array, the main criterion is the required number of trials [7].Mathematically, this is given as(12)N=A(L-1),
where N= Number of experiments,  A= Number of variables, and L= Number of levels.It was decided to conduct trials using the eight variables at seven levels (more numbers of levels means, more number of experiments, and better accuracy of results), that is,  10%, 25%, 40%, 55%, 70%, 85%, and 100% of the actual values of the variables. So, the number of trials required was found to be 49. The smallest standard orthogonal array that matches with this requirement isL49 (78), which will give a test matrix of 49 different configurations with various combinations of parameters and its physical values (eight parameters at seven levels). Three levels of noise variable (vehicle speed) are considered in the analysis and the test matrix used is shown in Table 1.Table 1
Design of Experiments for Pitch Response (L49 (78) OA).
Exp. no.Xu̇ZẇZq̇Mq̇ZwZqMwMqPitch angle (deg) forU =S/N ratio (dB)(%)(%)(%)(%)(%)(%)(%)(%)1.13 (m/s)1.59 (m/s)1.98 (m/s)11010101010101010-8.82-10.03-10.9330.03210252540557085100-33.43-34.46-34.9719.30310404055708510025-32.63-33.33-33.6719.57410555570851002540-7.82-7.98-8.0631.99510707085100254055-10.73-10.88-10.9729.28610858510025405570-35.01-37.61-39.0818.577101001002540557085-34.31-35.91-36.7318.9682510252525252525-17.12-18.56-19.4824.70925255585401010070-5.42-5.76-5.9634.861025408540101007055-53.19-58.33-60.9714.791125554010100705585-14.53-14.73-14.8326.661225701010070558540-28.52-29.16-29.4920.731325851007055854010-17.15-17.76-18.1125.0514251007055854010100-3.25-3.32-3.3539.61154010404040404040-21.04-22.14-22.7823.151640258570100551025-2.86-2.92-2.9540.721740407010055102585-10.93-11.33-11.5728.951840551005510258570-61.21-63.12-65.3613.981940705510258570100-42.29-45.13-46.6116.992040851025857010055-28.45-28.92-29.1520.8021401002585701005510-19.36-19.83-20.0924.08225510555555555555-23.01-23.79-24.2322.512355254010085257010-89.98-89.98-89.9810.922455401008525701040-6.82-7.45-7.8732.622555558525701040100-14.35-14.71-14.8626.692655702570104010085-65.01-67.02-69.2413.462755857010401008525-39.9-41.59-42.4317.6828551001040100852570-6.83-6.93-6.9933.20297010707070707070-24.09-24.66-24.9622.193070251055251004085-26.58-28.72-30.0120.913170405525100408510-89.98-89.98-89.9810.923270552510040851055-5.42-5.76-5.9634.863370701004085105525-16.63-16.93-17.0925.453470854085105525100-22.44-25.36-27.4121.9935701008510552510040-37.96-38.11-39.2318.30368510858585858585-24.73-25.15-25.3722.013785251001070402555-9.13-9.38-9.5130.593885401070402555100-28.02-29.41-30.1820.693985557040255510010-89.98-89.98-89.9810.924085704025551001070-4.45-4.64-4.7536.724185852555100107040-18.22-18.46-18.5924.6942851005510010704025-34.53-38.84-41.5918.314310010100100100100100100-25.13-25.46-25.6221.904410025702510855540-44.8-49.85-52.7616.154510040251085554070-12.3-12.53-12.6528.074610055108555407025-28.45-29.38-29.8620.684710070855540702510-13.46-14.2-14.6427.014810085554070251085-3.76-3.87-3.9438.2749100100407025108555-48.52-49.12-51.1216.09


## 5.1. Experimental Layout

The first step in the analysis is to choose the orthogonal array for trials. For selecting the orthogonal array, the main criterion is the required number of trials [7].Mathematically, this is given as(12)N=A(L-1),
where N= Number of experiments,  A= Number of variables, and L= Number of levels.It was decided to conduct trials using the eight variables at seven levels (more numbers of levels means, more number of experiments, and better accuracy of results), that is,  10%, 25%, 40%, 55%, 70%, 85%, and 100% of the actual values of the variables. So, the number of trials required was found to be 49. The smallest standard orthogonal array that matches with this requirement isL49 (78), which will give a test matrix of 49 different configurations with various combinations of parameters and its physical values (eight parameters at seven levels). Three levels of noise variable (vehicle speed) are considered in the analysis and the test matrix used is shown in Table 1.Table 1
Design of Experiments for Pitch Response (L49 (78) OA).
Exp. no.Xu̇ZẇZq̇Mq̇ZwZqMwMqPitch angle (deg) forU =S/N ratio (dB)(%)(%)(%)(%)(%)(%)(%)(%)1.13 (m/s)1.59 (m/s)1.98 (m/s)11010101010101010-8.82-10.03-10.9330.03210252540557085100-33.43-34.46-34.9719.30310404055708510025-32.63-33.33-33.6719.57410555570851002540-7.82-7.98-8.0631.99510707085100254055-10.73-10.88-10.9729.28610858510025405570-35.01-37.61-39.0818.577101001002540557085-34.31-35.91-36.7318.9682510252525252525-17.12-18.56-19.4824.70925255585401010070-5.42-5.76-5.9634.861025408540101007055-53.19-58.33-60.9714.791125554010100705585-14.53-14.73-14.8326.661225701010070558540-28.52-29.16-29.4920.731325851007055854010-17.15-17.76-18.1125.0514251007055854010100-3.25-3.32-3.3539.61154010404040404040-21.04-22.14-22.7823.151640258570100551025-2.86-2.92-2.9540.721740407010055102585-10.93-11.33-11.5728.951840551005510258570-61.21-63.12-65.3613.981940705510258570100-42.29-45.13-46.6116.992040851025857010055-28.45-28.92-29.1520.8021401002585701005510-19.36-19.83-20.0924.08225510555555555555-23.01-23.79-24.2322.512355254010085257010-89.98-89.98-89.9810.922455401008525701040-6.82-7.45-7.8732.622555558525701040100-14.35-14.71-14.8626.692655702570104010085-65.01-67.02-69.2413.462755857010401008525-39.9-41.59-42.4317.6828551001040100852570-6.83-6.93-6.9933.20297010707070707070-24.09-24.66-24.9622.193070251055251004085-26.58-28.72-30.0120.913170405525100408510-89.98-89.98-89.9810.923270552510040851055-5.42-5.76-5.9634.863370701004085105525-16.63-16.93-17.0925.453470854085105525100-22.44-25.36-27.4121.9935701008510552510040-37.96-38.11-39.2318.30368510858585858585-24.73-25.15-25.3722.013785251001070402555-9.13-9.38-9.5130.593885401070402555100-28.02-29.41-30.1820.693985557040255510010-89.98-89.98-89.9810.924085704025551001070-4.45-4.64-4.7536.724185852555100107040-18.22-18.46-18.5924.6942851005510010704025-34.53-38.84-41.5918.314310010100100100100100100-25.13-25.46-25.6221.904410025702510855540-44.8-49.85-52.7616.154510040251085554070-12.3-12.53-12.6528.074610055108555407025-28.45-29.38-29.8620.684710070855540702510-13.46-14.2-14.6427.014810085554070251085-3.76-3.87-3.9438.2749100100407025108555-48.52-49.12-51.1216.09


## 6. Results and Discussions

The simulation trials were conducted for the states mentioned in the OA. For each experiment, the pitch angle was recorded (refer to Table1) and the SNR (S/N ratio) (using smaller the better characteristics) was calculated. For example, the SNR for experiment 1 was calculated as(13)SNR=-10log(1n∑i=1nyi2)=-10log(13[(-8.82)2+(-10.03)2+(-10.93)2])+50=30.03dB.
(50 dB is added to avoid negative values for S/N ratio, to make comparisons easy).The plot of SNR for various parameters is shown in Figure4, which indicates a larger variation for Mw and Zw and higher sensitivity of these parameters.Figure 4
Response curves.In order to identify the optimal parameter combination and the contribution ratio of each parameter, statistical analysis of the results was carried out using Pareto-ANOVA analysis.For Pareto-ANOVA analysis, the sum of S/N ratios, sum of squares of differences, and contribution ratio were calculated as follows:(14)Sij=∑(SNR)ij,SSi=∑j=27(Si1-Sij)2,CRi=SSi∑SSix100,
where Sij = Sum of SNR of ith parameter at jth level, (SNR)ij = S/N ratio of ith parameter at jth level, SSi = Sum of squares of differences for ith parameter, and CRi = Percentage of contribution ratio.For example, the calculation of above values forMw is shown as follows:(15)SumoftheS/NratioforLevel1ofMw=30.03+39.61+40.72+32.62+34.86+36.72+38.27=252.83dB.(16)SumofsquaresofDifferencesofMw=(252.83-198.43)2+(252.83-171.46)2+(252.83-154.11)2+(252.83-129.22)2+(252.83-120.7)2+(252.83-139.82)2+(198.43-171.46)2+(198.43-154.11)2+(198.43-129.22)2+(198.43-120.7)2+(198.43-139.82)2+(171.46-154.11)2+(171.46-129.22)2+(171.46-120.7)2+(171.46-139.82)2+(154.11-129.22)2+(154.11-120.7)2+(154.11-139.82)2+(129.22-120.7)2+(129.22-139.82)2+(120.7-139.82)2=89949.91dB2.(17)Contributionratioin%ofMW=89949.91/(4489.1+1270.8+1406.7+4677.6+21419.2+4335+89949.9+8539.7)×100%=66.10%.
Table 2 shows the Pareto-ANOVA analysis details. It is found that the contribution ratio of Mw and Zw is significant compared to those of other parameters, showing higher sensitivity of these parameters. Also, this analysis gives the optimum combination of parameters to reduce the sensitivity.  For each parameter, the level which gives the highest SNR is the optimum level.Table 2
Pareto-ANOVA analysis.
ParametersSum of S/N ratio,SijSum of squares of differencesContribution ratio in %level 1level 2level 3level 4level 5level 6level 7Xu̇167.70186.40168.68157.08154.62163.92168.184489.0943.30Zẇ166.50173.45155.62165.76169.65167.06168.551270.7770.93Zq̇167.05169.15155.09173.84164.79168.10168.551406.6521.03Mq̇168.32154.93165.08168.29170.19185.53154.244677.6023.44Zw128.71140.81177.19171.52182.14178.84187.3721419.2315.74Zq186.76156.14156.99163.89166.88167.84168.074335.0053.19Mw252.83198.43171.46154.11129.22120.70139.8289949.9166.10Mq138.93167.11167.65168.93187.58169.23167.168539.6926.28From Table2, it can be seen that the optimum levels of parameters Xu̇,Zẇ,Zq̇,Mq̇,Zw,Zq,Mw, and Mq are 2, 2, 4, 6, 7, 1, 1, and 5, respectively. The corresponding physical values are 25%, 25%, 55%, 85%, 100%, 10%, 10%, and 70% of their original values.In order to show the effectiveness of this analysis in identifying optimum parameters to reduce the pitch angle variations of the AUV, simulations were carried out using the optimal values of the hydrodynamic parameters. The simulation results are summarized in Table3 and results are shown in Figures 5 and 6. It was observed that the pitch angle error has significantly been reduced by the use of optimized values. As a result, the errors in heave, surge, sway, and yaw were also reduced. Similarly, the closed loop control results show improved transient and steady state response (Figure 6). Further fine tuning of the parameter values may still reduce the errors to almost zero.Table 3
Comparative results.
ConditionXu̇ZẇZq̇Mq̇ZwZqMwMqPitch angle (deg) forU =(%)(%)(%)(%)(%)(%)(%)(%)1.13 (m/s)1.59 (m/s)Original values100100100100100100100100-25.13-25.46-25.62Optimized values25255585100101070-2.86-2.92-2.95Comparative results of AUV positions and orientations for open loop control.
(a)(b)(c)(d)(e)(f)Comparative results of positions and orientations for closed loop control.
(a)(b)(c)(d)(e)(f)The significance of the analysis is in the redesign of the AUV based on the sensitivity analysis. Once the most sensitive parameters were identified, designer can focus his attention on redesigning the physical features of the AUV to achieve the desired values for these parameters. For example, additional control planes in the aft or increased control plane areas will reduce the hydrodynamic moment (Mw) in the heave motion.


## 7. Conclusion

A systematic study on the sensitivity of various hydrodynamic parameters on the performance of an underwater robot was presented. Using statistical design techniques, critical parameters affecting the diving plane motion of the robot were identified and optimal combination of parameters to reduce the pitch angle error was determined. Simulation results have shown the effectiveness of the analysis in improving the dynamic performance. Further studies may include the robot motion in a three-dimensional plane and the sensitivity of environmental disturbance parameters as noise variables. Experimental validation of above findings will be taken up as soon as the prototype AUV is ready for trials.



---

*Source: 670340-2010-01-20.xml*