Next Article in Journal
A High-Efficient Reinforcement Learning Approach for Dexterous Manipulation
Next Article in Special Issue
How Free Swimming Fosters the Locomotion of a Purely Oscillating Fish-like Body
Previous Article in Journal
Effect of Platelet-Rich Plasma Addition on the Chemical Properties and Biological Activity of Calcium Sulfate Hemihydrate Bone Cement
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational Study of Stiffness-Tuning Strategies in Anguilliform Fish

1
School of Aerospace Engineering, Guizhou Institute of Technology, Guiyang 550003, China
2
School of Mechatronic Engineering, Guizhou University, Guiyang 550025, China
*
Author to whom correspondence should be addressed.
Biomimetics 2023, 8(2), 263; https://doi.org/10.3390/biomimetics8020263
Submission received: 21 May 2023 / Revised: 13 June 2023 / Accepted: 14 June 2023 / Published: 16 June 2023
(This article belongs to the Special Issue New Insights into Biological and Bioinspired Fluid Dynamics)

Abstract

:
Biological evidence demonstrates that fish can tune their body stiffness to improve thrust and efficiency during swimming locomotion. However, the stiffness-tuning strategies that maximize swimming speed or efficiency are still unclear. In the present study, a musculo-skeletal model of anguilliform fish is developed to study the properties of variable stiffness, in which the planar serial-parallel mechanism is used to model the body structure. The calcium ion model is adopted to simulate muscular activities and generate muscle force. Further, the relations among the forward speed, the swimming efficiency, and Young’s modulus of the fish body are investigated. The results show that for certain body stiffness, the swimming speed and efficiency are increased with the tail-beat frequency until reaching the maximum value and then decreased. The peak speed and efficiency are also increased with the amplitude of muscle actuation. Anguilliform fish tend to vary their body stiffness to improve the swimming speed and efficiency at a high tail-beat frequency or small amplitude of muscle actuation. Furthermore, the midline motions of anguilliform fish are analyzed by the complex orthogonal decomposition (COD) method, and the discussions of fish motions associated with the variable body stiffness and the tail-beat frequency are also presented. Overall, the optimal swimming performance of anguilliform fish benefits from the matching relationships among the muscle actuation, the body stiffness, and the tail-beat frequency.

1. Introduction

In the swimming cycles, most fish contract their muscles left and right rhythmically and propel forward by the undulations of the body and/or caudal fin. Now more and more studies show that swimming fish use muscular activities to adjust their body stiffness in real-time [1,2], and the property of various stiffness has significant influences on the thrust and/or swimming efficiency [3,4,5]. In early works, Wu [6] studied two-dimensional flexible hydrofoils and found that flexible hydrofoils could be more efficient than rigid ones at producing thrust. Several theoretical and experimental studies also confirmed that adding flexibility to a hydrofoil can improve thrust and/or efficiency [7,8,9]. Tytell et al. [10,11] developed both physical and computational models of lamprey fish and found that the bending stiffness of the fish body affected the wake structures behind the tail, as well as the swimming efficiency. McHenry et al. [12] also found that flexibility can be used to increase the forward speed of swimming fish. Additionally, several fish-like robots were developed to study the role of flexibility. For example, Jusufi et al. [13] used antagonistic pneumatic actuators to develop a robotic fish with variable body stiffness; White et al. [14] added flexible joints into a tuna-inspired robot to study the capability of high-frequency locomotion.
From the biomechanical aspect, flexible fish would maximize their body deformations when actuated at their resonant frequency because of the phenomenon of resonance [15]. In the experiments of robotic lamprey [16], the results suggested that the resonance played an important role in maximizing efficiency, with an increased stride length but only a very small expansion of wake energy. By contrast, several studies of a flexible foil [17,18] found that the peaks in amplitude did not always correspond to peaks in thrust and/or efficiency. For example, Quinn et al. [2] found that the efficiency of robotic fish may peak at just 0.4–0.6 times the resonant frequency. Therefore, someone believes that the thrust and/or efficiency of swimming fish is maximized by changing the aerodynamic variables, such as the camber, the angle of attack, and the undulatory motions. For instance, the flexibility of the fish body created a passive phase offset between the actuation force and its response, and this phase offset could reduce the effective angle of attack at the leading edge of the tail/fin, which affected the formation of leading-edge vortices (LEVs) and leads edge suction, as well as the swimming efficiency [19,20]. Meanwhile, flexibility also affects the temporal evolution of fish propulsive motions, which can reorient forces in ways that boost thrust and reduce drag [21]. Overall, the efforts to refine the role of flexibility are still ongoing, and there is no clear explanation about the stiffness-tuning strategies in fish locomotion.
Several dynamic models of swimming fish have been developed to investigate the influences of body stiffness. Williams et al. [22] used the viscoelasticity beam model to study the deformed patterns of the fish body, which was actuated by different types of muscle actuation. McMillen et al. [23,24] simplified anguilliform fish as an elastic rod and a revised muscle model was proposed to produce a more realistic actuation force. Moreover, Tytell et al. [1] built a neuromechanical model of lamprey fish to study the interactions between internal forces, body stiffness, and fluid environment. Porter et al. [25] developed a dynamic model to study the force transmission along the vertebral column of dogfish sharks.
By comparison, the dynamic model of swimming fish, proposed by McMillen et al. [24], is very suitable for studying the stiffness-tuning strategies of anguilliform fish. Different from the previous model, the redundantly actuated planar series-parallel structure is adopted to demonstrate the biological structure of the fish body, in which the stiffness can be adjusted as needed. We intend to expand the musculo-mechanical model to explore the influences of various body stiffness on swimming performance in the hope of gaining insight into the way it is performed in real fish. Some control variables, such as the strength and frequency of muscle activation, may be enlarged or narrowed as possible, and we will look for patterns that could give insight into the right matching relations among the muscle actuation, the swimming performance, and the variable body stiffness. For example, we will measure the swimming speeds at different frequencies for a given level of muscle activation or at different levels of muscle activation for a given body stiffness. Although how fish use active muscle tensioning to tune their body stiffness to improve thrust and/or efficiency in real-time cannot be fully explained at the moment, we try to throw light on this problem in the present study.
The stiffness-tuning strategies in anguilliform fish can be implemented in next-generation fish-like robots or underwater vehicles, with an immediate improvement in forward speed and propulsive efficiency. The stiffness-tuning mechanisms and automatic stiffness control based on various swimming conditions are also promising applications. Moreover, while the stiffness-tuning strategies are poised to improve fish-like robots, they would offer insights into new studies of live fish locomotion, with the requirements of new challenging experiments, such as invasive surgery with implanting transducers and electrodes into the fish body, vivo stiffness measurements of various fish species, etc. The rest of this paper is organized as follows. The dynamic model of anguilliform fish is presented in Section 2. The computational results are demonstrated in Section 3, and the influences of various body stiffness on swimming performance are analyzed in Section 4. The discussions and conclusions are presented in Section 5 and Section 6, respectively.

2. Dynamic model of Anguilliform Fish

2.1. Series-Parallel Structure of Fish Body

In nature, lamprey fish have elongated bodies with narrow margins of caudal fins, and they undulate like a flapping ribbon. In the present study, the geometry of the fish body is measured from a real lamprey fish [26]. The cross-section shape of the fish body is approximated as elliptic, and the major and minor axes are labeled by a x and b x , respectively.
a x = 0.08 L x x 2 ,   0 x < 0.04 L   0.04 L 0.03 L x 0.04 L 0.91 L 2 ,   0.04 L x 0.95 L L x / 5 ,   0.95 L < x L  
b x = 0.08 L 1 x 0.51 L 0.51 L 2
In Equations (1) and (2), L is the full length of the fish body. In general, the motions of anguilliform fish occur only in a horizontal plane, and a two-dimensional fish model with constant height is acceptable. As shown in Figure 1, a planar serial–parallel mechanism with variable stiffness is proposed to model the biological structure of anguilliform fish.
As displayed in Figure 2, the body of anguilliform fish is modeled by a chain of interconnected planar serial-parallel structure units, which consists of rigid links, two springs, two dampers, and a revolute joint. The revolute joints link each vertebra into a continuing spine, representing the flexible notochord. The muscles are attached to the rigid vertebrae, which are modeled by a parallel structure of spring and damper. When the muscles are elongated or contracted, the actuation is imposed along two sides of the fish body, and the whole serial–parallel structure deforms like a flexible fish body. Additionally, the geometry parameters of real fish are adopted to design the serial–parallel structure, and the muscle actuations are determined by the muscle length and the velocity of its shortening.
The working principle of a single rotational serial–parallel unit is shown in Figure 2. The top platform A1A3A2 is supported by the middle rigid leg OB3 and the elastic legs A1B1 and A2B2. The elastic legs A1B1 and A2B2 on both sides connect the rotating pairs of upper revolute joints A1 and A2 on the top platform and the rotating pairs of lower revolute joints B1 and B2 on the fixed platform B1B3B2. The rigid middle leg OB3 is fixed on the fixed platform. The top platform rotates around the rotating center with a single degree of freedom. The muscle actuations change the internal forces resulting from the elastic legs A1B1 and A2B2, and the internal forces balance each other in the closed mechanism. The ribs A1A3 (or B1B3) and A2A3 (or B2B3) have the same length. The longitudinal length h is the spatial distance between the adjacent serial–parallel units.
For the discrete fish model, the dynamic equations of i th serial–parallel unit are expressed as
m i a i = h W i + F i , J i ψ ¨ i = M i + L i × F i ,  
where m i is the mass of i th serial–parallel unit; a i is the acceleration of serial–parallel unit; h W i is the hydrodynamic force applied on the i th serial–parallel unit; F i is the inner force transmitted from the adjacent structure unit; L i is the arm of force F i . J i is the inertia moment of serial–parallel unit; M i denotes the bending moment generated by muscles.
In Equation (4), the hydrodynamic force W is calculated by Taylor’s resistance model [27], and the components in the normal direction n and the tangential direction τ can be written as:
W = W n n + W τ τ W n = a ρ f v n v n + 8 ρ f a μ f v n v n W τ = 2.7 2 ρ f a μ f v n v τ  
where v n and v τ represent the normal and tangential velocity along the surface of the fish body, respectively; ρ f is the fluid density; μ f is the fluid viscosity.

2.2. Force Analyses

In Figure 2, the bending fish body is modeled by the planar serial-parallel mechanism, in which the joints connecting each serial-parallel unit of length h is actuated by a pair of spring-dashpot-actuators in parallel, anchored to arms of length w that project normally from the links’ midpoints. These arms represent myosepta, the connective tissue layers to which the muscle fibers connect. The linear springs and dashpots represent the tissue’s passive viscoelasticity, described by the spring constant k and the damping coefficient γ . The actuation forces f L i and f R i are generated by the contractile muscles on the right and left sides of i th planar serial-parallel unit, respectively. Suppressing the dependence on i and denoting the relative extensions Δ R and Δ L , the total forces on the right and left sides can be written as
G R t = f R t + k Δ R + γ Δ R ˙
G L t = f L t + k Δ L + γ Δ L ˙
when Δ R > 0 , Δ L > 0 , the springs are in tension, hence generating contractile forces. The forces are applied at a distance w from the centerline of the fish body, so elementary trigonometry gives:
Δ R = C D ¯ h h = 2 d c o s α h h = c o s ψ i 2 1 2 w h s i n ψ i 2
Δ L = A B ¯ h h = 2 d c o s α h h = c o s ψ i 2 1 + 2 w h s i n ψ i 2
where ψ i = φ i + 1 φ i is the angle between neighboring links. Finally, the moment arms L R , L L to the joint along normal from the lines A1B1 and A2B2 on which the forces act are computed as
L R = d sin α = d cos β ψ i / 2 = w c o s ψ i 2 + h 2 s i n ψ i 2
L L = d sin α = d cos β + ψ i / 2 = w c o s ψ i 2 h 2 s i n ψ i 2
For the i th serial–parallel unit with small rotation angles, the resulting torque at the joint is given by
M i = G R i L R i G L i L L i f L i t f R i t w + 2 υ w 2 h 2 4 f L i t + f R i t φ i + 1 φ i h + 2 γ w 2 φ ˙ i + 1 φ ˙ i h
The forces   f L i and f R i are calculated by the calcium kinetics model (see Section 2.3). To approximate a uniform distribution of the muscle, we set w x = b x / 2 , b x is the half-width of the fish body, as described in Equation (2). According to the muscle model, the generated actuation force is a function of muscle fibers; therefore, we propose that the stiffness k and damping γ are proportional to the cross-sectional area of the fish body. Thus we set
k = a b k ¯ , γ = a b γ ¯
where k ¯ and γ ¯ are the stiffness and damping of muscle fibers, with the units N/m2 and N s/m2, respectively. In particular, Young’s modulus E is used to describe the body stiffness, which is written in terms of the spring stiffness as E = 2 k ¯ / π .

2.3. Muscle Activation Model

According to the study of Hill [28], each myotome is modeled as a contractile element (CE) in series with an elastic element (SE). The CE and SE are located in series and experience equal forces at steady-state. The active force is generated by the CE, while the positive force is exerted by the SE, which is modeled as a linear spring. In reality, the stretch of the SE is not instantaneous due to the activation of the CE. The actuation force f R , L is dependent upon both the muscle length and its shortening velocity. The force density is proportional to muscle activation, scaled by the amplitude factor D that represents the force per unit muscle cross-section generated by unit activation:
f R , L = χ R , L t D
In Equation (14), χ R , L t   refers to muscle activation, and the subscripts R and L indicate the right and left side, respectively. The muscle activation χ t , in turn, produces contractile forces in the muscle fibers. In the present study, the model of calcium kinetics is adopted to describe the muscle activation χ R , L t , and it has been revised by the experiments of single myotomes of fish muscle [23]. The parameters k 1 ~ k 4 of calcium, the kinetics model is the same as the data in [24], and the intrinsic shape determined by the model of calcium kinetics is illustrated in Figure 3. The results show that the actuation force agreed well with the experimental data, and it is very close to the muscle force in refs. [22,23]. The results also suggested that the calcium kinetics model can accurately simulate muscle activation and force generation that occur during fish swimming.
In anguilliform fish, the musculature consists of approximately 100 myotomes, which are activated sequentially from just behind the gills towards the tail, alternating on the left and right sides of the body [22]. During steady swimming, the relative timing between muscle activation and body curvature is shown in Figure 4. The dashed lines along the fish body show the regions of muscle activation. The solid bars represent the electromyogram (EMG) activity measured at four different points along the left side of the body, plotted against time, and the sinusoidal lines are a measure of the local curvature of the body, which gives a measure of muscle length. This relative timing in Figure 4 is illustrated as a changing phase delay between the beginning of muscle stimulation and shortening. The relative timing between activation and curvature is considered in the present muscle model, which allows a larger frequency of activation, different from that in McMillen et al. [24]. Additionally, we change the amplitude factor D to study the influences of muscle actuation.

2.4. Solution of Fish Dynamic Model

In the fish dynamic model, the body is divided into N sections, and a chain of serial-parallel units with each of length h , with mass m i at each pivot, is redundantly actuated by active muscle forces. If the entire fish body is actuated, 6 ( N 1 ) first-order ordinary differential equations (ODEs) describe the muscle forces in the N -link chain, and with the 3 N second order ODEs, they jointly determine the body dynamics. The Newton iterative scheme, the spatial and temporal discrete methods provided in reference [23] (Appendix (2)∼(4) in [23]), are adopted to solve the dynamic equations. The boundary conditions at the fish head and tail is free, i.e., the force and moment are zero. The flow chart to solve the finite-difference discretization equations of the fish model is presented in Figure 5. The algorithm of ode45 in the Matlab software (Matlab R2022a) is used to solve the numerical equations, and the time step is 0.001 s.

3. Computational Results of Fish Undulatory Motions

3.1. The Variables of Fish Model

In the fish model, the shape of anguilliform fish is described by Equation (2). The body length is 20 cm, and it is discrete into 40 segments. The passive properties of the fish body are described by the density ρ , Young’s modulus E and the viscosity μ . The water viscosity μ f is 10−3 Pa·s, and the density ρ f is 1 g/cm3. The density of fish bodies is as same as that of water. Real fish have complex networks of muscle, cartilage, bone, skin, and organs, and they are more complicated than any physical fish model. Therefore, it is necessary to simplify anguilliform fish as a rod with heterogeneous viscoelastic materials, which can be described by the constant Young’s modulus and viscosity coefficient.
In the study of eel fish [29], Young’s modulus is 0.39 MPa, and the viscosity coefficient is 11 kPa·s. In Ref. [30], Young’s modulus of hagfish is 0.30 MPa, and the viscosity coefficients are 2 kPa·s and 5.7 kPa·s when the tail-beat frequencies are 3 Hz and 1 Hz, respectively. Therefore, it is reasonable to set the initial Young’s modulus E as 0.1 MPa, and the initial viscosity μ as 3 kPa·s. The fish head, occupying 1/10 of the body length, is without activation muscles, while the remaining 9/10 of the fish body that is capable of activation works as the activation region. To study the influence of various body stiffness, the Young’s modulus is varied in the possible range, and the variables of the fish model are listed in Table 1. The nomenclature of all parameters in the fish model is listed in Appendix A.

3.2. Swimming Performance of Anguilliform Fish

In our simulations, anguilliform fish starts from rest and accelerates to a quasi-steady state with maximum forward velocity along a straight way. The lateral and forward velocities of the virtual fish are displayed in Figure 6. The solid red line is the forward speed, and the dashed curve is the averaged forward speed over the tail-beat period. The dotted line is the lateral speed at the mass center of the fish body, with a sine wave with an oscillating frequency of 1 Hz. The mean asymptotic speed (MAS) is about 0.60 BL/s (body lengths per second). Besides, we divide the fish body into 100 segments and compare the resulting forward speeds. The results show that the speed obtained from a body of 40 sections differs by less than 0.1% from that of 100 sections.
According to the large-amplitude elongated body theory [31], swimming efficiency can be calculated by
η = 1 2 1 + U V
where U is the forward speed, and V is the wave speed. In the present case, the wave speed is 1.05 BL/s, and thus the swimming efficiency is 78.6%. The successive shapes of anguilliform fish in one tail-beat cycle are shown in Figure 7. The simulations readily yield results that are qualitatively similar to the movements of real fish. Color intensity on the fish body indicates the absolute value of the active bending moment. The red region shows the propagation curve of the active moment crest, while the blue line shows the propagation curve of the positive moment crest. The red and blue dashed lines near two sides of the body indicate the positive and negative curvature regions, respectively.
During one tail-beat cycle, the tracings of swimming fish that approximate its body centerline are demonstrated in Figure 8. The oscillating amplitude of the fishtail is obviously larger than that of the fish head. The undulatory motions in Figure 8a result from the dynamic fish model, which is similar to the tracings in Figure 8b from real fish [32]. These results demonstrate that the propagating wave can naturally emerge in a forced planar serial-parallel structure with redundant muscle actuations, which can be used to study the influences of various body stiffness.

3.3. Complex Modal Characteristics of Body Motions

As shown in Figure 8, it is very difficult to distinguish the undulatory motions of fish bodies accurately. In our previous study [33], the COD (complex orthogonal decomposition) method is adopted to describe the differences in the undulatory motions, including the oscillating amplitude of the fish head, the tail’s attack angle, the wavelength. For the integrity of this study, the COD method is briefly introduced. At first, the midline motions of fish body y x , t are dispersed to finite elements y x m , t n , x m = m L / M m = 1 ,   2 ,   ,   M is different positions along the fish’s body, and t n = n T / N   n = 1 ,   2   N is different time instants in one tail-beat period T . Secondly, a column vector ym(t) is constructed as [ym(t1) ym(t2) ym(t3)…ym(tn)]T, and then the Hilbert transform is used to transform ym(t) into the complex signals zm(t):
zm(t) = ym(t) + iH(ym(t)), i2 = −1.
where H(ym(t)) is the Hilbert transform of ym(t). At last, a complex matrix Z M × N = [z1(tn), z2(tn), ……, zm(tn)]T is constructed, and the eigenvalues λ m   and eigenvectors w m of matrix Z M × N are solved. The traveling index ξ is defined by
ξ = 1 / cond ( real w 1 , imag w 1
where cond() is the operating symbol of the condition number; w 1 is the first mode of the eigenvectors w m ; real w 1 is the real part of w 1 , and imag w 1 is the imaginary part of w 1 . The traveling index is used to describe the independence between the real and imaginary parts. When the real and imaginary vectors are dependent, the motions are pure standing waves, and the traveling index is 0; when the real and imaginary vectors are completely independent, the motions are pure traveling waves, and the traveling index is 1.
In the COD method, the undulatory motions are decomposed into two parts: pure standing waves and pure traveling waves, and the traveling index is defined by the ratio between the parts of the standing wave and the traveling wave. The midline motions in Figure 8 are decomposed into two parts, pure standing waves, and pure traveling waves, as shown in Figure 9, and the traveling index is 0.672. The results indicate that the midline motions of anguilliform fish are mixed by standing waves and traveling waves.

4. Swimming Performance of Anguilliform Fish

4.1. Fish Body with Fixed Body Stiffness

Biological evidence showed that the range of muscle activation was 0.1% to 20% of maximum isometric strength [22], so the magnitude of muscle activation is changed by timing an actuation coefficient D , which is also varied from 0.01 to 0.2. In this section, the swimming performance of anguilliform fish is studied when the stiffness of the fish body is fixed, and four different Young’s moduli, 0.006 MPa, 0.06 MPa, 0.6 MPa, and 1.0 MPa, are selected in this section. The relations among the tail-beat frequency, the amplitude of muscle actuation, and the swimming speed are shown in Figure 10. The tail-beat frequency is increased from 0.5 Hz to 5 Hz in the interval of 0.5 Hz, and the results of swimming fish with fixed body stiffness are showed as below:
(1)
For fixed body stiffness, the forward speed is increased with the tail-beat frequency and then decreased from the maximum value slowly. The maximum speeds are varied with the actuation amplitude and the tail-beat frequency. For the cases in Figure 10a, Young’s modulus of the fish body is 0.006 MPa, when the actuation coefficients are 0.03 and 0.18, the maximum swimming speeds are 0.5 BL/s and 1.2 BL/s, and the corresponding tail-beat frequency is 1 Hz and 2.5 Hz, respectively. As a whole, for different fixed body stiffness, the maximum speeds are increased with the amplitudes of muscle activation.
(2)
For different fixed body stiffness in Figure 10, the maximum forward speed is also increased with Young’s modulus of the fish body. For example, when the amplitude of muscle activation D is 0.03, the maximum speeds are 0.5 BL/s in Figure 10a (E = 0.006 MPa, f = 1 Hz) and 1.3 BL/s in Figure 10b (E = 0.06 MPa, f = 3 Hz), respectively. It is notable that when Young’s modulus of the fish body is small, the amplitude of muscle actuation has a more significant influence on the swimming speed. The influence would be weak, with the increase of the muscle actuations and body stiffness, in particular, when f ≥ 3 Hz, E ≥ 0.06 MPa.
(3)
As to the muscle actuation amplitude, it has a complex relationship to the swimming speed, which depends on Young’s modulus of the fish body and the tail-beat frequency. In Figure 10a, when the tail-beat frequency is smaller than 1.5 Hz, the influence of muscle actuations on the swimming speed is weak, while the influence is strengthened with the increase of tail-beat frequency. By comparison, the influence of the amplitude of muscle actuations on the swimming speed in Figure 10c is strong when the tail-beat frequency is smaller than 3 Hz, while it becomes weak with the increase of tail-beat frequency. Another interesting finding in Figure 10c is that if the tail-beat frequency f ≥ 3 Hz, the swimming speed is increased with the amplitude of muscle actuations, but if f < 3 Hz, the swimming speed is decreased with the amplitude of muscle actuations. This phenomenon is also appeared in Figure 10a,b,d, although the critical frequency varies.
With different muscle actuations, the influences of tail-beat frequency on swimming efficiency are shown in Figure 11. Similar to the results in Figure 10, the swimming efficiency is also closely related to the tail-beat frequency, the body stiffness, and the amplitude of muscle actuation. For most cases in Figure 11, The swimming efficiency is firstly increased to the maximum and then decreased with the tail-beat frequency. The maximum swimming efficiency is also determined by the body stiffness, the muscle actuation, and the tail-beat frequency, and the changing pattern is similar to that of maximum swimming speed. But for the small actuation, i.e., D = 0.03, the swimming efficiency is decreased with the tail-beat frequency. As a whole, the swimming efficiency of anguilliform fish with small Young’s modulus (E = 0.006 MPa in Figure 11a) is smaller than that with large Young’s modulus in Figure 11b–d. Besides, when Young’s modulus of the fish body is larger than 0.06 MPa, the effects of muscle activation amplitude on the swimming efficiency is smaller, while the tail-beat frequency has a larger effect.
Combined the results in Figure 10 and Figure 11, for anguilliform fish with fixed body Young’s modulus, when the muscle activation amplitude is small (D = 0.03 N/m3 in the present study), the dynamic model of anguilliform fish is under-actuated and the forward speed and swimming efficiency are decreased with the tail-beat frequency. When the muscle activation is strengthened, the fish body is actuated appropriately, and as the increases of tail-beat frequency, the maximum efficiency and the maximum speed can be obtained. It is notable that the maximum efficiency is appeared at first, and then the maximum speed is followed. For example, when D = 0.05 N/m3 and E = 0.006 MPa, the maximum swimming efficiency is reached at f = 1 Hz, while the maximum forward speed is obtained at f = 1.5 Hz. Further, when the fish body is over-actuated, we find that with the increase of tail-beat frequency, both the forward speed and swimming efficiency cannot be improved anymore.

4.2. Fish Body with Variable Body Stiffness

4.2.1. Influences on the Swimming Performance

In this section, Young’s modulus of the fish body is ranged from 10 kPa to 1 MPa, and the influences of variable body stiffness on swimming speed are shown in Figure 12. In Figure 12a, the tail-beat frequency is 1 Hz. For the cases D = 0.03 and 0.05, the swimming speed is increased with Young’s modulus of the fish body. However, for other cases, D ≥ 0.07, and the swimming speed is decreased with Young’s modulus of the fish body. In Figure 12b, the tail-beat frequency is 2 Hz, and for all the cases with different actuation amplitudes, the swimming speed is increased with Young’s modulus of the fish body. Further, we find that: (1) if E 0.01   MPa , the swimming speed is proportional to the actuation amplitude; (2) if E > 0.1 MPa , the swimming speed is inversely proportional to the actuation amplitude (except for D = 0.03). In Figure 12c–e, the forward speed is increased with both the actuation amplitudes and Young’s modulus of the fish body. It is notable that under the same actuation, the forward speed of anguilliform fish is decreased with the tail-beat frequency. For example, for the cases with D = 0.15, when the tail-beat frequency is 3 Hz, 4 Hz, and 5 Hz, the swimming speeds are 1.5 BL/s, 1.3 BL/s, and 0.8 BL/s, respectively.
As shown in Figure 13, the swimming efficiency of anguilliform fish is also varied with Young’s modulus of the fish body. The changing patterns of swimming efficiency are similar to the results of swimming speed in Figure 12. For the tail-beat frequency is 1 Hz in Figure 13a, when D = 0.03 or 0.05, the swimming efficiency is increased with Young’s modulus of the fish body, reaching a steady state. However, for other actuation amplitudes, D 0.07, the swimming efficiency is decreased with Young’s modulus of the fish body. For the tail-beat frequency 3 Hz, 4 Hz, and 5 Hz, the swimming efficiency is increased with Young’s modulus and the muscle actuation, as shown in Figure 13c–e. The differences among these figures include: (1) with the fixed actuation amplitude and Young’s modulus of the fish body, the swimming efficiency is highest at f = 3 Hz and then decreases with the tail-beat frequency. (2) for the case of 3 Hz, the swimming efficiency is not increased with the actuation amplitude when D ≥ 0.11. This phenomenon also appeared in Figure 13d,e.

4.2.2. Influences on the Midline Motions

The midline motions of anguilliform fish are also varied with Young’s modulus of fish body and tail-beat frequency, and they are investigated by the COD method. Meanwhile, the swimming speed, efficiency, and traveling index are used to evaluate the propulsive abilities of anguilliform fish. The results are presented in Figure 14 and Figure 15.
In Figure 14, the amplitude of muscle force D is 0.01, and the tail-beat frequency is f = 3 Hz. The changing pattern of swimming speed can be divided into three zones: low-speed zone, speed-increasing zone, and high-speed zone. In each zone, the typical Young’s modulus of the fish body, 0.004 MPa, 0.02 MPa, and 0.4 MPa, are selected to demonstrate the influences of variable stiffness on the midline motions. The undulatory motions of swimming fish are decomposed by the COD method, which is evaluated by the traveling index, as described in Section 3.3. As shown in Figure 14, although Young’s modulus of the fish body is varied from 0.004 MPa to 0.4 MPa, the basic pattern of undulatory motions are similar, i.e., the oscillating amplitude of the fish head and the tail are larger than that in the middle part. For the case of E = 0.004 MPa, the fish model is under-actuated, and the envelope amplitude of midline motions is smaller than in other cases (E = 0.02 MPa and 0.4 MPa). As a result, both the swimming speed and the swimming efficiency are lower, and the traveling index is also small. For the cases of E = 0.02 MPa and E = 0.4 MPa, the envelope amplitudes of midline motions are large enough, but there exist some differences, such as: (1) the oscillating amplitude of fish head is larger in the case of E = 0.02 MPa; (2) when E = 0.02 MPa, the lowest point of the envelope curves is located at about 0.2 L from the fish head, while there is no lowest point in the case of E = 0.4 MPa. Therefore, it is reasonable to predict that the large Young’s modulus of the fish body is easily prone to small head oscillations, and the envelope curviness of midline motions has a significant influence on the swimming efficiency.
The traveling index of midline motions is also varied with the tail-beat frequency. As shown in Figure 15, the amplitude of muscle force D is 0.01, and Young’s modulus of the fish body is fixed at 0.02 MPa. Three typical tail-beat frequencies, 1 Hz, 3 Hz, and 5 Hz, are selected to present their effects on the midline motions. When the tail-beat frequency is 1 Hz, it also means that the muscle activation frequency is relatively slow (as same as 1 Hz), which results in a large envelope amplitude of undulatory motions. The swimming efficiency is 80%, and the swimming speed is 0.62 BL/s. By comparison, when the tail-beat frequency is 3 Hz, the swimming performance is similar to the case of 1 Hz, although the envelope amplitude becomes small. However, for the case of 5 Hz, the muscle activation frequency is very fast, and the maximum amplitude of muscle activation (as shown in Figure 3) cannot be reached. The muscle activations cannot be fully used due to lacking enough reaction time, and thus the swimming efficiency is very low. Meanwhile, the traveling index is also very small.

5. Discussions

5.1. The Properties of Various Body Stiffness

To demonstrate the properties of various body stiffness, two different situations of anguilliform fish are selected: (A) Fixed body stiffness, E = 0.002 MPa, (B) Various body stiffness, Young’s modulus of the fish body ranged from 0.002 MPa to 1 MPa. Except for the body stiffness, the tail-beat frequency and the muscle actuations are the same. In detail, the tail-beat frequency is ranged from 0.5 Hz to 5 Hz, and four muscle actuations with different amplitudes, D = 0.03, 0.07, 0.11, and 0.15, are chosen. The comparisons of swimming speed and swimming efficiency are shown in Figure 16 and Figure 17, respectively.
Compared the swimming speeds of anguilliform fish with or without various body stiffness in Figure 16, we find that the most remarkable changes in the speed increment are appeared at the muscle actuation D = 0.03, while the mean increments of swimming speed are minimum when D = 0.15. Under certain muscle actuations, the speed increments are high when the tail-beat frequency is larger than 2.0 Hz. One distinctive feature is that when the muscle actuation is small, the influences of variable body stiffness are more obvious at lower tail-beat frequency, for example, in the cases of 1.5 Hz in Figure 16a,b. These results suggest that anguilliform fish may phone to vary their body stiffness to improve their swimming speed when the tail-beat frequency is large and/or the amplitude of muscle actuation is small.
For anguilliform fish with various body stiffness, the increments of swimming efficiency are shown in Figure 17. Similar to the results in Figure 16, for certain muscle actuation, the increments of swimming efficiency are large when the tail-beat frequency is between 1.5 Hz to 3.5 Hz, while the gap between cases A and B becomes small with the continuing increase of the tail-beat frequency. In Figure 17b–d, it is notable that the influence of variable stiffness on swimming efficiency is rather weak when the tail-beat frequency is lower than 1 Hz. Therefore, it is more effective for anguilliform fish to improve their swimming efficiency by varying the body stiffness at a certain range of tail-beat frequency. Additionally, more work should be performed to study the relationship between swimming efficiency and variable body stiffness because there are many approaches to calculating swimming efficiency.

5.2. The Effects of Hydrodynamic Forces

The hydrodynamic forces of anguilliform fish are modeled by the Taylor resistance theory in the present study. In this section, two different cases, with or without considering the hydrodynamic force, are selected to study the influences of the fluid force on the midline motions. In these cases, the amplitude of muscle actuation D is 0.01, and the tail-beat frequency is 1 Hz. Figure 18 shows the midline motion of anguilliform fish in a frame in which the center of mass is fixed. Different Young’s modulus E = 0.01 MPa, 0.05 MPa, 0.09 MPa are shown in three rows, and the left and right panels show the midline motions of the fish body with or without fluid forces, respectively. As a whole, the small curvature amplitude yields smaller motions in the central portions of the fish body, while the tail experiences smaller changes in curvature and tends to act as a ‘paddle’. The results show that the oscillation amplitudes of fish head and tail are obviously smaller than those without water force, which is consistent with the biology observations of anguilliform fish in reference [23]. Additionally, the results also show that the inclusion of hydrodynamic forces has a significant effect on the envelope shape of midline motions, although Young’s modulus is varied from 0.001 MPa to 0.09 MPa.
Further, the inclusion of hydrodynamic force also contributes to enlarging the traveling index. In our previous study [33], we collected 17 samples of anguilliform fish to analyze the range of traveling index, and the biological range is about 0.74~0.90. In present simulations, the range of the traveling index is 0.44–0.68, which is lower than the biological data. Besides, we compare the midline motions produced from fish dynamic model with those from real fish, and find that anguilliform fish deform its entire body with larger magnitude in nature, and the oscillation amplitude of the fish head is minimum. Therefore, it is reasonable to predict that the large oscillation amplitude of fish heads is harmful to the propulsive performance of anguilliform fish, and it also has a notable impact on the traveling index of midline motions.
From the perspective of dynamic analyses, the swimming fish can be regarded as a deformed viscoelastic rod immersed in water, and the deformation patterns of the fish body are dominated by their passive viscoelastic properties. Based on the viscoelasticity beam theory, the solution of this dynamic system can be a pure standing wave when the viscosity is zero (the dynamic model is a proportional damping system), and it could be a fish-like wave (a mixed standing and traveling wave) if the dynamic model is a non-proportional damping system. The strong damping modifies the system impedance and leads to an imbalance between the incident and reflected waves. This phenomenon was also proved by a series of bio-locomotion experiments [34], in which the transition mechanism from standing wave to traveling wave was analyzed. The wave propagation transforming from a standing wave to a traveling wave can be properly evaluated by the traveling index. That is the reason why we choose the traveling index to study the inner relations between the dynamic parameters and the undulatory motions of swimming fish. Therefore, the deformed patterns of anguilliform fish are depended on the stiffness, viscosity, and inertial mass. In the present study, the viscosity of a fish body is 1.0 kPa·s, so the influences of body viscosity and the interactive mechanism of muscle activation should also be considered in further work. In addition, we predict that both the stiffness and viscosity of the fish body are adjusted as needed, and these adjustments may be helpful to get more realistic midline motions of swimming fish.

6. Conclusions

In the present study, the stiffness-tuning strategies of anguilliform fish are studied by developing a dynamic model of swimming fish, with considerations of the morphological parameters, the muscle actuations, and the hydrodynamic forces. The properties of variable body stiffness have been analyzed in terms of its complex relations to muscle actuation, forward speed, and swimming efficiency. Moreover, the midline motions of swimming fish are decomposed by the COD method, and the traveling index is used to evaluate the differences in undulatory motions. The main conclusions are summarized below:
(1)
The maximum swimming speed and efficiency of anguilliform fish can be obtained by matching the body stiffness, the tail-beat frequency, and the muscle actuations appropriately. In general, when the body stiffness of a fish body is fixed, the swimming speed is increased to the maximum value and then decreased with tail-beat frequency. The maximum speed is increased with the amplitudes of muscle actuation and Young’s modulus of the fish body. The changing patterns of swimming efficiency are similar to that of swimming speed.
(2)
For the stiffness-turning strategies, anguilliform fish tend to vary their body stiffness to improve the swimming speed and efficiency at high tail-beat frequency. However, the influences of variable stiffness on swimming speed and efficiency are weakened when continuing the increase of tail-beat frequency. As to the swimming efficiency, it is also effective to vary their body stiffness when anguilliform fish are activated by the small amplitude of muscle actuations.
(3)
The midline motions produced by the dynamic model are similar to those of anguilliform fish observed from nature, but the range of the traveling index is lower than the biological range. It is partly because of the large oscillations of fish heads, and we also find the inclusion of hydrodynamic force is useful to reduce the head oscillations.
Overall, this study contributes to our understanding of how anguilliform fish obtain super swimming performance by adjusting the body stiffness and the muscle actuations, and it is also helpful to offer a general guidance to design the robotic fish with variable stiffness.

Author Contributions

Conceptualization, Z.C.; methodology, Z.C. and X.Z.; formal analysis, Z.C. and X.Z.; resources, Z.C.; writing—original draft preparation, X.Z.; writing—review and editing, Z.C.; project administration, Z.C.; funding acquisition, Z.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (No. 12002097, No. 12262006), the Scientific Start-up Project of Guizhou Institute of Technology (XJGC20190956), and the Natural Science Foundation of Guizhou Province (No. ZK[2021]-266, No. [2017]5789-20).

Data Availability Statement

The numerical and experimental data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The nomenclature of fish model parameters.
SymbolMeaning
Parameters of physical fish model L Full length of fish body
f Tail-beat (or actuation ) frequency
E Young’s modulus of fish body
ρ Density of fish body
μ Viscosity coefficient of fish body
η Swimming efficiency
U Forward speed
V Wave speed
λ a c t Actuation wavelength
D Amplitude coefficient of muscle actuation
a x Half-thickness of fish body
b x Half-width of fish body
x Distance along the axis of fish body
Parameters of Taylor’s resistance model W Hydrodynamic force
W n Normal component of hydrodynamic force
W τ Tangential component of hydrodynamic force
n Normal direction
τ Tangential direction
v n Normal velocity along the surface of fish body
v τ Tangential velocity along the surface of fish body
ρ f Fluid density
μ f Fluid viscosity
Parameters of dynamic fish model with planar serial-parallel structure h Length of serial-parallel unit
N Segments number of fish body
m i Mass of the i th serial–parallel unit
a i Acceleration of the i th serial–parallel unit
W i Hydrodynamic force applied on the i th serial–parallel unit
ρ f Density of fluid flow
μ f Viscosity coefficient of fluid flow
F i Inner force transmitted from the adjacent structure unit
L i The arm of force F i .
J i Inertia moment of the i th serial–parallel unit
M i Bending moment of the i th serial–parallel unit
f R Actuation force on the right side
f L Actuation force on the left side
f R i Actuation force on the right side of i th planar serial-parallel unit
f L i Actuation force on the left side of i th planar serial-parallel unit
G R t Total force on the right side
G L t Total force on the left side
L R Moment arm to the joint on the right side
L L Moment arm to the joint on the left side
Δ R Relative extension Δ R on the right side
Δ L Relative extension Δ L on the left side
ψ i Relative angle between the neighboring links
φ i Rotation angle of i th serial–parallel unit
k Spring constant of the tissue’s viscoelasticity
γ Damping coefficient of the tissue’s viscoelasticity
k ¯ Stiffness of muscle fibers
γ ¯ Damping of muscle fibers
χ R t Muscle activation on the right side
χ L t Muscle activation on the left side
k 1 ~ k 4 Parameters of calcium kinetics model
Parameters of COD method y x , t Midline motions of fish body
x m Different positions along fish body
t n Different time instants in one tail-beat period
T Tail-beat period
ξ Traveling index
λ m Eigenvalues of matrix Z M × N
w m Eigenvectors of matrix Z M × N
H( )Hilbert transform
cond( )Operating symbol of condition number
w 1 First mode of the eigenvectors w m
real w 1 Real part of w 1
imag w 1 Imaginary part of w 1 .

References

  1. Tytell, E.D.; Hsu, C.Y.; Williams, T.L.; Cohen, A.H.; Fauci, L.J. Interactions between internal forces, body stiffness, and fluid environment in a neuromechanical model of lamprey swimming. Proc. Natl. Acad. Sci. USA 2010, 107, 19832–19837. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Quinn, D.; Lauder, G. Tunable stiffness in fish robotics: Mechanisms and advantages. Bioinspir. Biomim. 2021, 17, 011002. [Google Scholar] [CrossRef] [PubMed]
  3. Long, J.H.; Nipper, K.S. The importance of body stiffness in undulatory propulsion. Am. Zool. 1996, 36, 678–694. [Google Scholar] [CrossRef]
  4. Fish, F.E.; Lauder, G.V. Control surfaces of aquatic vertebrates: Active and passive design and function. J. Exp. Biol. 2017, 220, 4351–4363. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Quinn, D.B.; Lauder, G.V.; Smits Lucas, K.N.; Johnson, N.; Tirrell, G.; Colin, S.P.; Gemmell, B.J.; Dabiri, J.O.; Costello, J.H. Bending rules for animal propulsion. Nat. Commun. 2014, 5, 3293. [Google Scholar]
  6. Wu, T.Y.-T. Hydromechanics of swimming propulsion. Part 1. Swimming of a two-dimensional flexible plate at variable forward speeds in an inviscid fluid. J. Fluid Mech. 1971, 46, 337–355. [Google Scholar] [CrossRef] [Green Version]
  7. Hua, R.-N.; Zhu, L.; Lu, X.-Y. Locomotion of a flapping flexible plate. Phys. Fluids. 2013, 25, 121901. [Google Scholar] [CrossRef]
  8. Quinn, D.B.; Lauder, G.V.; Smits, A.J. Scaling the propulsive performance of heaving flexible panels. J. Fluid Mech. 2014, 738, 250–267. [Google Scholar] [CrossRef] [Green Version]
  9. Liu, H.; Taylor, B.; Curet, O.M. Fin ray stiffness and fin morphology control ribbon-fin-based propulsion. Soft Robot. 2017, 4, 103–116. [Google Scholar] [CrossRef]
  10. Tytell, E.D.; Leftwich, M.C.; Hsu, C.Y.; Griffith, B.E.; Cohen, A.H.; Smits, A.J.; Hamlet, C.; Fauci, L.J. Role of body stiffness in undulatory swimming: Insights from robotic and computational models. Phys. Rev. Fluids. 2016, 1, 073202. [Google Scholar] [CrossRef]
  11. Tytell, E.D.; Hsu, C.Y.; Fauci, L.J. The role of mechanical resonance in the neural control of swimming in fishes. Zoology 2014, 117, 48–56. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. McHenry, M.J.; Pell, C.A.; Long, J.H. Mechanical control of swimming speed: Stiffness and axial wave form in undulating fish models. J. Exp. Biol. 1995, 198, 293–305. [Google Scholar] [CrossRef] [PubMed]
  13. Wolf, Z.; Jusufi, A.; Vogt, D.M.; Lauder, G.V. Fish-like aquatic propulsion studied using a pneumatically-actuated soft-robotic model. Bioinspir. Biomim. 2020, 15, 46008. [Google Scholar] [CrossRef]
  14. White, C.H.; Lauder, G.V.; Bart-Smith, H. Tunabot Flex: A tuna-inspired robot with body flexibility improves high-performance swimming. Bioinspir. Biomim. 2021, 16, 26019. [Google Scholar] [CrossRef]
  15. Pabst, D.A. Springs in swimming animals. Am. Zool. 1996, 36, 723–735. [Google Scholar] [CrossRef]
  16. Leftwich, M.C.; Tytell, E.D.; Cohen, A.H.; Smits, A.J. Wake structures behind a swimming robotic lamprey with a passively flexible tail. J. Exp. Biol. 2012, 215, 416–425. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Floryan, D.; Rowley, C.W. Clarifying the relationship between efficiency and resonance for flexible inertial swimmers. J. Fluid Mech. 2018, 853, 271–300. [Google Scholar] [CrossRef] [Green Version]
  18. Jusufi, A.; Vogt, D.M.; Wood, R.J.; Lauder, G.V. Undulatory swimming performance and body stiffness modulation in a soft robotic fish-inspired physical model. Soft Robot. 2017, 4, 202–210. [Google Scholar] [CrossRef] [Green Version]
  19. Shoele, K.; Zhu, Q. Leading edge strengthening and the propulsion performance of flexible ray fins. J. Fluid Mech. 2012, 693, 402. [Google Scholar] [CrossRef]
  20. Eldredge, J.D.; Jones, A.R. Leading-edge vortices: Mechanics and modeling. Annu. Rev. Fluid Mech. 2019, 51, 75–104. [Google Scholar] [CrossRef]
  21. Low, K.H.; Chong, C.W.; Zhou, C. Performance study of a fish robot propelled by a flexible caudal fin. In Proceedings of the 2010 IEEE International Conference on Robotics and Automation, Anchorage, AK, USA, 3–7 May 2010; pp. 90–95. [Google Scholar]
  22. Williams, T.L.; Bowtell, G.; Curtin, N.A. Predicting force generation by lamprey muscle during applied sinusoidal movement using a simple dynamic model. J. Exp. Biol. 1998, 201, 869–875. [Google Scholar] [CrossRef] [PubMed]
  23. McMillen, T.; Holmes, P. An elastic rod model for anguilliform swimming. J. Math. Biol. 2006, 53, 843–886. [Google Scholar] [CrossRef] [PubMed]
  24. McMillen, T.; Williams, T.; Holmes, P.; Daniel, T. Nonlinear muscles, passive viscoelasticity and body taper conspire to create neuromechanical phase lags in anguilliform swimmers. PLoS Comput. Biol. 2008, 4, e1000157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Porter, M.E.; Ewoldt, R.H.; Long, J.H. Automatic control: The vertebral column of dogfish sharks behaves as a continuously variable transmission with smoothly shifting functions. J. Exp. Biol. 2016, 219, 2908–2919. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Kern, S.; Koumoutsakos, P. Simulations of optimized anguilliform swimming. J. Exp. Biol. 2006, 209, 4841. [Google Scholar] [CrossRef] [Green Version]
  27. Taylor, G. Analysis of the swimming of long and narrow animals. Proc. R. Soc. Lond. A 1952, 214, 158–183. [Google Scholar]
  28. Hill, A. The abrupt transition from rest to activity in muscle. Philos. Trans. R. Soc. Lond. B 1949, 136, 399–420. [Google Scholar]
  29. Long, J.H., Jr. Muscles, elastic energy, and the dynamics of body stiffness in swimming eels. Am. Zool. 1998, 38, 771–792. [Google Scholar] [CrossRef]
  30. Long, J.H., Jr.; Koob-Emunds, M.; Sinwell, B.; Koob, T.J. The notochord of hagfish Myxine glutinosa: Visco-elastic properties and mechanical functions during steady swimming. J. Exp. Biol. 2002, 205, 3819–3831. [Google Scholar] [CrossRef]
  31. Lighthill, M.J. Large-amplitude elongated-body theory of fish locomotion. P. Roy. Soc. B-Biol. Sci. 1971, 179, 125–138. [Google Scholar]
  32. Bowtell, G.; Williams, T. Anguilliform body dynamics: Modelling the interaction between muscle activation and body curvature. Philos. Trans. R. Soc. Lond. B 1991, 334, 385–390. [Google Scholar] [CrossRef] [PubMed]
  33. Cui, Z.; Yang, Z.X.; Shen, L.; Jiang, H. Complex modal analysis of the movements of swimming fish propelled by body and/or caudal fin. Wave Motion. 2018, 78, 83–97. [Google Scholar] [CrossRef]
  34. Ramananarivo, S.; Godoy-Diana, R.; Thiria, B. Propagating waves in bounded elastic media: Transition from standing waves to anguilliform kinematics. Euro. Phys. Lett. 2014, 105, 54003. [Google Scholar] [CrossRef]
Figure 1. Dynamic model of anguilliform fish described by the redundantly actuated planar serial-parallel structures.
Figure 1. Dynamic model of anguilliform fish described by the redundantly actuated planar serial-parallel structures.
Biomimetics 08 00263 g001
Figure 2. Schematic of the redundantly planar serial-parallel structures with muscle actuation.
Figure 2. Schematic of the redundantly planar serial-parallel structures with muscle actuation.
Biomimetics 08 00263 g002
Figure 3. The intrinsic shape of actuation force generated by the muscle model of calcium ion.
Figure 3. The intrinsic shape of actuation force generated by the muscle model of calcium ion.
Biomimetics 08 00263 g003
Figure 4. Sketch of the electromyograms (EMG) signal passage and the maximum convex (red lines) and concave (blue lines) curvature of bending body during swimming cycles.
Figure 4. Sketch of the electromyograms (EMG) signal passage and the maximum convex (red lines) and concave (blue lines) curvature of bending body during swimming cycles.
Biomimetics 08 00263 g004
Figure 5. Solution flow chart to solve the dynamic model of anguilliform fish.
Figure 5. Solution flow chart to solve the dynamic model of anguilliform fish.
Biomimetics 08 00263 g005
Figure 6. Time history of the lateral and forward velocities of anguilliform fish.
Figure 6. Time history of the lateral and forward velocities of anguilliform fish.
Biomimetics 08 00263 g006
Figure 7. Undulatory motions of swimming fish during one tail-beat cycle (Red region: the active moment in fish body; Blue region: the positive moment in fish body; Red dashed lines on one side of body: the positive curvature region; Blue dashed lines near one side of body: negative curvature region).
Figure 7. Undulatory motions of swimming fish during one tail-beat cycle (Red region: the active moment in fish body; Blue region: the positive moment in fish body; Red dashed lines on one side of body: the positive curvature region; Blue dashed lines near one side of body: negative curvature region).
Biomimetics 08 00263 g007
Figure 8. The undulatory motions of anguilliform fish in one tail-beat period. (a) Centerlines of the fish model. (b) Tracings from real lamprey fish (reproduced from [32]).
Figure 8. The undulatory motions of anguilliform fish in one tail-beat period. (a) Centerlines of the fish model. (b) Tracings from real lamprey fish (reproduced from [32]).
Biomimetics 08 00263 g008
Figure 9. The decomposed parts of fish midline motions: (a) The part of pure standing wave; (b) The part of pure traveling wave.
Figure 9. The decomposed parts of fish midline motions: (a) The part of pure standing wave; (b) The part of pure traveling wave.
Biomimetics 08 00263 g009
Figure 10. The influences of tail-beat frequency on swimming speed under different muscle actuation amplitudes. (a) E = 0.006 MPa (b) E = 0.06 MPa (c) E = 0.6 MPa (d) E = 1 MPa.
Figure 10. The influences of tail-beat frequency on swimming speed under different muscle actuation amplitudes. (a) E = 0.006 MPa (b) E = 0.06 MPa (c) E = 0.6 MPa (d) E = 1 MPa.
Biomimetics 08 00263 g010aBiomimetics 08 00263 g010b
Figure 11. The influences of tail-beat frequency on swimming efficiency under different muscle actuation amplitudes. (a) E = 0.006 MPa (b) E = 0.06 MPa (c) E = 0.6 MPa (d) E = 1 MPa.
Figure 11. The influences of tail-beat frequency on swimming efficiency under different muscle actuation amplitudes. (a) E = 0.006 MPa (b) E = 0.06 MPa (c) E = 0.6 MPa (d) E = 1 MPa.
Biomimetics 08 00263 g011aBiomimetics 08 00263 g011b
Figure 12. The influences of Young’s modulus on swimming speed under different muscle actuation amplitudes. (a) f = 1 Hz (b) f = 2 Hz (c) f = 3 Hz (d) f = 4 Hz (e) f = 5 Hz.
Figure 12. The influences of Young’s modulus on swimming speed under different muscle actuation amplitudes. (a) f = 1 Hz (b) f = 2 Hz (c) f = 3 Hz (d) f = 4 Hz (e) f = 5 Hz.
Biomimetics 08 00263 g012
Figure 13. The influences of Young’s modulus on swimming efficiency under different muscle actuation amplitudes. (a) f = 1 Hz (b) f = 2 Hz (c) f = 3 Hz (d) f = 4 Hz (e) f = 5 Hz.
Figure 13. The influences of Young’s modulus on swimming efficiency under different muscle actuation amplitudes. (a) f = 1 Hz (b) f = 2 Hz (c) f = 3 Hz (d) f = 4 Hz (e) f = 5 Hz.
Biomimetics 08 00263 g013aBiomimetics 08 00263 g013b
Figure 14. The traveling index of midline motions varied with Young’s modulus of fish body.
Figure 14. The traveling index of midline motions varied with Young’s modulus of fish body.
Biomimetics 08 00263 g014
Figure 15. The traveling index of midline motions varied with the tail-beat frequency.
Figure 15. The traveling index of midline motions varied with the tail-beat frequency.
Biomimetics 08 00263 g015
Figure 16. The comparisons of the swimming speed of anguilliform fish with fixed or variable body stiffness. (a) D = 0.03 (b) D = 0.07 (c) D = 0.11 (d) D = 0.15.
Figure 16. The comparisons of the swimming speed of anguilliform fish with fixed or variable body stiffness. (a) D = 0.03 (b) D = 0.07 (c) D = 0.11 (d) D = 0.15.
Biomimetics 08 00263 g016
Figure 17. The comparisons of the swimming efficiency of anguilliform fish with fixed or variable body stiffness. (a) D = 0.03 (b) D = 0.07 (c) D = 0.11 (d) D = 0.15.
Figure 17. The comparisons of the swimming efficiency of anguilliform fish with fixed or variable body stiffness. (a) D = 0.03 (b) D = 0.07 (c) D = 0.11 (d) D = 0.15.
Biomimetics 08 00263 g017
Figure 18. Successive snapshots of fish midline motions. The left and right panels show the motions with and without fluid forces.
Figure 18. Successive snapshots of fish midline motions. The left and right panels show the motions with and without fluid forces.
Biomimetics 08 00263 g018
Table 1. The variables of fish dynamic model.
Table 1. The variables of fish dynamic model.
ParametersSymbolsValues
Length of fish body L 20 cm
Body segment number N 40
Water density ρ f 1.0 × 103 kg/m3
Viscosity coefficients of water μ f 1.0 mPa·s
Density of fish body ρ 1.0 × 103 kg/m3
Young’s modulus of fish body E 0.1 MPa
Viscosity coefficient of fish body μ 3 kPa·s
Actuation frequency f 1 Hz
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cui, Z.; Zhang, X. Computational Study of Stiffness-Tuning Strategies in Anguilliform Fish. Biomimetics 2023, 8, 263. https://doi.org/10.3390/biomimetics8020263

AMA Style

Cui Z, Zhang X. Computational Study of Stiffness-Tuning Strategies in Anguilliform Fish. Biomimetics. 2023; 8(2):263. https://doi.org/10.3390/biomimetics8020263

Chicago/Turabian Style

Cui, Zuo, and Xuyao Zhang. 2023. "Computational Study of Stiffness-Tuning Strategies in Anguilliform Fish" Biomimetics 8, no. 2: 263. https://doi.org/10.3390/biomimetics8020263

APA Style

Cui, Z., & Zhang, X. (2023). Computational Study of Stiffness-Tuning Strategies in Anguilliform Fish. Biomimetics, 8(2), 263. https://doi.org/10.3390/biomimetics8020263

Article Metrics

Back to TopTop