Dynamic flight stability of a hovering model insect: lateral motion

The lateral dynamic flight stability of a hovering model insect （dronefly） was studied using the method of computational fluid dynamics to compute the stability derivatives and the techniques of eigenvalue and eigenvector analysis for solving the equations of motion. The main results are as following. （i） Three natural modes of motion were identified： one unstable slow divergence mode （mode 1）, one stable slow oscillatory mode （mode 2）, and one stable fast subsidence mode （mode 3）. Modes 1 and 2 mainly consist of a rotation about the horizontal longitudinal axis （x-axis） and a side translation; mode 3 mainly consists of a rotation about the x-axis and a rotation about the vertical axis. （ii） Approximate analytical expressions of the eigenvalues are derived, which give physical insight into the genesis of the natural modes of motion. （iii） For the unstable divergence mode, td, the time for initial disturbances to double, is about 9 times the wingbeat period （the longitudinal motion of the model insect was shown to be also unstable and td of the longitudinal unstable mode is about 14 times the wingbeat period）. Thus, although the flight is not dynamically stable, the instability does not grow very fast and the insect has enough time to control its wing motion to suppress the disturbances.     

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）.     

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).     

Control for going from hovering to small speed flight of a model insect

The longitudinal steady-state control for going from hovering to small speed flight of a model insect is studied, using the method of computational fluid dynamics to compute the aerodynamic derivatives and the techniques based on the linear theories of stability and control for determining the non-zero equilibrium points. Morphological and certain kinematical data of droneflies are used for the model insect. A change in the mean stroke angle （δФ） results in a horizontal forward or backward flight; a change in the stroke amplitude （δФ） or a equal change in the down- and upstroke angles of attack （δα1） results in a vertical climb or decent; a proper combination of δФ and δФ controls （or δФ and δα1 controls） can give a flight of any （small） speed in any desired direction.     

Lateral dynamic flight stability of hovering insects: theory vs. numerical simulation

In the present paper, the lateral dynamic flight stability properties of two hovering model insects are predicted by an approximate theory based on the averaged model, and computed by numerical simulation that solves the complete equations of motion coupled with the Navier-Stokes equations. Comparison between the theoretical and simulational results provides a test to the validity of the assumptions made in the theory. One of the insects is a model dronefly which has relatively high wingbeat frequency (164 Hz) and the other is a model hawkmoth which has relatively low wingbeat frequency (26 Hz). The following conclusion has been drawn. The theory based on the averaged model works well for the lateral motion of the dronefly. For the hawkmoth, relatively large quantitative differences exist between theory and simulation. This is because the lateral non-dimensional eigenvalues of the hawkmoth are not very small compared with the non-dimensional flapping frequency (the largest lateral non-dimensional eigenvalue is only about 10% smaller than the non-dimensional flapping frequency). Nevertheless, the theory can still correctly predict variational trends of the dynamic properties of the hawkmoth’s lateral motion.     

Asymptotic Technical Stability in Measure of the Nonlinear Controllable Motion of an Elastic Flight Vehicle

The asymptotic technical stability conditions for the nonlinear controllable motion of an elongated elastic flight vehicle are formulated. A longitudinal vertical flight at constant velocity is considered. The variation in the cross section, transverse strains, and vibrations of the vehicle are taken into account. The asymptotic technical stability criteria formulated depend on the key parameters of the controlled process, including the contribution of the weight due to curving of the vehicle axis to the transverse load and the aerodynamic force during a vertical flight     

Unsteady aerodynamic forces and power requirements of a bumblebee in forward flight

Aerodynamic forces and power requirements in forward flight in a bumblebee (Bombus terrestris) were studied using the method of computational fluid dynamics. Actual wing kinematic data of free flight were used in the study (the speed ranges from 0 m/s to 4.5 m/s; advance ratio ranges from 0–0.66). The bumblebee employs the delayed stall mechanism and the fast pitching-up rotation mechanism to produce vertical force and thrust. The leading-edge vortex does not shed in the translatory phase of the half-strokes and is much more concentrated than that of the fruit fly in a previous study. At hovering and low-speed flight, the vertical force is produced by both the half-strokes and is contributed by wing lift; at medium and high speeds, the vertical force is mainly produced during the downstroke and is contributed by both wing lift and wing drag. At all speeds the thrust is mainly produced in the upstroke and is contributed by wing drag. The power requirement at low to medium speeds is not very different from that of hovering and is relatively large at the highest speed (advance ratio 0.66), i.e. the power curve is J-shaped. Except at the highest flight speed, storing energy elastically can save power up to 20%–30%. At the highest speed, because of the large increase of aerodynamic torque and the slight decrease of inertial torque (due to the smaller stroke amplitude and stroke frequency used), the power requirement is dominated by aerodynamic power and the effect of elastic storage of energy on power requirement is limited.The project supported by the National Natural Science Foundation of China (10232010) and the National Aeronautic Science fund of China (03A51049)The English text was polished by Xing Zhang.     

Aeroelastic stability analysis considering a continuous flight envelope

This paper presents a new methodology to analyze aeroelastic stability in a continuous range of flight envelope with varying parameter of velocity and altitude. The focus of the paper is to demonstrate that linear matrix inequalities can be used to evaluate the aeroelastic stability in a region of flight envelope instead of a single point, like classical methods. The proposed methodology can also be used to study if a system remains stable during an arbitrary motion from one point to another in the flight envelope, i.e., when the problem becomes time-variant. The main idea is to represent the system as a polytopic differential inclusion system using rational function approximation to write the model in time domain. The theory is outlined and simulations are carried out on the benchmark AGARD 445.6 wing to demonstrate the method. The classical pk-method is used for comparing results and validating the approach. It is shown that this method is efficient to identify stability regions in the flight envelope.     

Stabilization control of a hovering model insect: lateral motion

Our previous study shows that the lateral disturbance motion of a model drone fly does not have inherent stability (passive stability),because of the existence of an unstable divergence mode.But drone flies are observed to fly stably.Constantly active control must be applied to stabilize the flight.In this study,we investigate the lateral stabilization control of the model drone fly.The method of computational fluid dynamics is used to compute the lateral control derivatives and the techniques of eigenvalue and eigenvector analysis and modal decomposition are used for solving the equations of motion.Controllability analysis shows that although inherently unstable,the lateral disturbance motion is controllable.By feeding back the state variables (i.e.lateral translation velocity,yaw rate,roll rate and roll angle,which can be measured by the sensory system of the insect) to produce anti-symmetrical changes in stroke amplitude and/or in angle of attack between the left and right wings,the motion can be stabilized,explaining why the drone flies can fly stably even if the flight is passively unstable.     

Stabilization control of a bumblebee in hovering and forward flight

Our previous study shows that the hovering and forward flight of a bumblebee do not have inherent stability （passive stability）. But the bumblebees are observed to fly stably. Stabilization control must have been applied. In this study, we investigate the longitudinal stabilization control of the bumblebee. The method of computational fluid dynamics is used to compute the control derivatives and the techniques of eigenvalue and eigenvector analysis and modal decomposition are used for solving the equations of motion. Controllability analysis shows that at all flight speeds considered, although inherently unstable, the flight is controllable. By feedbacking the state variables, i.e. vertical and horizontal velocities, pitching rate and pitch angle （which can be measured by the sensory system of the insect）, to produce changes in stroke angle and angle of attack of the wings, the flight can be stabilized, explaining why the bumblebees can fly stably even if they are passively unstable.     

Thermal convection in a cylindrical annulus under a combined effect of the radial and vertical gravity

The stability of the flow of a dielectric fluid confined in a cylindrical annulus submitted to a radial temperature gradient and a radial electric field is investigated theoretically and experimentally. The radial temperature gradient induces a vertical Archimedean buoyancy and a radial dielectrophoretic buoyancy. These two forces intervene simultaneously in the destabilization of the flow, leading to the occurrence of four types of modes depending on the relative intensity of these two buoyancies and on the fluid's properties: hydrodynamic and thermal modes that are axisymmetric and oscillatory, stationary columnar modes and electric modes which are stationary and non-axisymmetric modes. Experiments performed in a parabolic flight show the existence of non-axisymmetric modes that should be either columnar or helicoidal vortices.     

旋转三角形绳系卫星编队系统动态稳定性分析

本文研究了平动点附近的旋转三角形绳系卫星编队系统的动态稳定性问题.为了便于分析,假设系绳始终保持张紧状态,且不计系绳的质量.为了建立更加符合实际的数学模型,考虑了姿态运动和轨道运动的相互影响,即姿轨耦合,并在此基础上推导出了考虑姿轨耦合的非线性动力学方程.最后用数值方法模拟了在不同的旋转速度情况下系统的动态稳定性,模拟结果表明旋转速率越大系统稳定性越好.     

Dynamic Processes in a Spheroidal Seismic-Damping Mechanism

A mathematical model of a variable-structure system of solids with rolling friction and unilateral sliding friction bonds is described. The model can be applied to seismic isolation mechanisms. Conditions for transition between the possible structures of the equations of state are formulated. The behavior of the system with kinematically defined motion of the base is analyzed as an example__________Translated from Prikladnaya Mekhanika, Vol. 41, No. 1, pp. 119–127, January 2005.     

Modelling of the diffusion of fluid particles in turbulent flows approaching a circular cylinder by random Fourier modes and random flight

To investigate the diffusion of fluid particles around a cylinder in a turbulent flow, we have developed two new types of model for simulating the trajectory of particles:(1) a model combining random Fourier modes and random flight (RF); (2) a pure kinematic simulation (KS) by random Fourier modes. In model 1 the large-scale turbulence is simulated by a sum of random Fourier modes varying in space and time, and the small-scale random motion of particles is simply modelled by an Itô type of stochastic differential equation with a memory time comparable to the Lagrangian time scaleT _s ^L of the small-scale motion. In model 2, both large- and small-scale turbulence is simulated using random Fourier modes. The change of turbulence around the cylinder is modelled by rapid distortion theory (RDT), although the small-scale motion of particles in the RF model is simply assumed to keep the homogeneous random behaviour. These models give very similar and realistic trajectories showing rapid changes of direction due to the small-scale motion.     

Dynamic stability of a pipe conveying fluid with an uncertain computational model

This paper deals with the problem of a pipe conveying fluid of interest in several engineering applications, such as micro-systems or drill-string dynamics. The deterministic stability analysis developed by Paidoussis and Issid (1974) is extended to the case for which there are model uncertainties induced by modeling errors in the computational model. The aim of this work is twofold: (1) to propose a probabilistic model for the fluid–structure interaction considering modeling errors and (2) to analyze the stability and reliability of the stochastic system. The Euler–Bernoulli beam model is used to model the pipe and the plug flow model is used to take into account the internal flow in the pipe. The resulting differential equation is discretized by means of the finite element method and a reduced-order model is constructed from some eigenmodes of the beam. A probabilistic approach is used to model uncertainties in the fluid–structure interaction. The proposed strategy takes into account global uncertainties related to the noninertial coupled fluid forces (related to damping and stiffness). The resulting random eigenvalue problem is used to analyze flutter and divergence unstable modes of the system for different values of the dimensionless flow speed. The numerical results show the random response of the system for different levels of uncertainty, and the reliability of the system for different dimensionless speeds and levels of uncertainty.     

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.     

Technical stability of nonlinear states of an elastic vehicle in vertical flight

Sufficient conditions of technical stability of nonlinear dynamic states of extended elastic flying systems in controlled longitudinal vertical flight are obtained. In these flying systems, the effect of variation of their crosssectional area, transverse strains, and oscillations is taken into account. The formulated criteria of technical stability depend on the basic parameters of the process controlled, namely, on the increment of the transverse load due to the curvature of the axis of the system and aerodynamic forces during vertical flight.     

Stability domains of the vertical equilibrium state of a triple simple pendulum

The stability boundaries for the vertical equilibrium position of a triple simple pendulum subject to asymmetric follower force are analyzed. The effect of the upper spring and the magnitude of the follower force on the roots of the characteristic equation is studied Translated from Prikladnaya Mekhanika, Vol. 44, No. 10, pp. 122–134, October 2008.     

单桩竖向动力阻抗计算模型研究

运用土动力学和结构动力学原理,同时考虑桩周土的弱化效应和桩-土界面的相对滑移效应,利用数理方程方法分别求解单桩与桩周近场土域及远场土域的振动方程,建立了竖向荷载作用下单桩动力阻抗函数的简化计算方法,提出了一种改进的非线性动力Winkler模型,确定了模型中各物理元件的参数,进而通过对比分析验证了计算模型的合理性,从而为桩-土-上部结构耦合系统的非线性分析奠定了基础。