The control of flight force by a flapping wing: lift and drag production
Department of Integrative Biology, University of California, Berkeley, CA 94720, USA
*e-mail: sane{at}socrates.berkeley.edu
Accepted May 18, 2001
![]() |
Summary |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
(i) For a short, symmetrical wing flip, mean lift was optimized by a stroke amplitude of 180° and an angle of attack of 50°. At all stroke amplitudes, mean drag increased monotonically with increasing angle of attack. Translational quasi-steady predictions better matched the measured values at high stroke amplitude than at low stroke amplitude. This discrepancy was due to the increasing importance of rotational mechanisms in kinematic patterns with low stroke amplitude.
(ii) For a 180° stroke amplitude and a 45° angle of attack, lift was maximized by short-duration flips occurring just slightly in advance of stroke reversal. Symmetrical rotations produced similarly high performance. Wing rotation that occurred after stroke reversal, however, produced very low mean lift.
(iii) The production of aerodynamic forces was sensitive to changes in the magnitude of the wings deviation from the mean stroke plane (stroke deviation) as well as to the actual shape of the wing tip trajectory. However, in all examples, stroke deviation lowered aerodynamic performance relative to the no deviation case. This attenuation was due, in part, to a trade-off between lift and a radially directed component of total aerodynamic force. Thus, while we found no evidence that stroke deviation can augment lift, it nevertheless may be used to modulate forces on the two wings. Thus, insects might use such changes in wing kinematics during steering maneuvers to generate appropriate force moments.
(iv) While quasi-steady estimates failed to capture the time course of measured lift for nearly all kinematic patterns, they did predict with reasonable accuracy stroke-averaged values for the mean lift coefficient. However, quasi-steady estimates grossly underestimated the magnitude of the mean drag coefficient under all conditions. This discrepancy was due to the contribution of rotational effects that steady-state estimates do not capture. This result suggests that many prior estimates of mechanical power based on wing kinematics may have been grossly underestimated.
Key words: flapping flight, quasi-steady force, unsteady aerodynamics, fruit fly, Drosophila melanogaster, added mass, delayed stall, rotational circulation, wake capture.
![]() |
Introduction |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
In a few cases, researchers have attempted to capture the free-flight kinematics of maneuvering insects (Ennos, 1989; Ruppell, 1989). While such analyses are essential because they reveal what insects actually do with their wings when steering, free-flight studies are limited because it is not yet feasible to relate the changes in wing kinematics directly to changes in instantaneous aerodynamic forces. An alternative approach is to measure instantaneous forces on tethered insects (Cloupeau et al., 1979; Wilkin, 1990; Zanker, 1990b; Zanker and Götz, 1990; Dickinson and Götz, 1996). However, forces and stroke kinematics measured on tethered insects may not accurately represent those generated in free flight. Further, since tethered flight force transducers measure whole-body forces, it is not possible to resolve the instantaneous aerodynamic forces generated by individual wings. A third approach is to calculate the aerodynamic forces generated by arbitrary stroke kinematics using computational fluid dynamics (Liu et al., 1998; Wang, 2000). However, because of the critical role of unsteady mechanisms and three-dimensional flow structure in insect flight aerodynamics (Ellington et al., 1996; Dickinson et al., 1999), theoretical or numerical approaches have, as yet, offered only limited insight into the aerodynamics of steering.
Given the current limitations in studies of both real animals and numerical simulations, we have chosen to study the problem of maneuverability using a dynamically scaled model of a flapping insect. Aerodynamic models have proved valuable in the study of insect flight, particularly in the identification and analysis of unsteady aerodynamics (Bennett, 1977; Maxworthy, 1979; Spedding and Maxworthy, 1986; Dickinson and Götz, 1993; Ellington et al., 1996; Dickinson et al., 1999). In large part through the use of mechanical models, researchers have identified an array of mechanisms that collectively account for the elevated aerodynamic performance of flapping wings. These include the clap and fling (Spedding and Maxworthy, 1986), dynamic stall (Dickinson and Götz, 1993; Ellington et al., 1996), rotational lift (Bennett, 1970; Dickinson et al., 1999) and wake capture (Dickinson, 1994; Dickinson et al., 1999). Now that the various mechanisms responsible for the elevated aerodynamic performance of insect wings have been identified, it is possible to tackle the question of how animals manipulate such mechanisms to steer and maneuver.
In this study, we use a dynamically scaled mechanical model of Drosophila melanogaster to investigate how changes in wing kinematics affect the production of aerodynamic forces. In particular, we explore the influence of five behaviorally relevant kinematic parameters: stroke amplitude, angle of attack, the timing and duration of wing rotation and stroke plane deviation. We chose this particular set of parameters because fruit flies actively vary them during flight maneuvers (Götz et al., 1979; Zanker, 1990a; Dickinson et al., 1993; Lehmann and Dickinson, 1998). However, the goal of this project is not to replicate the precise kinematics of free flying insects per se, but rather to map aerodynamic forces within a broad parameter space that encompasses the variation seen among insects.
From the instantaneous force records, we calculate time-averaged aerodynamic force coefficients, lift-to-drag ratios and other measures of aerodynamic performance. The resultant data set is useful in identifying the kinematic parameters that most influence the magnitude and direction of aerodynamic forces generated by flapping wings. In a companion paper (S. P. Sane and M. H. Dickinson, in preparation), we will extend the analysis by considering the instantaneous and time-averaged force moments generated about the yaw, pitch and roll axes. The comprehensive parameter maps generated in these studies should be of help to biologists who wish to know the aerodynamic consequences of observed changes in wing kinematics as well as to engineers who wish to optimize the performance of small biomimetic flying robots. In addition, these data provide experimental validations for numerical simulations of the fluid motion around flapping wings.
![]() |
Materials and methods |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
Stroke kinematics
In the absence of wing deformation, the kinematics of the wings may be uniquely described by specifying the time course of three angles: stroke position, (t), angle of attack,
(t), and stroke deviation,
(t) (Fig.1B). In all experiments, the angular position of the wing within the stroke plane was described by a triangular waveform, which maintains a constant translational velocity throughout each half-stroke. The waveform was smoothed to minimize inertial accelerations during stroke reversal and to match more closely published stroke kinematics from a variety of insects (Ellington, 1984b; Zanker, 1990a). For smoothing, we filtered the triangular waveform using a zero-phase-delay low-pass two-pole Butterworth filter with a cut-off frequency equal to 10 times the stroke frequency of 0.17Hz. The peak-to-peak amplitude of the stroke angle waveform could be varied in each experiment. The angle of attack was described by a trapezoidal wave function, which maintained a constant angle of attack during each half-stroke and constant rotational velocity during stroke reversal. The shape of this waveform in each experiment was determined by setting the mid-stroke angles of attack during the upstroke and downstroke and by specifying the starting and stopping points for wing rotation. The resulting function was then smoothed using a low-pass filter with identical characteristics to that used for the stroke position waveform. We used two functions to describe stroke deviation: an oval pattern in which the wing tip deviated from the stroke plane according to a half-sine-wave per stroke period and a figure-of-eight pattern in which the stroke deviation varied as a full sine-wave. These patterns were chosen because they roughly approximate patterns described for a variety of insects (Ellington, 1984b; Zanker, 1990a).
To create the kinematic patterns used in this study, we varied any or all of six parameters: (i) the stroke amplitude, (ii) the mid-stroke angle of attack during upstroke and downstroke, (iii) the timing of wing rotation at dorsal and ventral reversal, (iv) the duration of the stroke reversal, (v) the shape of the wing tip trajectory (oval or figure-of-eight) and (vi) the angular deviation from the mean stroke plane during the upstroke and downstroke (Fig.1C). In most of the experiments, the deviation amplitude was set to zero, such that the wing tip remained within the stroke plane throughout the cycle. Under these conditions, the kinematics of the wing stroke were symmetrical such that the upstroke and downstroke were mirror images of one another. Only in trials using oval stroke deviations were the kinematics of the two strokes not identical. The frequency of the wing stroke (0.17Hz) remained constant in all experiments, as did the upstroke-to-downstroke duration ratio, which was fixed at 1. We constructed the kinematic patterns using a custom-designed MATLAB program (Mathworks) to convert the angular trajectories into a series of stepper motor commands.
Force measurements
Signals from the two-dimensional sensor were acquired using a National Instruments data-acquisition board (model BNC 2090) in a PC running custom-designed software written in MATLAB. Data were filtered on-line with an active four-pole Bessel filter with a cut-off frequency of 10Hz and off-line with a zero-phase-delay low-pass digital Butterworth filter with a cut-off frequency of 3Hz, which was 17.6 times the wing stroke frequency. Apart from increasing the high-frequency components resulting from motor jitter, increasing the cut-off frequency of the filter did not alter the time course of the force traces.
Each experiment consisted of one burst of four consecutive wing strokes following pre-programmed kinematics. The wing begins the first downstroke in still fluid, whereas during the subsequent strokes it moves through a wake created by the preceding strokes. As a result, the time course of forces generated during the first stroke is markedly different from those of subsequent strokes. For this reason, the data from the first stroke were excluded from this analysis, while those from the three subsequent strokes were averaged. Thus, each presented trace represents an average of three force records. After subtracting gravitational forces, the forces measured from the normal and parallel channels were transformed into lift, drag, thrust and radial components.
Added mass
The measured force at the wing base consists of gravitational, inertial and aerodynamic components. The gravitational contribution of the sensor and wing mass to the total force signal was easily calculated and subtracted from the measured force traces. The inertial components represent the acceleration forces on the mass of the sensor and wing as well as the added mass of the fluid around the wing. To examine the contribution of the inertial effects of the wing mass and sensor, we replaced the wing surface with an aerodynamically neutral inertial model of the wing. The aerodynamically neutral model was essentially a brass knob with the same mass and center of mass inside the oil as the Plexiglas wing. Because of its low surface area, the brass knob generated negligible aerodynamic forces compared with the Plexiglas wing. For any arbitrary kinematic pattern, the resulting force traces for the brass knob could be entirely accounted for by gravity. Thus, the inertial forces generated by flapping this brass model, and therefore the Plexiglas wing, were negligible and have been ignored. Compared with gravity and wing inertia, the non-circulatory forces due to added mass are more difficult to measure because the fluid acceleration induced by a moving wing changes dynamically as the wing rotates, decelerates or accelerates (Daniel, 1984).
To estimate the magnitude of added mass, we used an approximation derived for motions of an infinitesimally thin two-dimensional plate in an inviscid fluid (Sedov, 1965). Using blade element method, we adapted it to the case of a three-dimensional wing rotating around an axis located at one-quarter chord length from the leading edge. The force contribution normal to the wing surface due to the added mass inertia is given by:
|
where is the fluid density, R is the wing length,
is the mean chord length,
and
(
) are the non-dimensional
radial position along the wing and non-dimensional chord length,
respectively (for nomenclature, see Ellington, 1984a),
is the
angular position of the wing and
is the angle of attack. Using
equation 1, we calculated an estimate of added mass inertia for each
set of kinematics. As illustrated by the representative trace in
Fig.2, the absolute
contribution of added mass to net forces on the wing was quite small
in all cases. Further, by comparing equation 1 and equation 3 (see
below), it can be seen that the added mass forces
(
R2
2) scale in proportion to aerodynamic forces
(
R3
) for geometrically similar wings. Thus, for identical kinematics and
geometry, added mass will have the same physical effect on a model
wing as on the wing of a fly, provided that the Reynolds number is the
same. Because both added mass and aerodynamic contributions are
biologically relevant, we chose not to subtract the inertial estimates
from the presented force traces.
|
![]() |
where 0 is flip start and
is flip duration. As before, we determined the values of flip timing and flip duration that maximized mean lift within the 11x9 sets of trials. In the third set of experiments, we explored how aerodynamic forces varied with changes in deviation from the mean stroke plane using the values of stroke amplitude, angle of attack, flip timing and flip duration that maximize mean lift. In a set of 22 experiments, the total stroke deviation was varied from -50° (-25° during the downstroke, +25° during the upstroke) to +50° (+25° during the downstroke, -25° during the upstroke) in steps of 10°. The deviation followed either a half-sine per stroke or a full-sine per stroke time course, yielding oval and figure-of-eight trajectories, respectively (Fig.1B).
Data analysis
Since the conventions for lift and drag existing in current aerodynamic literature address non-flapping and primarily two-dimensional kinematic patterns, it is necessary to modify them slightly for three-dimensional motions. These modifications are based on the following two considerations. First, it is convenient to use a reference frame based on the insect body rather than its wing so that the measurements relate directly to free-flight studies. Second, the standard convention should apply when the kinematics reduce to two-dimensional motion. With these constraints in mind, we adopted the following convention: net aerodynamic force, defined as the total force on the wing, is resolved into three components: lift, drag and radial force. In hovering flight, lift must offset the gravitational force on the animals body mass. Hence, we define lift as the component of the net aerodynamic force perpendicular to the mean stroke plane of the wing regardless of its actual instantaneous trajectory. Since the mean stroke plane was horizontal in all experiments, lift is always the vertical force component. Drag is defined as the force component in the horizontal direction, opposing the wing movement. The radial component accounts for the remaining force component in the horizontal plane. For motions with no stroke deviation, these definitions reduce to the standard convention: lift is orthogonal to drag, and the radial component vanishes. With stroke deviation, the total normal pressure force consists of orthogonal vertical and radial components, each orthogonal to the drag vector in the horizontal direction.
From the forces on each wing, we calculated the corresponding mean force coefficients using an equation derived from blade element theory (Ellington, 1984c; Dickinson et al., 1999):
|
where is the magnitude of a specific force
component (lift, drag, radial, total) averaged over the stroke,
is stroke amplitude, n is wing beat frequency,
is the
mean non-dimensional angular velocity of the wing and
22(S)
is the non-dimensional second moment of wing area. The radial force
component changes sign when the wing crosses the mean stroke plane. As
a result, in the oval as well as the
figure-of-eight patterns, the mean radial coefficient
often averages to zero and is uninformative. For this reason, we base
our measurement of the average radial force coefficient on the
absolute values. Lift-to-drag ratios were calculated by dividing the
mean lift coefficient by the mean drag coefficient. Similarly,
radial-to-drag ratios were obtained by dividing the mean absolute
radial coefficients by the mean drag coefficients. To calculate the
ratio of mean lift to profile power, we estimated mean profile power,
, based on the
time-averaged product of instantaneous drag, D(t),
and instantaneous velocity, vwing(t):
|
where T is the stroke period. This calculation of mean
profile power ignores the power required to rotate the wing in
place. Mean lift was calculated as the time
average of instantaneous lift throughout the stroke.
Measures of the quasi-steady-state translational force
coefficients and
were derived from 180° sweeps of wing motion with fixed angles of
attack, as described elsewhere (Dickinson et al., 1999). The
equations that best fit measured translational force coefficients as
functions of angle of attack,
, for the model wing are
(Dickinson et al.,
1999):
![]() |
and
![]() |
These equations are used to generate quasi-steady translational estimates for comparison with measured values.
![]() |
Results |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
|
The influence of rotational circulation is most easily seen when the angle of attack during the upstroke and downstroke is zero (Fig.3G,H). Under these conditions, the vorticity generated during the translational phase of the stroke is minimal and, thus, the magnitude of the wake capture effect should be small. The influence of the wake is not entirely absent, however, because the process of rotation generates and sheds vorticity through which the wing must translate at the start of each stroke. Further, the vorticity created by rotation is particularly strong when the translational angle of attack is zero, because the wing must flip over by 180° during stroke reversal, making the angular velocity of the wing particularly large. At the end of the each stroke in Fig.3G,H, the influence of rotational circulation is manifest as a transient increase in lift and drag that exceeds the quasi-steady prediction. After stroke reversal, the continuing rotation of the wing generates a pressure force with opposite polarity, resulting in negative lift. The time course of this rotational effect is complicated by the presence of an added mass inertia and a modest amount of wake capture at the start of each stroke. The influence of these multiple mechanisms is manifest by the positive peak in lift immediately following stroke reversal due to wake capture, which is followed by the negative peak due to rotational circulation. The pattern of an early positive peak in lift followed by a later negative peak is seen throughout the traces in Fig.3.
The rest of the traces in Fig.3 illustrate the complex interactions among delayed stall, rotational circulation and wake capture that result from changes in angle of attack and stroke amplitude. At angles of attack of 30 and 50°, the wing generates lift throughout the stroke due to delayed stall (Fig.3CF). The influence of rotational lift is reduced as the angle of attack increases, however, because the wing flips over a smaller arc with lower angular velocity. This effect can be seen by comparing the relative magnitude of the force peaks at the end of each stroke in Fig.3B,D,F,H. In contrast, the influence of wake capture is greater at higher angles of attack because the vorticity shed into the wake at the end of the stroke is stronger. This effect can be seen by comparing the relative size of the force transients at the start of each stroke in the same panels. Thus, changing stroke amplitude and angle of attack has a complex but interpretable influence on the magnitude of the different unsteady mechanisms. The kinematics that optimize the aerodynamic performance of the wing will reflect these complex interactions. The maximum mean lift-to-drag ratio (0.8) occurred at an angle of attack of 30° and amplitude of 180°. The forces corresponding to this optimal condition are shown in Fig.3F.
To provide a more comprehensive picture of how force production
changes with kinematics, we calculated the mean lift, drag and net
force coefficients averaged throughout the stroke and plotted them for
all pairs of stroke amplitude and angle of attack values. Fig.4A,B depicts the mean total
aerodynamic force,
,
and force coefficient,
,
in pseudocolor maps. These two maps differ from one another because
the mean force coefficient is normalized with respect to the square of
stroke amplitude, which is a variable in these experiments. Thus,
while
increases with increasing stroke amplitude (Fig.4A),
decreases with increasing stroke amplitude (Fig.4B). Both parameters rise with increasing
angle of attack. The influence of stroke amplitude and angle of attack
on the mean lift coefficient
and the mean drag
coefficient
is shown in Fig.5A,C. For a fixed stroke amplitude,
exhibits a broad maximum ranging from 1.8 to 2.0 between angles of
attack of 40 and 50°. As expected,
rises monotonically with increasing angle of attack for any given
value of stroke amplitude. It is worth noting that the range of
values is much higher than has been previously reported for
Drosophila virilis wings under steady-state conditions
(Vogel, 1967) or
estimated on the basis of Reynolds number (CD
0.7;
Ellington,
1984c). Although much less pronounced than the dependence
on angle of attack, lift tends to rise, whereas drag falls, with
increasing stroke amplitude.
|
|
The effects of flip timing and flip duration
In all subsequent experiments, we set the stroke amplitude to
180° and the mid-stroke angle of attack to 45°, which are the
values that maximized lift production in the first set of
experiments. Next, we systematically varied the timing and duration of
wing rotation to examine their effects on force production. Sample
force traces selected from 99 pairs of flip start and flip duration
are shown in Fig.6AH, with the corresponding values
of flip timing, f, given in the upper left corner of
each panel. A comparison of Fig.6A,C,D illustrates that a long-duration
flip,
=0.5, can produce quite different forces
depending on when the flip occurs. If the flip occurs symmetrically
about the stroke reversal (
f=0, Fig.6A),
is quite large
(
=1.54)
and the time course is well approximated by the quasi-steady
predictions. An advanced flip (
f=-0.25,
Fig.6C) results in very
low mean lift
(
=0.36),
but produces a fairly prominent wake capture peak. In contrast, a
rotation after stroke reversal (
f=+0.25,
Fig.6D) results in mean
negative lift
(
=-0.28)
because of the adverse effects of rotational circulation following
stroke reversal.
|
The maps of mean force coefficients as a function of flip duration and flip start are shown in Fig.7A,C,E, with comparable quasi-steady, translational predictions shown in Fig.7B,D,F. Flip timing, the non-dimensional time when the mid-point of the flip occurs, is indicated by the inclined parallel lines on each graph. Both
and
are strongly influenced by flip timing and duration. For example, at a flip duration of 0.1,
varies with flip timing from as low as -1.5 to as high as +2. The comparable values of the quasi-steady translational estimate,
, also vary, but over a smaller range (from -1 to 1.6).
|
The effects of stroke plane deviation
Using kinematic values for stroke amplitude, angle of attack, flip duration and flip start that maximized lift production (=180°,
=45°,
=0.1 and
f=-0.05), we tested how forces vary with deviation from the mean stroke plane in a set of 22 experiments. The peak-to-peak magnitude of stroke deviation was varied from -50 to +50° in 10° steps for both the half-sine (oval) and full-sine (figure-of-eight) patterns. It is worth noting that, in the oval pattern, an upward deviation at the start of the downstroke requires a downward deviation at the start of the upstroke and vice versa. This is not the case for the figure-of-eight pattern, in which the two half-strokes are mirror images of one another.
Fig.8 shows a selection of force traces resulting from different patterns of stroke deviation. In general, the figure-of-eight pattern had a more profound influence on the magnitude and time course of force production than did the oval pattern. In both cases, however, the direction of stroke deviation at the start of each translational phase greatly influenced the magnitude of the force transient at the start of the stroke. For example, in the figure-of-eight pattern shown in Fig.8C, each stroke begins with an upward motion, and the lift and drag transients at the start of each stroke are quite small. In contrast, the comparable kinematic pattern that starts with a downward motion (Fig.8D) generates sizeable force peaks. A similar trend is seen in the oval patterns (Fig.8A,B). The upstroke in Fig.8A and the downstroke in Fig.8B, which both start with a downward motion, are marked by sizeable force peaks at the start of translation, whereas the strokes that begin with upward motion are not. This dependence of the early force transient on the direction of deviation is explained in part by an increase in the aerodynamic angle of attack caused by the downward motion of the wing. However, the measured force peaks are much greater than the quasi-steady estimates, which take into account this effect, suggesting that there is a substantial wake effect at the start of each stroke. The influence of the wake is stronger if the wing moves downwards, towards the descending vorticity of the previous stroke, than if it moves upwards, away from the wake. In all cases, the significance of radial forces increases with increasing deviation. As expected, like the lift and drag forces, the time course of radial forces is also dependent on the shape of the wing trajectory.
|
Fig.9 summarizes the effects of stroke deviation on the time-averaged measured forces and quasi-steady predictions. While the influence of stroke deviation on the time course of the aerodynamic forces is large, its impact on the mean force coefficients is surprisingly small. This indicates that the differences in the dynamics of force production noted in Fig.8 tend to average out over the stroke. For both the oval and figure-of-eight deviation trajectories (Fig.9A), the mean lift and drag coefficients decreased with increasing absolute deviation (Fig.9C,D). The changes in average performance for the oval deviation pattern should be symmetrical around zero deviation, since oval patterns with positive and negative deviations are mirror images of one another. Thus, the downstroke in Fig.8A should resemble the upstroke in Fig.8B, and the upstroke in Fig.8A should resemble the downstroke in Fig.8B. The asymmetry in these measurements results from the mechanical play in the gear mechanism of the robot. However, the asymmetry in the performance of the figure-of-eight patterns around zero deviation represents, in part, a real aerodynamic effect (Fig.9B,C,D). In this case, a positive deviation will result in a downward motion at the start of both the upstroke and downstroke, whereas a negative deviation indicates upward motion at the beginning of both strokes. Downward deviation should enhance wake capture, as described above. Values of
and
fall off faster with increasing positive deviation close to zero deviation. However, at large deviations, the coefficients for negative deviations are lower than the coefficients for positive deviations (Fig.9C,D). There is little effect on the
/
ratio because the influence of stroke deviation is nearly identical for both lift and drag. In contrast, because of the linear nature of sine functions at low angles,
/
ratios appear linear with the small range of stroke deviation in our experiments (Fig.9E).
|
Ratio of mean lift to mean profile power
Within the range of Reynolds numbers relevant for fruit flies, the
total mechanical power required to flap the wings is dominated by
profile power, , the cost to
overcome drag on the flapping wings (Lehmann and Dickinson,
1997). Fig. 10AC shows how
, estimated using equation 6,
varies with changes in the kinematics parameters. Fig.10DF shows the
corresponding ratios of mean lift to mean profile power for the same
kinematics. The
estimates vary
extensively, even within subregions of the parameter maps in which the
values of
are high enough to support flight. This result suggests that it may
be difficult to estimate mechanical power in free or tethered flight
solely on the basis of measures of stroke amplitude. In particular,
varies extensively with the timing
of wing rotation and the angle of attack, parameters that are revealed
only by extensive three-dimensional reconstructions of wing
kinematics.
|
![]() |
Discussion |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Comparison between the mechanical model and real flies
The maximum mean unsteady lift coefficients in this study are in
excellent agreement with measurements on tethered Drosophila
spp., which generated elevated flight force in response to optomotor
stimuli. Peak
values in tethered flight were 1.9 (Lehmann and Dickinson,
1998), which is precisely the same maximum value measured
on the model wing. Tethered flies generating just enough force to
support body weight produce a
value of 1.6. The tethered flight estimates relied on the assumption
that both up- and downstrokes generated lift, which appears to be
correct given the time course of forces generated when
Drosophila spp. kinematics are played through the model
(Dickinson et al.,
1999). Tethered Drosophila spp. flap with
morphologically maximum stroke amplitude when producing peak lift
(Lehmann and
Dickinson, 1998), which is consistent with the present
results. Unfortunately, apart from stroke amplitude, we do not have
adequate knowledge of other kinematic parameters during the peak
performance of real flies to compare them with the values that
maximized lift on the model.
As seen in Fig.5B, the maximum quasi-steady translational estimate of the mean lift coefficient,
was 1.6. Given that this value is sufficient to explain the forces required for a fruit fly to hover (Lehmann and Dickinson, 1998), it is tempting to claim that the quasi-steady estimates are sufficient to explain insect flight (Jensen, 1956). However, a simple comparison of time-averaged lift coefficients is not a robust test of the quasi-steady model. While the mean values might be similar, the time histories of measured forces and quasi-steady translational estimates differ greatly (Fig.3, Fig.6, Fig.8). Further, the precise time history of lift and drag is critical for calculating force moments and, thus, essential to considerations of stability and maneuverability. In addition, since many insects can fly while supporting forces nearly twice their body weight (Marden, 1987), any model must explain not only the forces required to hover, but also those required for maximal lift.
The inadequacies of the translational quasi-steady model are even more apparent when considering drag. In all experiments, we observed that the measured drag coefficients were greater in magnitude than those estimated from measured translational force coefficients. Also, the range of drag coefficients measured here (CD=06.5) is substantially higher than previously reported measurements from both real and model Drosophila virilis wings (CD=0.21; Vogel, 1967). The previous experimental values were based on steady-state measurements and thus excluded the contribution of rotational circulation and wake capture as well as added mass. However, even the range of steady-state coefficients in the previous study (0<CL<1, 0<CD< 1) (Vogel, 1967) is substantially lower than the quasi-steady estimates in this study (0<CL,t<1.91; 0.37<CD,t<3.47). At present, the reason for this discrepancy is not clear.
The relative importance of unsteady mechanisms
The relative contributions of delayed stall, rotational circulation and wake capture to total force production vary with the precise kinematics of the stroke. In general, the importance of delayed stall increases with stroke amplitude because the wing can integrate the influence of the leading edge vortex over a greater distance. In contrast, rotational effects become more important with lower stroke amplitude. For kinematic patterns mimicking those of hoverflies with a 60° stroke amplitude, rotational effects account for more than half the total force (Dickinson et al., 1999). Because of these differences, quasi-steady translational estimates should more closely resemble the measured values for kinematic patterns with large stroke amplitudes and deviate from measured values for those with lower stroke amplitudes. These predictions are borne out by the force traces shown in Fig.3 and the maps of
in Fig.5A,B. The quasi-steady translational model predicts maximum lift at an angle of attack of 45°. While this is close to the maximum value for the measured forces at high stroke amplitude, the measured lift maximum is shifted to higher angles of attack (60°) at low stroke amplitude.
There are two possible explanations for this shift. First, rotational circulation might make a greater contribution to mean lift at low stroke amplitude. However, this is unlikely to be true for symmetrical flips, for which rotational circulation enhances lift only before stroke reversal, but attenuates it after stroke reversal (Dickinson et al., 1999). Also, at stroke reversal, the wing rotates less for a 60° angle of attack than it does for a 45° angle of attack. Thus, the proportional contribution due to rotational lift is further minimized. Second, the shift in the lift maximum may reflect an increasing importance of wake capture at low stroke amplitude. The large contribution of wake capture can be easily seen in the forces generated by hoverfly-like patterns (Dickinson et al., 1999) as well as in the drag traces in Fig.3. The increasing importance of wake capture also explains the changes in drag coefficients at low stroke amplitude. While the translational quasi-steady model predicts that drag should be independent of amplitude (Fig.5D), the measured mean drag coefficient clearly increases at smaller stroke amplitudes.
Separating the effects of wing rotation from wake capture
Since rotational circulation, wake capture and added mass usually occur together during stroke reversal, it is often difficult to separate these effects. After estimating the effect of added mass on force traces before and after subtraction (Fig.2), we concluded that, at these Reynolds numbers, the magnitude of added mass is small compared with rotational circulation or wake capture. To isolate these two mechanisms, it is helpful to focus on kinematic patterns in which either the entire flip occurs prior to stroke reversal (Fig.6B,E,F) or the wing does not flip at all (Fig.3A,B). In the case of an advanced flip, the force peak that exceeds the translational estimate prior to stroke reversal may be attributed to rotational circulation, while the force peak after stroke reversal is due to wake capture. In the case of no flip, the large drag peak at the start of each stroke is due to wake capture, in which vorticity shed from the previous stroke elevates force by inducing an increase in flow velocity towards the wing (Dickinson, 1994; Dickinson et al., 1999). Because of the squared dependence of forces on relative velocity, even small changes in flow can cause a large elevation in force.
The influence of wake capture should be reduced when the wing translates at a 0° angle of attack during the previous stroke. Under these circumstances, the wing sheds less vorticity at the end of each stroke, since it did not create a leading edge vortex during translation. There is some vorticity present, however, due to rotational circulation at the end of the stroke. Nevertheless, in these cases (Fig.3G,H), wake capture is reduced, which partially isolates the influence of rotational circulation. The results indicate that a wing rotating early enhances lift, but that if the rotation continues after stroke reversal the force dips below the quasi-steady-state prediction because, after changing direction, rotational circulation will act to counter the circulation produced by translation.
By systematically varying the kinematic parameters, we determined that lift was optimized at a 45° angle of attack, 180° stroke amplitude, 10% flip duration and approximately -5% flip advance. Even subtle deviations from the optimal combination drastically decreased lift. A 5% delay in flip timing relative to the optimal conditions (Fig.6H) is less effective in generating rotational lift and decreases the average coefficient by as much as 25%. Small decreases in flip duration cause comparable changes. These results suggest that the control of flip timing and duration provides an easy and powerful way of modulating the forces on two wings, as is required to alter force moments during steering. Tethered flies advance wing rotation on the outside of a turn in response to a turning visual stimulus (Dickinson et al., 1993). Since early rotation enhances whereas delayed rotation attenuates lift, early wing rotation on the outside of a turn should produce a greater force on that side, thus creating an appropriate turning moment in the opposite direction. While flip advance is usually coupled with an increase in stroke amplitude, flies do possess the ability to decouple the modulation of rotation and amplitude (Dickinson et al., 1993). A more thorough analysis of steering must consider the influence of wing kinematics on force moments, which will be the topic of a separate paper (S. P. Sane and M. H. Dickinson, in preparation)
In nearly all the experiments presented here, the kinematics of the upstroke and downstroke were identical. Under these conditions, the contributions of the upstroke and downstroke to forces are equal, and the net force vector is oriented normal to the mean stroke plane. However, symmetrical wing strokes are probably more the exception in actual insects than the rule. Thus, further experiments must focus on asymmetric strokes that can cause the net force vector to be tilted relative to the stroke plane. Adjusting the symmetry of the two strokes is undoubtedly critical for the regulation of force moments in flight control.
The influence of stroke deviation
The wings of insects do not always beat back and forth within a flat stroke plane, but may exhibit large and complex deviations. Within studies of flies, the true form and function of these deviations have been subject to much debate (for a review, see Dickinson and Tu, 1997). While the existence of oval and figure-of-eight trajectories was often attributed to differences in experimental procedure (tethered versus free flight, thoracic versus abdominal tether, still air versus moving air), studies of neuromuscular control indicate that flies have the ability actively to modulate stroke deviation by altering the activity of steering muscles (Tu and Dickinson, 1994; Lehmann and Götz, 1996). The aerodynamic utility of these changes, however, is still not known.
We can divide the effects of stroke deviation on aerodynamic mechanisms into three categories: the introduction of radial force components, alteration of aerodynamic angle of attack and the efficacy of wake capture. The presence of one or more of these effects can substantially alter the dynamics of force generation even when the time-averaged net force produced by the wing remains largely unaltered (filled circles, Fig.9B). For three-dimensional flapping, the main influence of stroke deviation on flight forces is the presence of radial force components acting orthogonal to lift and drag (Fig.8AE). This force component, although ignored in much of the insect flight literature, is of considerable importance during large stroke deviations. Although the magnitude of the average net force vector remains nearly constant with changes in deviation, both lift and drag coefficients decrease with increasing positive and negative stroke deviations, while the radial component increases. In this way, deviation offers a versatile means of manipulating forces through simple geometric alterations in stroke trajectory.
Stroke deviation also influences force production by altering the angle of attack with respect to the wing path as well as its instantaneous velocity. Upward motion of the wing at the start of the stroke results in a decrease in the angle of attack, while downward motion causes an increase in the angle of attack. The converse is true at the end of the stroke. In addition, the velocity of the wing is higher, because the wing must cover a greater distance when it follows a curved trajectory for a given stroke frequency. This effect is clear in the quasi-steady prediction for the oval wingbeat patterns (Fig.8A,B). The predicted forces are higher at the start of strokes beginning with downward motion and at the end of strokes that finish with downward motion. However, a comparison of measured and predicted force traces indicates that the wing generates large force transients at the start of strokes beginning with downward motion that cannot be explained by the changes in aerodynamic angle of attack and wing velocity, which are incorporated in the quasi-steady model (Fig.8B,D). The most likely explanation for these large peaks is that wake capture is enhanced by downward motion and attenuated by upward motion at the start of the stroke. By moving downwards in the direction of the descending wake, the wing may intercept a stronger flow field than if it moved upwards. However, a rigorous test of this hypothesis will require a detailed visualization of the velocity distribution during stroke reversal.
The interaction between stroke deviation and wake capture is more extreme in the case of figure-of-eight patterns, in which the two half-strokes are mirror images of one another and can thus both begin with downward motion. The quasi-steady model predicts that mean lift should be maximal with a +20° deviation, which means that the wing begins each stroke by moving upwards. The enhanced performance of upward strokes is due to the fact that the angle of attack is high in the middle of the stroke when the wing velocity is greatest. The measured forces do not, however, exhibit this enhanced performance with small positive deviations. On the contrary, the net performance drops off more quickly with positive deviation (upward motion at start of strokes) than with negative deviation (downward motion at start of stroke). This result is again consistent with a more effective wake capture due to downward motion at the start of translation.
In spite of the large differences in instantaneous forces, the magnitudes of time-averaged forces in all cases are similar for the oval and figure-of-eight trajectories (Fig.9). For oval trajectories, the measured forces behave as predicted, with maximum lift and drag coefficients occurring at zero deviation and decreasing for higher absolute deviation. In contrast, the figure-of-eight trajectories behave quite differently from the corresponding quasi-steady estimates, varying in both magnitude and behavior. For quasi-steady estimates, the lift maximum occurs at a 20° positive deviation, whereas the drag maximum increases monotonically with increasing positive deviations (open red circles in Fig.9C,D respectively). In contrast, the measured values in both cases are maximal at zero deviation (filled blue and red circles in Fig.9C,D) and then decrease for higher absolute deviations.
Drag and efficiency
The values of mean profile drag
obtained in the current study are much higher than values in the
literature based either on steady-state measurements of wing profiles
(Vogel, 1967) or
calculations based on Reynolds number (Ellington, 1984d; Lehmann and Dickinson,
1997; Wakeling and Ellington,
1997). In particular, the equation:
![]() |
which represents the rough average of the Reynolds-number-dependence of a flat plate parallel and perpendicular to flow (Ellington, 1984c), may grossly underestimate the profile drag on a flapping wing. In the current study, this equation predicts a
value of approximately 0.7 for large stroke amplitudes, which is less than 25% of the measured values for kinematic patterns capable of sustaining flight
(
3). One obvious reason for this large discrepancy is that the process of wing rotation generates large instantaneous values of drag which, even though lasting only a small portion of the stroke, substantially increase the mean values. However, even the drag contributed by delayed stall during translation (
2) is higher than that predicted by equation 7. These large values for the mean drag coefficient have several important implications for flight energetics and muscle physiology. First, profile power, which scales linearly with the mean drag coefficient, may be several times greater than predicted in estimates using equation 7. Thus, under many conditions, profile power, and not induced power, may represent the dominant component in aerodynamic power. Further, estimates of total mechanical power based solely on induced power are likely grossly to underestimate flight costs. The elevated values of drag also call into question previous attempts to estimate elastic storage, which depend upon inertial power costs being high relative to aerodynamic costs (Dickinson and Lighton, 1995). If the ratio of aerodynamic cost to inertial cost is high, the energetic benefits of elastic storage are greatly reduced. Second, because previous estimates of mechanical power may be several times too low, calculations of muscle efficiency, based on the ratio of mechanical power to total metabolic power, may be substantially lower than actual values. Such estimates for the muscle efficiency of asynchronous flight muscle range from 8 to 12%. When corrected for the underestimate of drag, we find that previous estimates of muscle efficiency presented for fruit flies (Lehmann and Dickinson, 1997) should be adjusted to 2025%. If correct, these adjusted values help to solve a perplexing problem in insect flight physiology. Because Ca2+ cycling represents a sizeable fraction of the energetic cost in cyclic contractions, stretch-activated muscle should accrue an energetic savings and concomitant increase in mechanical efficiency. This prediction was not supported, however, by the previous estimates of efficiency. After correction with more accurate values of drag, the estimates of mechanical efficiency in asynchronous muscle are now consistent with the reduced cost of Ca2+ cycling.
In addition to indicating that the mean level of drag may be higher than previously expected, the results also show that the value of drag is quite sensitive to subtle changes in wing kinematics (Fig.5C, Fig.7C, Fig.9D). Mapping an estimate of profile power indicates that this sensitivity of the drag coefficient to wing kinematics translates into a sensitivity of flight cost as well (Fig.10A,B,C). This result again has consequences for estimates of flight costs and mechanical efficiency from free-flight kinematics based solely on measures of stroke frequency and amplitude (Ellington, 1984d). In particular, for a given stroke amplitude, profile power varies depending upon the timing and duration of wing rotation (Fig.10B). Thus, without knowledge of the precise stroke kinematics, estimates of free-flight energetics must be viewed with caution.
Concluding remarks
The data presented in this paper represent a comprehensive analysis of the effects that wing kinematics have on the generation of aerodynamic forces in flapping flight. The results reveal a complex system in which subtle alterations in stroke kinematics can have quite large effects on force production. In future studies, it will be an important challenge to determine how real insects use this sensitively to adjust the balance of forces and force moments as they steer and maneuver through their environment.
![]() |
List of symbols |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
![]() |
Acknowledgments |
---|
![]() |
References |
---|
![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
---|
Arbas, E. (1986). Control of hindlimb posture by wind-sensitive hairs and antennae during locust flight. J. Comp. Physiol. A 159, 849857.[Medline]
Bennett, L. (1970). Insect flight: lift and the rate of change of incidence. Science 167, 177179.
Bennett, L. (1977). Clap and fling aerodynamics an experimental evaluation. J. Exp. Biol. 69, 261272.
Cloupeau, M., Devillers, J. F. and Devezeaux, D. (1979). Direct measurements of instantaneous lift in desert locust; comparison with Jensens experiments on detached wings. J. Exp. Biol. 80, 115.
Daniel, T. L. (1984). Unsteady aspects of aquatic locomotion. Am. Zool. 24, 121134.
David, C. T. (1978). The relationship between body angle and flight speed in free flying Drosophila. Physiol. Ent. 3, 191195.
Dickinson, M. H. (1994). The effects of wing rotation on unsteady aerodynamic performance at low Reynolds numbers. J. Exp. Biol. 192, 179206.
Dickinson, M. H. and Götz, K. G. (1993). Unsteady aerodynamic performance of model wings at low Reynolds numbers. J. Exp. Biol. 174, 4564.
Dickinson, M. H. and Götz, K. G. (1996). The wake dynamics and flight forces of the fruit fly Drosophila melanogaster. J. Exp. Biol. 199, 20852104.
Dickinson, M. H., Lehmann, F.-O. and Götz, K. G. (1993). The active control of wing rotation by Drosophila. J. Exp. Biol. 182, 173189.
Dickinson, M. H., Lehmann, F.-O. and Sane, S. P. (1999). Wing rotation and the aerodynamic basis of insect flight. Science 284, 19541960.
Dickinson, M. H. and Lighton, J. R. B. (1995). Muscle efficiency and elastic storage in the flight motor of Drosophila. Science 128, 8789.
Dickinson, M. and Tu, M. (1997). The function of dipteran flight muscle. Comp. Biochem. Physiol. 116A, 223238.
Ellington, C. P. (1984a). The aerodynamics of hovering insect flight. II. Morphological parameters. Phil. Trans. R. Soc. Lond. B 305, 1740.
Ellington, C. P. (1984b). The aerodynamics of hovering insect flight. III. Kinematics. Phil. Trans. R. Soc. Lond. B 305, 4178.
Ellington, C. P. (1984c). The aerodynamics of hovering insect flight. IV. Aerodynamic mechanisms. Phil. Trans. R. Soc. Lond. B 305, 79113.
Ellington, C. P. (1984d). The aerodynamics of hovering insect flight. VI. Lift and power requirements. Phil. Trans. R. Soc. Lond. B 305, 145181.
Ellington, C., Vandenberg, C., Willmott, A. and Thomas, A. (1996). Leading-edge vortices in insect flight. Nature 384, 626630.
Ennos, A. R. (1989). The kinematics and aerodynamics of the free flight of some Diptera. J. Exp. Biol. 142, 4985.
Götz, K. G., Hengstenberg, B. and Biesinger, R. (1979). Optomotor control of wing beat and body posture in Drosophila. Biol. Cybernetics 35, 101112.
Jensen, M. (1956). Biology and physics of locust flight. III. The aerodynamics of locust flight. Phil. Trans. R. Soc. Lond. B 239, 511552.
Lehmann, F.-O. and Dickinson, M. H. (1997). The changes in power requirements and muscle efficiency during elevated force production in the fruit fly Drosophila melanogaster. J. Exp. Biol. 200, 11331143.
Lehmann, F. and Dickinson, M. (1998). The control of wing kinematics and flight forces in fruit flies (Drosophila spp.). J. Exp. Biol. 201, 385401.
Lehmann, F. and Götz, K. (1996). Activation phase ensures kinematic efficacy in flight-steering muscles of Drosophila melanogaster. J. Comp. Physiol. A 179, 311322.[Medline]
Liu, H., Ellington, C. P., Kawachi, K., VandenBerg, C. and Willmott, A. P. (1998). A computational fluid dynamic study of hawkmoth hovering. J. Exp. Biol. 201, 461477.
Lorez, M. (1995). Neural control of hindleg steering in flight in the locust. J. Exp. Biol. 198, 869875.
Marden, J. H. (1987). Maximum lift production during takeoff in flying animals. J. Exp. Biol. 130, 235258.
Maxworthy, T. (1979). Experiments on the Weis-Fogh mechanism of lift generation by insects in hovering flight. Part 1. Dynamics of the fling. J. Fluid Mech. 93, 4763.
May, M. L. and Hoy, R. R. (1990). Leg-induced steering in flying crickets. J. Exp. Biol. 151, 485488.
Robertson, R. and Johnson, A. (1993). Collision avoidance of flying locusts: steering torques and behaviour. J. Exp. Biol. 183, 3560.
Ruppell, G. (1989). Kinematic analysis of symmetrical flight manoeuvres of Odonata. J. Exp. Biol. 144, 1343.
Sedov, L. I. (1965). Two-Dimensional Problems in Hydrodynamics and Aerodynamics, pp. 2030. New York: Interscience Publishers.
Spedding, G. R. (1993). On the significance of unsteady effects in the aerodynamic performance of flying animals. Contemp. Math. 141, 401419.
Spedding, G. R. and Maxworthy, T. (1986). The generation of circulation and lift in a rigid two-dimensional fling. J. Fluid Mech. 165, 247272.
Tu, M. S. and Dickinson, M. H. (1994). Modulation of negative work output from a steering muscle of the blowfly Calliphora vicina. J. Exp. Biol. 192, 207224.
Vogel, S. (1967). Flight in Drosophila. III. Aerodynamic characteristics of fly wings and wing models. J. Exp. Biol. 46, 431443.
Wakeling, J. E. and Ellington, C. P. (1997). Dragonfly flight. III. Lift and power requirements. J. Exp. Biol. 200, 583600.
Wang, J. (2000). Vortex shedding and frequency selection in flapping flight. J. Fluid Mech. 410, 323341.
Wilkin, P. J. (1990). The instantaneous force on a desert locust, Schistocerca gregaria (Orthoptera: Acrididae), flying in a wing tunnel. J. Kansas Ent. Soc. 63, 316328.
Wortmann, M. and Zarnack, W. (1993). Wing movements and lift regulation in the flight of desert locusts. J. Exp. Biol. 182, 5769.
Zanker, J. M. (1988). How does lateral abdomen deflection contribute to flight control of Drosophila melanogaster. J. Comp. Physiol. A 162, 581588.
Zanker, J. M. (1990a). The wing beat of Drosophila melanogaster. I. Kinematics. Phil. Trans. R. Soc. Lond. B 327, 118.
Zanker, J. M. (1990b). The wing beat of Drosophila melanogaster. III. Control. Phil. Trans. R. Soc. Lond. B 327, 4564.
Zanker, J. M. and Götz, K. G. (1990). The wing beat of Drosophila melanogaster. II. Dynamics. Phil. Trans. R. Soc. Lond. B 327, 1944.