全机跨音速大迎角绕流的数值分析

应用有限体积法求解Ｅｕｌｅｒ方程计算了真实飞机外形的跨音速大迎角绕流流场及空气动力特性。本方法能自动捕获大后掠机翼前缘脱体涡，侧缘涡，机身体涡以及它们间的相互干扰。计算的机翼压力分布及全机气动力系数与实验值符合良好。表明了本文所采用的计算网格生成技术及流场计算方法是有效、可行的。     

伸缩翼飞机纵向飞行动力学研究

研究了可变翼展飞机在变翼展过程中的气动特性及纵向飞行动响应。首先基于多刚体系统动力学方法建立了可伸缩机翼飞机的动力学模型,进而获得了纵向飞行动力学方程;然后分析了飞机质心、质量矩、惯性矩随翼展的变化规律;最后通过实时气动力计算和纵向飞行动力学方程联立求解,获得了不同翼展下的飞机升力、阻力、力矩及纵向飞行动响应。结果表明:在相同迎角条件下,大翼展飞机的升力性能与小翼展飞机的升力性能相比有较大的优势;未加控制时机翼收缩导致飞机阻力和升力减小,飞行速度增大;为维持变翼过程的稳定飞行,需要对俯仰运动施加控制。     

飞机翼身干扰流动分离机理研究

采用CFD(计算流体力学)方法,结合风洞试验结果,对飞机翼身干扰流动分离机理进行了研究。分析了机翼根部马蹄涡、边角涡的内部结构和形成过程;研究了机身和机翼的压力分布及附面层,进而总结给出了飞机翼根流动分离机理;同时也验证了CFD方法对翼身干扰计算的准确性。研究结果表明:机翼前缘对气流的阻挡作用形成较高的逆压梯度,使气流产生回流形成马蹄涡,马蹄涡影响翼根表面附面层,使机翼机身附面层掺混失稳,向下游流动发展形成边角涡。     

机翼-机身-短舱组合体的气动力计算

引言包括机翼,机身和短舱的飞机气动力计算是一个大型课题。在亚声速和超声速范围内,有限基本解是公认的解理想位势流的最好的方法。由于对气动力计算的要求越来越高,对 ...     

三角翼动态大迎角气动力特性数值分析研究

采用数值计算方法,对三角翼从 0°上仰至 90°的动态流场和气动力特性进行了计算,并对俯仰角速度对三角翼流场和气动力特性的影响进行了计算分析。给出了三角翼纵向动态情况下的气动力系数变化,特别是大迎角横侧力矩系数的变化特征。结果表明,随着机翼俯仰角速度的提高,前缘分离涡破裂位置相对滞后,机翼升力和阻力系数明显增加,机翼抵抗旋涡非对称破裂的能力明显增强。     

具有迟滞非线性的高超声速机翼颤振分析

本文采用两自由度的二元机翼模型,研究高超声速机翼由于气动弹性引起的机翼颤振问题.考虑了由于机翼连接部位的松弛和摩擦引起的机翼迟滞非线性特性的影响,采用三阶活塞理论给出高超声速机翼的非线性气动力和气动力矩.通过数值模拟,获得系统的时域响应曲线和Poincare图,分析发现,随着系统参数的变化,二元机翼会出现极限环、分岔等复杂的动力学行为,并发现迟滞非线性参数对系统极限环幅值、分岔和混沌特性有较大影响.     

非定常俯抑振荡下的横向喷流数值模拟

采用高精度格式数值求解RANS方程,研究了定常状态下横向喷流流场,压力分布计算结果与实验结果基本吻合,并捕捉到喷流干扰流场中多种流动结构.在非定常计算过程中,飞行器的振动引起了法向力和俯仰力矩系数的相位滞后,推力放大因子随俯仰角周期变化.飞行器振动过程中,喷流流场的动态气动特性与稳态喷流有明显的区别,因此在利用横向喷流对飞行器进行姿态控制时,应该考虑由于飞行器姿态的变化对横向喷流所产生的非定常影响问题.     

旋转流场下模型大幅俯仰运动实验研究

在中航气动院FL-8低速风洞中,采用单自由度振荡机构进行了旋转流场下大幅俯仰运动的气动特性实验研究。模型在绕风速轴连续旋转的同时,进行给定频率和振幅绕体轴的俯仰振荡运动。测量了模型的动态气动特性,着重分析了不同运动参数对模型气动特性的影响。实验结果表明,旋转速度的存在使大振幅俯仰振荡实验中的滚转力矩和偏航力矩产生了明显的迟滞特性,但对俯仰力矩和法向力的迟滞特性影响不大。     

基于相交不稳定性的尾流自消散机翼实验研究

飞机尾流是复杂的流动现象,相关控制的研究常采用简化模型,抓住主要矛盾进行尾流不稳定性的学术探索. 采用结构化矩形机翼模型,通过添加扰流片来模拟襟翼的一种作动方式,引入一对与主翼涡反向的小涡,以期诱发尾涡的瑞利-路德维希相交不稳定性. 改变模型在水槽中的拖曳速度以及机翼攻角,采用粒子图像速度场仪定量研究单主翼尾涡发展特性以及双涡相互作用特性. 研究表明,未添加扰流片时,尾涡环量在45 个翼展内相对于初始环量衰减了10%;而添加了扰流片的实验中,在较好的实验参数组合情况下,主翼尾涡环量较初始环量降低35%~45%. 结果表明添加适当扰流片产生的反向小涡能诱发与主翼尾涡的相交不稳定性,在尾流涡系中引入自消散机制,加速机翼尾涡的消散过程,达到提早消弱尾涡的目的.      

弹性机翼阵风响应数值计算方法

建立了求解弹性机翼阵风响应的计算方法.在计算中,通过采用数值方法求解三维非定常Euler方程来获得气动特性;采用模态叠加的方法考虑弹性影响,实现了流体力学和弹性力学的耦合计算.通过对刚性机翼在攻角突然增大的阵风作用下的响应历程计算和二维NLR7301翼型的极限环振荡计算,对计算方法进行了验证.此后在"1-cos"阵风响应的计算中考虑弹性效应影响,先是只考虑了结构变形的前三个基本模态,弹性机翼气动力响应的计算结果与刚性机翼的响应计算结果有比较大的区别,弹性机翼阵风响应的升力峰值低于刚性机翼,这与文献中的结果是一致的.最后在计算中考虑了高阶弹性模态,计算结果表明:考虑高阶模态后,机翼气动力计算结果的总体变化趋势与只考虑前三个模态时基本一致,但结果中出现了高频的波动,波动的频率与高阶模态本身的频率有关.     

Validation of actuator disc circulation distribution for unsteady virtual blades model

The actuator disc method is an engineering approach to reduce computer resources in computational fluid dynamics (CFD) simulations of helicopter rotors or aeroplane propellers. Implementation of an actuator disc based on rotor circulation distribution allows for approximations to be made while reproducing the blade tip vortices. Radial circulation distributions can be formulated according to the nonuniform Heyson-Katzoff “typical load” in hover. In forward flight, the nonuniform disk models include “azimuthal” sin and cos terms to reproduce the blade cyclic motion. The azimuthal circulation distribution for a forward flight mode corresponds to trimmed conditions for the disk rolling and pitching moments. The amplitude of the cos harmonic is analysed and compared here with presented in references data and CFD simulations results.     

Force production and asymmetric deformation of a flexible flapping wing in forward flight

Insect wings usually are flexible and deform significantly under the combined inertial and aerodynamic load. To study the effect of wing flexibility on both lift and thrust production in forward flight, a two-dimensional numerical simulation is employed to compute the fluid–structure interaction of an elastic wing section translating in an inclined stroke plane while pitching around its leading ledge. The effects of the wing stiffness, mass ratio, stroke plane angle, and flight speed are considered. The results show that the passive pitching due to wing deformation can significantly increase thrust while either maintaining lift at the same level or increasing it simultaneously. Another important finding is that even though the wing structure and actuation kinematics are symmetric, chordwise deformation of the wing shows a larger magnitude during upstroke than during downstroke. The asymmetry is more pronounced when the wing has a low mass ratio so that the fluid-induced deformation is significant. Such an aerodynamic cause may serve as an additional mechanism for the asymmetric deformation pattern observed in real insects.     

Control of flight forces and moments by flapping wings of model bumblebee

The control of flight forces and moments by flapping wings of a model bumblebee is studied using the method of computational fluid dynamics.Hovering flight is taken as the reference flight:Wing kinematic parameters are varied with respect to their values at hovering flight.Moments about(and forces along)x,y,z axes that pass the center of mass are computed.Changing stroke amplitude(or wingbeat frequency)mainly produces a vertical force.Changing mean stroke angle mainly produces a pitch moment.Changing wing angle of attack,when down-and upstrokes have equal change,mainly produces a vertical force,while when down-and upstrokes have opposite changes,mainly produces a horizontal force and a pitch moment.Changing wing rotation timing,when dorsal and ventral rotations have the same timing,mainly produces a vertical force,while when dorsal and ventral rotations have opposite timings,mainly produces a pitch moment and a horizontal force.Changing rotation duration has very small effect on forces and moments.Anti-symmetrically changing stroke amplitude(or wingbeat frequency)of the contralateral wings mainly produces a roll moment.Anti-symmetrically changing angles of attack of the contralateral wings,when down-and upstrokes have equal change,mainly produces a roll moment,while when down-and upstrokes have opposite changes,mainly produces a yaw moment.Anti-symmetrically changing wing rotation timing of the contralateral wings,when dorsal and ventral rotations have the same timing,mainly produces a roll moment and a side force,while when dorsal and ventral rotations have opposite timings,mainly produces a yaw moment.Vertical force and moments about the three axes can be separately controlled by separate kinematic variables.A very fast rotation can be achieved with moderate changes in wing kinematics.     

Aeroelastic response of a coupled rotor/fuselage system in hovering and forward flight

Summary The aeroelastic response analysis of a coupled rotor/fuselage system is approached by iterative solution of the blade aeroelastic response in the non-inertial reference frame fixed on the hub, and the periodic response of the fuselage in the inertial reference frame. A model of the coupled system hinged with the flap and lag hinges, the pitching bearing which may not coincide with the hinges, and the sweeping-blade configuration is established. The moderate-deflection beam theory and the two-dimensional quasi-steady aerodynamic model are employed to model the aeroelastic blade, all the kinetic and inertial factors are taken into account in a unified manner. A five-nodes, 15-DOFs pre-twisted nonuniform beam element is developed for the discretization of the blade, three rigid-body-motion DOFs are introduced for the motion of the hinges and the bearing. The Hamilton's principle is employed to evaluate the equation of motion of the blade. The derived nonlinear ordinary differential equations with time-dependent periodic coefficients are solved by a modified quasi-linearization method, which is developed for the higher DOF periodic system. The resulting periodic forces and moments exerted on the fuselage by all the blades are evaluated every time, when the converged nonlinear periodic response of the blade is obtained under the consideration of the equilibrium of the blades. The fuselage structure is simplified to be a beam structure, the governing equation is established in the inertial reference frame and a two-nodes beam element is used to discretize the flexible fuselage. The periodic response of the fuselage is solved by a simple shooting method. The iteration of the rotor/fuselage response is continued, until the aeroelastic responses of the blade and the fuselage converge simultaneously. Both the hovering and the forward flight states can be considered. The results of a computed numerical example by the developed program are presented to verify in practice the economy of the modeling as well as the reliability and efficiency of the corresponding solving methods. Received 4 May 1998; accepted 11 August 1998     

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.     

Analysis on the nonlinear response of cracked rotor in hover flight

The motion equations for a Jeffcott rotor in hover flight are derived. A periodically sampled peak-to-peak value diagram is used for characterizing and distinguishing different types of nonlinear responses in hovering state. The nonlinear responses become more apparent when the rotor is running above the critical speed in flat flight. There are three ways for rotor responses going to chaos, namely through quasi-periodic, intermittence, or period-3 bifurcation to chaos. The hover flight might suppress some nonlinear responses. However, the position of axis center might obviously deflect, leading to either nonlinear response or peak-to-peak value jump near the fraction frequency of swing critical speed.     

A framework for CFD analysis of helicopter rotors in hover and forward flight

A framework is described and demonstrated for CFD analysis of helicopter rotors in hover and forward flight. Starting from the Navier–Stokes equations, the paper describes the periodic rotor blade motions required to trim the rotor in forward flight (blade flapping, blade lead‐lag and blade pitching) as well as the required mesh deformation. Throughout, the rotor blades are assumed to be rigid and the rotor to be fully articulated with separate hinges for each blade. The employed method allows for rotors with different numbers of blades and with various rotor hub layouts to be analysed. This method is then combined with a novel grid deformation strategy which preserves the quality of multi‐block structured, body‐fitted grids around the blades. The coupling of the CFD method with a rotor trimming approach is also described and implemented. The complete framework is validated for hovering and forward flying rotors and comparisons are made against available experimental data. Finally, suggestions for further development are put forward. For all cases, results were in good agreement with experiments and rapid convergence has been obtained. Comparisons between the present grid deformation method and transfinite interpolation were made highlighting the advantages of the current approach. Copyright © 2006 John Wiley & Sons, Ltd.     

两种椭圆翼型高速气动特性实验研究

新概念旋转机翼飞机的主机翼既能高速旋转作为旋翼,又可锁定作为固定翼,所以只能使用特殊的前后对称翼型。针对主机翼翼型的这一特殊要求,对16％相对厚度,相对弯度分别为0％和3％的两种椭圆翼型的高速气动特性进行了风洞实验研究,试验分别在中国空气动力研究发展中心FL-21风洞和荷兰代尔夫特大学TST-27风洞进行,采用表面测压和尾排型阻测量技术。试验结果的对比分析表明,有弯度椭圆翼型的升力和力矩特性优于无弯度椭圆翼型,而阻力特性和最大升阻比劣于无弯度椭圆翼型。试验结果为旋转机翼飞机主机翼翼型的选取提供了参考。     

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.     

Multiple time scale analysis of aircraft longitudinal dynamics with aerodynamic vectoring

The paper analyzes the longitudinal dynamics of sustained high angle-of-attack flights performed by a small agile Unmanned Aerial Vehicle (UAV) equipped with an aerodynamic vectoring feature. This unique feature is achieved through variable-incidence wing that allows the angle of incidence of the wing with respect to fuselage to be decoupled. The aerodynamics of the UAV is estimated through wind-tunnel experiments. Subsequently, the dynamics of the UAV motion are studied for trim conditions across the complete spectrum of velocity between hover and forward cruise flight. The short-period eigenvalue migration, as the trimmed flight conditions are changed, studied, and the peculiarities observed from this linear analysis are highlighted. Further, the Multiple Time Scales method in conjunction with bifurcation theory is used to obtain the approximate solutions of the multiple degree-of-freedom nonlinear equations of motion. The explicit analytical results obtained are useful to identify the key parameters affecting the dynamics and stability of the aircraft??s longitudinal motion. The conditions leading to limit cycle responses in the vicinity of the stall regime are also highlighted.