Studies of Store-Induced Limit-Cycle Oscillations Using a Model with Full System Nonlinearities

The authors investigate limit-cycle oscillations of a wing/store configuration. Unlike typical aeroelastic studies that are based upon a linearized form of the governing equations, herein full system nonlinearities are retained, and include transonic flow effects, coupled responses from the structure, and store-related kinematics and dynamics. Unsteady aerodynamic loads are modeled with the equations from transonic small disturbance theory. The structural dynamics for the cantilevered wing are modeled by the nonlinear equations of motion for a beam. The effects of general store-placement are modeled by the nonlinear equations of motion related to the position-induced nonlinear kinematics. Chordwise deformations of the wing surface, as well as pylon and store flexibility, are assumed negligible. Nonlinear responses are studied by examining bifurcation and related response characteristics using direct simulation. Particular attention is given to cases for which large-time, time-dependent behavior is dependent on initial conditions, as observed for some configurations in flight test. Comparisons of results in which selective nonlinearities are excluded indicate that the accurate prediction of nonlinear responses such as limit cycle oscillations (LCOs) may depend upon consideration of all nonlinearities related to the full system.     

Virtual flight simulation of a dual rotor micro air vehicle

In this paper, a new computational method is developed based on computational fluid dynamics (CFD) coupled with rigid body dynamics (RBD) and flight control law in an in-house programmed source code. The CFD solver is established based on momentum source method, preconditioning method, lower–upper symmetric Gauss–Seidel iteration method, and moving overset grid method. Two-equation shear–stress transport k ? ω turbulence model is employed to close the governing equations. Third-order Adams prediction-correction method is used to couple CFD and RBD in the inner iteration. The wing-rock motion of the delta wing is simulated to validate the capability of the computational method for virtual flight simulation. Finally, the developed computational method is employed to simulate the longitudinal virtual flight of a dual rotor micro air vehicle (MAV). Results show that the computational method can simulate the virtual flight of the dual rotor MAV.     

Dynamic flight stability of a model dronefly in vertical flight

The dynamic flight stability of a model dronefly in hovering and upward flight is studied.The method of computational fluid dynamics is used to compute the stability derivatives and the techniques of eigenvalue and eigenvector used to solve the equations of motion.The major finding is as following.Hovering flight of the model dronefly is unstable because of the existence of an unstable longitudinal and an unstable lateral natural mode of motion.Upward flight of the insect is also unstable,and the instability increases as the upward flight speed increases.Inertial force generated by the upward flight velocity coupled with the disturbance in pitching angular velocity is responsible for the enhancement of the instability.     

Computational study of a pitching bio-inspired corrugated airfoil

The phenomenon of insect flight has been of scientific interest for many years and is more recently inspiring modern engineering devices such as Micro Aerial Vehicles (MAVs). Insect flight is characterized by unsteady fluid dynamics at low Reynolds numbers. The importance of viscous effects to the successful flapping flight of insects has been identified and with the current state of computing and Computational Fluid Dynamics (CFD) these effects can now be studied in detail. The present work attempts to simplify this complex phenomenon by considering symmetric oscillating rotational motion of a wing (pitching). What is of interest in this study is how the shape of a corrugated idealized insect wing affects the performance and flow characteristics around the pitching wing. Two dimensional CFD on an oscillating wing has been performed and reported. Measurements were taken to ensure the accuracy of the computational solution and the results validated against experimental PIV results. A range of frequencies and rotational amplitudes have been investigated. Lift and drag coefficients have been analyzed for all cases to quantify the effects of unsteady flow features on the performance of the oscillating wing. It was found that the wing shape used in this study resulted in the viscous features formed on the top of the wing exhibiting high sensitivity to the oscillating conditions and these influenced the performance of the wing. The flow features formed on the bottom of the wing remained similar throughout the cases tested. In the pitching regime this wing profile did not perform as well as published results for smooth airfoils in terms of thrust and propulsion efficiency. However this may be due to reduced frequency effects becoming important at our high pitching amplitude which need to be investigated further. There may be other oscillatory regimes that more accurately represent flapping flight in which the corrugated foil outperforms a smooth counterpart but these are yet to be investigated. Further research in this area may help answer the question as to how evolutionarily significant other benefits of a corrugated wing, such as being light and strong, are compared to its aerodynamic properties, the present results seem to favor the former.     

Aeroelastic study of flexible flapping wings by a coupled lattice Boltzmann-finite element approach with immersed boundary method

In this paper, the behavior of two-dimensional symmetric flapping wings moving in a viscous fluid is investigated. Harmonic motion is applied to idealize flying organisms with flexible wings and extensive testing is carried out to investigate the resultant flight behavior related to the ability to take-off or accelerate the flapping wing system away from a starting location. Special attention is paid to analyze the effect of the main mechanical parameters, as well as the effect of lateral wind on flight performances. Moreover, aiming to investigate the possible benefits of flying in flocks, a couple of synchronously flapping wings is considered in addition to the single arrangement. The numerical simulations are performed by solving the fluid–structure interaction problem through a strongly coupled partitioned approach. Fluid dynamics are modeled at the mesoscopic scale by the lattice Boltzmann method. The resulting macroscopic quantities are derived, as usual, based on the statistical molecular-level interpretation.Wings are modeled by geometrically nonlinear, elastic beam finite elements and structure dynamics is solved by the time discontinuous Galerkin method. Fluid–structure interface conditions are handled using the immersed boundary method. The resultant numerical approach combines simplicity and high computational efficiency. A Monte Carlo simulation strategy is employed to characterize the flight behavior subjected to lateral wind. Various scenarios are discussed.     

Effect of the aerodynamic-induced parametric excitation on the longitudinal stability of hovering MAVs/insects

Longitudinal flight dynamics of hovering micro-air-vehicles and insects is considered. The natural oscillatory flapping motion of the wing leads to time-periodic stability derivatives (aerodynamic loads due to perturbation in the body-motion variables). Hence these terms play the role of parametric excitation on the flight dynamics. The main objective of this work is to assess the effects of these aerodynamic-induced parametric excitation terms that are neglected by the averaging analysis. The method of multiple scales is used to determine a second-order uniform expansion for the response of the time-periodic system at hand. The proposed approach is applied to the hovering flight dynamics of five insects that cover a wide range of operating frequency ratios to assess the applicability of the averaging analysis.     

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 controllability perspective of dynamic soaring

Dynamic soaring is an exquisite flying technique to acquire energy from the atmospheric wind shear. In this study, a geometric nonlinear controllability analysis of an unmanned aerial vehicle (UAV) under dynamic soaring conditions is performed. To achieve such an objective, the state-of-the-art mathematical tools of nonlinear controllability are summarized and presented to an aeronautical engineering audience. The dynamic soaring optimal control problem is then formulated and solved numerically. The controllability of the UAV along the optimal soaring trajectory is analyzed. More importantly, the geometric nonlinear controllability characteristics of generic flight dynamics are analyzed in the presence and absence of wind shear to provide a controllability explanation for the role of wind shear in the physics of dynamic soaring flight. It is found that the wind shear is instrumental in ensuring controllability as it allows the UAV attitude controls (pitch and roll) to play the role of thrust in controlling the flight path angle. The presented analysis represents a controllability-based mathematical proof for the energetics of flight physics.     

昆虫飞行的空气动力学

昆虫是最早出现、数量最多和体积最小的飞行者. 它们能悬停、跃升、急停、快速加速和转弯, 飞行技巧十分高超. 由于尺寸小, 因而翅膀的相对速度很小, 从而进行上述飞行所需的升力系数很大. 但昆虫翅膀的雷诺数又很低. 它们是如何在低雷诺数下产生高升力的, 是流体力学和生物学工作者都十分关心的问题. 近年来这一领域有了许多研究进展. 该文对这些进展进行综述, 并对今后工作提一些建议. 因2005 年前的工作已在几篇综述文章有了详细介绍, 该文主要介绍2005 年以来的工作. 首先简述昆虫翅的拍动运动及昆虫绕流的基本方程和相似参数; 然后对2005 年之前的工作做一简要回顾. 之后介绍2005 年后的进展, 依次为: 运动学观测; 前缘涡; 翅膀柔性变形及皱褶的影响; 拍动翅的尾涡结构; 翼/身、左右翅气动干扰及地面效应; 微小昆虫; 蝴蝶与蜻蜓; 机动飞行. 最后为对今后工作的建议.      

Extracting micro air vehicles aerodynamic forces and coefficients in free flight using visual motion tracking techniques

This paper describes a methodology to extract aerial vehicles’ aerodynamic characteristics from visually tracked trajectory data. The technique is being developed to study the aerodynamics of centimeter-scale aircraft and develop flight simulation models. Centimeter-scale aircraft remains a largely unstudied domain of aerodynamics, for which traditional techniques like wind tunnels and computational fluid dynamics have not yet been fully adapted and validated. The methodology takes advantage of recent progress in commercial, vision-based, motion-tracking systems. This system dispenses from on-board navigation sensors and enables indoor flight testing under controlled atmospheric conditions. Given the configuration of retro-reflective markers affixed onto the aerial vehicle, the vehicle’s six degrees-of-freedom motion can be determined in real time. Under disturbance-free conditions, the aerodynamic forces and moments can be determined from the vehicle’s inertial acceleration, and furthermore, for a fixed-wing vehicle, the aerodynamic angles can be plotted from the vehicle’s kinematics. By combining this information, we can determine the temporal evolution of the aerodynamic coefficients, as they change throughout a trajectory. An attractive feature of this technique is that trajectories are not limited to equilibrium conditions but can include non-equilibrium, maneuvering flight. Whereas in traditional wind-tunnel experiments, the operating conditions are set by the experimenter, here, the aerodynamic conditions are driven by the vehicle’s own dynamics. As a result, this methodology could be useful for characterizing the unsteady aerodynamics effects and their coupling with the aircraft flight dynamics, providing insight into aerodynamic phenomena taking place at centimeter scale flight.     

Near wake vortex dynamics of a hovering hawkmoth

Numerical investigation of vortex dynamics in near wake of a hovering hawkmoth and hovering aerodynamics is conducted to support the development of a biology-inspired dynamic flight simulator for flapping wingbased micro air vehicles. Realistic wing-body morphologies and kinematics are adopted in the numerical simulations. The computed results show 3D mechanisms of vortical flow structures in hawkmoth-like hovering. A horseshoe-shaped primary vortex is observed to wrap around each wing during the early down- and upstroke; the horseshoe-shaped vortex subsequently grows into a doughnut-shaped vortex ring with an intense jet-flow present in its core, forming a downwash. The doughnut-shaped vortex rings of the wing pair eventu- ally break up into two circular vortex rings as they propagate downstream in the wake. The aerodynamic yawing and rolling torques are canceled out due to the symmetric wing kinematics even though the aerodynamic pitching torque shows significant variation with time. On the other hand, the time- varying the aerodynamics pitching torque could make the body a longitudinal oscillation over one flapping cycle.     

三角翼大迎角流场结构和气动力特性的数值分析研究

采用数值计算方法对亚音速三角翼纵向及带有小侧滑情况下的流场结构和气动力特性进行了计算。文中给出了三角翼大迎角纵向情况下气动力、机翼前缘分离涡轴线位置和旋涡破裂位置随迎角的变化规律，以及带有横侧小扰动和小侧滑情况下流场结构的非对称性对气动力的影响。计算结果表明与实验结果符合较好。     

Effects of tip perturbation and wing locations on rolling oscillation induced by forebody vortices

The wing rock motion is frequently suffered by a wing-body configuration with low swept wing at high angle of attack. It is found from our experimental study that the tip perturbation and wing longitudinal locations affect significantly the wing rock motion of a wing-body. The natural tip perturbation would make the wing rock motion of a nondeterministic nature and an artificial mini-tip perturbation would make the wing rock motion deterministic. The artificial tip perturbation would, as its circumferential location is varied, generate three different types of motion patterns： （1） limit cycle oscillation, （2） irregular oscillation, （3） equilibrium state with tiny oscillation. The amplitude of rolling oscillation corresponding to the limit cycle oscillatory motion pattern is decreased with the wing location shifting downstream along the body axis.     

Saturation-based actuation for flapping MAVs in hovering and forward flight

Stringent weight and size constraints on flapping-wing microair-vehicles dictate minimal actuation. Unfortunately, hovering and forward flight require different wing motions and, as such, independent actuators. Therefore, either a hovering or a forward-flight requirement should be included in the mission and design statements of a flapping-wing microair-vehicle. This work proposes a design for an actuation mechanism that would provide the required kinematics in each flight condition using only one actuator. The idea is to exploit the nonlinear dynamics of the flapping wing to induce the saturation phenomenon. One physical spring in the plunging direction is needed along with a feedback of the plunging angle into the control torque of the actuator in the back and forth flapping direction. By detuning the feedback gains away from the saturation requirement, we obtain the flapping kinematics required for hovering. In contrast, tuning the feedback gains to induce the saturation phenomenon transfers the motion into the plunging direction. Moreover, the actuating torque (in the back and forth flapping direction) would then provide a direct control over the amplitude of the plunging motion, while the amplitude of the actuated flapping motion saturates and does not change as the amplitude of the actuating torque increases.     

Dynamic flight stability of hovering insects

The equations of motion of an insect with flapping wings are derived and then simplified to that of a flying body using the “rigid body” assumption. On the basis of the simplified equations of motion, the longitudinal dynamic flight stability of four insects (hoverfly, cranefly, dronefly and hawkmoth) in hovering flight is studied (the mass of the insects ranging from 11 to 1,648 mg and wingbeat frequency from 26 to 157 Hz). The method of computational fluid dynamics is used to compute the aerodynamic derivatives and the techniques of eigenvalue and eigenvector analysis are used to solve the equations of motion. The validity of the “rigid body” assumption is tested and how differences in size and wing kinematics influence the applicability of the “rigid body” assumption is investigated. The primary findings are: (1) For insects considered in the present study and those with relatively high wingbeat frequency (hoverfly, drone fly and bumblebee), the “rigid body” assumption is reasonable, and for those with relatively low wingbeat frequency (cranefly and howkmoth), the applicability of the “rigid body” assumption is questionable. (2) The same three natural modes of motion as those reported recently for a bumblebee are identified, i.e., one unstable oscillatory mode, one stable fast subsidence mode and one stable slow subsidence mode. (3) Approximate analytical expressions of the eigenvalues, which give physical insight into the genesis of the natural modes of motion, are derived. The expressions identify the speed derivative M _u (pitching moment produced by unit horizontal speed) as the primary source of the unstable oscillatory mode and the stable fast subsidence mode and Z _w (vertical force produced by unit vertical speed) as the primary source of the stable slow subsidence mode. The project supported by the National Natural Science Foundation of China (10232010 and 10472008).     

Dynamic flight stability of a bumblebee in forward flight

The longitudinal dynamic flight stability of a bumblebee in forward flight is studied. The method of computational fluid dynamics is used to compute the aerodynamic derivatives and the techniques of eigenvalue and eigenvector analysis are employed for solving the equations of motion. The primary findings are as the following. The forward flight of the bumblebee is not dynamically stable due to the existence of one （or two） unstable or approximately neutrally stable natural modes of motion. At hovering to medium flight speed [flight speed Ue = （0-3.5）m s^-1; advance ratio J = 0-0.44], the flight is weakly unstable or approximately neutrally stable; at high speed （Ue = 4.5 m s^-1; J = 0.57）, the flight becomes strongly unstable （initial disturbance double its value in only 3.5 wingbeats）.     

折叠翼飞行器的动力学建模与稳定控制

折叠翼飞行器在变形过程中,其动力学模型呈现多刚体、多自由度和强非线性特点,同时气动力/力矩、压心、质心和转动惯量等参数也会大幅度变化,严重影响飞行稳定性.由此,本论文将对飞行器的多刚体动力学建模与变形稳定控制进行研究.基于凯恩方法建立了折叠翼飞行器的多刚体动力学模型,并从中得到了变形所产生的附加力和力矩表达式.通过气动计算拟合出气动参数与折叠角之间的函数关系,由此分析了不同折叠角速度下飞行器的纵向动态特性,结果表明,折叠翼飞行器变形过程中速度、高度和俯仰角均会发生变化,飞行器无法保持稳定飞行.为此提出了一种基于自抗扰理论的飞行器变形过程中的稳定控制方法.将折叠翼飞行器纵向非线性动力学模型中存在的非线性项、耦合项以及参数时变项都视为系统内外总扰动,利用扩张状态观测器对总扰动进行实时估计和补偿,针对补偿后的系统设计PD控制器,实现了速度通道和高度通道的解耦控制.通过Lyapunov稳定性原理证明了系统的稳定性,并进行数学仿真验证.仿真结果表明,基于自抗扰理论设计的稳定控制器能够解决飞行器变形所带来的强非线性和参数时变等问题,保证飞行器的高精度稳定控制.     

DYNAMIC MODELING AND STABILITY CONTROL OF FOLDING WING AIRCRAFT 1)

折叠翼飞行器在变形过程中,其动力学模型呈现多刚体、多自由度和强非线性特点,同时气动力/力矩、压心、质心和转动惯量等参数也会大幅度变化,严重影响飞行稳定性. 由此,本论文将对飞行器的多刚体动力学建模与变形稳定控制进行研究.基于凯恩方法建立了折叠翼飞行器的多刚体动力学模型,并从中得到了变形所产生的附加力和力矩表达式.通过气动计算拟合出气动参数与折叠角之间的函数关系,由此分析了不同折叠角速度下飞行器的纵向动态特性, 结果表明,折叠翼飞行器变形过程中速度、高度和俯仰角均会发生变化,飞行器无法保持稳定飞行.为此提出了一种基于自抗扰理论的飞行器变形过程中的稳定控制方法.将折叠翼飞行器纵向非线性动力学模型中存在的非线性项、耦合项以及参数时变项都视为系统内外总扰动,利用扩张状态观测器对总扰动进行实时估计和补偿, 针对补偿后的系统设计PD控制器,实现了速度通道和高度通道的解耦控制.通过Lyapunov稳定性原理证明了系统的稳定性, 并进行数学仿真验证. 仿真结果表明,基于自抗扰理论设计的稳定控制器能够解决飞行器变形所带来的强非线性和参数时变等问题,保证飞行器的高精度稳定控制.     

Time-resolved wake structure and kinematics of bat flight

We present synchronized time-resolved measurements of the wing kinematics and wake velocities for a medium sized bat, Cynopterus brachyotis, flying at low-medium speed in a closed-return wind tunnel. Measurements of the motion of the body and wing joints, as well as the resultant wake velocities in the Trefftz plane are recorded at 200 Hz (approximately 28–31 measurements per wing beat). Circulation profiles are found to be quite repeatable although variations in the flight profile are visible in the wake vortex structures. The circulation has almost constant strength over the middle half of the wing beat (defined according the vertical motion of the wrist, beginning with the downstroke). A strong streamwise vortex is observed to be shed from the wingtip, growing in strength during the downstroke, and persisting during much of the upstroke. At relatively low flight speeds (4.3 m/s), a closed vortex structure behind the bat is postulated.