On mathematical modelling of insect flight dynamics in the context of micro air vehicles.

We discuss some aspects of mathematical modelling relevant to the dynamics of insect flight in the context of insect-like flapping-wing micro air vehicles (MAVs). MAVs are small flying vehicles developed to reconnoître in confined spaces. This requires power-efficient, highly-manoeuvrable, low-speed flight with stable hover. All of these attributes are present in insect flight and hence the focus on reproducing the functionality of insect flight by engineering means. Empirical research on insect flight dynamics is limited by experimental difficulties. Force and moment measurements require tethering the animal whose behaviour may then differ from free flight. The measurements are made when the insect actively tries to control its flight, so that its open-loop dynamics cannot be observed. Finally, investigation of the sensory-motor system responsible for flight is even more challenging. Despite these difficulties, much empirical progress has been made recently. Further progress, especially in the context of MAVs, can be achieved by the complementary information derived from appropriate mathematical modelling. The focus here is on a means of computing the data not easily available from experiments and also on making mathematical predictions to suggest new experiments. We consider two aspects of mathematical modelling for insect flight dynamics. The first one is theoretical (computational), as opposed to empirical, generation of the aerodynamic data required for the six-degrees-of-freedom equations of motion. For this purpose we first explain insect wing kinematics and the salient features of the corresponding flow. In this context, we show that aerodynamic modelling is a feasible option for certain flight regimes, focusing on a successful example of modelling hover. Such modelling progresses from the first principles of fluid mechanics, but relies on simplifications justified by the known flow phenomenology and/or geometric and kinematic symmetries. This is relevant to six types of fundamental manoeuvres, which we define as those flight conditions for which only one component of the translational and rotational body velocities is nonzero and constant. The second aspect of mathematical modelling for insect flight dynamics addressed here deals with the periodic character of the aerodynamic force and moment production. This leads to consideration of the types of solutions of nonlinear equations forced by nonlinear oscillations. In particular, the mechanism of synchronization seems relevant and should be investigated further.

Nonlinear time-periodic models of the longitudinal flight dynamics of desert locusts Schistocerca gregaria.

Previous studies of insect flight control have been statistical in approach, simply correlating wing kinematics with body kinematics or force production. Kinematics and forces are linked by Newtonian mechanics, so adopting a dynamics-based approach is necessary if we are to place the study of insect flight on its proper physical footing. Here we develop semi-empirical models of the longitudinal flight dynamics of desert locusts Schistocerca gregaria. We use instantaneous force-moment measurements from individual locusts to parametrize the nonlinear rigid body equations of motion. Since the instantaneous forces are approximately periodic, we represent them using Fourier series, which are embedded in the equations of motion to give a nonlinear time-periodic (NLTP) model. This is a proper mathematical generalization of an earlier linear-time invariant (LTI) model of locust flight dynamics, developed using previously published time-averaged versions of the instantaneous force recordings. We perform various numerical simulations, within the fitted range of the model, and across the range of body angles used by free-flying locusts, to explore the likely behaviour of the locusts upon release from the tether. Solutions of the NLTP models are compared with solutions of the nonlinear time-invariant (NLTI) models to which they reduce when the periodic terms are dropped. Both sets of models are unstable and therefore fail to explain locust flight stability fully. Nevertheless, whereas the measured forces include statistically significant harmonic content up to about the eighth harmonic, the simulated flight trajectories display no harmonic content above the fundamental forcing frequency. Hence, manoeuvre control in locusts will not directly reflect subtle changes in the higher harmonics of the wing beat, but must operate on a coarser time-scale. A state-space analysis of the NLTP models reveals orbital trajectories that are impossible to capture in the LTI and NLTI models, and inspires the hypothesis that asymptotic orbital stability is the proper definition of stability in flapping flight. Manoeuvre control on the scale of more than one wing beat would then consist in exciting transients from one asymptotically stable orbit to another. We summarize these hypotheses by proposing a limit-cycle analogy for flapping flight control and suggest experiments for verification of the limit-cycle control analogy hypothesis.

Insect-like flapping wing mechanism based on a double spherical Scotch yoke.

We describe the rationale, concept, design and implementation of a fixed-motion (non-adjustable) mechanism for insect-like flapping wing micro air vehicles in hover, inspired by two-winged flies (Diptera). This spatial (as opposed to planar) mechanism is based on the novel idea of a double spherical Scotch yoke. The mechanism was constructed for two main purposes: (i) as a test bed for aeromechanical research on hover in flapping flight, and (ii) as a precursor design for a future flapping wing micro air vehicle. Insects fly by oscillating (plunging) and rotating (pitching) their wings through large angles, while sweeping them forwards and backwards. During this motion the wing tip approximately traces a "figure-of-eight" or a "banana" and the wing changes the angle of attack (pitching) significantly. The kinematic and aerodynamic data from free-flying insects are sparse and uncertain, and it is not clear what aerodynamic consequences different wing motions have. Since acquiring the necessary kinematic and dynamic data from biological experiments remains a challenge, a synthetic, controlled study of insect-like flapping is not only of engineering value, but also of biological relevance. Micro air vehicles are defined as flying vehicles approximately 150 mm in size (hand-held), weighing 50-100g, and are developed to reconnoitre in confined spaces (inside buildings, tunnels, etc.). For this application, insect-like flapping wings are an attractive solution and hence the need to realize the functionality of insect flight by engineering means. Since the semi-span of the insect wing is constant, the kinematics are spatial; in fact, an approximate figure-of-eight/banana is traced on a sphere. Hence a natural mechanism implementing such kinematics should be (i) spherical and (ii) generate mathematically convenient curves expressing the figure-of-eight/banana shape. The double spherical Scotch yoke design has property (i) by definition and achieves (ii) by tracing spherical Lissajous curves.

On aerodynamic modelling of an insect-like flapping wing in hover for micro air vehicles.

This theoretical paper discusses recent advances in the fluid dynamics of insect and micro air vehicle (MAV) flight and considers theoretical analyses necessary for their future development. The main purpose is to propose a new conceptual framework and, within this framework, two analytic approaches to aerodynamic modelling of an insect-like flapping wing in hover in the context of MAVs. The motion involved is periodic and is composed of two half-cycles (downstroke and upstroke) which, in hover, are mirror images of each other. The downstroke begins with the wing in the uppermost and rearmost position and then sweeps forward while pitching up and plunging down. At the end of the half-cycle, the wing flips, so that the leading edge points backwards and the wing's lower surface becomes its upper side. The upstroke then follows by mirroring the downstroke kinematics and executing them in the opposite direction. Phenomenologically, the interpretation of the flow dynamics involved, and adopted here, is based on recent experimental evidence obtained by biologists from insect flight and related mechanical models. It is assumed that the flow is incompressible, has low Reynolds number and is laminar, and that two factors dominate: (i) forces generated by the bound leading-edge vortex, which models flow separation; and (ii) forces due to the attached part of the flow generated by the periodic pitching, plunging and sweeping. The first of these resembles the analogous phenomenon observed on sharp-edged delta wings and is treated as such. The second contribution is similar to the unsteady aerodynamics of attached flow on helicopter rotor blades and is interpreted accordingly. Theoretically, the fluid dynamic description is based on: (i) the superposition of the unsteady contributions of wing pitching, plunging and sweeping; and (ii) adding corrections due to the bound leading-edge vortex and wake distortion. Viscosity is accounted for indirectly by imposing the Kutta condition on the trailing edge and including the influence of the vortical structure on the leading edge. Mathematically, two analytic approaches are proposed. The first derives all the quantities of interest from the notion of circulation and leads to tractable integral equations. This is an application of the von Kármán-Sears unsteady wing theory and its nonlinear extensions due to McCune and Tavares; the latter can account for the bound leading-edge vortex and wake distortion. The second approach uses the velocity potential as the central concept and leads to relatively simple ordinary differential equations. It is a combination of two techniques: (i) unsteady aerodynamic modelling of attached flow on helicopter rotor blades; and (ii) Polhamus's leading-edge suction analogy. The first of these involves both frequency-domain (Theodorsen style) and time-domain (indicial) methods, including the effects of wing sweeping and returning wake. The second is a nonlinear correction accounting for the bound leading-edge vortex. Connections of the proposed framework with control engineering and aeroelasticity are pointed out.