Abstract
Animal propulsors such as wings and fins bend during motion and these bending patterns are believed to contribute to the high efficiency of animal movements compared with those of man-made designs. However, efforts to implement flexible designs have been met with contradictory performance results. Consequently, there is no clear understanding of the role played by propulsor flexibility or, more fundamentally, how flexible propulsors should be designed for optimal performance. Here we demonstrate that during steady-state motion by a wide range of animals, from fruit flies to humpback whales, operating in either air or water, natural propulsors bend in similar ways within a highly predictable range of characteristic motions. By providing empirical design criteria derived from natural propulsors that have convergently arrived at a limited design space, these results provide a new framework from which to understand and design flexible propulsors.
You have full access to this article via your institution.
Similar content being viewed by others
Introduction
Thrust production by propulsive structures is a fundamental requirement for animal movement through air or water. Models of animal propulsion often emulate patterns that have been established for hydrofoils and airfoils1,2,3,4,5 and, as with these engineering-based models, animal propulsor models often assume a rigid structure that oscillates through the surrounding fluid with varying types of heaving and pitching motions6,7. However, animal propulsors are not rigid, but instead, typically bend during motion (Fig. 1). Attempts to evaluate the role of propulsor bending on thrust production have been met with conflicting results. For example, although some empirical8,9,10, modelling11,12,13 and computational studies14,15,16,17 have indicated that flexible propulsors may permit higher thrust production than rigid counterparts, other evidence indicates that flexibility may sometimes lessen or prevent thrust production10,18,19. Evaluation of these results has been complicated by a lack of common bending criteria to guide experimental design. Consequently, bending has been projected to occur both evenly over the full length17,20 and only over a limited portion21,22 of the propulsor. Most frequently, flexibility has been evaluated in terms of material properties such as elastic modulus10,17 or flexural stiffness19,23,24 of the propulsor rather than actual bending patterns. Experimental testing of the effects of bending on thrust production have been influenced by the inherent control possible with human-engineered propulsors so that most work on flexibility has been based on specifically engineered propulsors in laboratory conditions. Rarely have bending patterns of living animal propulsors, which have evolved over millions of years and among multiple taxonomic groups, been measured to inform experimental manipulations.
Here we describe, based on measurements of bending patterns across diverse animal taxa, universal patterns of animal propulsor bending, which transcend differences in fluid medium, animal size or phylogenetic background.
Results
Patterns of propulsor tip bending
Comparison of video sequences (see Supplementary Data 1 and 2) for 59 animal species demonstrates clear, replicable patterns of spanwise propulsor bending during steady motion, which are similar over a broad range of animal sizes, fluid media and taxonomic groups. The position of propulsor bending is documented as the flexion ratio (length from propulsor base to flexion point of bending relative to the total propulsor length). In addition, the maximum extent of propulsor bending is measured as the maximum flexion angle (referred to simply as flexion angle). Both flexion ratio and angle vary significantly between individual species (analysis of variance (ANOVA), n=59, df=58, F=6.7 for ratio, 8.8 for angle, P<0.001 for both variables) but only within constrained ranges for both flexion ratio and flexion angle (Fig. 2). The two variables form a discrete set of combinations (Fig. 3), or morphospace25,26,27, within the range of potentially available combinations. Although the ranges of both variables are strongly limited, there is no significant linear relationship between the two measures of propulsor bending (simple linear regression, N=59, df=1, 57, P=0.46, r2=0.009), indicating that each variable may be selected independently. The relative consistency of these two bending measures across such a wide array of animal propulsors suggests strong conservation of bending geometry during steady motion among animal taxa operating within inertially dominated fluid flows.
Statistical analysis integrating phylogenetic signal
Individual species (two to four individuals of each species, see Supplementary Data 2) were comparatively robust for their propulsor flexion traits. Replicate cycles of propulsor motion for particular individuals (2–15 cycles for each individual of a species) do not contribute significantly to variance of either flexion ratio (ANOVA, n=557, df=11, P=0.99) or angle (ANOVA, n=574, df=14, P=0.28). When replicates for each individual are grouped, individuals of a species do not contribute significantly to variations in flexion ratio (ANOVA, n=157, df=3, P=0.38) or angle (ANOVA, n=157, df=3, P=0.51). Central tendencies of the flexion variables are relatively insensitive to small (<30%) errors introduced into the data set by the generation of random errors within the measurements (Supplementary Fig. 1). Owing to the low variation found between propulsive cycles of a single individual and between individuals of the same species, values for individuals of the same species were combined and analysis focused on variations in bending patterns between species.
The distributions of interspecific flexion variables were generally normal within species categories (size, fluid medium and taxonomic groupings). With one exception (flexion angle for the group of animals a metre or more in length), the distributions of flexion ratio and angle values for species categories did not deviate significantly from normal distributions (Supplementary Table 1). Similarly, variances within species categories did not deviate significantly from homogeneity (Supplementary Table 2). The similarity in flexion ratio and angle distributions within species categories therefore satisfied the essential assumptions for comparison between the categories using ANOVA methods.
We considered the potential role of phylogenetic signal in bending patterns among animal groups because phylogenetic relationships may contribute to pattern formation for a range of variables among related taxa. Most broadly, phylogenetic signal is the tendency for evolutionarily related organisms to resemble each other without identification of the mechanism that might cause such resemblances. Phylogenetic signal has been detected at significant levels in comparative studies of animal traits, particularly when more than 20 species are compared28 and consequently might be expected within the flexion data. This would have an impact on our analysis if similarities of bending patterns between animal groups are confounded by phylogenetic relatedness between groups. In this case, rather than each species representing a separate and independent set of observations, species that were phylogenetically closer might be expected to have more similar patterns and would not be statistically independent, as assumed during statistical hypothesis testing. Blomberg’s28 K values above 1.0 demonstrate strong phylogenetic relatedness for the variable considered whereas values below 1.0 indicate that species resemble each other less than expected based solely on a Brownian motion model of trait evolution along the phylogenetic tree. Phylogenetic signal was relatively low (<1.0) but significantly present within the data for both flexion variables (flexion ratio: K=0.56, P=0.004; flexion angle: K=0.67, P<0.001, Supplementary Table 3). The presence of significant phylogenetic signal within the flexion data necessitated a statistical approach (phylogenetic ANOVA29) that incorporates this information.
Phylogenetic ANOVA indicated that differences in flexion traits between categorical variables (size class, fluid medium and taxonomic group) were not significant when considering the phylogenetic signal within the data (all P>0.1, Supplementary Table 4). Graphic analysis of means and confidence intervals (Fig. 4), illustrated differences in flexion variables that would be considered small, but significant, using a conventional one-way ANOVA approach. However, inclusion of genetic distance information (Supplementary Fig. 2, Supplementary Data 3) within a phylogenetic ANOVA demonstrated that when the genetic relatedness of the species within the data was considered, the differences between major categorical variations were not significant (Supplementary Table 4).
We include a conventional ANOVA analysis that does not consider genetic distances because, although the phylogenetic ANOVA is the appropriate hypothesis testing method for a situation involving significant phylogenetic signal, comparison of phylogenetic and conventional ANOVA approaches provides insight into patterns within the data. Graphical evaluation using 95% confidence intervals demonstrated a very limited morphospace for both flexion variables (Figs 2, 3, 4) and the phylogenetic ANOVA demonstrates that the minor variations between grouping variables can be explained by phylogenetic proximity of the taxa being compared rather than some important influence of size, medium or structural differences between the taxa. However, the utility of the conventional ANOVA lies with its implication that, within the constrained design space for propulsor flexion, small absolute differences are most likely to be found among groups with comparatively large phylogenetic distances. For example, propulsor flexion angle is, on average, lower for insects than for several other animal taxa (Fig. 4c, Supplementary Table 5). Insects are a comparatively basal group among those we considered (Supplementary Fig. 2). Similarly, within the vertebrates, fish are more likely to have lower flexion ratios than other taxa (Fig. 4c, Supplementary Table 5). The differences indicated by the conventional ANOVA do not result in large variations within the overall design spaces of the flexion variables (Fig. 4c) but reflect the phylogenetic distances over which these lineages have evolved during convergence on a restricted design space for propulsor bending.
Discussion
The remarkable consistency of animal bending patterns across animal sizes, fluid media or taxonomic groups leads to some unexpected results. For example, it is not coincidental that, during steady motion, the wings of an airborne monarch butterfly bend in a proportionately similar location and to a similar maximum angle as the tail of a swimming bottlenose dolphin or the wing-foot of the molluscan pteropod, Clione antarctica (Fig. 2). These are examples of common bending patterns that remain consistent across groups with exoskeletons (insects), endoskeletons (birds, bats, fish and cetaceans) or hydrostatic skeletons (molluscs). The propulsors may be actuated (birds, bats, fish, cetaceans and molluscs) or completely passive (insects). The material compositions of these propulsors vary dramatically (for example, chitin, feathers and bones) among taxa.
These data suggest that, rather than being determined by structural designs or material compositions, bending patterns of animal propulsors are constrained to a relatively narrow range of kinematic criteria. In turn, this suggests a key distinction between material properties (for example, elastic modulus, flexural stiffness, material resonance and joint patterns) and systems requirements (bending kinematics) of animal propulsors. Material properties are necessary conditions for flexible propulsors, but the universality of animal bending patterns demonstrates that the necessary material and structural requirements for propulsor bending have been resolved in many ways by many different animal taxa. Most importantly, the universality of the bending patterns demonstrates that a limited set of system requirements, primarily bending kinematics, are fundamentally important for flexible biological propulsion.
What factor(s) drive natural selection to converge on highly constrained bending kinematics across such a wide range of animal groups? Energetic efficiency of thrust production probably provides both advantages and constraints on bending designs. Propulsor margin bending can increase vortex circulation and generate pressure gradients22 around propulsors, which results in thrust enhancement of a factor >102 relative to rigid propulsors21. The energy savings accrued through maximizing energetic efficiency during thrust production may strongly favour evolutionary convergence on designs that produce such thrust advantages. Vortex circulation, thrust production and efficiency vary with flexibility, suggesting that the flexibility of animal propulsors may be tuned to maximize thrust production and efficiency30,31,32,33 by controlling vorticity associated with propulsor bending. For example, although some bending is favourable, extensive bending can lessen thrust production efficiency18 by disrupting vortex organization along the span of a propulsor10. Similarly, specific spatial and temporal alignments of adjacent opposite-spin vortices at the tips of propulsor margins are crucial for production of fluid jets produced at vortex interfaces34,35. As propulsor geometries and kinematics control spatial and temporal distributions of vorticity36,37, propulsor geometries and kinematics are key determinants of thrust production in fluids. Even subtle bending kinematics influence alignment of vortices, which determine the energetic efficiency of propulsion38, and may be subject to strong natural selection. Small variations in bending patterns among taxonomic groups resemble model predictions that weak, but significant, phylogenetic signal may reflect strong stabilizing selection within constrained boundaries acting on phylogenetically disparate species groups39. The bending pattern morphospace demonstrates a degree of design latitude, but the generally narrow range of natural propulsor flexion patterns indicates strong selection during the evolution of animal propulsors towards a tightly constrained range of bending kinematics. The wide diversity of propulsor material and actuation traits following these design rules demonstrate that a diversity of properties have been modified by different animal taxa operating in different fluids to reach similar geometric arrangements of propulsor bending.
The narrow range of observed natural propulsor geometries parallels other documented kinematic patterns of natural propulsors. Steady propulsion has been found to conform to a Strouhal number (St=fA/U, where (f) stroke frequency and (A) amplitude are divided by (U) forward speed) range of ~0.2–0.4 for animals moving within inertially dominated fluid flows40. Like propulsor bending patterns, St conservation occurs uniformly across wide animal length scales, fluid media and taxonomic groups40,41. Also, similar to propulsor bending patterns, St conservation is believed to permit efficient generation and control of vortices, which results in maximum thrust efficiency per propulsor energy input42,43. Constrained propulsor bending and Strouhal number conservation may be related. However, further exploration of the specific fluid interactions underlying these universal patterns is necessary to clarify our understanding of animal propulsion and provide essential design rules for creating successful biomimetic propulsors. Human-engineered propulsors may then benefit from patterns developed over millions of years of evolutionary selection within the animal kingdom.
Methods
Data sources and quantification of propulsor flexion
Video sequences of animals (N=59 species) in steady-state swimming or flight were collected from various online sources (for example, YouTube, Vimeo and so on) using YouTube Downloader freeware (DVDVideoSoft). Internet addresses from which videos were accessed are given in Supplementary Data 1 and the complete data set, including all replicate sample values, is given in Supplementary Data 2. Animals were selected based on their use of oscillatory motions of their propulsors. Animals with highly undulatory motions such as eels and animals augmenting their forward motion with other appendages (such as the walking limbs of a reptile) were not considered. While there are several lineage-dependent patterns of propulsor bending associated with internal body structure (for example, joint structures of birds, bats, fish and mammals), our goal was to characterize bends solely due to propulsor flexibility, so we used the flexion point of the bend closest to the propulsor tip for all animal lineages. The selected animals represented divergent taxa (bats, birds, insects, fish, cetaceans and molluscs), and the propulsors comprised many different materials (skin membrane, feathers, chitin, bone and so on). Animal size range included ~0.005 m Limacina helicina, a pteropod mollusc (Re~40), up to ~8.4 m orca whale (Orcinus orca) (Re~107). Animals were selected based on availability of video, such that at least two (range 2–4) individuals of the same species were represented, with each of these individuals undergoing multiple (average 3, range 2–15) propulsive cycles (Supplementary Data 2). Rapid acceleration events often included transient, extensive propulsor bending and were not included in the analysis. Instead, sequences were chosen to represent steady-state conditions in which propulsor motions were relatively even between consecutive propulsor movement cycles.
Two variables were measured to characterize the kinematics of the propulsor margin: the location (flexion ratio) and the magnitude (flexion angle) of bending (Fig. 2). These features were measured on video frames using the application ImageJ (NIH). Flexion angle was defined as the angle of the margin’s bend away from the central propulsor axis. To obtain this angle, the obtuse angle between the flexion point and distal propulsor tip was measured using ImageJ. This obtuse angle was the supplement of the acute flexion angle (illustrated with examples for each major taxonomic group in Fig. 5). The flexion ratio was defined as the length along the propulsor from the propulsor base to the flexion point relative to the total length of the propulsor. High flexion ratios (close to 1.0) indicated a flexion point located more distally along the propulsor (Figs 2 and 5).
Measurements throughout a particular propulsion cycle required conversion of video sequences into a series of images representing successive stages of propulsor motion using Adobe Premiere Pro CS3 software. Alternate frames were used for each image series, except where propulsor oscillations occurred so quickly that all frames were used to provide adequate resolution of motions. The flexion angle was measured on each still image extracted from a video sequence. As the magnitude of the flexion angle varied throughout the propulsive cycle, only the maximum amplitude of flexion angle in each cycle was considered during further analyses. Flexion ratio was measured multiple times per propulsive cycle (average 3.5, s.d. 1.4) and the maximum flexion angle was measured over multiple propulsive cycles (average 2.9 cycles, s.d. 2.2) for each replicate individual of a species. The full range of bending variation over a pulsation cycle was evaluated by taking measurements throughout pulsation cycles (average 9.9, s.d. 7.5 measurements per cycle).
Propulsor tip bending occurred over a range of propulsor orientations relative to oncoming fluid flow. Therefore, both the streamwise tip bending of birds, bats, insects and molluscs, and the spanwise bending of fish and cetaceans were measured (Fig. 5). Accurate measurement of the flexion angle optimally required an orthogonal view of the propulsor and this was a key criterion for video sequence selection. For fish, this meant that we used only views looking directly upwards or downwards on the fish. For close-up sequences of flying animals, only directly orthogonal views were used. With distant views, commonly used for birds, only sequences in which the bird was flying towards or away from the camera were selected. The flexion data set does not have completely parallel flexion ratio and angle data—not all video sequences had identical numbers of replicates. This is because image sequences were not always suitable for collection of both types of data in the same scene. Consequently, sample sizes for each species were not uniform (Supplementary Data 2).
Owing to the diverse sources of videos used for flexion measurements, we evaluated the sensitivity of flexion measurements to measurement variations such as might accompany deviations from orthogonality. Random corruption of the original flexion data demonstrated that central tendencies of the flexion variables were relatively robust. Although variance around mean values increased, average values for flexion ration and angle remained stable when random errors of 10, 20 and 30% were introduced into the original data set (Supplementary Fig. 1). This pattern indicates that small (<30%) errors introduced into the data set during variable measurements or by minor deviations from orthogonality of images would be unlikely to alter central patterns within the data set.
Phylogenetic signal evaluation
A phylogenetic hypothesis of relationships among taxa and species used for propulsor flexion analysis was generated based on mitochondrial cytochrome c genetic data available through GenBank (Supplementary Data 3) using MEGA44. The resulting phylogenetic tree (Supplementary Fig. 2) was consistent with contemporary phylogenetic hypotheses describing animal evolution45,46. Phylogenetic distances from that tree were then used to calculate Blomberg’s K and its level of significance28 for both experimental variables.
The presence of significant phylogenetic signal influenced the choice of ANOVA used to compare flexion variables between animal groups. Significant phylogenetic signal within flexion traits (Supplementary Table 3) indicated that average values for phenotypes did not represent independent or random data points. This is true within a restricted lineage as well as within a larger grouping involving multiple animal lineages. Consequently, degrees of freedom may be inflated and significance levels derived from conventional tests may become unreliable when there is significant phylogenetic signal present28. However, phylogenetic signal can be included within ANOVA procedures to permit comparative analyses that account for evolutionary relatedness of traits among clades. Consequently, we used a phylogenetic ANOVA29 to analyse patterns of flexion kinematics among animal groups.
Additional information
How to cite this article: Lucas, K. N. et al. Bending rules for animal propulsion. Nat. Commun. 5:3293 doi: 10.1038/ncomms4293 (2014).
References
Lighthill, J. Mathematical Biofluiddynamics Society for Industrial and Applied Mathematics (1975).
Pennycuick, C. J. in Avian Biology (eds D. S. Farner & J. R. King) 1–75 (Academic Press, 1975).
Blake, R. Fish Locomotion Cambridge Univ. Press (1983).
Videler, J. J. Fish Swimming Chapman & Hall London (1993).
Alexander, R. M. Principles of Animal Locomotion Princeton University Press (2003).
Triantafyllou, M. S., Triantafyllou, G. S. & Gopalkrishnan, R. Wake mechanics for thrust generation in oscillating foils. Phys. Fluids 3, 2835–2837 (1991).
Anderson, J. M., Streitlien, K., Barrett, D. S. & Triantafyllou, M. S. Oscillating foils of high propulsive efficiency. J. Fluid Mech. 360, 41–72 (1998).
Bozkurttas, M., Dong, H., Mittal, R., Maden, K. P. & Lauder, G. V. Hydrodynamic Performance of Deformable Fish Fins and Flapping Foils 44th AIAA Aerospace Sciences Meeting (AIIA Paper, Reno, Nevada, 2008).
Toomey, J. & Eldredge, J. D. Numerical and experimental study of the fluid dynamics of a flapping wing with low order flexibility. Phys. Fluids 20, 073603 (2008).
Heathcote, S., Wang, Z. & Gursul, I. Effect of spanwise flexibility on flapping wing propulsion. J. Fluid Struct. 24, 183–199 (2008).
Mittal, R., Dong, H., Bozkurttas, M., Lauder, G. B. & Madden, P. Locomotion with flexible propulsors: II. Computational modeling of pectoral fin swimming in sunfish. Bioinsp. Biomim. 1, S35–S41 (2006).
Alben, S. Optimal flexibility of a flapping appendage in an inviscid fluid. J. Fluid Mech. 614, 355–380 (2008).
Michelin, S. & Smith, S. G. L. Resonance and propulsion performance of a heaving flexible wing. Phys. Fluids 21, 071902 (2009).
Vanella, M., Fitzgerald, T., Preidikman, S., Balaras, E. & Balachandran, B. Influence of flexibility on the aerodynamic performance of a hovering wing. J. Exp. Biol. 212, 96–106 (2009).
Young, J., Walker, S. M., Bomphrey, R. J., Taylor, G. K. & Thomas, A. L. R. Details of insect wing design and deformation enhance aerodynamic function and flight efficiency. Science 325, 1549–1552 (2009).
Le, T. Q., Ko, J. H., Byun, D., Park, S. H. & Park, H. C. Effect of chord flexure on aerodynamic performance of a flapping wing. J. Bionic Eng. 7, 87–94 (2010).
Kang, C.-K., Aono, H., Cesnik, C. E. S. & Shyy, W. Effects of flexibility on the aerodynamic performance of flapping wings. J. Fluid Mech. 689, 32–74 (2011).
Liu, P. & Bose, N. Propulsive performance from oscillating propulsors with spanwise flexibility. Proc. R. Soc. Lond. A 453, 1763–1770 (1997).
Tangorra, J. L., Lauder, G. V., Hunter, I. W., Mittal, R., Madden, P. G. A. & Bozkurttas, M. The effect of fin ray flexural rigidity on the propulsive forces generated by a biorobotic fish pectoral fin. J. Exp. Biol. 213, 4043–4054 (2010).
Miao, J. M. & Ho, M. H. Effect of flexure on aerodynamic propulsive efficiency of flapping flexible airfoil. J. Fluid. Struct. 22, 401–419 (2006).
Villanueva, A., Smith, C. & Priya, S. A biomimetic robotic jellyfish (Robojelly) actuated by shape memory alloy composite actuators. Bioinspir. Biomim. 6, 036004 (2011).
Colin, S. P. et al. Biomimetic and live medusae reveal the mechanistic advantages of a flexible bell margin. PLoS One 7, e48909 (2012).
Combes, S. A. & Daniel, T. L. Shape, flapping and flexion: wing and fin design for forward flight. J. Exp. Biol. 204, 2073–2085 (2001).
Combes, S. A. & Daniel, T. L. Flexural stiffness in insect wings II. Spatial distribution and dynamic wing bending. J. Exp. Biol. 206, 2989–2997 (2003).
Raup, D. M. & Michelson, A. Theoretical morphology of the coiled shell. Science 166, 994–995 (1965).
McGhee, G. R. Theoretical Morphology: The Concept and Its Applications Columbia University Press (1999).
Costello, J. H., Colin, S. P. & Dabiri, J. O. Medusan morphospace: phylogenetic constraints, biomechanical solutions, and ecological consequences. Invert. Biol. 127, 265–290 (2008).
Blomberg, S. P., Garland, T. Jr & Ives, A. R. Testing for phylogenetic signal in comparative data: behavioral traits are more labile. Evolution 57, 717–745 (2003).
Garland, T. Jr, Dickerman, A. W., Janis, C. M. & Jones, J. A. Phylogenetic analysis of covariance by computer simulation. Syst. Biol. 42, 265–292 (1993).
Ramananarivo, S., Godoy-Diana, R. & Thiria, B. Rather than resonance, flapping wing flyers may play on aerodynamics to improve performance. Proc. Natl Acad. Sci. USA 108, 5964–5969 (2011).
Wilson, M. M. & Eldredge, J. D. Performance improvement through passive mechanics in jellyfish-inspired swimming. Int. J. Non-Linear Mech. 46, 557–567 (2011).
Kang, C.-K. & Shyy, W. Scaling law and enhancement of lift generation of an insect-size hovering flexible wing. R. Soc. Interface 10, 20130361.
Zheng, L., Hedrick, T. L. & Mittal, R. Time-varying wing-twist improves aerodynamic efficiency of forward flight in butterflies. PLoS One 8, e53060 (2013).
Ahlborn, B., Harper, D., Blake, R., Ahlborn, D. & Cam, M. Fish without footprints. J. Theor. Biol. 148, 521–533 (1991).
Gopalkrishnan, R., Triantafyllou, M. S., Triantafyllou, G. S. & Barrett, D. S. Active vorticity control in a shear flow using a flapping foil. J. Fluid Mech. 274, 1–21 (1994).
Colin, S. P., Costello, J. H., Katija, K., Seymour, J. & Kiefer, K. Propulsion in Cubomedusae: mechanisms and utility. PLoS One 8, e56393 (2013).
Videler, J., Müller, U. K. & Stamhuis, E. J. Aquatic vertebrate locomotion: wakes from body waves. J. Exp. Biol. 202, 3423–3430 (1999).
Gemmell, B. J. et al. Passive energy recapture in jellyfish contributes to propulsive advantage over other metazoans. Proc. Natl Acad. Sci. 110, 17904–17909 (2013).
Revell, L. J., Harmon, L. J. & Collar, D. C. Phylogenetic signal, evolutionary process and rate. Syst. Biol. 57, 591–601 (2008).
Triantafyllou, M. S., Triantafyllou, G. S. & Yue, D. K. P. Hydrodynamics of fish swimming. Annu. Rev. Fluid Mech. 32, 33–53 (2000).
Taylor, G. K., Nudds, R. L. & Thomas, A. L. R. Flying and swimming animals cruise at a Strouhal number tuned for high power efficiency. Nature 425, 707–711 (2003).
Triantafyllou, M. S., Hover, F. S., Techet, A. H. & Yue, D. K. P. Review of hydrodynamic scaling laws in aquatic locomotion and fishlike swimming. Appl. Mech. Rev. 58, 226–237 (2005).
Wolfgang, M. J., Anderson, J. M., Grosenbaugh, M. A., Yue, D. K. P. & Triantafyllou, M. S. Near-body flow dynamics in swimming fish. J. Exp. Biol. 202, 2303 (1999).
Tamura, K. et al. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 28, 2731–2739 (2011).
Stuart, G. W., Moffett, K. & Leader, J. J. A comprehensive vertebrate phylogeny using vector representations of protein sequences from whole genomes. Mol. Biol. Evol. 19, 554–562 (2001).
Dunn, C. W. et al. Broad phylogenomic sampling improves resolution of the animal tree of life. Nature 452, 745–749 (2008).
Kembel, S. W. et al. Picante: R tools for integrating phylogenies and ecology. Bioinformatics 26, 1463–1464 (2010).
Harmon, L. J., Weir, J. T., Brock, C. D., Glor, R. E. & Challenger, W. GEIGER: investigating evolutionary radiations. Bioinformatics 24, 129–131 (2008).
R Core Team. R: A Language and Environment for Statistical Computing R Foundation for statistical computing: Vienna, Austria, ISBN 3-900051-07-0 (2012).
Acknowledgements
This work was supported by an ONR MURI award to J.H.C. (N000140810654) and National Science Foundation Grant OCE 0623508 to J.H.C., S.P.C. and J.O.D. We thank C. Syslo and K. Kiefer for assistance with video collection and A. Shomberg for statistical advice. We also thank photographers for use of their images in Fig. 1: (a) E. Schipul ( http://www.flickr.com/photos/eschipul), (b) M. Baird (Michael L. Baird, flickr.bairdphotos.com), (d) J. Fabrice, ( fabrice@marinemegafauna.org), (e) D. Wrobel, ( http://www.wrobelphoto.com), (f) E. Van 3000 from in Flanders fields—Belgiquistan—United Tribes of Europe (CC-BY-SA-2.0). We also thank photographers for use of their images in Fig. 5: G. Humes, The US Fish and Wildlife Service, B. Nyman, B. Weinstein, OpenCage and E. Van 3000.
Author information
Authors and Affiliations
Contributions
K.N.L., N.J., E.C. and G.T. collected the data; K.N.L., S.P.C., B.J.G., J.O.D., J.H.C. and W.T.B. analysed the data; and K.N.L. and J.H.C. wrote the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Supplementary information
Supplementary Figures and Tables
Supplementary Figures 1-2 and Supplementary Tables 1-5 (PDF 494 kb)
Supplementary Data 1
Video source locations and access dates. The source types refer to video sequences or size information about the species. Time information (minutes, seconds, playback frame number) are listed for the start location of each individual measured within a species. (XLS 123 kb)
Supplementary Data 2
Flexion data derived from video sequences. Medium refers to fluid medium (air or water), Size Class refers to the size groupings listed in Fig. 2, Group refers to the taxonomic group of the species. Abbreviations: Ind - Individual, Rep - Replicate. Individuals represent individual members of a species; replicates were replicated cycles for an individual of a species. (XLS 112 kb)
Supplementary Data 3
Accession information for GenBank (National Center for Biotechnology Information, US NIH) cytochrome c genetic data used to construct a phylogenetic tree for species used in the flexion propulsor data set. Surrogate species were used when target species data was absent. The order of surrogate species was based on phylogenetic proximity and chosen in rank of genus, family, infraclass, subclass, class and order until a surrogate species for the target species was found. (XLS 25 kb)
Rights and permissions
About this article
Cite this article
Lucas, K., Johnson, N., Beaulieu, W. et al. Bending rules for animal propulsion. Nat Commun 5, 3293 (2014). https://doi.org/10.1038/ncomms4293
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/ncomms4293
This article is cited by
-
Design and build of small-scale magnetic soft-bodied robots with multimodal locomotion
Nature Protocols (2024)
-
Escape performance in the cyclopoid copepod Oithona davisae
Scientific Reports (2024)
-
Numerical study and parameter prediction of a full-free flexible flament in wakes of flapping wing
Acta Mechanica Sinica (2024)
-
Numerical study of the correlation between fish school arrangement and propulsive performance
Artificial Life and Robotics (2023)
-
The locomotory apparatus and paraxial swimming in fossil and living marine reptiles: comparing Nothosauroidea, Plesiosauria, and Chelonioidea
PalZ (2021)