悬停飞昆虫拍动翼低雷诺数高升力气动机理研究

采用数值方法对果蝇翼悬停飞行拍动翼问题进行了模拟,介绍了前人生物观察和动态比例模型得来的简化拍动翼运动规律。果蝇翼模型平面形状采用Dickinson动态比例试验采用的形状,计算的无量纲参数根据Weis-Fogh,和Vogel果蝇捆绑飞行的数据得到。分别数值模拟了三种运动模态:超前模式(advanced mode),对称模式(symmetrical mode)和滞后模式(delayed mode)。数值模拟的结果同Dickinson的试验结果和孙茂的数值模拟结果进行了比较,吻合很好,表明本文的计算结果是合理的。根据数值模拟结果,结合气动力系数和流场结构进行分析,研究了果蝇悬停飞行获得高升力的流动机理。     

压电致动器扑翼结构动力学仿真

采用传统扑动机构的微型扑翼,气动效率低、扑动能耗大。采用压电致动器的微型扑翼通过压电材料进行致动使机翼产生上下扑动,有效地将扑动机构和机翼两个主要系统进行集成,不仅节省重量,同时它具有任务变形自适应能力强、气动效率高和扑动能耗小的特点。通过PCL语言建立采用压电致动器的扑翼有限元模型,结合同尺寸扑翼气动力试验数据,进行采用压电致动器扑翼结构仿真。利用仿真结果研制采用压电致动器扑翼原理样机,研究表明,采用压电致动器后扑翼扑动频率得到明显提高,压电片可有效控制机翼的弯扭变形,有助于提高扑翼的气动效率。     

基于CFD/CSD方法的蜻蜓柔性翼气动特性分析

给出了一种基于计算流体力学/计算结构力学（CFD/CSD）的双向流固耦合方法.通过交替数字二叉树（ADT）搜索技术识别流固网格之间的宿主-受体关系.采用局部插值算法完成两套网格系统之间的数据交换,并使用Delaunay图映射方法来完成气动网格的移动.将自编的非线性结构有限元程序、接口程序与南京航空航天大学（NUAA）微型飞行器中心的流体计算程序3D2MUFS相连接,应用于蜻蜓柔性翼拍动飞行的气动计算中.计算结果表明:柔性变形使得蜻蜓翼的时均举力系数从0.31提高到0.53,时均推力系数从0.07提高到0.13,证实了柔性变形能改善扑翼的气动性能.      

扑翼机翼气动力和惯性力对翼杆结构变形研究

通过有限元方法探究了机翼上气动力和惯性力对结构变形的影响.研究表明扑翼机翼产生的气动力和惯性力大小在同一量级.机翼惯性力对机翼结构变形具有不可忽视的影响.当机翼处于水平位置时机翼气动力对结构变形贡献最大,同时合力产生的最大变形也出现在机翼水平位置.仿生扑翼机翼结构变形研究为高升力、低能耗扑翼机翼的设计提供了理论指导.      

三维机翼的静气动弹性特性数值模拟研究

发展了一种计算流体动力学（CFD）和计算结构动力学（CSD）的耦合计算方法,对三维机翼的静气动弹性进行了数值模拟研究。采用三维欧拉方程为控制方程基于直角网格计算气动力,并耦合结构静平衡方程进行静气动弹性数值模拟,设计了CFD/CSD耦合计算的数据交换的方式。以M6机翼作为算例,进行了机翼静气动弹性的数值模拟,计算结果表明所发展的三维机翼静气动弹性数值模拟方法是合理可行的,并可作为机翼结构优化设计和考虑结构弹性变形影响的气动外形优化设计的基础。     

过渡状态的倾转旋翼气动特性计算分析(英文)

针对特殊的旋翼倾转运动,建立了一个过渡状态的旋翼非定常气动力数值计算方法.为正确模拟旋翼倾转运动使桨叶受到的附加惯性力及哥氏力作用,重新推导了旋翼倾转时的桨叶气动力模型和挥舞运动方程;为了适合于旋翼倾转时的入流和气动力计算,入流模型中考虑了倾转运动引起的旋翼尾迹弯曲影响.应用建立的方法,首先进行了旋翼配平计算,以验证计算模型,并给出了倾转旋翼的操纵量.然后,着重计算了旋翼在倾转过渡时的拉力.俯仰和滚转力矩随倾转角的变化,分析了倾转飞行时的前飞速度.倾转时间等对旋翼气动力的影响,得出了一些新的结论.     

一种非对称折叠扑翼的风洞试验与数值模拟

为了解鸟类翅膀折叠运动的作用,对一个专利中的折叠扑翼机构进行了数学建模,外翼的折叠运动由非定常过程中的气动力、弹性恢复力矩和惯性力决定.利用风洞试验和数值模拟两种方式对该折叠扑翼模型进行了研究.数值模拟和风洞试验结果表明:相对于非折叠扑翼,折叠翼能够有效提高平均升力;在一定范围内增加内翼扑动频率可以提高折叠翼平均升力系...     

扑翼飞行器动力系统建模方法

为快速评估扑翼飞行器的航时,便于针对不同扑动翼进行动力系统设计与优化,逐步减少实物验证与试飞,加快扑翼飞行器的研制,基于实验数据参数辨识的方法建立了包含直流无刷电机、电调（ESC）、锂电池和扑动机构等扑翼飞行器动力系统组件的动态模型,其中电机模型相对误差小于10%,锂电池动态模型相对误差小于6%;提出了一种基于风洞试验气动数据和功率数据的扑动轴瞬时气动载荷半经验高精度建模方法,解决了气动载荷测量较为困难的问题,模型确定系数大于0.89;集成以上模型后的扑翼飞行器仿真系统还包含扑动翼周期平均气动模型、平尾气动模型和纵向控制模型,确保仿真在动态配平状态下进行,可进行全任务剖面航时仿真,航时仿真与实际试飞结果相比误差小于3%。集成的扑翼飞行器仿真系统采用模块化建模思想,各模型参数独立可调,能进一步应用于扑翼飞行器多学科优化等研究。     

基于CFD方法对微型扑翼翼型设计的研究

基于CFD方法从扑翼设计的角度出发对扑翼的外形进行深入的研究.使用基于广义无限插值理论的代数网格生成方法开发适合扑翼运动的三维结构N-S动网格,并应用有限体积法结合双时间推进技术求解预处理后的低速三维非定常N-S方程,完成了对低速扑翼绕流的数值模拟.通过与文献算例的比较验证了本文方法的可靠性和正确性.在此基础上,对文献[8]中介绍的一种混合翼型的相关参数变化对扑翼气动特性的影响开展了深入的研究.     

微型扑翼飞行器非定常运动对平尾的影响

以西北工业大学自行研制的微型扑翼飞行器ASN211为研究对象,利用其简化的二维扑翼及平尾串列翼模型进行了非定常数值模拟,分析了扑翼俯仰运动及沉浮运动对平尾气动性能的影响。在数值模拟模块中,模型的俯仰运动及沉浮运动由动网格技术实现。通过计算流体力学(CFD)软件Fluent对此非定常流场进行数值计算,重点研究了扑翼非定常运动尾流对平尾气动效率的影响。定常状态与非定常时均条件下平尾升力曲线的对比分析表明,扑翼的非定常运动能够增大平尾的失速迎角及最大升力系数,因而使平尾的失速特性得到改善。     

基于直接力控制的阵风响应及阵风减缓研究

考察某无人机模型的阵风气动响应特性和基于直接力操纵的阵风减缓效果.在求解非定常Euler方程时引入“网格速度”方法模拟阵风边界条件,通过动态嵌套网格技术实现舵面运动.首先对NACA0006翼型迎角阶跃型阵风的气动力响应进行了验证计算,计算结果与理论结果、文献计算结果吻合良好.在此基础上对某无人机模型在迎角阶跃型、One-minus-cosine型阵风作用下的气动力响应过程进行了数值模拟,并分析评估了不同舵面运动方式对飞机气动性能的影响.最终研究比较了不同舵面运动方式下阵风气动响应的减缓效果.结果表明:通过合理设计舵面运动,可以有效达到抑制阵风引起的非定常气动干扰的目的.     

微型扑翼飞行器扑翼/尾翼气动干扰的数值研究

以二维刚性约束条件下的微型扑翼飞行器模型为研究对象,在动网格技术基础上,应用非定常数值分析手段对比分析了单翼/纵列翼布局的气动性能,深入研究了纵列翼缩减频率、扑翼—尾翼无量纲水平间距、来流攻角对其气动性能的影响.结果表明:①纵列翼尾翼对扑翼产生正效应干扰,相对于单翼布局,扑翼—尾翼无量纲水平间距为0.5倍翼型弦长时的纵列翼布局的推力系数和推进效率分别增加28.7%和5.7%;②缩减频率是影响推力的关键参数,随着缩减频率的增加,脱落涡的强度增加,推力系数增大.对于单翼、纵列翼两种布局模式,当缩减频率在1.0附近时推进效率达到最优;③对于纵列翼布局,在扑翼—尾翼无量纲水平间距为1.1倍翼型弦长时推进效率达到峰值;④在0°~20°来流攻角变化范围内,随着来流攻角的增加,升力系数增加,推力系数减小,当来流攻角大于9°时,两种布局的推力均为负值.      

微型扑翼飞行器的气动建模分析与试验

用计算流体力学的数值模拟方法研究了微扑翼飞行器的扑翼飞行的非定常空气动力学问题。在对昆虫扑翼飞行运动的仿生模拟基础上 ,对实际可飞的微扑翼飞行器的扑翼运动建立了三维翼型的运动学与空气动力学模型。利用任意拉格朗日欧拉 ( ALE)有限元方法求解出 N-S方程的数值解 ,证明简单扑翼布局所提供的升力足以克服微扑翼飞行器本身的重力使其飞行。在此基础上 ,分别计算并分析了拍动幅值、俯仰幅度以及扑翼频率等各种扑翼参数对升力的影响。最后 ,探索性的扑翼风洞试验与飞行试验结果在一定程度上验证了文中计算方法的可行性      

基于空间RURS机构的三维仿生扑翼机构设计与分析

为实现仿昆虫翼尖的空间“8”字型运动轨迹，设计了一种基于空间revolute-universal-revolute-spherical(RURS)四杆机构的扑翼机构，通过单自由度驱动即可输出三维的空间“8”字轨迹。运用Denavit-Hartenberg参数法建立了空间四杆机构的运动学模型，利用遗传算法对机构进行了优化，得到了利于扑翼飞行的机构参数。基于该空间四杆机构的优化结果，建立了一种微型的扑翼机构虚拟样机，通过ADAMS仿真得到其输出运动并验证了运动学理论计算的正确性。所设计的扑翼机构扑动幅度达到149.8°，扭转角度达到29.9°，且“8”字型扑动规律与昆虫翅膀的运动更为相近。扑翼机构的最大尺寸不超过5.8cm，仿真发现的时间非对称扑动对气动性能有一定提升，对于微型化、轻质化、高效化扑翼飞行器的研究具有重要的参考价值。      

扑翼升力特性的非定常涡格法计算研究

针对西北工业大学航空学院微型飞行器研究室研制的PY-1型扑翼机,建立了扑翼的运动模型,并采用非定常涡格法模拟了该扑翼的升力特性,分析了迎角、飞行速度对升力特性的影响。计算结果与实际飞行试验吻合,验证了非定常涡格法模拟扑翼非定常气动特性的正确性和有效性。研究结果可为扑翼机的气动设计提供参考。     

Correlation analysis of combined and separated effects of wing deformation and support system in the CAE-AVM study

The Chinese Aeronautical Establishment (CAE) Aerodynamic Validation Model (AVM) is a dual-purpose test geometry dedicated to verify the aerodynamic performance of a conceptual intercontinental jet aircraft and to provide a dataset for CFD software validation. To this end, a scaled model of the AVM was tested in the High-Speed Tunnel (HST) of the German-Dutch Wind-tunnels (DNW) with special test consideration and instrumentation. For complementary analysis of experimental results, specific CAE-AVM geometries are analyzed using a CAE in-house CFD code. The specific geometries consist of a baseline aircraft, an aircraft with a deformed wing shape, and an aircraft with both a deformed wing shape and a representation of the model support system used in the wind tunnel. Detailed analysis of numerical and experimental results is presented; both the combined and individual attributions of wing deformation and support system interference on wing pressure distributions and longitudinal aerodynamic characteristics are summarized.     

Nonlinear dynamics of a flapping rotary wing: Modeling and optimal wing kinematic analysis

The analysis of the passive rotation feature of a micro Flapping Rotary Wing(FRW)applicable for Micro Air Vehicle(MAV) design is presented in this paper. The dynamics of the wing and its influence on aerodynamic performance of FRW is studied at low Reynolds number(~10~3).The FRW is modeled as a simplified system of three rigid bodies: a rotary base with two flapping wings. The multibody dynamic theory is employed to derive the motion equations for FRW. A quasi-steady aerodynamic model is utilized for the calculation of the aerodynamic forces and moments. The dynamic motion process and the effects of the kinematics of wings on the dynamic rotational equilibrium of FWR and the aerodynamic performances are studied. The results show that the passive rotation motion of the wings is a continuous dynamic process which converges into an equilibrium rotary velocity due to the interaction between aerodynamic thrust, drag force and wing inertia. This causes a unique dynamic time-lag phenomena of lift generation for FRW, unlike the normal flapping wing flight vehicle driven by its own motor to actively rotate its wings. The analysis also shows that in order to acquire a high positive lift generation with high power efficiency and small dynamic time-lag, a relative high mid-up stroke angle within 7–15° and low mid-down stroke angle within -40° to -35° are necessary. The results provide a quantified guidance for design option of FRW together with the optimal kinematics of motion according to flight performance requirement.     

蝉翅翼的运动研究与三维数值分析

以蝉为研究对象,首先通过分析蝉翅翼在高速摄像机下的运动影像获得了其运动方程,然后建立了蝉翅翼的三维模型,并研究了基于Fluent软件的蝉翅翼周围流场数值模拟方法,最后模拟分析了三维蝉翅翼模型在不同运动参数(拍动频率、拍动幅值和最大扭转角度)下的平均升力系数和平均阻力系数.研究结果表明最大扭转角度对蝉翅翼气动力影响较大,因此在仿蝉扑翼飞行器的设计与操控中必须考虑扭转运动.      

带后缘小翼的旋翼振动载荷计算

建立了考虑弹性桨叶、刚性小翼的旋翼气动弹性分析模型和旋翼载荷计算方法.以广义质量和广义力的形式描述小翼惯性力和气动力对系统的影响,以非定常/动态失速模型计算剖面气动力,结合基于实验数据修正的组合气动模型计算带小翼部分的剖面气动力,集成大变形桨叶模型考虑弹性变形的非线性,以力积分法计算桨叶剖面振动载荷.通过计算分析与实验结果相比较,验证了建立的气动弹性模型和载荷计算方法.结果表明:建立的桨叶结构模型精度很高,气弹模型能够准确预测旋翼的振动载荷,挥舞弯矩平均误差控制在9.1%,使用修正的小翼气动模型能有效提高小翼运动时桨叶振动载荷的计算精度.      

基于CFD/CSD耦合方法的旋翼气动弹性载荷计算分析

为提高直升机前飞状态下旋翼非定常气动弹性载荷的预估精度,在旋翼气动弹性综合分析方法中引入旋翼CFD模块,建立了一套基于CFD/CSD松耦合分析的计算方法和程序。为高效解决流固耦合方法中由于桨叶挥舞、扭转等弹性变形带来的旋翼贴体网格变形问题,采用基于代数变换方法的网格变形技术,桨叶运动变形量和旋翼气动力信息通过流固交接面传递。旋翼流场分析方法中,主控方程采用耦合S-A湍流模型的Navier-Stokes方程,围绕旋翼流场的网格采用结构嵌套网格方法生成,无黏通量计算采用Roe格式,时间推进采用双时间法。旋翼结构分析中,考虑旋翼配平,基于Hamilton变分原理和20自由度Timoshenko梁模型求解弹性旋翼非线性运动方程。分别对CSD和CFD方法进行验证,在此基础上,计算了SA349/2旋翼桨叶在前飞状态下的非定常气动力、挥舞弯矩、摆振弯矩和扭转力矩,并与飞行测试数据进行了对比。计算表明:CFD/CSD耦合方法可以显著提高旋翼非定常气动弹性载荷的分析精度,精确捕捉桨叶表面压强峰值、激波位置等,表明本文发展的旋翼CFD/CSD耦合方法可以有效地运用到旋翼气动弹性载荷的预测分析中。