SMA驱动变体机翼后缘结构力学分析

对基于形状记忆合金（Shape memory alloy,SMA）丝驱动器的多关节变体机翼后缘模型,在指定飞行工况下结构的整体极限偏转角度是最重要的设计指标。本文根据SMA驱动器的力学本构模型,建立了在控制回复状态下的SMA-弹簧模型。基于N S方程编制了2D流场求解器,并使用网格重构技术建立了后缘偏转时的多套流场网格,通过将后缘结构所受的气动载荷力矩对应施加于后缘结构的各个关节,在ANSYS静力分析模块中对后缘的偏转过程进行了数值仿真,考察了变体后缘的整体极限偏转角度。为验证仿真结果,对变体后缘结构进行了模拟加载实验。实验结果表明,本文使用的数值计算方法能够较为准确地预测变体后缘的整体极限偏转角度。     

一种新PSO混合算法在直升机配平中的应用

直升机配平计算是动力学分析的基础,其实质是求解高维复杂的非线性方程组。针对经典算法与智能算法的特点与不足,提出了一种求解非线性方程组的新粒子群方法。在粒子群(Particle swarm optimization, PSO)算法的基础上,根据模拟退火(Simulated annealing, SA)思想,引入了嵌入式LM (Levenberg-marquardt)优化 算子。该方法充分发挥了3种算法的优势,克服了LM算法初值敏感性,PSO算法易陷入局部极值等问题。通过UH-60A直升机实例配平计算,验证了本文算例模型的准确性。在此基础 上,针对某一前飞状态下的配平算例,在收敛可靠性和计算效率上通过与其他算法进行对比,表明该算法具有可靠的收敛性和较高的计算效率,进一步验证了该算法在配平问题上的可信度与实用性,为直升机飞行动力学问题的处理提供了一种新的有效方法。     

风力机翼型阵风响应计算与分析

基于网格速度方法,对阵风条件下风力机翼型运动的气动特性进行数值计算。中心格式有限体积法、双时间推进法求解预处理后的非定常N-S方程。计算结果表明:在弦向阵风作用下,翼型的运动能有效降低响应波动峰值,即弦向阵风对运动影响较小;而在法向阵风作用下,气动响应波动很大。计算翼型在不同参数阵风作用下的响应过程,阵风的频率变化对响应影响不大。弦向阵风的幅值增大,叶片运动对阵风响应峰值减缓效果增加;法向阵风的幅值增大1倍,响应峰值也增加1倍。     

高阶精度DG/FV混合方法在二维粘性流动模拟中的推广

DG/FV混合方法因其具有紧致性、易于推广至高阶及相比同阶DGM计算量、存储量小等优点,已成功应用于一维/二维标量方程和Euler方程的求解。在此基础上,将该方法推广于二维三角形/矩形混合网格上的NavierStokes方程数值模拟,将格式形式精度提高至4~5阶。物理量的空间重构及离散使用DG/FV混合重构方法;无粘通量计算采用Roe格式;粘性通量计算采用BR2格式;时间方向离散采用高阶显式R-K方法或隐式方法。利用该方法计算了有解析解的Couette流动问题以验证几种格式的数值精度阶,并计算了层流平板流动和定常、非定常圆柱绕流问题等经典算例。计算结果表明DG/FV混合方法达到了设计的精度阶,在较粗的网格上亦能得到高精度的计算结果;定性分析和数值结果表明相比同阶DG方法单步计算量减少约40%。     

抛物化 Navier-Stokes 方程(PNS)和高氏 PNS 理论及其应用的评述

抛物化 NS 方程得到广泛应用,已经成为工业标准气动计算的基础。现有的八种抛物化 NS 方程有不同的名称,方程中粘性项的形式略有不同,其中的 PNS 和薄层（TL）NS 方程应用最多。但是这些方程都具有类似的数学性质,例如,当流向方向上马赫数大于1时,他们都是抛物型方程,可以采用空间推进算法（SMA）进行求解。与采用时间推进算法求解的 NS 方程或雷诺平均（RA）NS 方程相比,PNS-SMA 计算降低了空间的维数,节省了大量的存储空间和 CPU 计算时间。PNS-SMA 算法也获得了巨大的进展。但是,早期 PNS 研究在理论上是相当模糊的,高智在1990年提出的粘性／无粘干扰剪切流理论（ISF）弥补了这一不足。ISF 理论概括了 PNS 方程所能描述的基本流动,提出了其流动的运动规律及数学定义式,所导出的 ISF 方程组也属于 PNS 方程的一种。为了不增加新的名称,我们将 ISF 方程组也称为高氏 PNS 理论和方程组。这一理论在 NS 方程和RANS 方程的计算中均有重要的应用。例如,计算最优坐标系的选择以减少伪扩散,网格尺度选择及局部网格加密设计以捕捉高超声速流动中物体表面热流等的急剧变化,壁面压力边界条件的选择以及由高 PNS 导出的壁面判据来进行 NS 和 RANS 近壁数值解可信度评估。本文评述了一些初步的应用,进一步的应用和综合 PNS-SMA,RANS-SMA 以及 PSE-SMA 计算值得深入研究,这里 PSE指抛物化稳定性方程。     

DLR-F6/FX2B翼身组合体构型高阶精度数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS)和剪切应力输运(SST)两方程湍流模型,开展了DLR-F6和DLR-F6_FX2B 2种翼身组合体构型的高阶精度数值模拟,计算外形来自AIAA第三届阻力预测研讨会。主要目的是确认WCNS模拟跨声速典型运输机构型和预测局部构型变化引起的气动特性变化量的能力。在固定升力系数条件下,采用粗、中、细3套网格开展了网格收敛性研究,从气动力系数、压力系数分布、表面流态等方面研究了网格规模对DLR-F6和DLR-F6_FX2B翼身组合体数值模拟结果的影响;采用中等网格开展了来流迎角对2种翼身组合体气动特性的影响研究。通过与National Transonic Facility(NTF)的试验结果和CFL3D的计算结果对比,表明采用高阶精度计算方法得到了网格收敛的数值模拟结果,较好地模拟了DLR-F6翼身组合体局部修型引起的微小气动特性变化和翼身结合部流动特性的差异。     

Aerodynamic design optimization of nacelle/pylon position on an aircraft

The arbitrary space-shape free form deformation (FFD) method developed in this paper is based on non-uniform rational B-splines (NURBS) basis function and used for the integral parameterization of nacelle-pylon geometry. The multi-block structured grid deformation technique is established by Delaunay graph mapping method. The optimization objects of aerodynamic characteristics are evaluated by solving NavierStokes equations on the basis of multi-block structured grid. The advanced particle swarm optimization (PSO) is utilized as search algorithm, which com-bines the Kriging model as surrogate model during optimization. The optimization system is used for optimizing the nacelle location of DLR-F6 wing-body-pylon-nacelle. The results indicate that the aerodynamic interference between the parts is significantly reduced. The optimization design system established in this paper has extensive applications and engineering value.     

人员离机轨迹数值模拟及快速预测研究

离机轨迹研究在分析外挂物分离、导弹投放以及人员空中脱离过程等方面有着重要的参考意义。运用动态嵌套网格技术,通过CFD方法来模拟人员空中离机后的运动轨迹。提出一套基于本征正交分解(POD)方法的人员离机轨迹快速预测算法。将CFD数值模拟得出的离机轨迹作为样本,以飞行马赫数和攻角双参数变化为例,详细介绍POD方法预测离机轨迹的实现步骤。结果表明:POD方法相对于CFD方法可以大大减少计算时间,并且POD预测结果与CFD方法得到的结果吻合较好,为离机轨迹预测在工程应用上提供了一种快速简便的方法。     

一种基于共轭方程法求解黏性反问题的简化方法

 对于给定压力分布的黏性气动反问题,考虑到壁面微小扰动造成的压力变化主要由势流作用引起,因此可以简化用以获得目标函数对设计变量敏感性导数的共轭方程。将黏性流场插值到粗网格中作为彻体力模型方程的解,则其相应的共轭方程将以简单的源项取代原方程中复杂的黏性项。由于在粗网格中求解,网格数减少,同时收敛速度加快,简化的共轭方程计算时间可以减少到黏性方程的十分之一。典型的算例结果表明,对于附着的边界层流动简化方法计算得到的敏感性导数具有较高的精度,能够有效完成反问题设计且减少总的计算耗时。     

气动声学计算中一种声扰动方程的改进

研究了气动声学计算中的一种改进声扰动方程(IAPE)。采用三种不同的声传播模型,数值分析了声在均匀流以及剪切流中的传播问题。与线化欧拉方程(LEE)相比,除了个别位置预测的压力峰值较小以外,声扰动方程(APE)几乎能够预测到类似的物理特性;另外,APE方程具有无需求解密度方程,计算成本相对较少等优点。为弥补APE方程预测压力峰值较小的缺点,通过分析APE方程,这类方程中需要加入剪切流与声波相互作用的源项。为此,对APE方程进行改进,加入了一种声源项,计算结果表明,改进的APE方程(IAPE)能够和线化欧拉方程的计算结果保持一致。