Kinematics Modeling and Simulation of a Bionic Fish Tail System Based on Linear Hypocycloid
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.
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 . 85% of fish are swimming in BCF mode for power supply, supplemented by MPF mode to keep bodies balance, retreating, hovering, and turning movement .
Some researchers devote themselves to reveal kinematics and hydrodynamics of live fish. Lightill  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 . Vo et al.  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.  designed a planar carangiform robot fish with motion control algorithms to obtain experimental trajectory tracking results. Esposito et al.  presented a robotic fish caudal fin with six individually moveable fin rays based on sunfish tail. Fay et al.  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.  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
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
Therefore, the swing angle in (9) could be simplified as
(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 . 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) 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.
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.
This research is financially supported by National Nature Science Foundation of China (51205173) and by Jiangsu Government Scholarship for overseas studies.
J. J. Videler, Fish Swimming, Springer, Dordrecht, The Netherlands, 1993.View at: Publisher Site
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