Effects of unsteady deformation of flapping wing on its aerodynamic forces

Effects of unsteady deformation of a flapping model insect wing on its aerodynamic force production are studied by solving the Navier-Stokes equations on a dynamically deforming grid.Aerodynamic forces on the flapping wing are not much affected by considerable twist,but affected by camber deformation.The effect of combined camber and twist deformation is similar to that of camber deformation.With a deformation of 6% camber and 20°twist(typical values observed for wings of many insects),lift is increased bv 10%～20%and lift-to-drag ratio by around 10%compared with the case of a rigid flat-plate wing.As a result.the deformation can increase the maximum lift coefficient of an insect.and reduce its power requirement for flight.For example,for a hovering bumblebee with dynamically deforming wings(6?mber and 20°twist),aerodynamic power required is reduced by about 16%compared with the case of rigid wings.     

Numerical simulation of flapping‐wing insect hovering flight at unsteady flow

A computational fluid dynamics (CFD) analysis was conducted to study the unsteady aerodynamics of a virtual flying bumblebee during hovering flight. The integrated geometry of bumblebee was established to define the shape of a three‐dimensional virtual bumblebee model with beating its wings, accurately mimicking the three‐dimensional movements of wings during hovering flight. The kinematics data of wings documented from the measurement to the bumblebee in normal hovering flight aided by the high‐speed video. The Navier–Stokes equations are solved numerically. The solution provides the flow and pressure fields, from which the aerodynamic forces and vorticity wake structure are obtained. Insights into the unsteady aerodynamic force generation process are gained from the force and flow‐structure information. The CFD analysis has established an overall understanding of the viscous and unsteady flow around the virtual flying bumblebee and of the time course of instantaneous force production, which reveals that hovering flight is dominated by the unsteady aerodynamics of both the instantaneous dynamics and also the past history of the wing. A coherent leading‐edge vortex with axial flow and the attached wingtip vortex and trailing edge vortex were detected. The leading edge vortex, wing tip vortex and trailing edge vortex, which caused by the pressure difference between the upper and the lower surface of wings. The axial flow, which include the spanwise flow and chordwise flow, is derived from the spanwise pressure gradient and chordwise pressure gradient, will stabilize the vortex and gives it a characteristic spiral conical shape. Copyright © 2006 John Wiley & Sons, Ltd.     

Experimental investigation of the effect of chordwise flexibility on the aerodynamics of flapping wings in hovering flight

Ornithopters or mechanical birds produce aerodynamic lift and thrust through the flapping motion of their wings. Here, we use an experimental apparatus to investigate the effects of a wing's twisting stiffness on the generated thrust force and the power required at different flapping frequencies. A flapping wing system and an experimental set-up were designed to measure the unsteady aerodynamic and inertial forces, power usage and angular speed of the flapping wing motion. A data acquisition system was set-up to record important data with the appropriate sampling frequency. The aerodynamic performance of the vehicle under hovering (i.e., no wind) conditions was investigated. The lift and thrust that were produced were measured for different flapping frequencies and for various wings with different chordwise flexibilities. The results show the manner in which the elastic deformation and inertial flapping forces affect the dynamical behavior of the wing. It is shown that the generalization of the actuator disk theory is, at most, only valid for rigid wings, and for flexible wings, the power P varies by a power of about 1.0  of the thrust T. This aerodynamic information can also be used as benchmark data for unsteady flow solvers.     

Lift and power requirements of hovering insect flight

Lift and power requirements for hovering flight of eight species of insects are studied by solving the Navier-Stokes equation numerically. The solution provides velocity and pressure fields, from which unsteady aerodynamic forces and moments are obtained. The inertial torque of wing mass are computed analytically. The wing length of the insects ranges from 2 mm (fruit fly) to 52 mm (hawkmoth); Reynolds numbers Re (based on mean flapping speed and mean chord length) ranges from 75 to 3850. The primary findings are shown in the following: (1) Either small (R = 2mm, Re = 75), medium (R ≈ 10 mm, Re ≈ 500) or large (R ≈ 50 mm, Re ≈ 4 000) insects mainly employ the same high-lift mechanism, delayed stall, to produce lift in hovering flight. The midstroke angle of attack needed to produce a mean lift equal to the insect weight is approximately in the range of 25° to 45°, which is approximately in agreement with observation. (2) For the small insect (fruit fly) and for the medium and large insects with relatively small wingbeat frequency (cranefly, ladybird and hawkmoth), the specific power ranges from 18 to 39W·kg^-1 , the major part of the power is due to aerodynamic force, and the elastic storage of negative work does not change the specific power greatly. However for medium and large insects with relatively large wingbeat frequency (hover fly, dronefly, honey bee and bumble bee), the specific power ranges from 39 to 61 W·kg^-1 , the major part of the power is due to wing inertia, and the elastic storage of negative work can decrease the specific power by approximately 33%. (3) For the case of power being mainly contributed by aerodynamic force (fruit fly, cranefly, ladybird and hawkmoth), the specific power is proportional to the product of the wingbeat frequency, the stroke amplitude, the wing length and the drag-to-lift ratio. For the case of power being mainly contributed by wing inertia (hoverfly, dronefly, honey bee and bumble bee), the specific power (without elastic storage) is proportional to the product of the cubic of wingbeat frequency, the square of the stroke amplitude, the square of the wing length and the ratio of wing mass to insect mass.     

A novel MAV with treadmill motion of wing

The Magnus effect is well known phenomena for producing high lift values from spinning symmetrical geometries such as cylinders, spheres, or disks. But, the Magnus force may also be produced by treadmill motion of aerodynamic bodies. To accomplish this, the skin of aerodynamic bodies may circulate with a constant circumferential speed. Here, a novel wing with treadmill motion of skin is introduced which may generate lift at zero air speeds. The new wing may lead to micro aerial vehicle configurations for vertical take-off or landing. To prove the concept, the NACA0015 aerofoil section with circulating skin is computationally investigated. Two cases of stationary air and moving air are studied. It is observed that lift can be generated in stationary air although drag force is also high. For moving air, the lift and drag forces may be adopted between the incidence angles 20° to 25° where lift can posses high values and drag can remain moderate.     

Experimental study on interaction of aerodynamics with flexible wings of flapping vehicles in hovering and cruise flight

Flapping wings are promising lift and thrust generators, especially for very low Reynolds numbers. To investigate aeroelastic effects of flexible wings (specifically, wing’s twisting stiffness) on hovering and cruising aerodynamic performance, a flapping-wing system and an experimental setup were designed and built. This system measures the unsteady aerodynamic and inertial forces, power usage, and angular speed of the flapping wing motion for different flapping frequencies and for various wings with different chordwise flexibility. Aerodynamic performance of the vehicle for both no wind (hovering) and cruise condition was investigated. Results show how elastic deformations caused by interaction of inertial and aerodynamic forces with the flexible structure may affect specific power consumption. This information was used here to find a more suitable structural design. The best selected design in our tests performs up to 30% better than others (i.e., less energy consumption for the same lift or thrust generation). This measured aerodynamic information could also be used as a benchmarking data for unsteady flow solvers.     

Measurement of Thrust and Lift Forces Associated With Drag of Compliant Flapping Wing for Micro Air Vehicles Using a New Test Stand Design

Compliant wing designs have the potential of improving flapping wing Micro-Air Vehicles (MAVs). Designing compliant wings requires a detailed understanding of the effect of compliance on the generation of thrust and lift forces. The low force and high-frequency measurements associated with these forces necessitated a new versatile test stand design that uses a 250 g load cell along with a rigid linear air bearing to minimize friction and the dynamic behavior of the test stand while isolating only the stationary thrust or lift force associated with drag generated by the wing. Moreover, this stand is relatively inexpensive and hence can be easily utilized by wing designers to optimize the wing compliance and shape. The frequency response of the wing is accurately resolved, along with wing compliance on the thrust and lift profiles. The effects of the thrust and lift force generated as a function of flapping frequency were also determined. A semi-empirical aerodynamic model of the thrust and lift generated by the flapping wing MAV on the new test stand was developed and used to evaluate the measurements. This model accounted for the drag force and the effects of the wing compliance. There was good correlation between the model predictions and experimental measurements. Also, the increase in average thrust due to increased wing compliance was experimentally quantified for the first time using the new test stand. Thus, our measurements for the first time reveal the detrimental influence of excessive compliance on drag forces during high frequency operation. In addition, we were also able to observe the useful effect of compliance on the generation of extra thrust at the beginning and end of upstrokes and downstrokes of the flapping motion.     

Aerodynamic mechanism of forces generated by twisting model-wing in bat flapping flight

The aerodynamic mechanism of the bat wing membrane Mong the lateral border of its body is studied. The twist-morphing that alters the angle of attack （AOA） along the span-wise direction is observed widely during bat flapping flight. An assumption is made that the linearly distributed AOA is along the span-wise direction. The plate with the aspect ratio of 3 is used to model a bat wing. A three-dimensional （3D） unsteady panel method is used to predict the aerodynamic forces generated by the flapping plate with leading edge separation. It is found that, relative to the rigid wing flapping, twisting motion can increase the averaged lift by as much as 25% and produce thrust instead of drag. Furthermore, the aerodynamic forces （lift/drag） generated by a twisting plate-wing are similar to those of a pitching rigid-wing, meaning that the twisting in bat flight has the same function as the supination/pronation motion in insect flight.     

A flow control mechanism in wing flapping with stroke asymmetry during insect forward flight

A theoretical modeling approach as well as an unsteady analytical method is used to study aerodynamic characteristics of wing flapping with asymmetric stroke-cycles in connection with an oblique stroke plane during insect forward flight. It is revealed that the aerodynamic asymmetry between the downstroke and the upstroke due to stroke-asymmetrical flapping is a key to understand the flow physics of generation and modulation of the lift and the thrust. Predicted results for examples of given kinematics validate more specifically some viewpoints that the wing lift is more easily produced when the forward speed is higher and the thrust is harder, and the lift and the thrust are generated mainly during downstroke and upstroke, respectively. The effects of three controlling parameters, i.e. the angles of tilted stroke plane, the different downstroke duration ratios, and the different angles of attack in both down- and up-stroke, are further discussed. It is found that larger oblique angles of stroke planes generate larger thrust but smaller lift; larger downstroke duration ratios lead to larger thrust, while making little change in lift and input aerodynamic power; and again, a small increase of the angle of attack in downstroke or upstroke may cause remarkable changes in aerodynamic performance in the relevant stroke.The project supported by the National Natural Science Foundation of China (10072066, 90305009) and the Chinese Academy of Sciences (KJCX-SW-L04, KJCX2-SW-L2)The English text was polished by Ron Marshall.     

The aerodynamic forces and pressure distribution of a revolving pigeon wing

The aerodynamic forces acting on a revolving dried pigeon wing and a flat card replica were measured with a propeller rig, effectively simulating a wing in continual downstroke. Two methods were adopted: direct measurement of the reaction vertical force and torque via a forceplate, and a map of the pressures along and across the wing measured with differential pressure sensors. Wings were tested at Reynolds numbers up to 108,000, typical for slow-flying pigeons, and considerably above previous similar measurements applied to insect and hummingbird wing and wing models. The pigeon wing out-performed the flat card replica, reaching lift coefficients of 1.64 compared with 1.44. Both real and model wings achieved much higher maximum lift coefficients, and at much higher geometric angles of attack (43°), than would be expected from wings tested in a windtunnel simulating translating flight. It therefore appears that some high-lift mechanisms, possibly analogous to those of slow-flying insects, may be available for birds flapping with wings at high angles of attack. The net magnitude and orientation of aerodynamic forces acting on a revolving pigeon wing can be determined from the differential pressure maps with a moderate degree of precision. With increasing angle of attack, variability in the pressure signals suddenly increases at an angle of attack between 33° and 38°, close to the angle of highest vertical force coefficient or lift coefficient; stall appears to be delayed compared with measurements from wings in windtunnels.     

昆虫飞行的高升力机理

对近年来关于昆虫产生非定常高升力的研究进行了综述和归纳.这方面的工作对生物学研究和微型飞行器等微型机械的仿生设计有重要意义.研究表明:果蝇等昆虫翅膀的拍动运动可产生很大的非定常升力,其平均值是定常值的2~3倍,足够平衡昆虫的重量,并有较大的富余用于机动飞行;产生高升力有三个因素:一是拍动开始阶段翅的快速加速运动,二是拍动中的不失速机制,三是拍动结束阶段翅的快速上仰运动.人们从能耗的角度考察了这些非定常高升力机制的正确性和可行性.当作悬停飞行的果蝇用以上机制产生平衡其重量的升力时,其比功率(支持单位身体质量所需的功率)约为29W/kg, 生化/机械效率约为17%. 这些值与人们基于对昆虫肌肉力学特性的研究所预估的值接近.果蝇前飞时,其比功率随速度变化的曲线是一J形曲线,而不是象飞机或鸟的那样是一U形曲线;这与人们基于昆虫新陈代谢率的测量数据所推断的结果一致.对于蜻蜒等(功能上)有前、后两对翅膀的昆虫,有以下初步结果:翅的下拍主要产生升力,上挥主要产生推力;下拍时的平均升力系数可达2~3,十分大,上挥时的平均推力系数可达1~2, 也很大,它们主要由非定常效应产生;前、后翅的相互干扰并未起增大升力和推力的作用,反而有一定的不利作用.      

Aerodynamic Performance of a Gliding Swallowtail Butterfly Wing Model

In the present study, we perform a wind-tunnel experiment to investigate the aerodynamic performance of a gliding swallowtail-butterfly wing model having a low aspect ratio. The drag, lift and pitching moment are directly measured using a 6-axis force/torque sensor. The lift coefficient increases rapidly at attack angles less than 10° and then slowly at larger attack angles. The lift coefficient does not fall off rapidly even at quite high angles of attack, showing the characteristics of low-aspect-ratio wings. On the other hand, the drag coefficient increases more rapidly at higher angles of attack due to the increase in the effective area responsible for the drag. The maximum lift-to-drag ratio of the present modeled swallowtail butterfly wing is larger than those of wings of fruitfly and bumblebee, and even comparable to those of wings of birds such as the petrel and starling. From the measurement of pitching moment, we show that the modeled swallowtail butterfly wing has a longitudinal static stability. Flow visualization shows that the flow separated from the leading edge reattaches on the wing surface at α < 15°, forming a small separation bubble, and full separation occurs at α ≥ 15°. On the other hand, strong wing-tip vortices are observed in the wake at α ≥ 5° and they are an important source of the lift as well as the main reason for broad stall. Finally, in the absence of long hind-wing tails, the lift and longitudinal static stability are reduced, indicating that the hind-wing tails play an important role in enhancing the aerodynamic performance.     

蜻蜓非均匀柔性前翼扑动飞行的气动性能研究

为了探究柔性对于蜻蜓前翼在扑动向前飞行时的气动性能, 本文根据蜻蜓前翼的实际参数建立蜻蜓前翼模型, 提出了两种柔性分布方式即均匀柔性分布和沿蜻蜓前翼弦向的变柔性分布. 本文通过STAR-CCM+软件, 首先采用重叠网格和双向流固耦合技术, 用于实现蜻蜓前翼的扑动流固耦合, 其次通过改变蜻蜓前翼固体区域的杨氏模量函数从而实现蜻蜓前翼的两种不同柔性分布. 结果表明, 在均匀柔性分布条件下, 柔性翼在杨氏模量较小时的升力系数和阻力系数曲线的变化规律滞后于刚性翼半周期并且给飞行增加阻力, 但是随着杨氏模量的逐渐增加即柔性逐渐减小, 蜻蜓前翼受到的阻力减小, 获得的推力增加且推力给予蜻蜓前飞的动量增量、加速度以及时均推力系数先增加后减小. 在合理的非均匀柔性分布条件下, 柔性翼显著提高推力系数峰值和时均推力系数, 在扑动前飞时, 给予蜻蜓前翼较大的动量增量以及加速度. 两种柔性分布方式的蜻蜓前翼与刚性翼对比之下, 蜻蜓前翼在柔性为非均匀柔性分布时可以获得更好的气动性能.      

Experimental investigation on aerodynamic performance of a flapping wing vehicle in forward flight

The aerodynamic performance of a flexible membrane flapping wing has been investigated here. For this purpose, a flapping-wing system and an experimental set-up were designed to measure the unsteady aerodynamic forces of the flapping wing motion. A one-component force balance was set up to record the temporal variations of aerodynamic forces. The flapping wing was studied in a large low-speed wind tunnel. The lift and thrust of this mechanism were measured for different flapping frequencies, angles of attack and for various wind tunnel velocities. Results indicate that the thrust increases with the flapping frequency. An increase in the wind tunnel speed and flow angle of attack leads to reduction in the thrust value and increases the lift component. The aerodynamic and performance parameters were nondimensionalized. Appropriate models were introduced which show its aerodynamic performance and may be used in the design process and also optimization of the flapping wing.     

昆虫拍翼方式的非定常流动物理再探讨

基于提出的理论模化方法来探讨昆虫拍翼方式的非定常流动物理. 以悬停飞行为 例,通过对拍翼运动的分析,不仅解释了昆虫利用高频拍翼的方式为何能够克服低雷诺数带 来的气动局限性(St \gg 1/Re),而且还指出高升力产生和调节的3个流动 控制因素：(1) 由于拍翼的变速运动即时引起了流体动力响应,这种附加惯性效应 可产生瞬时的高升力; (2) 保持前缘涡不脱离翼面有助于减少升力的下降; (3) 增大后缘涡的强度并加速其脱离后缘能够有效地提高升力.     

On the aerodynamic characteristics of hovering rigid and flexible hawkmoth-like wings

Insect wings are subjected to fluid, inertia and gravitational forces during flapping flight. Owing to their limited rigidity, they bent under the influence of these forces. Numerical study by Hamamoto et al. (Adv Robot 21(1–2):1–21, 2007) showed that a flexible wing is able to generate almost as much lift as a rigid wing during flapping. In this paper, we take a closer look at the relationship between wing flexibility (or stiffness) and aerodynamic force generation in flapping hovering flight. The experimental study was conducted in two stages. The first stage consisted of detailed force measurement and flow visualization of a rigid hawkmoth-like wing undergoing hovering hawkmoth flapping motion and simple harmonic flapping motion, with the aim of establishing a benchmark database for the second stage, which involved hawkmoth-like wing of different flexibility performing the same flapping motions. Hawkmoth motion was conducted at Re = 7,254 and reduced frequency of 0.26, while simple harmonic flapping motion at Re = 7,800 and 11,700, and reduced frequency of 0.25. Results show that aerodynamic force generation on the rigid wing is governed primarily by the combined effect of wing acceleration and leading edge vortex generated on the upper surface of the wing, while the remnants of the wake vortices generated from the previous stroke play only a minor role. Our results from the flexible wing study, while generally supportive of the finding by Hamamoto et al. (Adv Robot 21(1–2):1–21, 2007), also reveal the existence of a critical stiffness constant, below which lift coefficient deteriorates significantly. This finding suggests that although using flexible wing in micro air vehicle application may be beneficial in term of lightweight, too much flexibility can lead to deterioration in flapping performance in terms of aerodynamic force generation. The results further show that wings with stiffness constant above the critical value can deliver mean lift coefficient almost the same as a rigid wing when executing hawkmoth motion, but lower than the rigid wing when performing a simple harmonic motion. In all cases studied (7,800 ≤ Re ≤ 11,700), the Reynolds number does not alter the force generation significantly.     

Aerodynamics and mechanisms of elementary morphing models for flapping wing in forward flight of bat

Large active wing deformation is a significant way to generate high aerodynamic forces required in bat's flapping flight. Besides the twisting, elementary morphing models of a bat wing are proposed, including wing-bending in the spanwise direction, wing-cambering in the chordwise direction, and wing area-changing. A plate of aspect ratio 3 is used to model a bat wing, and a three-dimensional unsteady panel method is used to predict the aerodynamic forces. It is found that the cambering model has great positive influence on the lift, followed by the area-changing model and then the bending model. Further study indicates that the vortex control is a main mechanism to produce high aerodynamic forces. The mechanisms of aerodynamic force enhancement are asymmetry of the cambered wing and amplification effects of wing area-changing and wing bending. Lift and thrust are generated mainly during downstroke, and they are almost negligible during upstroke by the integrated morphing model-wing.     

Design,manufacturing, and flight testing of a fixed wing micro air vehicle with Zimmerman planform

An optimized and comprehensive method is used to design and manufacture a fixed wing micro air vehicle (MAV) with Zimmerman planform. The design process includes four stages which are the specification of the flight mission, determination of the best aspect ratio, identification of the optimum wing loading and thrust loading values, and estimation of the weight of the structural components of the MAV. To do this, various statistical and analytical methods are utilized. Based on an aerodynamic analysis, the results show that an optimum aspect ratio that maximizes the performance of the Zimmerman MAV for a well-defined cruise speed is determined. Considering six possible flights, a constraint analysis is performed and an optimum wing loading value is determined. It is shown that the computational method is beneficial to determine the exact masses for the structural components including the wing, fuselage, and vertical tail. Using the 3D panel method, the determination of the shape of the reflexed airfoil for the MAV is successfully done by minimizing the drag force and the angle of attack to use less powerful motor and avoid any stall effect, respectively. A stability analysis is then performed to check the safe flight of the designed vehicle. During test flight, the results show that the designed Zimmerman MAV satisfies the pre-defined specification. The final characteristics of the manufactured MAV are: wingspan of 44 cm, weight of 450 g, aspect ratio of 1.51, cruise speed of 20 m/s, and flight endurance of 20 min.     

Characterization of vortical structures and loads based on time-resolved PIV for asymmetric hovering flapping flight

Flight agility, resistance to gusts, capability to hover coupled with a low noise generation might have been some of the reasons why insects are among the oldest species observed in nature. Biologists and aerodynamicists focused on analyzing such flight performances for diverse purposes: understanding the essence of flapping wings aerodynamics and applying this wing concept to the development of micro-air vehicles (MAVs). In order to put into evidence the fundamentally non-linear unsteady mechanisms responsible for the amount of lift generated by a flapping wing (Dickinson et al. in Science 284:1954–1960, 1999), experimental and numerical studies were carried out on typical insect model wings and kinematics. On the other hand, in the recent context of MAVs development, it is of particular interest to study simplified non-biological flapping configurations which could lead to lift and/or efficiency enhancement. In this paper, we propose a parametrical study of a NACA0012 profile undergoing asymmetric hovering flapping motions at Reynolds 1000. On the contrary to normal hovering, which has been widely studied as being the most common configuration observed in the world of insects, asymmetric hovering is characterized by an inclined stroke plane. Besides the fact that the vertical force is hence a combination of both lift and drag (Wang in J Exp Biol 207:1137–1150, 2004), the specificity of such motions resides in the vortex dynamics which present distinct behaviours, whether the upstroke angle of attack leads to a partially attached or a strong separated flow, giving more or less importance to the wake capture phenomenon. A direct consequence of the previous remarks relies on the enhancement of aerodynamic efficiency with asymmetry. If several studies reported results based on the asymmetric flapping motion of dragonfly, only few works concentrated on parametrizing asymmetric motions (e.g. Wang in Phys Rev Lett 85:2216–2219, 2000). The present study relies on TR-PIV measurements which allow determination of the vorticity fields and provide a basis to evaluate the resulting unsteady forces through the momemtum equation approach.