基于自适应非结构嵌套网格的旋翼流场模拟

针对直升机旋翼CFD仿真的复杂性,提出了改进的适合于格心格式求解器的非结构嵌套网格算法。采用自适应网格技术在旋翼流场仿真的整个过程中进行网格的自适应加密和疏化操作,以更好地捕捉桨尖涡等流动细节。对于频繁的自适应过程中产生的大量重复点和无用点,采用了高效的交替数字树算法(Alternating digital tree,ADT)和标记-删除-移动算法(Mark,delete,move,MDM)进行删除,节约了不必要的存储。针对格心格式的求解器,采用了基于梯度的网格间插值方式,简化了网格间数值传递的复杂性,同时不降低求解器的精度。对CaradonnaTung旋翼悬停算例和HLISHAPE 7A旋翼悬停算例进行了模拟验证,计算值与实验值吻合,表明本文建立的方法具有良好的鲁棒性和有效性。最后,与未采用自适应时求解器对桨尖涡的捕捉效果进行了对比,结果表明本文所采用的方法可以明显地提高求解器对桨尖涡的捕捉。     

基于弹簧系统网格变形方法的旋翼气弹耦合分析

在旋翼气动弹性耦合（CFD/CSD）分析中引入弹簧系统网格变形方法,建立了一套适合于旋翼气动载荷分析的CFD/CSD耦合方法。为了解决CFD/CSD耦合中关键的网格变形问题,旋翼桨叶贴体网格变形采用基于“ball-vertex” 弹簧系统的动态网格方法,通过添加冗余约束,避免了畸形网格单元的产生。旋翼流场计算采用基于Navier-Stokes(N-S)方程的CFD模块,对基于运动嵌套网格 的空间流场进行求解,湍流模型采用B-L模型。结构分析采用基于中等变形梁理论的CSD模块,基于Hamilton变分原理建立旋翼桨叶动力学方程。首先对振荡NACA0012翼型的流场进行了求解,验证了网格变形模块和CFD模块的有效性,然后采用UH-60A直升机旋翼作为算例对结构动力学模块进行数值验证。在此基础上,计算了UH-60A直升机旋翼桨叶在前飞状态下的非定常气动载荷,并与飞行测试数据进行了对比。计算结果表明,文中的弹簧系统网格变形方 法可以有效地用于旋翼CFD/CSD耦合计算分析,提高了旋翼气弹载荷的预测精度。     

结冰风洞试验段水滴分布特性分析

喷雾系统是结冰风洞的主要组成部分,在结冰风洞试验段直接进行不同状态粒子分布特性的测量,会耗费巨大的成本。为此,本文提出一种采用实验测试和数值计算相结合的手段研究结冰风洞试验段水滴分布特性的方法。通过搭建独立的喷雾粒子试验系统,得到喷嘴出口处的粒子分布特性,在此基础上,采用数值方法计算不同水滴在风洞内的运动及传质传热过程,得到不同水滴蒸发之后的直径,进而获得试验段粒子的分布特性。采用该方法对典型雾化状态下3m×2m结冰风洞试验段粒子分布特性进行了研究,对比了空气湿度的影响。研究发现:(1)喷嘴出口处的初始喷雾粒子与试验段的水滴均保持近似正态的分布,试验段的水滴平均直径(MVD)与初始MVD接近,蒸发不能引起明显的MVD变化;(2)虽然空气湿度越小,水滴蒸发量越大,但湿度为100%时试验段水滴的 MVD比湿度为70%时小。研究成果为结冰风洞喷雾系统设计和调试提供了较好的技术基础。     

基于独立桨距控制的电控旋翼主动振动控制

电控旋翼利用桨叶后缘襟翼偏转通过气弹作用带动桨叶变距,从而实现对旋翼的控制。但由于其系统的复杂性,桨叶间易存在扭转刚度和质量不相似,从而引起较严重的旋翼振动。针对该情况,提出了基于独立桨距控制的主动振动控制方法,并以某原理性电控旋翼为算例进行了数值仿真,验证了所提出控制方法对电控旋翼由扭转刚度和质量不相似引起振动的减振有效性,最佳减振水平可达90%。     

共轴刚性旋翼气动干扰数值计算方法

建立了一个适用于共轴刚性旋翼气动特性分析的数值模拟方法。该方法采用任意拉格朗日欧拉方法(Arbitrary Lagrange Euler,ALE)描述的可压缩Navier-Stokes(N-S)方程求解流场,采用低数值耗散的Roe格式进行空间离散;使用多重嵌套网格方法以模拟双旋翼的运动。针对共轴刚性旋翼配平,引入"差量修正"策略解决了传统配平中雅克比矩阵计算复杂的问题。首先,对Harrington-2共轴双旋翼的悬停气动性能进行了计算,然后,对某2 m直径共轴双旋翼的悬停及前飞状态进行了计算,并与试验值进行了对比。结果表明:在典型状态下拉力系数的计算结果与试验值误差在3%以内,扭矩系数的计算结果与试验值误差基本在5%以内;所采用的数值计算方法对旋翼涡尾迹特征具有较高的捕捉精度,可以有效模拟共轴刚性旋翼悬停和小速度前飞下的复杂流场及其细节特征。     

针对旋翼流场模拟的自适应网格算法

针对基于格心格式求解器的旋翼流场模拟,提出了相应的自适应笛卡尔网格的数据存储结构及自适应算法。给出了相应的单元处理策略,简化了对自适应笛卡尔网格的处理;对于频繁的自适应加密过程中产生的大量重复点,采用了高效的交替数字树算法(Alternating digital tree,ADT)予以删除;对于自适应疏化过程中产生的大量无用点,提出了标记-删除-移动(Mark,delete,move,MDM)算法予以快速地删除,减少了不必要的计算资源消耗。对CaradonnaTung旋翼在不同悬停状态下进行了模拟验证,对比了压力分布系数与桨尖涡位置。之后对HELISHAPE 7A旋翼在前飞情况下进行了模拟验证,计算值与实验值吻合。此外,求解器对桨尖涡的捕捉效果得到了明显的提高,表明本文方法具有良好的有效性与鲁棒性。     

气动除冰类飞机结冰风洞试验适航审定技术

民用飞机为获得型号合格证,应按照有关结冰适航规章条款进行结冰适航验证,如何解读适航条款要求并制定有效的符合性验证流程,是进行适航合格审定的关键。本文以Y12F气动除冰飞机结冰风洞试验的实际工程为例,以相关结冰适航文件为基础,结合国际最新的飞机结冰研究成果,研究并总结了目标试验状态设定、设备选择、试验状态的等效转换、模型研制、风洞试验等适航审定要求和技术。构建的结冰风洞试验适航审定方法,有效地指导完成了Y12F飞机除冰系统的适航验证工作,完成中国民用航空局（CAAC）和美国联邦航空管理局（FAA）同时审查,为获得CAAC和FAA型号合格证奠定良好基础。     

气动除冰飞机结冰风洞试验技术

民用飞机为获得型号合格证,应按照有关结冰适航规章条款进行结冰适航验证,结冰风洞试验是获取临界冰形的有效途径。本文以Y12F飞机结冰风洞试验实际工程过程为例,总结与分析了气动除冰飞机结冰风洞试验的模型设计、气动除冰套模拟、试验状态转换、试验流程、冰形测量等关键技术,并给出代表性试验结果。试验结果表明:在典型最大结冰条件下,除冰套工作正常,除冰套循环工作期间正常除冰;除冰套工作间歇,机翼前缘结冰表现为前缘形成较为光滑、厚度约为5~6mm的冰帽,上翼面产生1道高度为3~4mm的冰脊,下翼面形成3道2~4mm的冰脊。     

四旋翼倾转飞行器操纵冗余设计

建立了四旋翼倾转飞行器旋翼、机翼、机身的非线性气动模型和飞行动力学模型,并在Matlab/Simulink环境下搭建了仿真模型,根据仿真模型计算各个操纵面的操纵效率,确定了不同状态下各个通道需要用到的操纵方式,结合飞行动力学模型和操纵方式,在纵向通道上进行了全包线的配平计算,表明了操纵策略的可行性。     

基于合成射流的旋翼翼型动态失速控制研究

针对直升机旋翼工作环境下来流速度和迎角(Angle of attack,AoA)耦合引起的动态失速问题,建立了基于合成射流的旋翼动态失速控制的数值模拟方法。采用运动嵌套网格方法,通过对翼型的平移和旋转实现变来流速度-变迎角的耦合。以积分形式的雷诺平均N-S方程为主控方程,空间离散使用Roe格式,时间离散为隐式LU-SGS方法,以OA209翼型为研究对象,在翼型上表面放置合成射流激振器,开展了射流位置、动量系数、无量纲频率以及偏角等参数对轻度失速、深度失速下翼型动态失速控制的研究。研究发现,轻度失速下,射流位置靠近气流分离点时(20%c附近,c为翼型弦长),对逆压梯度引起的轻度失速控制效果最佳。深度失速下气流分离点虽在5%c之前,但射流位于前缘分离泡后端(10%c附近)时控制效果较好。大迎角需要较大的动量系数才能有效控制。射流频率对涡结构的尺寸和数量会产生一定影响,能改变气动特性波动幅度。较小的射流偏角对轻度失速的控制更有效,而深度失速则需要较大的偏角。