Research Article | Open Access

# Kinematics Modeling and Simulation of a Bionic Fish Tail System Based on Linear Hypocycloid

**Academic Editor:**Jan Harm Koolstra

#### Abstract

Kinematics and simulation study on a two-joint linear hypocycloid tail driving system composed of a special planetary gear system and a linkage mechanism are conducted in this paper. First, the composition and working principle of the linear hypocycloid tail transmission system are introduced and analyzed. Second, the kinematics study on the transmission mechanism is conducted with graphical method of vector equation. The relationships between the caudal peduncle stroke, the tail fin swing angle, and the phase difference with structure parameters are studied, and further optimization of structure sizes (i.e., linkage length, sun gear’s diameter, the intersection angle between planet gears, etc.) is developed. At last, simulation and comparative study on a biofish in sample parameters with a live fish of Carp is conducted in MATLAB. The study would serve for underwater vehicles thruster design and its mechanism.

#### 1. Introduction

Bionic propulsion device inspired by fish swimming skills to replace traditional underwater devices has caught much attention of biologists and engineers all over the world. Compared to traditional screw propellers, bionic fish propulsion has its unique advantages in high efficiency, low noise, and great mobility [1–3]. Previous investigations have shown that fish swimming in BCF mode can obtain a larger propulsion force during escape and prey, but fish swimming in MPF mode can obtain higher stability and maneuverability [4]. 85% of fish are swimming in BCF mode for power supply, supplemented by MPF mode to keep bodies balance, retreating, hovering, and turning movement [1].

Some researchers devote themselves to reveal kinematics and hydrodynamics of live fish. Lightill [5] put forward the elongated-body theory where the movements of any horizontal section of caudal fin, with yaw angle fluctuating in phase with its velocity of lateral translation, were studied for different positions of the yawing axis. In addition, the proposed theory was extended to the large-amplitude elongated-body theory so that a prediction of instantaneous reactive force between fish and water was achieved for fish motions of arbitrary amplitude [6]. Vo et al. [7] proposed an analytical optimization approach which can guarantee the maximum propulsive velocity of fish robot in the given parametric conditions. Other researchers proposed various bionic vehicles to mimic the swimming of live fish. Morgansen et al. [8] designed a planar carangiform robot fish with motion control algorithms to obtain experimental trajectory tracking results. Esposito et al. [9] presented a robotic fish caudal fin with six individually moveable fin rays based on sunfish tail. Fay et al. [10] developed a wireless aquatic bionic fish with a wireless video camera, a controller, and polypyrrole actuators to detect and analyze pollutants in natural waters. Yun et al. [11] applied a special waving caudal fin with vertical phase differences to reduce reaction torque and to improve bionic fish’s velocity and stability.

The purpose of this paper is to design a linear hypocycloid driving mechanism, which has advantages of combining speed reducer with transformation mechanism, and adjustable phase difference between caudal peduncle and tail fin. In this paper, kinematics of the driving system was analyzed comprehensively, and structural parameters optimization is developed for mimicking a real fish tail’s oscillating motion, which is verified by a further comparative study with a live Carp.

#### 2. Structure and Working Principle

The bionic fish tail system based on linear hypocycloid is composed of a special planetary gear train with the reference circle of the planetary gear whose reference radius is half of the sun gear and a plane linkage which can form a variable triangle motion relation, as shown in Figure 1(a). The proposed special planetary gear train involves a type planetary carrier, two planetary gears installed at the planetary carrier in two parallel planes, and the sun gear. The plane linkage consists of two links: one end of two links is connected to rotate with a fixed point located at the reference circle of a planetary gear, respectively, and the other end of two links is connected to rotate at a certain point.

**(a) 3D structure**

**(b) The working mechanism**

The working mechanism of the bionic fish tail system based on linear hypocycloid is shown in Figure 1(b). When the two planetary gears 3 and 4 meshed with the sun gear 2 under the driving of the planetary carrier 1 at a certain speed, Point or Point located at the reference circle of the planetary gear 4 or 3 will be reciprocating along the connected line of the point and the sun gear’s central point. In case of Point and Point designed at one same diameter line of the sun gear 2, both Point and Point would be reciprocating in the diameter line with a phase difference decided by shape’s angle of planetary carrier 1. With the reciprocating motion of Point and Point with a certain phase difference, a motion triangle in the plane linkage mechanism is formed. As a result, the tail fin link 6 will rotate around Point and make reciprocating movement with Point along the diameter line . Obviously, the relative length between and , the length of or , shape’s angle of planetary carrier, and the diameter of the sun gear will be important parameters for the driving system, which will be further discussed in this paper.

#### 3. Kinematics on the Tail Driving System

##### 3.1. Foundation of Coordinate System

Three coordinate systems are employed for the tail driving system as shown in Figure 2. The first fixed coordinate system is - with the centre of the sun gear 2 as the origin Point , axis points at the horizontal direction, and axis is upward. With the planetary gear 4 meshing with the sun gear, the second moving coordinate system - is connected to the connection line of Point and Point , the origin point is central point of the sun gear, is along the line of , and is perpendicular to the line of .

##### 3.2. Kinematics of the Tail Driving System

###### 3.2.1. Kinematics of the Planetary Gear Train

Kinematic relation of the linear hypocycloid planetary gear train can be deduced easily based on relative kinematics. When the planetary gear meshed with the sun gear at pitch Point shown in Figure 2, in coordinate system, velocity vector equation could be written as

Here, is the relative velocity of Point to Point , is the relative velocity of Point to Point , is the absolute velocity of the planetary gear 4 at the meshing Point . When the sun gear is fixed, . Then (1) could be replaced as

Here, is the angular velocity of planetary gear 4 in its revolution, and is the angular velocity of planetary gear 4 in its rotating motion. , and is the reference radius of planetary gear 4.

Therefore, , which shows and with the same magnitude but in opposite direction, as shown in Figure 2.

###### 3.2.2. The Working Mechanism of Linear Hypocycloid

For the proposed hypocycloidal gear train in this paper, the track of arbitrary point located at the reference circle of planetary gear should be one certain diameter line of the sun gear. The specific working mechanism of linear hypocycloid is shown in Figure 3. In the moving coordinate system , an arbitrary point located at the reference circle of planetary gear 4 is marked as Point , and the position of would be supposed to move to a new spot marked with after planetary motion with any time . Based on the closed vector triangle , equations would be gained as follows:

Here, is the angle of and , and is the angle between and . , where is the reference diameter of the sun gear 2, and with the proposed relation .

As a result, must be equal to zero, so (4) could be simplified as

Equation (4) shows that Point is reciprocating along the diameter line in harmonic motion, and its stroke is the diameter of the sun gear.

If Points and were selected at the same diameter line of the sun gear but in different planetary gears installed at planetary carrier, Point should also do reciprocating motion along the same diameter line of sun gear with a certain phase difference , and the phase difference is decided by shape’s angle of the planetary carrier. Therefore, the motion equation of Point could be described as

The proposed and were deduced in the moving coordinate system . If putting and into the fix coordinate system , the equations of Point and Point could be deduced as

Here, is the angle between the axis and the axis.

###### 3.2.3. The Working Mechanism of the Motion Triangle

To simplify kinematic analysis, the link mechanism in the driving system would be replaced with an equivalent mechanism by removing the planetary gear train, shown in Figure 4. In the equivalent mechanism, link 5, tail fin link 6, slider 7, and slider 8 are connected at Points , , and to form the equivalent mechanism. Here, link 5 would rotate around Point of slider 8, the tail fin link 6 would rotate around Point of slider 7, and the tail fin link 6 and link 5 are connected to rotate at Point . Slider 7 and slider 8 do the same reciprocating motion with Points and , respectively, so the motion triangle still remained the same as the driving system. As a result, the tail fin link 6 would gain a composite motion of reciprocating with slider 7 along line and oscillating around Point .

Based on vector triangles and , kinematics model of the equivalent mechanism would be established as

Here, is the swing angle of the tail fin link 6, is the swing angle of link 5, and are rod lengths of the tail fin link 6 and link 5, respectively, is the motion position of slider 7, is the motion position of Point in axis, and is the motion position of Point in axis.

The instantaneous position distance of the two sliders could be described as

Based on (7) and (8), displacement equations of the tail fin link 6 with composite motion of reciprocating and oscillating could be described as

With derivation of (9), velocity equations of the tail fin link could be deduced as

Here, is the angular velocity of the tail fin link 6, is the velocity of slider 7, and is velocity difference between slider 7 and slider 8.

With derivation of (10), acceleration equations of the tail fin link could be deduced as

Here, is the angular accelerated velocity of the tail fin link 6, is the accelerated velocity of Point along , is the accelerated velocity of slider 7, and is the accelerated velocities difference of two sliders.

#### 4. Optimal Design on Structural Parameters

It is obvious that specific parameters such as the phase difference, rod length, relative length of two rods, and reference diameter of the sun gear will directly or indirectly affect the behavior of tail link 6. In this chapter, we are focused on developing optimal parameters to make the motion triangle which existed in the whole cycle avoid some extreme situations and make the tail fin’s behavior mimic real fish’s caudal fin. In order to simplify the problem, the influence of friction and gravity is supposed to be ignored.

##### 4.1. Rod Length Relation of the Planar Linkage

Based on the motion triangle , the side length relation could be written as

Simultaneous (12) and (8), the rod length must be satisfied with the following equation:

Therefore, the swing angle in (9) could be simplified as

Based on (14), the length principle with rad/s, , and mm is verified in MATLAB, as shown in Figure 5(a).

**(a) The swing angle -time curve with different rod length relations**

**(b) The extreme positions with**

**(c) The swing scope and motion trail with**

From Figure 5(a), when the rod lengths of the tail link 6 and the link 5 are unequal, the swing angle of the tail link 6 would vary irregularly and discontinuously, two sudden change points with peaking at 90° and −90° as shown in Figure 5(b). The two extreme positions at horizontal direction would cause destruction of the mechanism in its weak link joint. Only when = , as shown in Figure 5(c), the swing angle could vary regularly and smoothly with a sinusoidal motion in the whole cycle, and the trajectories of Point are two sine waves which are symmetric about axis.

##### 4.2. The Relation of Structural Parameters

The rod length of the tail link is also decided by stroke value in reciprocating, swing amplitude of the tail link, and phase difference with , and the rod length can be described by

With parameters mm and swing amplitude = 30°, 60°, and 90°, respectively, the relation between the rod length and phase difference is shown in Figure 6. With parameters mm and phase difference = 30°, 60°, and 90°, respectively, the relationship between the rod length and swing amplitude is shown in Figure 7.

From Figure 6, if the swing amplitude is fixed, the rod length can be adjusted to be shorter by decreasing the phase difference to gain more compact structure, and the rod length would peak when phase difference reaches 180°. The smaller the swing amplitude becomes, the faster the growth rate of the rod length would tend to grow. From Figure 7, the rod length could be adjusted to be shorter by decreasing the swing amplitude with a fixed phase difference, and the rod length will be the minimal length when the swing amplitude reaches 90°.

Except for the phase difference and the swing amplitude, the rod length is still determined by the reference diameter of the sun gear. With and , , and , respectively, the rod length varies with the reference diameter of the sun gear shown in Figure 8. The rod length increases linearly with the reference diameter of the sun gear, and the increase of the rod length would slow down with increasing the swing amplitude.

#### 5. Design for an Application Example of the Driving System

##### 5.1. A Sample of the Tail Driving System

A sample of the tail driving system was designed with parameters mm, , rad/s, and mm. The motion equation of the tail fin link could be specific as

##### 5.2. Parameters of the Live Fish

The researchers of National University of Singapore have observed a real Carp with a length of 190 mm by PIV [12]. In their works, they selected four feature points for estimating Carp’s joint angles but focused on kinematics studies of Point located on Carp’s peduncle and Point located on its tail end, as shown in Figure 9. The motion of Point and the swing angle of Point were collected by videotaping the motion of fish at 60 frames per second using a video recording system. From the empirical observation in “Cruise” swimming, Carp continued to move in a nearly straight line at a constant speed, and the caudal fin flapped periodically. The trajectory of the caudal fin’s swing angle is shown in Figure 10(a) with dotted line, and the trajectory of caudal peduncle’s reciprocating is shown in Figure 10(b) with dotted line.

**(a)**

**(b)**

**(a) Comparative analysis on the swing angle of the tail fin**

**(b) Comparative analysis on the reciprocating motion of the caudal peduncle**

##### 5.3. Simulation and Comparative Study

From Figure 10, the symmetry harmonic movement of the driving system coincided with the live fish. As Figures 10(a) and 10(b) have shown, the swing angle of the caudal fin in our designed driving system is similar to the live fish, and the swing amplitude could reach 45°. The reciprocating motion of the caudal peduncle is also similar to the live fish, and the stroke could reach 18 cm. When the caudal peduncle is situated at the limiting positions ( cm or −9 cm), the swing angle of the tail link will be zero. When the caudal peduncle is in the balance position with , the swing angle of tail fin reaches the maximum: that is, °. Above all, the phase difference between the motion of caudal peduncle and tail fin’s swing is about 90° and the phase difference of two planetary gears is 180°; that is, the phase difference of the actual output from the mechanism equals the half of input phase difference.

#### 6. Conclusion

Based on linear hypocycloid, the driving system composed of a planetary gear train and a linkage has been developed for the tail transmission of two-joint bionic fish. The model and kinematics analysis of tail driving system were deduced by vector graphic method in this paper. The optimization of structure sizes was comprehensively studied to improve the kinematics performance. The simulation and comparative study on a bionic fish with a live fish were conducted so as to testify the feasibility of the driving system.

The results of structure size optimization show that the two rod lengths are equal to realize the tail oscillating continuously in cycle. In addition, the diameter of the sun gear , the phase difference of two planetary gears, and the tail swing amplitude together affect the rod length . The rod length will increase with the growth of and but decrease with the growth of .

With the optimized structure, simulation and comparative study with a sample size ( mm, mm, , rad/s, and mm) in MATLAB have been conducted with experimental results of a live Carp to verify the feasibility of the driving system. These studies will work for future experiment study and the development of the mechanism design in underwater propulsion.

#### Conflict of Interests

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

#### Acknowledgments

This research is financially supported by National Nature Science Foundation of China (51205173) and by Jiangsu Government Scholarship for overseas studies.

#### References

- L. Wang, M. Xu, B. Liu, K. H. Low, J. Yang, and S. Zhang, “A three-dimensional kinematics analysis of a Koi Carp pectoral fin by digital image processing,”
*Journal of Bionic Engineering*, vol. 10, no. 2, pp. 210–221, 2013. View at: Publisher Site | Google Scholar - J. J. Videler,
*Fish Swimming*, Springer, Dordrecht, The Netherlands, 1993. View at: Publisher Site - Y. Zhang, J. He, and G. Zhang, “Measurement on morphology and kinematics of crucian vertebral joints,”
*Journal of Bionic Engineering*, vol. 8, no. 1, pp. 10–17, 2011. View at: Publisher Site | Google Scholar - J. L. Tangorra, S. N. Davidson, I. W. Hunter et al., “The development of a biologically inspired propulsor for unmanned underwater vehicles,”
*IEEE Journal of Oceanic Engineering*, vol. 32, no. 3, pp. 533–550, 2007. View at: Publisher Site | Google Scholar - M. J. Lightill, “Aquatic animal propulsion of high hydromechanical efficiency,”
*Journal of Fluid Mechanics*, vol. 44, no. 2, pp. 265–301, 1970. View at: Publisher Site | Google Scholar - M. J. Lightill, “Large-amplitude elongated-body theory of fish locomotion,”
*Proceedings of the Royal Society of London B: Biological Sciences*, vol. 179, pp. 125–138, 1971. View at: Publisher Site | Google Scholar - T. Q. Vo, H. S. Kim, and B. R. Lee, “Propulsive velocity optimization of 3-joint fish robot using genetic-hill climbing algorithm,”
*Journal of Bionic Engineering*, vol. 6, no. 4, pp. 415–429, 2009. View at: Publisher Site | Google Scholar - K. A. Morgansen, V. Duindam, R. J. Mason, J. W. Burdick, and R. M. Murray, “Nonlinear control methods for planar carangiform robot fish locomotion,” in
*Proceedings of the IEEE International Conference on Robotics and Automation (ICRA '01)*, pp. 427–434, May 2001. View at: Google Scholar - C. J. Esposito, J. L. Tangorra, B. E. Flammang, and G. V. Lauder, “A robotic fish caudal fin: effects of stiffness and motor program on locomotor performance,”
*The Journal of Experimental Biology*, vol. 215, no. 1, pp. 56–67, 2012. View at: Publisher Site | Google Scholar - C. Fay, K.-T. Lau, S. Beirne et al., “Wireless aquatic navigator for detection and analysis (WANDA),”
*Sensors & Actuators B*, vol. 150, no. 1, pp. 425–435, 2010. View at: Publisher Site | Google Scholar - D. Yun, K.-S. Kim, S. Kim, J. Kyung, and S. Lee, “Actuation of a robotic fish caudal fin for low reaction torque,”
*Review of Scientific Instruments*, vol. 82, no. 7, Article ID 075114, 2011. View at: Publisher Site | Google Scholar - Q. Ren, J. Xu, L. Fan, and X. Niu, “A GIM-based biomimetic learning approach for motion generation of a multi-joint robotic fish,”
*Journal of Bionic Engineering*, vol. 10, no. 4, pp. 423–433, 2013. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2015 Shu-yan Wang 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.