Next Article in Journal
Experimental Study on Extreme Hydrodynamic Loading on Pipelines Part 2: Induced Force Analysis
Next Article in Special Issue
Received Signal Strength Indication (RSSI) of 2.4 GHz and 5 GHz Wireless Local Area Network Systems Projected over Land and Sea for Near-Shore Maritime Robot Operations
Previous Article in Journal
Validation of the Hazard and Vulnerability Analysis of Coastal Erosion in the Caribbean and Pacific Coast of Colombia
Previous Article in Special Issue
Second Path Planning for Unmanned Surface Vehicle Considering the Constraint of Motion Performance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling Analysis and Simulation of Viscous Hydrodynamic Model of Single-DOF Manipulator

School of Mechanical Engineering, Hebei University of Technology, Tianjin 300132, China
*
Author to whom correspondence should be addressed.
J. Mar. Sci. Eng. 2019, 7(8), 261; https://doi.org/10.3390/jmse7080261
Submission received: 15 July 2019 / Revised: 7 August 2019 / Accepted: 7 August 2019 / Published: 9 August 2019
(This article belongs to the Special Issue Intelligent Marine Robotics Modelling, Simulation and Applications)

Abstract

:
Hydrodynamic modeling is the basis of the precise control research of underwater manipulators. Viscous hydrodynamics, an important part of the hydrodynamic model, directly affects the accuracy of the dynamic model and the control model of the manipulator. Considering the limited research on viscous hydrodynamics of underwater manipulators and the difficulty in measuring viscous hydrodynamic coefficients, the viscous hydrodynamic model in the form of Taylor expansion is analyzed and established. Through carrying out simulation calculations, curve fitting and regression analysis, positional derivatives, rotational derivatives, and coupling derivatives in the viscous hydrodynamic model, are determined. This model provides a crucial theoretical foundation and reference data for subsequent related research.

1. Introduction

With the rapid development of the marine industry, the application of robots mounted on underwater vehicles has become increasingly widespread. Underwater manipulators are mainly employed in the acquisition and exploration, etc. Their underwater performance can be determined by their speed and accuracy. Underwater manipulators are usually subjected to large forces and moments during their operation. These forces and moments mainly include gravity, inertial hydrodynamics, and viscous hydrodynamics. There are more solutions to the determination of gravity and inertial hydrodynamics. Therefore, the viscous hydrodynamic model has become a research focus that needs to be broken through because of its large number of coefficients and the difficulty of measurement.
References [1,2,3] verified the accuracy of fluid simulation calculations using software, such as FLUENT, by comparing the simulated calculation values of hydrodynamics with the measured values of engineering methods. References [4,5,6] studied the drag and additional mass force of the manipulator underwater and obtained the inertial hydrodynamic model by the CFD (Computational Fluid Dynamics) simulation, but the viscous hydrodynamics were not considered as key targets in the study. References [7,8,9,10] numerically calculated the viscous hydrodynamic coefficients of submarine and ship models with complex shapes under the specified navigation conditions and further improved the handling performance during navigation. Reference [11] determined the viscous kinetic coefficients of the manipulator through aerodynamic experiments and carried out a sea trial.
At present, research on viscous hydrodynamics mostly focus on specific ship models and underwater vehicles, and less on manipulators loaded on them. Therefore, a numerical calculation method of the viscous hydrodynamic coefficient based on single-DOF (Degree of Freedom) manipulators is put forward. The method is based on ANSYS Fluent. UDF (user-defined function) and dynamic mesh which are used to simulate the rotational motion of the manipulator. The inlet and outlet conditions, as well as boundary conditions of the fluid domain, are changed. Also, the viscous hydrodynamic model could be obtained by curve fitting and regression analysis.

2. Mathematical Model and Theoretical Analysis

2.1. Mathematical Modeling of Viscous Hydrodynamics

Since its surrounding flow field is changed during the operation of the underwater manipulator; the arm is subjected to the reaction of the water body caused by force generated by the water body. The hydrodynamics of the moving manipulator in water include two parts: Inertial hydrodynamics caused by acceleration and viscous hydrodynamics caused by friction. Inertial hydrodynamics is generally determined by empirical formulas, while coefficients of viscous hydrodynamics are complex and difficult to measure.
This paper is intended to use a single-DOF manipulator model mounted on a fixed base. The arm is a homogeneous lightweight rod with a circular cross-section, 50 mm in diameter (d), 500 mm in arm length (l), at length to diameter ratio of 10, belonging to the slender pole. The establishment of its coordinate system is shown in Figure 1.
When the manipulator performs the constant motion, the influence of acceleration and angular acceleration on the motion variable could be ignored based on the “slow-motion” assumption. So the mere considerations for viscous hydrodynamics are the effects of the velocity and angular velocity in the motion variable. It can be represented as: U = [ u   v   w   p   q   r ] T , where u stands for the velocity in the OX direction, v for the velocity in the OY direction, w for the velocity in the OZ direction, p for the angular velocity of the rotation around the OX axis, q for the angular velocity of the rotation around the OY axis, and r for the angular velocity of the rotation around the OZ axis.
The viscous hydrodynamics can be expressed as a multivariate function F = f ( u , v , w , p , q , r ) , and the six-dimensional component of viscous hydrodynamics F = [ X   Y   Z   K   M   N ] T could also be displayed in the above functional form, the direction of which is shown in Figure 2.
The viscous hydrodynamic F performs Taylor expansion in the form of Equation (1):
f ( U ) = f ( U k ) + ( U U k ) T f ( U k ) + 1 2 ! ( U U k ) T H ( U k ) ( U U k ) + ο n
where U k = [ u 0   v 0   w 0   p 0   q 0   r 0 ] T , U k is the initial constant.
H ( U k ) =   [ 2 f ( U k ) u 2 2 f ( U k ) u v 2 f ( U k ) u r 2 f ( U k ) u v 2 f ( U k ) v 2 2 f ( U k ) v r 2 f ( U k ) u r 2 f ( U k ) v r 2 f ( U k ) r 2 ]
Considering the left-right and front-back symmetry and the motion limit of the manipulator, several coefficients in the second-order expansions are zero, and the remaining others are non-negligible hydrodynamic coefficients, as in Equation (2):
{ X = X uu u 2 + X vv v 2 + X ww w 2 + X rr r 2 + X vr vr Y = Y v v + Y r r + Y v | v | v | v | + Y r | r | r | r | + Y v | r | v | r | Z = Z w w + Z w | w | w | w | + Z vv v 2 + Z rr r 2 + Z vr vr K = K v v + K v | v | v + K r r + K r | r | r | r | M = M w w + M w | w | w | w | + M vv v 2 + M rr r 2 + M vr vr N = N v v + N r r + N v | v | v | v | + N r | r | r | r | + N v | r | v | r |
Viscous hydrodynamics of the manipulator could be calculated in FLUENT. The regression analysis of the calculated viscous hydrodynamics and corresponding velocities could be performed to obtain unknown coefficients in Equation (2) in MATLAB. The first derivative coefficients associated only with the linear velocity (u, v, w) are the position derivatives, and the first derivative coefficients related to the angular velocity (p, q, r) are the rotational derivatives. The coefficients caused by joint changes of two or more parameters are the coupling derivatives.

2.2. Control Equation

To analyze hydrodynamics of underwater manipulators, it is generally assumed that the fluid is isothermal and incompressible, also as a constant flow that magnitude and direction do not change with time. The continuity equation (Equation (3)) and the Navier–Stokes equation (Equation (4)) serve as the two basic equations necessary to solve the flow problem of viscous fluid. These equations are generally described in the form of partial differentials:
u = 0
u t + ( u ) u = f 1 ρ p + ν 2 u
The form of expression of the time-averaged continuity equation does not change. Instead, the tensor of the Reynolds stress is added to the formula after the N–S equation time-averaged, which leads to the closure problem of the original equation. The Reynolds stress is about 10 2 Pa, which is indicative of the turbulent flow and cannot be directly ignored. Therefore, it is necessary to introduce a proper turbulence model to make a modeling description of the increased Reynolds stress.
For the incompressible isothermal turbulent water environment, the basic equations of turbulence include the DNS equations (direct numerical simulation), LES equations (large eddy simulation), and RANS equations (Renault time average). The two equations (DNS equations and LES equations) are of limited use due to the requirements of a large number of computational grids. At present, most engineering calculations adopt the RANS equation to solve closed equations formed by introducing the turbulence model and, thus, obtain the time-average value of the turbulent elements. To solve the viscous hydrodynamics under steady-state, the appropriate turbulence model is the key to the numerical simulations in this paper.

3. Calculation Method Analysis

3.1. Calculation Domain Establishment

Generally speaking, the larger the size of the calculation domain is, the closer it is to the real working condition. The downside is that it will increase the amount of calculation and prolong the calculation period. If the calculation domain is too small, the boundary conditions and calculation results are difficult to match the real working conditions. Therefore, it is very important to choose the size of the calculation domain reasonably.
Based on previous experience and multiple numerical practices [6,12,13], this paper establishes the domain as the computational domain, shown in Figure 3b. The specific sizes are as follows:
  • Front boundary: 1.5 H
  • Back boundary: 2 H
  • Side boundary: 1.5 H
where, H stands for the sum of the height of the arm and the base.
The motion parameters of the simulation are shown in Figure 3, α stands for the rotation angle around the OZ axis, β for the angle of the arm to the flow direction, V for the inlet flow velocity of the domain, and the linear velocity components of the manipulator are as shown in Equation (5):
{ u = Vcos β cos α v = Vcos β sin α w = Vsin β

3.2. Meshing

The calculation of CFD requires well-distributed grids. Usually, grids are classified into structured and unstructured grids. The unstructured grid means that internal points in the mesh area do not have the same adjacent unit, with no regular topology, not arranged in layers, and the distribution of mesh nodes is arbitrary. Therefore, it is more flexible than structured grid. Unstructured grids can be optimized by using certain criteria in the process of their generation. Ultimately, they can be displayed as high-quality meshes, which are suitable for complex geometry, easy to control grid size, and node density. Moreover, the adoption of random data structures is conducive to mesh adaptation. It is difficult to model the structured mesh due to the shape of the manipulator, instead, the unstructured mesh is easy to combine with the dynamic mesh technology. Thus, the unstructured mesh is adopted in the research.
We compare basic mesh to dense mesh for the sake of grid independence verification. Basic mesh in the vicinity of the manipulator is shown in Figure 4a and dense mesh is shown in Figure 4b. The number of elements and nodes are shown in Table 1.
The rotation process of the manipulator is realized by UDF and dynamic mesh technology in Fluent. Dynamic meshing is performed by using the spring smoothing model, which approximates connecting lines between the grid nodes as springs, and the position of nodes after smoothing is obtained by calculating force equilibrium equations between them. In the calculation process, meshes with a large aberration rate or huge change in size are grouped together to re-divide the partial meshes.

3.3. The Solution of Position Derivatives, Rotation Derivatives, and Coupling Derivatives

Solving the unknown coefficients in Equation (2) is the key to the research task, where Yv, Nv, Zw, and Mw are positional derivatives, Yr and Nr are rotational derivatives, and the remaining unknown coefficients are coupled derivatives.
The positional derivatives could be obtained by simulating the low-speed wind tunnel test, and the rotational derivatives are obtained by measuring the viscous hydrodynamic of the model at different rotational angular velocities. The number of coupled derivatives is large, and the viscous hydrodynamic subjected to the rotational manipulator is measured when β is not 0. After extensive experiments, the viscous hydrodynamic coefficients were obtained by least-square regression analysis.
The calculation domain inlet is set as the velocity boundary, the outlet as the free-flow condition, and the calculation domain wall as the non-slip fixed wall.
According to the research of the turbulence model in the reference [13], the standard k–ω model has many advantages, such as good numerical stability, accurate solution of pressure gradient, low Reynolds number influence, compressibility effect, and shear flow diffusion. It has better adaptability in calculating the flow problem of the boundary layer separation. It is one of the most widely used turbulence models for the viscous hydrodynamic solution. Therefore, the standard k–ω model is used as the turbulence model in this paper. The equation of the kinetic energy k and the turbulent frequency ω are as shown in Equation (6). The determined model parameters are shown in Table 2, according to the reference [13]:
{ ρ k t + ( ρ U k ) = [ ( μ + μ t σ k ) k ] + P k β ρ k ω ρ ω t + ( ρ U ω ) = [ ( μ + μ t σ ω ) ω ] + α t ω k P k β t ρ ω 2

4. Result Analysis

Comparing the viscous hydrodynamics of the manipulator measured by adopting the basic and dense mesh, respectively, we conclude that the differences are in the range from 1.7% to 4.3%. The differences are sufficiently low and negligible.
Computations are executed in a 64-bit processor consist of CPU (Intel Core i5-8400 @ 2.80 GHz) and 8 GB accessible memory, and take 37 h with at least 40 iterations per time step.
Using the size of calculation domain set in Section 3.1, the flow velocity of the calculation domain is kept constant, the angle β between the manipulator and the incoming flow direction is adjusted, the viscous hydrodynamics at different angles are calculated, and then the viscous hydrodynamic position derivative is obtained by linear analysis. The position derivative calculation contents are shown in Table 3.
Since the different degrees of α and β cause the manipulator to have different linear velocities in the OX, OY, and OZ directions, viscous hydrodynamics of the arm is measured in the horizontal plane XOZ and the vertical plane XOY, respectively. The values of the vertical force Y and the pitching moment N measured at the different linear velocity v in the OY direction are shown in Table 4.
The position derivatives Yv and Nv are the first derivative coefficients of the linear velocity v, so the least squares curve fitting of the vertical force Y and the pitching moment N in viscous hydrodynamics to the linear velocity v is performed, as shown in Figure 5 and Figure 6. The derivative value at the median 0 points of the line velocity v is taken as the position derivative.
Same as above, the values of the lateral force Z and the yaw moment M measured at the different linear velocity w in the OZ direction are shown in Table 5.
The position derivatives Zw and Mw are the first derivative coefficients of the linear velocity w, so the least squares curve fitting of the viscous hydrodynamic lateral force Z and the yawing moment M to the linear velocity w is performed, as shown in Figure 7 and Figure 8. The derivative at 0 point of the median line velocity w is taken as the position derivative.
The positional derivatives could be obtained as shown in Table 6 below:
Since the manipulator studied in this paper is a single DOF, only the viscous hydrodynamics are measured when it rotates around the OZ axis. UDF is applied to adjust the rotational velocity of the arm. The measurement conditions of the rotation derivatives are shown in Table 7.
This paper calculates the hydrodynamics of the rotating underwater manipulator when the flow is stationary. The arm rotates in the XOY plane. The angular velocity r is adjusted according to Table 7, and the instantaneous viscous hydrodynamics of the arm during the rotation from −10° to 10° around the OZ axis are measured. The measured vertical force Y and pitching moment N are shown in Table 8.
The rotational derivatives Yr and Nr are the first derivative coefficients of the angular velocity r, similar to the solution method of the position derivatives. The least squares curve fitting of the vertical force Y and the pitching moment N in viscous hydrodynamics to the angular velocity r is performed, as shown in Figure 9 and Figure 10. The derivative at the median angular velocity r of 1.25 rad/s is taken as the rotation derivative.
The rotational derivatives could be obtained as shown in Table 9 below:
The coupled derivative is the second-order viscous hydrodynamic coefficient of the arm under complex motion conditions. Through the transient iterative calculation of the motion of the manipulator at different α and β angles and various angular velocities, the viscous hydrodynamics of each motion state are measured, and many coupling derivatives are obtained by least-squares regression analysis. The coupling derivative solution conditions are shown in Table 10 below.
The coupled derivatives are obtained by regression analysis, and the obtained coupled derivatives are processed without dimensioning. The values of the parameters are as shown in Table 11.

5. Conclusions

The hydrodynamics of underwater manipulators during operation are complex and difficult to predict. It is analyzed that components of the hydrodynamics include inertial hydrodynamics caused by acceleration and viscous hydrodynamics caused by friction. This paper takes viscous hydrodynamics as the research target.
According to the research on dynamics of AUV (Autonomous Underwater Vehicle), ROV(Remote Operated Vehicle) and ships in other references, the viscous hydrodynamic model of the manipulator is analyzed in the form of Taylor expansion, and the viscous hydrodynamics are measured by using a single-DOF manipulator to simulate the underwater motion, and the viscous hydrodynamic coefficients among the model are calculated via regression analysis. An accurate viscous hydrodynamic model is obtained to predict viscous hydrodynamics of manipulators during operation at any attitude.
This model is the basis for the feedforward control and is helpful to study control stability of underwater manipulators. We believe that although the simulations in this paper were performed for single-DOF manipulators, the modeling method may be extended for the manipulator with multi-DOF and more complicated shapes as is in the case of real underwater manipulators.

Author Contributions

M.Z. Provided supervision and funding support for the project. X.L. conducted formal analysis, designed all the experiments, and subsequently drafted the manuscript. Y.T. conceived the original ideas, provided supervision to the project and reviewed the writing.

Funding

This research was supported by the National Key Research and Development Program of China (Grant No. 2018YFB1309401), and the Hebei Provincial Higher Education Science and Technology Research Project (Grant No. QN2018090).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Stern, F.; Agdrup, K.; Kim, S.Y.; Hochbaum, A.C.; Rhee, K.P.; Quadvlieg, F.; Perdon, P.; Hino, T.; Broglia, R.; Gorski, J. Experience from SIMMAN 2008—The First Workshop on Verification and Validation of Ship Maneuvering Simulation Methods. J. Ship Res. 2011, 55, 135–147. [Google Scholar]
  2. Larsson, L. Numerical Ship Hydrodynamics; Springer: Dordrecht, The Netherlands, 2013. [Google Scholar]
  3. Irwin, R.P.; Chauvet, C. Quantifying Hydrodynamic Coefficients of Complex Structures; IEEE: New York, NY, USA, 2007; pp. 1341–1345. [Google Scholar]
  4. Mclain, T.W.; Rock, S.M. Development and Experimental Validation of an Underwater Manipulator Hydrodynamic Model. Int. J. Robot. Res. 1998, 17, 748–759. [Google Scholar] [CrossRef] [Green Version]
  5. Pazmiño, R.S.; Cena, C.E.G.; Arocha, C.A.; Santonja, R.A. Experiences and results from designing and developing a 6 DoF underwater parallel robot. Robot. Auton. Syst. 2011, 59, 101–112. [Google Scholar] [CrossRef]
  6. Kołodziejczyk, W. Preliminary Study of Hydrodynamic Load on an Underwater Robotic Manipulator. J. Autom. Mob. Robot. Intell. Syst. 2015, 9, 11–17. [Google Scholar] [CrossRef]
  7. Arabshahi, A.; Beddhu, M.; Briley, W.R.; Chen, J.P.; Gaither, A.; Gaither, K.; Janus, J.M.; Jiang, M.; Marcum, D.; Mcginley, J. A Perspective on Naval Hydrodynamic Flow Simulations. In Proceedings of the Symposium on Naval Hydrodynamics, Washington, DC, USA, 9–14 August 1998. [Google Scholar]
  8. Ohmori, T.; Fujino, M.; Miyata, H. A study on flow field around full ship forms in maneuvering motion. J. Mar. Sci. Technol. 1998, 3, 22–29. [Google Scholar] [CrossRef]
  9. Tahara, Y.; Longo, J.; Stern, F. Comparison of CFD and EFD for the Series 60 CB = 0.6 in steady drift motion. J. Mar. Sci. Technol. 2002, 7, 17–30. [Google Scholar] [CrossRef]
  10. Ueno, M. Hydrodynamic derivatives and motion response of a submersible surface ship in unbounded water. Ocean Eng. 2010, 37, 879–890. [Google Scholar] [CrossRef]
  11. Filaretov, V.F.; Konoplin, A.J.; Getman, A.V. Experimental Determination of the Viscous Friction Coefficients for Calculation of the Force Impacts on the Moving Links of the Underwater Manipulators. Mekhatronika Avtom. Upr. 2015, 16, 738–743. [Google Scholar] [CrossRef]
  12. Kolodziejczyk, W. The method of determination of transient hydrodynamic coefficients for a single DOF underwater manipulator. Ocean Eng. 2018, 153, 122–131. [Google Scholar] [CrossRef]
  13. Zhiqiang, H.U.; Ruiwen, Y.I.; Lin, Y.; Haitao, G.U.; Wang, C. Numerical calculation methods for hydrodynamics of unmanned underwater vehicles based on body-fixed coordinate frames. Chin. Sci. Bull. 2013, 58, 55. [Google Scholar]
Figure 1. Manipulator model and coordinate system.
Figure 1. Manipulator model and coordinate system.
Jmse 07 00261 g001
Figure 2. Six-dimensional component of viscous hydrodynamic F.
Figure 2. Six-dimensional component of viscous hydrodynamic F.
Jmse 07 00261 g002
Figure 3. Schematic diagram of the calculation domain and α / β .
Figure 3. Schematic diagram of the calculation domain and α / β .
Jmse 07 00261 g003
Figure 4. Meshing of manipulator section.
Figure 4. Meshing of manipulator section.
Jmse 07 00261 g004
Figure 5. Fitted curve of vertical force Y and line velocity v.
Figure 5. Fitted curve of vertical force Y and line velocity v.
Jmse 07 00261 g005
Figure 6. Fitted curve of pitching moment N and linear velocity v.
Figure 6. Fitted curve of pitching moment N and linear velocity v.
Jmse 07 00261 g006
Figure 7. Fitted curve of lateral force Z and linear velocity w.
Figure 7. Fitted curve of lateral force Z and linear velocity w.
Jmse 07 00261 g007
Figure 8. Fitted curve of yaw moment M and line speed w.
Figure 8. Fitted curve of yaw moment M and line speed w.
Jmse 07 00261 g008
Figure 9. Fitted curve of vertical force Y and angular velocity r.
Figure 9. Fitted curve of vertical force Y and angular velocity r.
Jmse 07 00261 g009
Figure 10. Fitted curve of pitching moment N and angular velocity r.
Figure 10. Fitted curve of pitching moment N and angular velocity r.
Jmse 07 00261 g010
Table 1. Parameters of basic and dense mesh.
Table 1. Parameters of basic and dense mesh.
Maximum LayersGrowth RateElementsNodes
Basic Mesh51.2636,457114,559
Dense Mesh51.21,429,231401,609
Table 2. Standard k–ω model parameters.
Table 2. Standard k–ω model parameters.
σ k σ ω α t β t β
2.02.05/90.0750.09
Table 3. Solving cases of position derivatives.
Table 3. Solving cases of position derivatives.
α 0°, ±2°, ±5°, ±7°, ±10°
β 0°, ±2°, ±4°, ±6°, ±8°, ±10°
r 0 rad/s
V 1 m/s
Table 4. Measurements of vertical force Y and pitching moment N at different linear velocity v.
Table 4. Measurements of vertical force Y and pitching moment N at different linear velocity v.
v (m/s)Y (N)N (N·m)
−0.1714−0.03730.0677
−0.1200−0.02350.0673
−0.0860−0.01350.0670
−0.03400.00030.0660
00.01040.0658
0.03400.01900.0652
0.08600.03310.0649
0.12000.04200.0640
0.17140.05450.0628
Table 5. Measurement of the lateral force Z and the yaw moment M at the different linear velocity w.
Table 5. Measurement of the lateral force Z and the yaw moment M at the different linear velocity w.
w (m/s)Z (N)M (e−5·N·m)
−0.1740−0.0364−4.8199
−0.1390−0.0280−2.5920
−0.1050−0.01841.5570
−0.0700−0.01002.0860
−0.0350−0.00252.1628
00.00664.1640
0.03500.01655.4699
0.07000.02508.9650
0.10500.03289.3461
0.13900.043012.3140
0.17400.052616.7242
Table 6. Position derivatives of the manipulator in viscous hydrodynamics.
Table 6. Position derivatives of the manipulator in viscous hydrodynamics.
Y v N v Z w M w
1.4352−0.14721.35200.0117
Table 7. Solving cases of rotational derivatives.
Table 7. Solving cases of rotational derivatives.
α −10°–10°
β
r0.5, 0.75, 1, 1.25, 1.5, 1.75, 2.0 (rad/s)
V0 m/s
Table 8. Measurements of vertical force Y and pitching moment N at different angular velocities r.
Table 8. Measurements of vertical force Y and pitching moment N at different angular velocities r.
r (rad/s)Y (N)N (N·m)
0.50.00330.0444
0.750.01130.0731
1.00.01820.1023
1.250.03060.1462
1.50.03990.1750
1.750.04920.2181
2.00.05920.2613
Table 9. Rotation derivatives of the manipulator in viscous hydrodynamics.
Table 9. Rotation derivatives of the manipulator in viscous hydrodynamics.
Y r N r
0.40433.0869
Table 10. Solving cases of coupling derivatives.
Table 10. Solving cases of coupling derivatives.
α 0°, ±2°, ±5°, ±7°, ±10°
β 0°, ±2°, ±4°, ±6°, ±8°, ±10°
r0.5 rad/s, 1 rad/s
V1 m/s
Table 11. Coupling derivatives of the manipulator in viscous hydrodynamics.
Table 11. Coupling derivatives of the manipulator in viscous hydrodynamics.
X uu −0.6960 X vv −3.8896
X ww −74.4904 X rr 44.8128
X vr −275.0688 Y v | v | 1.172
Y r | r | −6.2688 Y v | r | 31.6064
Z w | w | 13.0336 Z vv 2.1088
Z rr −0.0512 Z vr −0.1264
K v −0.0011 K v | v | −0.0160
K r 0.1515 K r | r | 0.0256
M w | w | −1.0816 M vv 0.0144
M rr 0.0384 M vr −0.2144
N v | v | 32.9184 N r | r | −25.0688
N v | r | 130.4128

Share and Cite

MDPI and ACS Style

Zhang, M.; Liu, X.; Tian, Y. Modeling Analysis and Simulation of Viscous Hydrodynamic Model of Single-DOF Manipulator. J. Mar. Sci. Eng. 2019, 7, 261. https://doi.org/10.3390/jmse7080261

AMA Style

Zhang M, Liu X, Tian Y. Modeling Analysis and Simulation of Viscous Hydrodynamic Model of Single-DOF Manipulator. Journal of Marine Science and Engineering. 2019; 7(8):261. https://doi.org/10.3390/jmse7080261

Chicago/Turabian Style

Zhang, Minglu, Xiaoyu Liu, and Ying Tian. 2019. "Modeling Analysis and Simulation of Viscous Hydrodynamic Model of Single-DOF Manipulator" Journal of Marine Science and Engineering 7, no. 8: 261. https://doi.org/10.3390/jmse7080261

APA Style

Zhang, M., Liu, X., & Tian, Y. (2019). Modeling Analysis and Simulation of Viscous Hydrodynamic Model of Single-DOF Manipulator. Journal of Marine Science and Engineering, 7(8), 261. https://doi.org/10.3390/jmse7080261

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop