Faculty of Mechanical Engineering, Curtin University, 98009 Miri, Sarawak, Malaysia
Received Date: February 17, 2017; Accepted Date: March 22, 2017; Published Date: March 28, 2017
Citation: Yi DSL, Al-Qrimli HF (2017) Identification of Hydrodynamics Coefficient of Underwater Vehicle Using Free Decay Pendulum Method. J Powder Metall Min 6: 161. doi:10.4172/2168-9806.1000161
Copyright: © 2017 Yi DSL, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Visit for more related articles at Journal of Powder Metallurgy & Mining
In this research paper, free decay method is proposed to identify the hydrodynamics coefficient of ROVs. This method is preferable for smaller-scale testing as it is a smaller and simpler method. In order to carry the experiment, an experimental rig is designed and constructed. The experimental rig is used to create an underwater environment which is important in the identification of the hydrodynamics coefficient. First, the equations of motion are formulated that can be used to describe the motion of pendulum given a set of parameters which included measured parameters such as mass and length of pendulum and also the unknown hydrodynamic coefficient which are added mass and the damping coefficients. Then, experiments are conducted where the displacement of the pendulum is recorded over time by using camera. The desired parameters are found by finding the combination of parameters such that the mathematical model fits the experimental data. Through this research project, it was found that the results obtained from this method are not reliable as claimed by a few authors. Although the data processing procedure was a success, the results obtained vary widely with different initial condition. This research project shows that in the current state, this method is not a good method in evaluating the hydrodynamic coefficient of the ROVs. Thus, an alternative method should be proposed or improvements should be made to produce a better result.
Computational dynamics; ROVs; CAD design
In the present study, Remotely Operated Vehicles (ROVs) are underwater unmanned robots that have been widely used in the ocean activities. ROVs are able to carry out tasks such as navigation, control and guidance for the deep ocean exploration which are beyond the capability of human beings. Thus, the identification of hydrodynamics coefficient of ROVs becomes a great concern . Haidar F AL-Qrimli used many methods to identify the hydrodynamic parameters of ROVs. One of the most common method is Planar Motion Mechanism (PMM)/tow tank test. This method was introduced by authors [2-4] as a routine dynamics testing. PMM is capable of measuring the forces and moments in 6DOF which allows a complete identification. However, it is highly time-consuming and expensive to perform the test. According to Morrision and Yoerger , system identification method is also recommended to be used in the study of underwater vehicles. System identification method is the method of using onboard sensors and thrusters in the process of identification. Although there are many methods that can be used to identify the coefficient, free decay pendulum method is proposed to identify the hydrodynamics coefficient of ROVs because it is claimed by authors [2,3,6] that it is a repeatable and reliable method. As presented by Tan et al. , drag coefficient of any given models can be estimated by using fourth order regression method. By using CFD modelling, the model was subjected to pre-conditioning phase. The simulation was then run to obtain the drag profile in real-time. In their studies, it stated that the accuracy of the model is based on the frequency of trial runs and also the efficiency of the model’s sensors. It also stated that basic model of a cylinder was used to exhibit the accuracy and robustness of the proposed model. Thus, there is no proper experimental data to prove the validation of the data obtained. Further investigation such as identifying other fitting methods was also needed as stated in this paper to enable an accurate prediction. The aim of this project is to determine the hydrodynamic coefficient of ROV model through the free decay pendulum method.
The equations of motion are formulated from Newton’s second law. As the equation of motion can be expressed as: ƩF=ma, as it can be seen in Figure 1a.
Where ƩF is the total sum of forces on the body, ‘m’ is the mass and ‘a’ is the acceleration. In the pendulum experiment the considered forces include: hydrodynamic drag, buoyancy and weight. In this section the individual forces are discussed in detail leading to the formulation of the final equation of motion.
From the free body diagram,
The equation of drag force:
Substituting equation (2) and (3) into (1), the total equation of motion:
The free decay pendulum experiment is a method of determining the hydrodynamic parameters of a submerged body. First the equations of motion are formulated that can describe the motion of the pendulum given a set of parameters which include measured parameters such as mass or pendulum length and also the unknown hydrodynamic coefficients which are added mass and the damping coefficients. Experiments are conducted where the displacement of the pendulum is recorded over time. The desired coefficients are found by finding the combination of parameters such that the mathematical model fits the experimental data.
Experimental rig is one of the most important phases in this project. As the rig wasn’t an available facility to have due to the cost, a test tank was required to conduct the experiments. The tank is designed to be portable to allow us to move the tank easily from place to place. The external tank structure is made of joinable steel tubules with a PVC coated polyester fabric to contain the water. It is therefore disassemblable. The pendulum is mounted on a beam that spans the distance as it can be seen in Figure 1b.
During the experiment, the pendulum is lifted to the starting position. The camera is started to begin the recording and at the same time the pendulum is released to decay in the water. After a short amount of time, the camera is stopped as the pendulum comes to at rest. The video recorded is processed by using image processing method to obtain the data and least squares fitting is applied to obtain hydrodynamic coefficients. The ROV can be seen inside the rig when its fill of water in Figure 2 using a GoPro camera to record the experiment. The images of the motion of pendulum are required for analyzing the hydrodynamics coefficient (Figure 2).
The experiment is carried out and a camera is used to record the video. Images of the video are obtained by splitting the video in frames. After that, the images are processed by using MATLAB program. In this experiment, the motion of the pendulum is represented by the coordinates of a dot. In order to find the coordinate of dots, circle detection is required. By using MATLAB software, image processing method is used to process the image of motion. First, the image is read by using loop function to repeat the image reading process. Next, the contrast of the image is changed to enable the program to identify the circles easier. The image is then cropped to minimize the error during image processing. The circle is detected by using function in MATLAB. Angle is determined by using MATLAB Atan2 function. The angle obtained is then converted into displacement by using Equation 3. Lastly, regression method is applied to compare the theoretical and experimental data. Least squares regression is a technique to match analytical models with experimental data. In order to do this, the algorithm estimates the combination of parameters that reduce the error between the data set and the model in the least squares sense. In this experiment, this method is used to compute the estimated added mass and drag coefficients.
From the total equation of motion there are few parameters which can be calculated or measured without carrying out the experiment. These known parameters are simplified in Table 1. This volume corresponds to a buoyancy force of 24.43N. Buoyancy acts in the opposite direction of the weight so their effects may be combined to simplify analysis. The experiment is carried out multiple times with varying initial displacements in both directions to have a better analysis of the experiment. All the data obtained are analyzed and shown in this section (Table 1).
|Volume (V)||0.0025 m3|
|Density of water (ρ)||1000 kg/m3|
|Gravitational force (g)||9.81 m/s2|
|Length of pendulum (L)||0.76 m|
|Frontal Area (A)||0.0322 m2|
|Buoyancy (B)||24.43 N|
Table 1: Summary of known parameters and its values.
In order to ease the experiment, few assumptions were made such that the model is rigid and full submerged once in water. Water is assumed to be an ideal fluid that is incompressible, in viscid and irrational. The rod is assumed to be a point mass and frictional forces of the rod and pivot are negligible. The experiment results data is carried out multiple times with varying initial displacements in both directions to have a better analysis of the experiment. All the data obtained are analyzed and shown in this section. The result obtained by using initial position of 75°, the video of the motion is analyzed using Matlab software is shown in Figure 3. By using initial position of 20°, the video of the motion is analyzed using Matlab software. The result obtained is shown in Figure 4. In addition Figures 5-7 shows the initial positions of -15°, 5° and position of approximately -5° respectively. From observing the graphs obtained, it can be concluded that the data with larger initial position do not fit well with the model. This can also be proven by the values of RMS error calculated which also show that the data with larger initial position have higher value of error. However, the RMS error values, shown in Table 2, for all the data are relatively small which only range from 0.005 to 0.001. This shows that the data extracted from the free decay pendulum method can actually fit well with the model derived (Table 2 and Figures 3-7).
|Initial position (°)||RMS error|
Table 2: RMS error values.
From Table 3, the added mass obtained for initial position 75, 20, -15, 5 and -5 are 4.2577 × 10-14 kg, 19.266 kg, 3.0988 kg, 15.098 kg and 19.169 kg respectively while the drag coefficient calculated are 0.457, 2.999, 6.897 × 10-13, 6.903 × 10-13 and 6.896 × 10-13 respectively. By comparing these values, we can observe that the estimated values are not consistent and vary widely depending on the starting angle. The tests run at 5 and -5 degrees show the most consistency, this indicates that the output is highly sensitive to initial conditions (Table 3).
|Initial position (°)||Ma (kg)||CdA||Cd|
|75||4.2577 × 10-14||1.4719 × 10-2||0.457|
|20||19.266||9.6585 × 10-2||2.999|
|-15||3.0988||2.2208 × 10-14||6.897 × 10-13|
|5||15.098||2.2228 × 10-14||6.903 x10-13|
|-5||19.169||2.2204 × 10-14||6.896 x 10-13|
Table 3: Summary of experimental data.
This research project uses a free decay test to identify the added mass and the drag coefficient of the openROV robot. The data processing procedure was a success as it was able to extract displacement-time data from a video recording. Through the analysis of the experimental data, curve fitting routine was proven successful as it was able to match the data with the model well. This is proven by visual inspection of the graphs and also the low RMS values obtained. However, the results obtained shows that the parameter output is highly sensitive to initial conditions as a result, the estimated parameters vary widely. Thus, from this research project, it shows that the free decay method is not a good method in finding the hydrodynamic coefficient and added mass. The recommendation would be by improving the analytical model of the system which may include acceleration effects such as Basset force, time varying effects such as flow development and also others such as friction in bearings and pendulum arm.
I would like to express the deepest appreciation to Dr. Haidar Al-Qrimli and Mr. Zachary MacChesney for guiding and helping me throughout this research project. Besides, I would like to thank to the lab assistant for allowing me to conduct my experiment in the lab.