Next Article in Journal
A Generalization of the Hausdorff Dimension Theorem for Deterministic Fractals
Next Article in Special Issue
Nonlinearly Preconditioned FETI Solver for Substructured Formulations of Nonlinear Problems
Previous Article in Journal
Time-Dependent Reliability Analysis of Plate-Stiffened Prismatic Pressure Vessel with Corrosion
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phi-Bonacci Butterfly Stroke Numbers to Assess Self-Similarity in Elite Swimmers

by
Cristiano Maria Verrelli
1,*,
Cristian Romagnoli
2,3,
Roxanne Jackson
4,
Ivo Ferretti
5,
Giuseppe Annino
3,6,7 and
Vincenzo Bonaiuto
3
1
Department of the Electronic Engineering, University of Rome Tor Vergata, via del Politecnico 1, 00133 Rome, Italy
2
Department for Life Quality Studies, University of Bologna, 47900 Rimini, Italy
3
Sport Engineering Lab, Department of Industrial Engineering, University of Rome Tor Vergata, 00133 Rome, Italy
4
Institute of Automation Technology, Otto-von-Guericke University Magdeburg, 39106 Magdeburg, Germany
5
Biomechanical and Video-Analysis Area for the National Teams of Federazione Italiana Nuoto (FIN), 00100 Rome, Italy
6
Department of Medicine Systems, University of Rome Tor Vergata, 00133 Rome, Italy
7
Centro di Biomedicina Spaziale, University of Rome Tor Vergata, 00133 Rome, Italy
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(13), 1545; https://doi.org/10.3390/math9131545
Submission received: 21 May 2021 / Revised: 13 June 2021 / Accepted: 22 June 2021 / Published: 1 July 2021
(This article belongs to the Special Issue Advanced Numerical Methods in Computational Solid Mechanics)

Abstract

:
A harmonically self-similar temporal partition, which turns out to be subtly exhibited by elite swimmers at middle distance pace, is formally defined for one of the most technically advanced swimming strokes—the butterfly. This partition relies on the generalized Fibonacci sequence and the golden ratio. Quantitative indices, named ϕ -bonacci butterfly stroke numbers, are proposed to assess such an aforementioned hidden time-harmonic and self-similar structure. An experimental validation on seven international-level swimmers and two national-level swimmers was included. The results of this paper accordingly extend the previous findings in the literature regarding human walking and running at a comfortable speed and front crawl swimming strokes at a middle/long distance pace.

1. Introduction

Fibonacci numbers play crucial roles in combinatorial mathematics and elementary number theory. Although such numbers have been investigated for centuries, they continue to intrigue mathematicians and researchers in many areas of human endeavor (as can be seen in [1] and the references therein), while providing new tools for expanding the frontiers of mathematical study. In respect to this, the golden ratio  ϕ 1.61803 , namely the positive solution to the equation x 2 = 1 + x , appears in some very fundamental relationships involving numbers, with one of the most basic occurrences of the golden ratio involving the use of two seed values and a simple Fibonacci-like additive recursion relationship.
Walking and running are human gait modes exhibiting different mechanics and energetics. A double support phase, i.e., both limbs are in ground contact, identifies the walking gait, whereas a double float phase, i.e., no limb is in ground contact, identifies the running gait. However, both physiological (symmetric and recursive) human walking and running are characterized, from a temporal point of view, by only four specific time intervals, associated with the durations of gait cycle, swing, stance and double support (double float) phases. More precisely, physiological symmetric walking (running) is classically recognized to exhibit:
  • A stance (swing) duration, t ST ( t SW ), being close to 62% of gait cycle duration t GC ;
  • A swing (stance) duration being close to 38% of gait cycle duration;
  • A double support (double float) duration, t DS ( t DF ), being consequently close to 24% of gait cycle duration.
As recently formally recognized in [2], the above sequence { 0.24 , 0.38 , 0.62 , 1 } can be viewed as a slight approximation of the special sequence in walking:
Σ : { t DS / t GC , t SW / t GC , t ST / t GC , 1 } = { 1 / ϕ 3 , 1 / ϕ 2 , 1 / ϕ , 1 }
and its conceptual counterpart in running:
Σ : { t DF / t GC , t ST / t GC , t SW / t GC , 1 } = { 1 / ϕ 3 , 1 / ϕ 2 , 1 / ϕ , 1 } .
The above sequence Σ is nothing but a generalized four-length Fibonacci sequence [3] (with the real numbers 1 / ϕ 3 , 1 / ϕ 2 as seeds), generally defined, for any values γ and δ , by ( ψ is the negative solution to x 2 = 1 + x ):
U n = γ φ n + δ ψ n
satisfying the recurrence:
U n = γ φ n 1 + δ ψ n 1 + γ φ n 2 + δ ψ n 2 = U n 1 + U n 2 .
Regarding Σ : (i) the general Fibonacci sequence structure is rooted in the following duration constraints of symmetric walking (and running, respectively): t DS + t SW = t ST (respectively, t DF + t ST = t SW ), t SW + t ST = t GC ; (ii) the specific ϕ -dependent (temporally harmonic) self-similar structure relies on the special chain of ratios t SW / t DS = t ST / t SW = t GC / t ST (respectively, t ST / t DF = t SW / t ST = t GC / t SW ) all being equal to ϕ . Indeed, 1 / ϕ 3 0.23608 , 1 / ϕ 2 0.38198 , 1 / ϕ 0.61804 . In other words, conditions in (i)–(ii) apparently induce a repeated self-proportional partition (namely self-similarity), in accordance with the fact that two quantities are in the golden ratio if their ratio is the same as the ratio of their sum to the larger of the two quantities. This way, walking and running implicitly involve a fractal nature, in which the structure of the larger scale resembles the structure of the sub-unit and in which one of the simplest ways of transformation, i.e., the new domain composed of two previous ones, is highlighted.
Very recent experimental and theoretical analyses—inspired by the aforementioned cyclic human movements in walking and running (see also [4,5]—have found harmonic structures to even appear in front crawl swimming [6] at a middle/long distance pace, at which the presence of redundant movements is reduced via in-phase synchronization with the induced waves in the water. With respect to this, the recent paper [6] not only provides a mathematical framework and experimental consistency for recognizing preliminary evidence, for elite front-crawl swimmers swimming at a middle/long distance pace, a recovery phase duration that is very close to the fundamental unit 3 u f (where u f is equal to the duration of the front crawl stroke divided by 12) but it also illustrates the existence of harmonic structures—in both their simple and enhanced versions—that elite swimmers seek to reproduce in order to improve their performance. We also refer the reader to [7,8] for a different scenario involving the dimensionless Strouhal number.
A more technically advanced swimming stroke is the butterfly. Simultaneous strokes, such as breaststroke and butterfly, are in fact considered to be highly technical due to the complex coordination of the arm and leg actions [9]. In particular, in butterfly, the out-of-water arms recovery is facilitated by the leg undulation. To be effective, however, the kick must appear as a consequence of a body wave-like cefalo–caudal undulation motion [10]. Having a glide time with the arms extended forward at the top of the stroke is certainly a strategy used to reduce the energetic cost (metabolic power/velocity) for long-distance swims [11]. This way, the head, trunk and upper limbs are profiled in a streamlined position in order to glide, and consequently, provide some rest at each stroke during a long butterfly swim (some aquatic animals such as giant cormorants, penguins, and dolphins improve the metabolic efficiency of swimming by adopting locomotion patterns with alternating periods of propulsion and gliding [10]). However, this is not effective for achieving a higher stroke rate [10] and for avoiding great instantaneous velocity fluctuations. On the other hand, when velocity and stroke rate increase, coordination becomes closer to an in-phase mode [12,13], just like in human locomotion and in quadrupedal coordination. Indeed, the significant skill effect in [10] indicates that elite swimmers—whom are the object of analysis in this paper—have stronger arm/leg synchronization than the sub-elite swimmers: elite swimmers adopt a shorter glide to overcome great forward resistance and generate higher forces during the arm pull; sub-elite swimmers often compensate for coordination mistakes by applying greater force. With respect to this, we also refer the reader to [14]—for devices performing propulsion analysis in swimming; [15]—for a quantitative evaluation of phases of turns during competition; and [16]—for the role of the hip movement in the stroke mechanics.
The aim of this paper was to extend, for the first time in the literature—to the best of our knowledge—the aforementioned findings regarding human walking and running at comfortable speed and front crawl swimming strokes at a middle/long distance pace: a harmonically self-similar temporal partition, which relies on the generalized Fibonacci sequence and the golden ratio, is formally defined for the highly complex and upper and lower-limbs-coordinated butterfly stroke. Quantitative indices, named ϕ -bonacci butterfly stroke numbers, are accordingly proposed to assess such a hidden time-harmonic and self-similar structure being subtly exhibited by elite swimmers at middle distance pace.

2. Generalized Fibonacci Sequences in Butterfly Stroke

This section presents the main result of the paper.

2.1. Butterfly Stroke Phases

Different arm and leg phases (see Figure 1) can be identified in a butterfly stroke S with duration t S , here characterized by two leg undulations for one arm stroke [10,12,17]. They are in order, for the arms:
(i)
Entry and catch phase (between the entry of the hands into the water and the beginning of their backward movement), with duration t EC ;
(ii)
Pull phase (between the beginning of the backward movement of the hands and their entry into the plane vertical to the shoulders), with duration t PL ;
(iii)
Push phase (between the positioning of the hands below the shoulders and their exit from the water), with duration t PS ;
(iv)
Recovery phase (between the arrival of the hands at the water level and their subsequent entry into the water), with duration t RE ,
and for the legs:
(i)
Downward phase 1 (between the high and low break-even points– at first occurrence–of the feet during the first undulation), with duration t K 1 ;
(ii)
Upward phase 1 (between the low and high break-even points–at first occurrence–of the feet during the first undulation), with duration t U 1 ;
(iii)
Downward phase 2 (between the high and low break-even points–at first occurrence–of the feet during the second undulation), with duration t K 2 ;
(iv)
Upward phase 2 (between the low and high break-even points–at first occurrence–of the feet during the second undulation), with duration t U 2 .
The inter-limb coordination of such phases is governed by the time delays (with a positive or negative sign):
  • T 1 , between the start of Entry and catch phase and the start of Downward phase 1;
  • T 2 , between the start of Pull phase and the start of Upward phase 1;
  • T 3 , between the start of Push phase, and the start of Downward phase 2;
  • T 4 , between the start of Recovery phase, and the start of Upward phase 2.
As shown in [10], elite swimmers show a smaller total time gap, taken as the sum—expressed as the percentage of a complete stroke—of the absolute values of T 1 , T 2 , T 3 , T 4 and used to assess the effectiveness of the global arm–leg coordination.

2.2. Self-Similarity and Enhanced Self-Similarity in Butterfly Stroke

We first define the repetitive butterfly stroke, in which T 1 of the subsequent stroke equals T 1 of the stroke under investigation, so that the arm stroke (namely from the start of Entry and catch phase to the start of the subsequent Entry and catch phase) equals the leg stroke (namely, from the start of Downward phase 1 to the start of the subsequent Downward phase 1). A repetitive butterfly stroke that is characterized by its values ( T 1 / t S = a [ % ] , T 2 / t S = b [ % ] , T 3 / t S = c [ % ] , T 4 / t S = d [ % ] ) is denoted hereafter by ( a , b , c , d ) -repetitive butterfly stroke. We then define an ( a , b , c , d ) -repetitive butterfly stroke that is kick-to-kick temporal symmetric by imposing the equality of the Upward phase 1 and Upward phase 2 durations:
t U 1 = t U 2 + T 1 .
Furthermore, a highly coordinated  ( a , b , c , d ) -kick-to-kick temporal symmetric repetitive butterfly stroke is an ( a , b , c , d ) -kick-to-kick temporal symmetric repetitive butterfly stroke that exhibits a relatively small positive value for the normalized total delay:
N T D = a 2 + b 2 + c 2 + d 2 [ % ] .
In light of [6], a self-similar partition in highly coordinated ( a , b , c , d ) -kick-to-kick temporal symmetric repetitive butterfly strokes is presented hereafter.
Theorem 1. 
( F 4 ):In a highly coordinated ( a , b , c , d ) -kick-to-kick temporal symmetric repetitive butterfly stroke, the sequence:
F 4 : { t K 1 + t K 2 , t U 1 , t K 1 + t K 2 + t U 1 , t S }
is a generalized Fibonacci sequence of length 4. If t U 1 / ( t K 1 + t K 2 ) = ϕ , then the chain of inequalities:
t K 1 + t K 2 + t U 1 t U 1 = t S t K 1 + t K 2 + t U 1 = ϕ
holds and the sequence F 4 in (7) has a self-similar structure, with t K 1 + t K 2 + t U 1 61.804 % of t S , t U 1 38.198 % of t S and t K 1 + t K 2 23.608 % of t S .
Proof. 
As in [6], the sequence F 4 is represented through the discrete-time system:
y ( k + 2 ) = y ( k + 1 ) + y ( k ) , k = 0 , 1 ,
with y ( 0 ) = t K 1 + t K 2 , y ( 1 ) = t U 1 and y ( 2 ) = t S t U 2 , y ( 3 ) = t S . Then, it suffices to write the solution to its state-space description:
ξ ( k + 1 ) = M ξ ( k ) , k = 0 , 1
as
ξ ( 1 ) = M ξ ( 0 ) ,
ξ ( 2 ) = M 2 ξ ( 0 ) ,
with ξ ( l ) being the vector [ y ( l ) , y ( l + 1 ) ] T , l = 0 , 1 , 2 , and M = [ 0 , 1 ; 1 , 1 ] with distinct eigenvalues ϕ , 1 ϕ and orthogonal eigenvectors ( c ϕ = ( 1 + ϕ 2 ) 1 / 2 ):
v ϕ = c ϕ [ 1 , ϕ ] T , v ( 1 ϕ ) = c ϕ [ ϕ , 1 ] T .
The following theorem characterizes a stronger self-similarity property, referred to as enhanced self-similarity. It relies on the following constraint (When t K 2 > t K 1 , it suffices to switch t K 1 and t K 2 in the C -constraint, in the sequence of Corollary 1 (including Remark ( F 7 )), in the corresponding index (24), and in the fourth comment Section 4):
C : t U 1 = t K 1 + 2 t K 2
Theorem 2. 
( F 6 ):In a highly coordinated ( a , b , c , d ) -kick-to-kick temporal symmetric repetitive butterfly stroke under C , the sequence (7) enforced with { t K 1 , t K 2 } to the left, meaning that:
F 6 : { t K 1 , t K 2 , t K 1 + t K 2 , t U 1 , t K 1 + t K 2 + t U 1 , t S } ,
is a generalized Fibonacci sequence of length 6. The equality t U 1 / ( t K 1 + t K 2 ) = ϕ makes the sequence F 6 possess an enhanced self-similar structure, with t K 1 + t K 2 + t U 1 61.804 % of t S , t U 1 38.198 % of t S , t K 1 + t K 2 23.608 % of t S , t K 2 14.591 % of t S , t K 1 9.0175 % of t S .
Proof. 
Sequence F 6 can be represented by the discrete-time system:
y ( k + 2 ) = y ( k + 1 ) + y ( k ) , k = 0 , 1 , 2 , 3 ,
with initial conditions y ( 0 ) = t K 1 , y ( 1 ) = t K 2 and y ( 2 ) = t K 1 + t K 2 , y ( 3 ) = t U 1 , y ( 4 ) = t K 1 + t K 2 + t U 1 , y ( 5 ) = t S . The previous analysis carried out in the proof of Theorem 1 can be extended to cover the generalized 6-length Fibonacci sequence case. In particular, the state-space representation of the system above reads:
ξ ( k + 1 ) = M ξ ( k ) , k = 0 , 1 , 2 , 3 ,
with ξ ( l ) = [ y ( l ) , y ( l + 1 ) ] T , l = 0 , 1 , 2 , 3 , 4 , and M as in the proof of Theorem 1. We now explicitly write:
y ( 2 ) y ( 3 ) = ϕ 2 β v ϕ + ( 1 ϕ ) 2 α v ( 1 ϕ )
y ( 3 ) y ( 4 ) = ϕ 3 β v ϕ + ( 1 ϕ ) 3 α v ( 1 ϕ )
y ( 4 ) y ( 5 ) = ϕ 4 β v ϕ + ( 1 ϕ ) 4 α v ( 1 ϕ ) ,
in terms of:
α = ξ ( 0 ) , v ( 1 ϕ ) = ( ϕ y ( 0 ) y ( 1 ) ) c ϕ
β = ξ ( 0 ) , v ϕ = ( y ( 0 ) + ϕ y ( 1 ) ) c ϕ = ( 1 + ϕ 2 ) c ϕ y ( 0 ) .
When α is equal to zero, i.e., the initial vector ξ ( 0 ) has no components along the v ( 1 ϕ ) - direction, y ( 1 ) / y ( 0 ) = ϕ holds and y ( 2 ) / y ( 1 ) = ϕ , y ( 3 ) / y ( 2 ) = ϕ , y ( 4 ) / y ( 3 ) = ϕ , y ( 5 ) / y ( 4 ) = ϕ hold too.
Remark 1. 
( F 7 ):In the very special case of T 1 = t K 2 t K 1 for the highly coordinated ( a , b , c , d ) -kick-to-kick temporal symmetric repetitive butterfly stroke, the sequence (15), once enforced with T 1 to the left, becomes a generalized Fibonacci sequence of length 7, with the equality t U 1 / ( t K 1 + t K 2 ) = ϕ making the sequence F 7 possess a strongly enhanced self-similar structure; T 1 5.5735 % of t S is additionally obtained.

2.3. Quantitative Measures of Self-Similarity and Enhanced Self-Similarities

Two indices I f , 4 , I f , 6 —named ϕ-bonacci butterfly stroke number and enhanced ϕ-bonacci butterfly stroke number—can be naturally introduced, in order to quantitatively assess self-similarity or enhanced self-similarities of highly coordinated ( a , b , c , d ) -kick-to-kick temporal symmetric repetitive butterfly strokes. They are, in order:
I f , 4 = 100 [ ( t K 1 + t K 2 + t U 1 ) / t S 0.62 2 + t U 1 / t S 0.38 2 + ( t K 1 + t K 2 ) / t S 0.24 2 ] 1 / 2
I f , 6 = 100 [ ( t K 1 + t K 2 + t U 1 ) / t S 0.62 2 + t U 1 / t S 0.38 2 + ( t K 1 + t K 2 ) / t S 0.24 2 + t K 2 / t S 0.14 2 + t K 1 / t S 0.10 2 ] 1 / 2 ,
where the values 0.62 , 0.38 , 0.24 , 0.14 , 0.10 are used to approximate 0.61804 , 0.38198 , 0.23608 , 0.14591 , 0.090175 of Theorem 1 and Theorem 2, respectively. The smaller such indices are, the stronger the self-similarity is. The just self-similarity of Theorem 1 may clearly lead to a non-zero value for I f , 6 that tends towards zero when the self-similarity tends to turn into the enhanced self-similarity.

3. Experimental Analysis

The feasibility of the preceding analysis is here illustrated by the dedicated analysis of butterfly stroke training sessions for: (i) seven international-level swimmers, namely IL1 (male, 31 y, 190 cm, 80 kg), IL2 (female, 31 y, 170 cm, 65 kg), IL3 (female, 27 y, 168 cm, 58 kg), IL4 (male, 20 y, 196 cm, 80 kg), IL5 (male, 19 y, 180 cm, 73 kg), IL6 (male, 19 y, 193 cm, 85 kg), IL7 (female, 27 y, 173 cm, 66 kg); (ii) two national-level swimmers, namely, NL1 (female, 19 y, 166 cm, 55 kg), NL2 (male, 25 y, 185 cm, 95 kg), all of them swimming at their own middle distance pace. In particular, the above international-level swimmers and national-level swimmers are reported, within the corresponding sets, in order of physical shape (measured as race performance capabilities) at the moment of data acquisition. While the international-level swimmers IL1–IL7 compete at major international events on a regular basis and hold national/international records, the national-level swimmers NL1–NL2 are national medalists who compete on a regular basis of major national events (with NL1 being close to becoming an international-level swimmer). The analysis was performed by using high frame rate videos (100 for IL1, NL2; 120 fps for IL2–IL7, NL1) of stable strokes via the 2D BioMovie ERGO system at http://www.infolabmedia.eu/ (accessed on 7 June 2021).

3.1. Phase Durations and Interlimb Coordination

Phase and delay durations for all the swimmers IL1–IL7, NL1–NL2 are reported in Table 1, Table 2, Table 3 and Table 4: ( a , b , c , d ) -(almost) kick-to-kick temporal symmetric repetitive butterfly strokes S (under constraints C ) are exhibited (Kick-to-kick temporal symmetry appears to be almost verified for IL1–IL7, NL1–NL2, with the modulus of the difference between t U 1 and t U 2 + T 1 belonging to the set [ 5 , 67 ] ms. Constraint C even appears to be almost verified for IL1–IL7, NL1–NL2, with the modulus of the difference between t U 1 and t K 1 + 2 t K 2 belonging to the set [ 1 , 81 ] ms.), characterized by the ( a , b , c , d ) - and NTD-values reported in Table 3 and Table 4. In accordance with the evidence of [10] on elite swimmers, all of such relatively small NTDs define highly coordinated strokes. All nine swimmers’ strokes exhibit a negative delay T 2 (leading to a lag time in glide position) and a positive delay T 3 . On the other hand, in contrast to the international-level swimmers IL1, IL3–IL6, the international-level swimmers IL2, IL7 and the national-level swimmers NL1, NL2 are characterized by a small negative superposition of two contradictory actions ( T 1 < 0 ). Furthermore, T 4 = 0 for NL1, whereas T 4 is negative for IL1–IL7 and NL2.

3.2. Self-Similarity Analysis

The aggregate phase percentage values for the international-level swimmers IL1–IL7 and the national-level swimmers NL1–NL2 are reported in Table 5 and Table 6, along with the values for the indices I f , 4 and I f , 6 in (23) and (24). Comments are in order:
  • Rather small values are obtained for IL1–IL7, with IL1’s one being the smallest, owing to the strict closeness of the corresponding phase percentage values to 62 % , 38 % , 24 % , 14 % , and 10 % ;
  • Relatively large reductions in the self-similarity and enhanced self-similarity magnitudes (especially of the latter) appear for the national-level swimmers NL1–NL2 when compared to the international-level swimmers IL1–IL7;
  • The I f , 4 - and I f , 6 - values turn out to reproduce the order of physical shape within the two swimmers’ set;
  • IL5 even presents an a = T 1 / t S -value ( 4.71 % ) that is close to the one ( 5.57 % ) characterizing the strongly enhanced self-similar structure of Remark ( F 7 );
  • The slight differences in the phase durations of Table 1 and Table 2 (compare, for instance, IL2 to NL1, or IL4 to IL5, or IL2 to IL3), which lead to the differences in self-similarity magnitudes of Table 5 and Table 6, have been successfully identified via the high frame rate analysis used in this paper, with the self-similarity information complementing the delay partition values of Table 3 and Table 4;
  • Larger percentage reductions in enhanced self-similarity (with respect to self-similarity) are exhibited by IL2, IL3, IL5—when compared to IL1, IL4, IL5, IL6, NL1, NL2—so that the I f , 6 - values for IL2–IL3 and IL5 tend to thicken (more than the I f , 4 -ones) towards the I f , 6 - values for IL4 and IL6, respectively.

4. Discussion

By following up the direction highlighting the connections of human walking and running gaits and front crawl swimming strokes with the generalized Fibonacci sequences and golden ratio, this paper aimed to provide an answer to the question regarding the existence of subtle and unveiled self-similar structures at the root of the butterfly stroke of elite swimmers, for whom—as in dolphins and fish [18]—energy accrued by raising the upper body is transmitted along the body through wave-like movements [19]. Harmonically characterizing butterfly strokes through self-similarity—in both its simple and enhanced versions—and the ϕ -bonacci butterfly stroke numbers constitutes the original contribution of this paper, as a relevant application of computational number theory involving generalized Fibonacci sequences.
When compared to the front crawl scenario of [6] (concerning the arm stroke):
  • Theorem 1 and Theorem 2 provide partition constraints that regard the leg stroke, with the inter-delay composition linking the arm phase partition with the leg phase partition;
  • The four-length generalized Fibonacci sequence F 4 of Theorem 1 is generated by t K 1 + t K 2 , which plays the role of the recovery phase duration in the front crawl of [6] and of the double support (double float) duration in walking (running) of [2];
  • The durations t U 1 and t U 2 + T 1 play the role of t F a and t F b in the front crawl stroke in [6] and of the right and left swing (right and left stance) durations in the asymmetric walking (running) of [2,5], with the involved equality between the durations of such phases simply transposing the swing (stance)-symmetry constraint of symmetric walking (running);
  • In light of the twelve-tone equal temperament-based interpretation, constraint C imposes that the stroke duration t S of an ( a , b , c , d ) -kick-to-kick temporal symmetric repetitive butterfly stroke equals the sum of the durations of eight disjoint sub-phases, three among them with duration t K 1 and five among them with duration t K 2 : according to [6], notes D4, E4, G4, C5 correspond, in the (suspended and restored) C-variant Cadd2-chord, to the frequencies: f 3 ( D 4 ) = f 0 2 3 12 = 293.7 Hz, f 5 ( E 4 ) = f 0 2 5 12 = 329.6 Hz, f 8 ( G 4 ) = f 0 2 8 12 = 391.96 Hz, f 13 ( C 5 ) = f 0 2 13 12 = 2 f 1 ( C 4 ) = 523.2 Hz ( f 0 = 246.9175 Hz), with n = 3 , 5 , 8 , 13 corresponding to elements of a Fibonacci sequence and with the ratios 5 / 3 1.667 , 8 / 5 = 1.6 , 13 / 8 = 1.625 of consecutive elements of such sequence quickly getting close to ϕ ;
  • The above results—again occurring in elite swimmers—confirm that, in contrast to walking, a precise swimming technique (coming from a considerable amount of practice and instruction repeating specifically precise movements for a long enough time, while producing rhythmic motor patterns through the interrelationship between cortical input, central pattern generator (CPG), and sensory feedback) is relevantly involved (the intra-cyclic variation of the horizontal velocity of the hip is larger in non-expert swimmers than their expert counterparts [20]) as recalled by [6], the fractal dimension is 1.8 1.9 for highly qualified expert swimmers [20], whereas it is 1.1 1.4 for on-land walking [21].

5. Practical Implications

Self-similarity, as defined in this paper, may become a reference point in the formulation of advanced training programs (on the basis of (23), (24)) for elite swimmers, while: (i) going beyond the evaluation of the inter-limb synchronicity level, classically used to assess the performance of high-level butterfly stroke swimmers; and (ii) not affecting the peculiar stroke duration of the swimmer. It is clear that improvements of the muscle strength remain valid as well, but at a different level of analysis. The findings of this paper might also provide comparative (quantitative) information about: (i) the physical recovery level of the swimmer after a hard and intense workout; (ii) the conditioning of the rhythmic neural patterns of the swimmer during the contingent-moment-based swimming actions; and (iii) the swimmers’ improvements over time and the necessary correction actions to improve the performance. Indeed, at a high-level, at which inter-individual coordination variability reduces to a narrow range of movement variants [20], each small change that a practitioner may perform might lead to a relevant change in the dynamic interaction of the systems components [22].

6. Strengths and Limits

The derivations of this paper seem to provide meaningful information about both the swimmer’s level and physical shape, which complement the analysis of the normalized total delay. The results of this paper also seek to open new research directions aiming to connect computational number theory (involving generalized Fibonacci sequence) to human neural system modulating principles at the root of highly effective cyclic rhythmic patterns of swimming. Anyway, the results of the present study should be interpreted with caution, owing to its specific limits, such as: (i) the small size of the samples enrolled in the experiments; and (ii) the dependence on high-frame rate analyses. Many aspects should be further investigated in future studies. One of them relies on the fact that the present analysis is focused on the proportions among the stroke phases, limiting this type of analysis to temporal features. Furthermore, it could be certainly interesting to put self-similarity in relationship with the role played by different sensory information in constructing and maintaining the aforementioned stroke harmony.

7. Forecast

It is worth noticing how Figure 3 of [23] (200 m-pace, high-level male swimmer) seems to respect such a self-similar partition in breaststroke, with the arm recovery and the out-sweep phases in place of t F a , t F b and the arm in-sweep phase in place of the front crawl stroke recovery phase. On the other hand, we recognized the same self-similar partition in the fourth stroke of the start phase (which accounts for a greater percentage of the race success [24]) in elite sprint kayaking contests, once the sequence generated by min { t Ext , t Pull } , t Imm + max { t Ext , t Pull } is considered under the experimentally verified constraint t Rec ( 2 t Imm + t Pull + t Ext + max { t Ext , t Pull } ) , in terms of the paddle immersion, pull, paddle extraction and recovery phases of the stroke. Analogously, the reader can appreciate in Table 2 and Figure 2 of [25] how, in the forward swing of the forehand tennis stroke (cross court and down the line situations) of elite tennis players, the same self-similar partition can be found owing to: (i) a temporal distance between the instants of maximum racquet rotation angle and maximum hip linear velocity that is almost equal to the temporal distance between the instants of maximum trunk angular velocity and ball impact; (ii) instants of maximum wrist linear velocity and maximum trunk angular velocity that occur at about 24% and 38% of the temporal distance duration between the ball impact and the maximum racquet rotation angle, with these proportions not being exhibited by high-level tennis players. All of these scenarios seem to suggest that harmonic properties in top-level athletes—possessing talent, natural aptitudes or skills— are an expression of a harmonic rhythm generated by CPG, in which the repetitive proportions of phases act as an attractor for a smooth, efficient, graceful and melodic flow of movements in which different limbs, joints and muscles are controlled at the same time and with a high level of coordination.

8. Conclusions

A temporal partition for the butterfly stroke that is harmonically characterized by Fibonacci sequence-based self-similarity—in both its simple and enhanced versions—was presented and discussed. High frame rate video-based experiments on nine elite swimmers (seven of international level and two of national level) were included. The practical implications, strengths and limits of the analysis were discussed.

Author Contributions

Conceptualization: C.M.V., I.F., V.B., G.A.; Methodology: C.M.V., V.B., G.A., R.J.; Software and Resources: C.M.V., V.B., C.R., I.F.; Formal Analysis: C.M.V.; Validation, Investigation: C.M.V., I.F., V.B., C.R., G.A.; Writing—Original Draft: C.M.V.; Writing—Review and Editing: C.M.V., V.B., R.J., G.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

The Internal Research Board (University of Rome Tor Vergata) examined the documentation about this paper and recognized that: it is in accordance with ethical principles for medical research involving human subjects of Declaration of Helsinki; it just assesses motor skills during standard practice routines, while involving no invasive technique.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

No publicly archived datasets were analyzed or generated during the study.

Acknowledgments

The first author is indebted to Marcello Molinari, Technical Director of “New Country Club Frascati” and “Federazione Italiana Tennis (FIT)”—National Master, for drawing attention to the phase partitions of the forehand tennis stroke.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, L.; Wang, X. The power sums involving fibonacci polynomials and their applications. Symmetry 2019, 11, 635. [Google Scholar] [CrossRef] [Green Version]
  2. Verrelli, C.M.; Iosa, M.; Roselli, P.; Pisani, A.; Giannini, F.; Saggio, G. Generalized finite-length Fibonacci sequences in healthy and pathological human walking: Gait recursivity, asymmetry, consistency, self-similarity, and variability. In Frontiers in Human Neuroscience; Special Issue: Rhythmic Patterns in Neuroscience and Human Physiology; to be published.
  3. Horadam, A.F. A generalized Fibonacci sequence. Am. Math. Mon. 1961, 68, 455–459. [Google Scholar] [CrossRef]
  4. Iosa, M.; Fusco, A.; Marchetti, F.; Morone, G.; Caltagirone, C.; Paolucci, S.; Peppe, A. The golden ratio of gait harmony: Repetitive proportions of repetitive gait phases. Biomed Res. Int. 2013, 2013, 918642. [Google Scholar] [CrossRef] [Green Version]
  5. Marino, R.; Verrelli, C.M.; Gnucci, M. Synchronicity rectangle for temporal gait analysis: Application to Parkinson’s Disease. Biomed. Signal Process. Control 2020, 62, 102156. [Google Scholar] [CrossRef]
  6. Verrelli, C.M.; Romagnoli, C.; Jackson, R.R.; Ferretti, I.; Annino, G.; Bonaiuto, V. Front crawl stroke in swimming: Ratios of phase durations and self-similarity. J. Biomech. 2021, 118, 110267. [Google Scholar] [CrossRef] [PubMed]
  7. Taylor, G.K.; Nudds, R.L.; Thomas, A.L.R. Flying and swimming animals cruise at a Strouhal numer tuned for high power efficiency. Lett. Nat. 2003, 425, 707–711. [Google Scholar] [CrossRef] [PubMed]
  8. Eloy, C. Optimal Strouhal number for swimming animals. J. Fluids Struct. 2012, 30, 215–218. [Google Scholar] [CrossRef] [Green Version]
  9. Barbosa, T.M.; Keskinen, K.L.; Fernandes, R.; Colaço, P.; Lima, A.B.; Vilas-Boas, J.P. Energy cost and intracyclic variation of the velocity of the centre of mass in butterfly stroke. Eur. J. Appl. Physiol. 2005, 93, 519–523. [Google Scholar] [CrossRef]
  10. Seifert, L.; Boulesteix, L.; Chollet, D.; Vilas-Boas, J.P. Differences in spatial-temporal parameters and arm–leg coordination in butterfly stroke as a function of race pace, skill and gender. Hum. Mov. Sci. 2008, 27, 96–111. [Google Scholar] [CrossRef]
  11. Di Prampero, P.E.; Pendergast, D.R.; Wilson, D.W.; Rennie, D.W. Blood Lactic Acid Concentrations in High Velocity Swimming. In Swimming Medicine IV, Proceedings of the Fourth International Congress on Swimming Medicine, Stockholm, Sweden; Eriksson, B., Furberg, B., Eds.; University Park Press: Baltimore, MD, USA, 1978; pp. 249–261. [Google Scholar]
  12. Chollet, D.; Seifert, L.; Boulesteix, L.; Carter, M. Arm–leg coordination in elite butterfly swimmers. Int. J. Sport Med. 2006, 27, 322–329. [Google Scholar] [CrossRef]
  13. Seifert, L.; Delignieres, D.; Chollet, D. Effect of expertise on butterfly stroke coordination. J. Sport. Sci. 2007, 25, 131–141. [Google Scholar] [CrossRef] [PubMed]
  14. Lanotte, N.; Annino, G.; Bifaretti, S.; Gatta, G.; Romagnoli, C.; Salvucci, A.; Bonaiuto, V. A new device for propulsion analysis in swimming. Proceedings 2018, 2, 285. [Google Scholar] [CrossRef] [Green Version]
  15. Tourny-Chollet, C.; Chollet, D.; Hogie, S.; Papparodopoulos, C. Kinematic analysis of butterfly turns of international and national swimmers. J. Sport. Sci. 2002, 20, 383–390. [Google Scholar] [CrossRef]
  16. Mason, B.R.; Tong, Z.; Richards, R.J. A biomechanical analysis of the butterfly stroke. Excel 1991, 7, 6–11. [Google Scholar]
  17. Chollet, D.; Chalies, S.; Chatard, J.C. A new index of coordination for the crawl: Description and usefulness. Int. J. Sport. Med. 1999, 21, 54–59. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Ungerechts, B.E. A comparison of the movements of the rear parts of dolphins and butterfly swimmers. In Biomechanics and Medicine in Swimming; Hollander, A.P., Huijing, P.A., de Groot, G., Eds.; Human Kinetics: Champaign, IL, USA, 1983; pp. 215–221. [Google Scholar]
  19. Sanders, R.H.; Cappaert, J.M.; Devlin, R.K. Wave Characteristics of butterfly swimming. J. Biomech. 1995, 28, 9–16. [Google Scholar] [CrossRef]
  20. Barbosa, T.M.; Goh, W.X.; Morais, J.E.; Costa, M.J. Variation of linear and nonlinear parameters in the swim strokes according to the level of expertise. Mot. Control 2017, 21, 312–326. [Google Scholar] [CrossRef]
  21. Schiffman, J.M.; Chelidze, D.; Adams, A.; Segala, D.B.; Hasselquist, L. Nonlinear analysis of gait kinematics to track changes in oxygen consumption in prolonged load carriage walking: A pilot study. J. Biomech. 2009, 42, 2196–2199. [Google Scholar] [CrossRef] [PubMed]
  22. Barbosa, T.M.; Goh, W.X.; Morais, J.E.; Costa, M.J. Comparison of classical kinematics, entropy, and fractal properties as measures of complexity of the motor system in swimming. Front. Psychol. 2016, 7, 1–7. [Google Scholar] [CrossRef] [Green Version]
  23. Takagi, H.; Sugimoto, S.; Nishijima, N.; Wilson, B. Differences in stroke phases, arm-leg coordination and velocity fluctuation due to event, gender and performance level in breaststroke. Sport. Biomech. 2004, 3, 15–27. [Google Scholar] [CrossRef] [PubMed]
  24. Ualí, I.; Herrero, A.J.; Garatachea, N.; Marín, P.J.; Alvear-Ordenes, I.; García-López, D. Maximal strength on different resistance training rowing exercises predicts start phase performance in elite kayakers. J. Strenght Cond. Res. 2012, 26, 941–946. [Google Scholar] [CrossRef] [PubMed]
  25. Landlinger, J.; Lindinger, S.; Stöggl, T.; Wagner, H.; Müller, E. Key factors and timing patterns in the tennis forehand of different skill levels. J. Sport. Sci. Med. 2010, 9, 643–651. [Google Scholar]
Figure 1. Butterfly stroke: arm and leg phases (case of T 2 , T 3 , T 4 > 0 ) according to [12].
Figure 1. Butterfly stroke: arm and leg phases (case of T 2 , T 3 , T 4 > 0 ) according to [12].
Mathematics 09 01545 g001
Table 1. Experimental data for swimmers IL1–IL7 and NL1–NL2: arm phase durations (in seconds).
Table 1. Experimental data for swimmers IL1–IL7 and NL1–NL2: arm phase durations (in seconds).
t EC t PL t PS t RE t S
IL1 0.220 0.220 0.270 0.440 1.150
IL2 0.233 0.183 0.316 0.333 1.065
IL3 0.308 0.158 0.266 0.333 1.065
IL4 0.375 0.200 0.258 0.407 1.240
IL5 0.291 0.291 0.250 0.400 1.232
IL6 0.300 0.283 0.333 0.367 1.283
IL7 0.234 0.175 0.358 0.417 1.184
NL1 0.266 0.166 0.308 0.326 1.066
NL2 0.220 0.150 0.200 0.330 0.900
Table 2. Experimental data for swimmers IL1–IL7 and NL1–NL2: leg phase durations (in seconds).
Table 2. Experimental data for swimmers IL1–IL7 and NL1–NL2: leg phase durations (in seconds).
t K 1 t K 2 t U 1 t U 2 t S
IL1 0.120 0.160 0.420 0.450 1.150
IL2 0.108 0.166 0.400 0.466 1.065
IL3 0.125 0.150 0.391 0.383 1.065
IL4 0.133 0.191 0.441 0.450 1.240
IL5 0.120 0.141 0.483 0.430 1.232
IL6 0.150 0.166 0.517 0.417 1.283
IL7 0.150 0.158 0.467 0.475 1.184
NL1 0.108 0.183 0.417 0.408 1.066
NL2 0.120 0.140 0.330 0.330 0.900
Table 3. Delays and related values for swimmers IL1–IL7 and NL1–NL2: delay durations (in seconds) and NTD.
Table 3. Delays and related values for swimmers IL1–IL7 and NL1–NL2: delay durations (in seconds) and NTD.
T 1 T 2 T 3 T 4 NTD
IL1    0.000 0.100 0.100 0.010 12.320 %
IL2 0.075 0.200 0.017 0.133 23.680 %
IL3    0.016 0.167 0.066 0.050 17.570 %
IL4    0.025 0.217 0.024 0.043 18.060 %
IL5    0.058 0.113 0.079 0.030 12.380 %
IL6    0.033 0.117 0.117 0.050 13.720 %
IL7 0.066 0.150 0.142 0.058 18.960 %
NL1 0.050 0.208 0.043 0.082 21.860 %
NL2 0.020 0.120 0.060     0.000 15.070 %
Table 4. Delays and related values for swimmers IL1–IL7 and NL1–NL2: ( a , b , c , d ) -values.
Table 4. Delays and related values for swimmers IL1–IL7 and NL1–NL2: ( a , b , c , d ) -values.
abcd
IL1    0.00 % 8.69 % 8.69 % 0.87 %
IL2 7.04 % 18.78 % 1.60 % 12.49 %
IL3    1.50 % 15.68 % 6.20 % 4.69 %
IL4    2.02 % 17.50 % 1.94 % 3.47 %
IL5    4.71 % 9.17 % 6.41 % 2.44 %
IL6    2.57 % 9.12 % 9.12 % 3.90 %
IL7 5.57 % 12.67 % 11.99 % 4.90 %
NL1 4.69 % 19.51 % 4.03 % 7.69 %
NL2 2.22 % 13.33 % 6.67 %     0.00 %
Table 5. Data analysis for swimmers IL1–IL7 and NL1–NL2: aggregate phases (percentage values).
Table 5. Data analysis for swimmers IL1–IL7 and NL1–NL2: aggregate phases (percentage values).
( t K 1 + t K 2 + t U 1 ) / t S t U 1 / t S ( t K 1 + t K 2 ) / t S t K 2 / t S t K 1 / t S
IL1 60.87 36.52 24.35 13.91 10.43
IL2 63.29 37.56 25.73 15.59 10.14
IL3 62.54 36.71 25.82 14.08 11.74
IL4 61.69 35.56 26.13 15.40 10.73
IL5 60.39 39.20 21.19 11.44 9.74
IL6 64.93 40.30 24.63 12.94 11.69
IL7 65.46 39.44 26.01 13.34 12.67
NL1 66.42 39.12 27.30 17.17 10.13
NL2 65.56 36.67 28.89 15.56 13.33
Table 6. Data analysis for swimmers IL1–IL7 and NL1–NL2: values of indices (23) and (24).
Table 6. Data analysis for swimmers IL1–IL7 and NL1–NL2: values of indices (23) and (24).
I f , 4 I f , 6
IL1 1.89 1.94
IL2 2.20 2.72
IL3 2.30 2.88
IL4 3.25 3.61
IL5 3.45 4.31
IL6 3.78 4.27
IL7 4.25 5.06
NL1 5.63 6.46
NL2 6.19 7.20
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Verrelli, C.M.; Romagnoli, C.; Jackson, R.; Ferretti, I.; Annino, G.; Bonaiuto, V. Phi-Bonacci Butterfly Stroke Numbers to Assess Self-Similarity in Elite Swimmers. Mathematics 2021, 9, 1545. https://doi.org/10.3390/math9131545

AMA Style

Verrelli CM, Romagnoli C, Jackson R, Ferretti I, Annino G, Bonaiuto V. Phi-Bonacci Butterfly Stroke Numbers to Assess Self-Similarity in Elite Swimmers. Mathematics. 2021; 9(13):1545. https://doi.org/10.3390/math9131545

Chicago/Turabian Style

Verrelli, Cristiano Maria, Cristian Romagnoli, Roxanne Jackson, Ivo Ferretti, Giuseppe Annino, and Vincenzo Bonaiuto. 2021. "Phi-Bonacci Butterfly Stroke Numbers to Assess Self-Similarity in Elite Swimmers" Mathematics 9, no. 13: 1545. https://doi.org/10.3390/math9131545

APA Style

Verrelli, C. M., Romagnoli, C., Jackson, R., Ferretti, I., Annino, G., & Bonaiuto, V. (2021). Phi-Bonacci Butterfly Stroke Numbers to Assess Self-Similarity in Elite Swimmers. Mathematics, 9(13), 1545. https://doi.org/10.3390/math9131545

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

Article Metrics

Back to TopTop