基于CFD/CSD的非线性气动弹性分析方法

提出了一种基于计算流体力学/计算结构动力学(CFD/CSD)的非线性气动弹性分析方法,并应用于切尖三角翼的非线性颤振和极限环振荡(LCO)研究。该方法将非线性有限元(FEM)和CFD计算相结合,并辅以高精度的界面插值,能够分析结构和气动非线性共存的气动弹性问题。结构部分以四边形平板壳元为基础,采用更新的拉格朗日(UL)方法分析结构大变形引起的几何非线性问题。气动部分以Navier-Stokes(N-S)方程作为控制方程,采用CFD方法计算跨声速气动力。机翼的非线性颤振计算表明了方法的有效性。最后应用该方法研究了切尖三角翼的LCO现象,其计算精度明显优于已有结果。     

三角翼和双三角翼前缘涡及其卷吸作用的计算

本文提出一种计算模型,用以计算亚声速大迎角下前缘分离机翼的流动和气动特性。通过在涡轴上分布线涡/线汇,这种计算模型包括了对前缘自由涡面、涡核及涡的卷吸作用的模拟。由于它是在非线性离散涡法的基础上发展起来的,因而具有计算过程简捷的特点。对三角翼及双三角翼气动特性的计算表明,计算值与实验值符合得相当好。计算还表明,在计算中计及和不计及涡的卷吸作用能引起计算载荷相当大的变化。     

压敏漆测量三角翼模型气动载荷的超声速风洞实验

采用压敏漆在超声速风洞中测量三角翼模型气动载荷,得到了应用压敏漆测压与天平测力相结合获取的三角翼模型气动载荷对比实验结果。简述了压敏漆原理、标定、试验设置和数字图像数据处理方法。由于尺寸限制,模型没有开设测压孔。对三角翼模型进行了数值计算,对试验和计算得到的结果进行了比较分析。     

基于径向基函数的网格变形及非线性气动弹性时域仿真研究

为开展非线性气动弹性研究,基于非线性结构有限元软件NASTRAN和自主研制的多块结构化计算流体力学(CFD)求解器,开发了一套基于计算流体力学/计算结构动力学(CFD/CSD)耦合求解方法的气动弹性时域仿真程序.该程序采用径向基函数(RBF)交换两套求解器之间的数据并进行网格变形.为提高RBF方法的效率,构造了基于多次插值的空间待插值点精简算法.在多次插值过程中,每次插值的对象为上次插值的误差,并同时限制插值区域,以此实现了空间待插值网格数的精简.数个网格变形的算例表明该方法可支持大变形运动,并且具有较高的计算效率.采用此程序开展了AGARD 445.6机翼颤振计算、大展弦比机翼的静气动弹性计算与切尖三角翼极限环振荡(LCO)现象的动气动弹性仿真,结果揭示了当机翼展弦比较大或者响应幅值较大时,结构非线性对于气动弹性有显著影响.     

非定常自由来流对静态三角翼气动特性的影响

在南航非定常风洞中,运用动态测力、测压和流动显示技术,详细研究了非定常自由来流对静态三角翼气动特性的影响和三角翼背风面空间流场结构的变化.研究结果表明,在不同攻角下,随来流速度的脉动三角翼气动特性产生的变化不同.非定常自由来流对静态三角翼气动特性产生的影响,主要是由于来流风速的变化对三角翼上翼面的流动结构产生的影响所造成,特别是在静态失速攻角前后,这种影响最为明显,它使原先翼面上的破碎涡流变成了集中涡流.     

三角翼DBD等离子体流动控制研究进展

现代战机采用较多的三角翼，在大迎角绕流时存在前缘涡破裂等气动问题。作为新型主动流动控制技术，等离子体激励频带宽、响应快、结构简单、便于闭环控制，在解决三角翼气动问题上具有潜力。回顾了介质阻挡放电（DBD）等离子体气动激励的基本原理，及其用于三角翼前缘涡控制的研究进展。从来流条件、几何构型、激励参数等方面分析了DBD等离子体激励对流动控制效果的影响规律；结合不同激励频率下流场演化特性，分析了流动控制机理。最后，从理论研究和工程应用的角度，对三角翼前缘涡控制的发展进行总结展望。     

基于仿真分析的压铆力预测模型

铆接是飞机结构机械连接的主要方式,压铆力作为影响铆接质量的重要因素,传统的计算方法通常是建立在体积不变的假设之上,且不考虑钉杆材料被压入钉孔部分的体积,导致模型误差较大。为此,依据仿真分析结果和铆钉材料的流动趋势,引入体积缩减系数来描述钉杆被压入钉孔部分的体积,建立压铆力预测模型,并与已有的试验数据进行对比,结果表明计算值与试验值的一致性较好,该模型可以用来预测铆接过程中压铆力的大小。     

含鸟撞变形叶片的压气机气动性能数值模拟

为研究受到鸟撞前后压气机气动性能的变化，提出了鸟撞叶片结构-气动分析几何模型的转化流程，基于NASA Rotor37 转子得到了鸟撞变形叶片几何模型，分别建立了鸟撞前后的全通道气动性能CFD计算分析模型，在设计转速下开展了全3维黏性 流场数值模拟，并与Rotor37转子部件气动性能试验数据进行了对比分析。结果表明：模拟结果与试验结果非常接近，证明了该数 值模拟方法有效；鸟撞后叶片变形区域攻角增大导致的局部气流分离及并发的气流低速流动的耦合是转子气动性能恶化与转子 进入失稳工况的主要原因，含鸟撞变形叶片的转子压比、效率等气动性能参数明显降低，稳定工作边界明显缩小。所发展的气动 性能数值模拟方法与流程可有效地预估含变形叶片的压气机稳态气动性能。     

考虑机翼几何非线性的气动弹性建模与分析

大展弦比柔性机翼结构重量轻、气动效率高,广泛应用于高空长航时无人机(UAVs)。飞行过程中,这类机翼在气动力作用下发生大变形,线性结构模型不再适用,需要建立考虑几何大变形的结构模型。采用牛顿力学方法推导了考虑结构几何非线性的机翼结构动力学模型,该方法推导过程简洁、物理意义明确,可以与Hodges基于哈密顿原理的推导方法相互补充,相互验证。为了能够更准确地求解大展弦比柔性机翼的非定常气动力,建立了能够考虑机翼三维效应且适用于机翼空间大变形的非定常气动力模型。基于建立的非线性结构模型和非定常气动力模型,采用松耦合方法建立了非线性气动弹性模型,并通过算例验证了气弹模型的准确性。研究结果表明,大展弦比柔性机翼颤振速度对来流迎角和机翼的展长均较为敏感;当来流速度大于颤振速度时,由于几何非线性,机翼振动并未发散而是形成稳定的极限环振荡(LCO);随着来流速度进一步增加,机翼再次穿过临界稳定点,由不稳定系统变为稳定系统,直到随着速度的增加系统再次达到临界稳定状态。     

纳秒脉冲等离子体激励控制小后掠三角翼低速绕流试验

为探索纳秒脉冲介质阻挡放电(NS DBD)对小后掠尖前缘三角翼的流动控制效果和作用机理,进行NS DBD用于改善其气动特性的测力试验和流动显示试验。当来流速度分别为30m/s和45m/s时,测力试验结果表明位于机翼前缘的NS DBD能很好地改善三角翼大迎角气动特性,其中来流速度为45m/s时最大升力系数提高了18.3%;研究了脉冲激励频率对流动控制效果的影响规律,最佳的无量纲激励频率F+≈1~2。在来流速度为20m/s时,采用粒子图像测速仪(PIV)研究了不同迎角下激励前后机翼背风面流场,表明NS DBD可改善上翼面旋涡结构,使分离涡附体并得到加强。基于试验结果,认为NS DBD进行三角翼前缘涡控制的机理是激励诱导分离剪切层周期性产生附体的分离涡,从而维持了上翼面大迎角时的涡升力。     

大展弦比复合材料机翼失速颤振分析

研究了大展弦比复合材料机翼在较大迎角状态下的失速颤振特性,探讨了结构几何非线性和由复合材料剪裁产生的刚度耦合效果对机翼失速颤振特性的影响.首先,将复合材料机翼建模为转角和位移均可为有限值的非线性薄壁单闭室截面Euler梁,并在综合考虑结构几何非线性、气动非线性和材料各向异性对机翼运动状态的影响的基础上,建立机翼的运动微分方程.然后,使用小扰动分析的方法得到机翼在平衡位置附近的振动方程,采用ONERA半经验的非定常失速气动力模型,获得机翼在平衡位置附近的非线性失速颤振分析方程.最后,利用谐波平衡法求解并判定机翼颤振稳定性.通过算例,首先验证了算法的正确性,然后研究了几何非线性对失速颤振的影响,并讨论不同的复合材料铺层方式导致机翼失速特性的改变.     

大后掠翼前缘涡对其颤振特性的影响

大迎角三角翼的前缘涡不仅可以改善其气动力特性,也会显著影响机翼的气动弹性特性.运用基于Euler方程的非定常气动力降阶模型(ROM)方法,耦合结构运动方程,在状态空间内建立了气动弹性分析模型,研究了70°削尖三角翼的大迎角颤振特性.研究结果显示前缘涡对该机翼颤振特性的影响不可忽略.颤振速度随迎角的增加而大幅降低,迎角α=20°时的颤振速度比α=0°时降低了22%.发现了颤振特性随迎角变化时出现的不连续现象,并揭示了该现象是由于系统颤振分支随着静态迎角的增加发生转移所致.     

Limit cycle oscillation control of wing with static output feedback control method

This paper presents a static output feedback controller (SOFC) for aeroelastic control of a cantilevered rectangular wing in low subsonic flow. For this purpose, an optimal formulation of this control method is developed, and a solution method is proposed for the related matrix equation. This optimal solution is obtained by solving combined Lyapunov and Riccati equations. At first these equations are transformed into a set of nonlinear algebraic equations and then are solved with iterative Newton–Raphson?s method. This solution method is applicable to the full state feedback case. The controller is designed to extend flutter boundary and suppress limit cycle oscillation (LCO) of a low aspect ratio rectangular nonlinear structural wing. This structural nonlinearity is given by Von Karman plate theory. Both full and reduced order aerodynamic models are examined based on the modified vortex lattice theory. Results show combination of SOFC with reduced order aerodynamic model would be an effective choice for aeroelastic stabilization, and this controller has a very comparable result with linear quadratic regulator (LQR).     

Control law design for transonic aeroservoelasticity

Existing computational transonic aeroservoelastic researches focus on directly coupling the structural dynamic equations, CFD solver and servo system in time domain, study the effect of the given feedback control laws on the responses of the aeroelastic system. These works have not involved the design of the flutter active control law. The non-linearity of transonic flow brings great difficulties to aeroservoelastic analysis and design. Recent research of the unsteady aerodynamic reduced order models (ROM) based on CFD provides a challenging approach for transonic aeroservoelastic analysis and design. Coupling the structural state equations with the aerodynamic state equations of the wing and the control surface based on the ROM, we construct a transonic aeroservoelastic model in state-space. Then the sub-optimal control method based on output feedback is used to design the flutter suppressing law. The study first demonstrates the open loop of the Benchmark Active Controls Technology (BACT) wing. The computational results of the CFD direct simulation method and the ROM analysis method are both agree well with the experimental data. Then both the closed loop time responses and the flutter results by ROM technique are compared with those of numerical aeroservoelastic simulation based on Euler codes to validate the correctness of the design method of the control law and aeroservoelastic analysis method. An increase of up to 20% of the speed index can be achieved by the control law designed by sub-optimal control method for this model.     

大后掠翼前缘涡对其颤振特性的影响

大迎角三角翼的前缘涡不仅可以改善其气动力特性,也会显著影响机翼的气动弹性特性。运用基于Euler方程的非定常气动力降阶模型（ROM）方法,耦合结构运动方程,在状态空间内建立了气动弹性分析模型,研究了70°削尖三角翼的大迎角颤振特性。研究结果显示前缘涡对该机翼颤振特性的影响不可忽略。颤振速度随迎角的增加而大幅降低,迎角α=20°时的颤振速度比α=0°时降低了22％。发现了颤振特性随迎角变化时出现的不连续现象,并揭示了该现象是由于系统颤振分支随着静态迎角的增加发生转移所致。     

双三角翼布局歼击教练飞机的概念优化设计

双三角翼气动布局比三角翼飞机具有更好的大攻角空气动力特性。引入了评估教练机训练效能的作战分析法 ,研究了双三角机翼布局飞机空气动力特性的工程计算途径以及飞机性能指标的确定方法。以训练效能作为目标函数并选取机翼平面形状的几何参数为设计变量 ,采用多变量数值寻优方法 ,在战术技术指标及相关几何约束条件下 ,对某高级教练机的双三角机翼气动布局方案进行了优化选择。算例表明最优方案不仅比原准方案具有更高的训练效能 ,还改善了结构的受力情况 ,与工程实践吻合。     

大尖削机翼非定常跨音速绕流的数值解法(英文)

本文给出一种计算飞机机翼上定常和非定常跨音速气动力的数值解法,使用了一种特殊设计的坐标变换,是用时间精确交替方向隐式(ADI)有限差分算法来求解非定常跨音速修正三元小扰动位势方程。给出了F5战斗机机翼的数值结果并与XTRAN3S,ATRAN3S及试验结果进行了比较,表明本方法是有效的和经济的。