基于N-S方程的剪刀式尾桨前飞状态气动力计算研究

建立了一个适用于计算剪刀式尾桨前飞状态气动力的数值方法。在该方法中,采用惯性坐标系下的非定常Navier-Stokes方程作为主控方程,使用低数值耗散的Roe格式进行空间离散,并应用高效的隐式LU-SGS格式进行时间推进;采用双时间法以模拟尾桨前飞时的非定常现象;使用运动嵌套网格方法以计入尾桨桨叶的旋转运动。应用建立的方法,首先对有试验数据可作对比的Lynx直升机尾桨悬停状态气动性能进行了计算,验证了所建立方法对尾桨气动性能计算的有效性;同时,针对Helishape 7A旋翼进行了前飞状态气动载荷计算,验证了该方法对旋翼(尾桨)前飞气动载荷计算的有效性。最后,对剪刀式尾桨前飞状态进行了气动力计算,并分析了构型参数对剪刀式尾桨气动力的影响。计算结果表明,前飞状态时,相同总距下剪刀式尾桨的平均拉力要大于常规尾桨,且剪刀式尾桨在90°-180°-270°方位角内的桨-涡干扰现象较常规尾桨严重。     

旋翼干扰下的直升机尾桨气动噪声计算研究

针对旋翼干扰下的直升机尾桨气动噪声计算问题,将CFD方法与Farassat 1A(F1A)公式相结合,提出了一种适用于旋翼干扰下的尾桨气动噪声数值计算方法。首先,应用所提出的方法进行了噪声的算例验证,并与参考文献结果进行了对比;然后,着重针对悬停和前飞状态旋翼干扰下的尾桨噪声及气动特性进行了计算分析,并与孤立尾桨状态进行了对比。计算分析结果表明:悬停状态下,旋翼干扰下的尾桨气动噪声水平显著大于孤立尾桨;而在前飞状态,由于旋翼尾迹对尾桨桨盘平面产生更强的气动干扰,这一现象会更加明显;同时,在旋翼干扰作用下,尾桨噪声的主传播方向也会发生明显变化。     

四旋翼无人机流场及气动干扰数值模拟研究

主要研究了一般的四旋翼无人机在悬停和前飞状态下的流场数值模拟方法和前飞状态下的气动干扰。以普通"十"字构型布局的四旋翼无人机为例,建立了基于CFD的网格分区、旋转网格和雷诺时均方程等高精度算法的流场数值模拟方法,并进行了试验验证。以此为基础,进一步对四旋翼无人机在不同前飞状态下的气动干扰进行数值模拟,以加深理解四旋翼飞行器的空气动力学特性,可为其总体气动布局设计和飞控设计提供气动力变化规律支持。     

直升机旋翼对尾桨非定常气动载荷的影响

悬停和侧滑状态的直升机主旋翼桨尖涡将穿透尾桨桨尖平面,由此导致尾桨非定常气动载荷发生明显变化。为更准确地模拟由主旋翼/尾桨干扰产生的尾桨非定常气动载荷变化,通过在面元压力项中增加由旋翼桨尖涡诱导的时变项,体现旋翼桨尖涡速度和几何时变对桨叶非定常压力的影响,同时采用涡面镜像法修正涡粒子法的黏性项,确保桨叶附近区域旋翼涡量守恒,建立旋翼尾迹对尾桨叶的非定常气动干扰模型,并耦合面元/黏性涡粒子法,构建直升机主旋翼/尾桨干扰下的尾桨非定常气动载荷分析方法。通过计算AH-1G旋翼桨叶非定常气动载荷特性,并与实验测量值、计算流体力学(CFD)计算结果对比,验证本文非定常气动干扰模型的有效性。随后基于NASA ROBIN(Rotor Body Interaction)模型分析悬停、侧风和60°右侧滑状态主旋翼对尾桨非定常气动载荷的影响,分析表明主旋翼尾迹对尾桨非定常气动载荷影响显著。悬停状态的主旋翼/尾桨干扰导致尾桨拉力平均值下降、非定常气动载荷显著增加;左侧风状态,主旋翼/尾桨干扰削弱尾桨"涡环"程度,显著增加尾桨拉力和非定常气动载荷;60°右侧滑状态,主旋翼/尾桨干扰导致尾桨拉力损失最大,且在低速侧滑状态出现尾桨拉力"迅速恢复"现象,尾桨非定常气动载荷幅值迅速增加。     

基于嵌套网格方法的直升机尾梁边条气动干扰数值模拟研究

采用基于运动嵌套网格技术的计算流体力学方法开展旋翼/机身气动干扰计算分析,研究了45°、90°、120°三个安装角度尾梁边条对直升机气动特性的影响。数值分析结果表明,悬停状态尾梁边条可为机身提供侧向力、增加机身扭矩,能够帮助尾桨卸载。其中,120°尾梁边条悬停状态增加机身扭矩效果最明显。且安装120°尾梁边条后直升机在前飞状态的气动特性基本没有影响。     

基于Navier-Stokes方程的旋翼前飞状态气动特性数值模拟研究

建立了一个适用于旋翼前飞状态非定常流场和气动载荷计算的数值方法。在该方法中,控制方程为惯性坐标系下的三维非定常Navier-Stokes方程,空间方向上将三阶逆风格式(MUSCL)与通量差分分裂方法相结合。为较好地模拟流动分离,采用了一方程的Spalart-Allmaras湍流模型,时间方向则采用双时间法推进求解。为计入桨叶的旋转、周期性挥舞和变距运动,采用了运动嵌套网格方法。此外,为了更真实地反映旋翼桨叶的实际运动,发展了一个旋翼配平分析模型及求解方法。应用所建立的方法,首先计算了具有先进气动外形的四片桨叶的UH-60A直升机旋翼悬停状态的流场及气动性能,验证了该方法对悬停状态旋翼气动性能数值模拟的有效性。在此基础上,着重对Caradonna模型旋翼无升力前飞状态和AH-1G直升机旋翼桨-涡干扰前飞状态的气动特性(表面压强分布、扭矩、气动载荷)进行了计算,数值结果和试验值吻合良好,表明本文建立的CFD方法对旋翼前飞状态流场的模拟和气动载荷的计算是有效的。     

直升机旋翼/尾桨/垂尾气动干扰计算研究

建立了一个基于计算流体力学(CFD)技术的直升机旋翼/尾桨/垂尾气动干扰分析方法。在该方法中,选取Navier-Stokes方程为控制方程,使用二阶迎风的Roe格式进行空间离散,并选取隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式进行时间推进,湍流模型为B-L(Baldwin-Lomax)模型;为了实现旋翼、尾桨和垂尾网格之间的流场信息交换,采用运动嵌套网格方法。应用所建立的方法,对Helishape 7A旋翼和Lynx直升机尾桨进行了算例计算,并与试验数据进行对比,验证了方法的正确性。着重针对旋翼/尾桨干扰特性进行了计算,并进一步计入垂尾的干扰,对垂尾/尾桨干扰以及旋翼/尾桨/垂尾干扰特性进行了研究,分析了旋翼、尾桨和垂尾相互干扰的规律。结果表明:对于不同的垂尾/尾桨构型,阻塞面积越大,对应的尾桨拉力也较大,但尾桨和垂尾获得的净拉力却减小,且不同阻塞面积下,推力式构型尾桨比拉力式尾桨具有更大的净拉力;然而,在前飞过程中,直升机垂尾对旋翼与尾桨干扰的影响很小。     

应用多块对接结构网格方法的直升机涵道尾桨气动特性分析

基于湍流Navier-Stokes方程,建立了一套直升机涵道尾桨流场及气动特性的CFD分析方法.在该方法中,针对直升机涵道尾桨流场的特点以及动量源方法对网格系统的要求,提出了一种多块对接网格生成方法,分块贴体网格采用求解Poisson方程的方法生成.通过结合涵道尾桨的运动方式、几何特征及气动特征,建立一个包含动量源项的N-S方程的涵道尾桨流场计算方法和迭代流程;为较好地捕捉涵道壁附近存在的旋涡及分离现象,采用了一方程S-A湍流模型.通过对ONERA M6机翼的绕流特性的计算分析,以及对涵道尾桨在悬停状态下的拉力、诱导速度、载荷分布等的计算分析,验证了该CFD方法的有效性.在此基础上,采用该方法进行直升机悬停、侧飞和前飞状态下涵道尾桨流场与气动特性的数值分析,得到了关于涵道尾桨流场和气动特性的一些有意义的结论.     

旋翼气动力数值模拟准确度影响因素计算

针对旋翼气动力CFD计算时的一些常见影响因素的选择问题,进行了数值计算讨论.计算求解器选择自研代码,该代码采用基于共享内存(OMP)的并行化方法,在嵌套网格系统下求解Navier-Stokes方程,空间方向上使用2阶迎风格式,时间方向采用双时间方法.讨论的问题包括背景网格数目对悬停桨尖涡捕捉、对前飞状态激波捕捉和气动力的影响,远场动量理论修正作用对计算结果影响等.对于讨论的不同网格数目情况下,计算的悬停状态旋翼拉力差异在6%左右,扭矩差异2.6%左右,旋翼前飞阻力和扭矩最大差异为2.58%和3.28%.是否采用动量修正旋翼拉力相差3%~5%.通过计算,给出了针对以上因素计算结果差异的范围,可以为旋翼气动力计算工作提供一定的参考.      

剪刀式尾桨悬停状态气动力及噪声特性计算研究

建立了一个基于计算流体力学(CFD) /FW-H(Ffowcs Williams-Hawkings)方程的预测剪刀式尾桨悬停状态气动性能和噪声特性的分析方法.该方法首先采用CFD方法对尾桨流场进行求解,并应用嵌套网格技术对流场空间进行离散.控制方程采用非惯性坐标系下的Navier-Stokes方程,空间方向采用二阶迎风格式(Roe格式)进行求解,时间方向采用隐式LU-SGS(Lower-Upper Symmetric Gauss-Seidel)格式进行推进.在此基础上,采用FW-H方程将尾桨噪声声压扰动传播至远场,以获得尾桨的噪声特性.应用该方法对两种剪刀式尾桨构型(“L”构型和“U”构型)进行了计算研究,对比分析了剪刀式尾桨在气动力和噪声方面与常规尾桨的差别,以及两个重要构型参数(剪刀角和轴向间距)对剪刀式尾桨气动力和噪声的影响规律.计算结果表明,构型参数对剪刀式尾桨气动力和噪声特性影响很大,合理地选择构型参数可以降低尾桨噪声水平.     

倾转四旋翼飞行器垂直飞行状态气动特性

综合采用基于滑移网格技术的计算流体力学(CFD)方法与悬停状态气动干扰试验方法,对倾转四旋翼(QTR)飞行器垂直飞行状态的流场进行模拟与试验,研究飞行器垂直飞行状态气动特性以及部分参数对气动特性的影响.结果表明:倾转四旋翼飞行器在垂直飞行状态,前后旋翼之间干扰不明显,但旋翼与机翼的干扰明显;旋翼旋向对旋翼与机翼的干扰不...     

剪刀式尾桨前飞状态气动噪声特性计算

基于计算流体力学（CFD）建立了适用于剪刀式尾桨的流场计算模型，采用嵌套网格方法模拟尾桨运动，采用双时间方法进行时间推进。针对不同的尾桨构型，采用高效配平策略“差量法”将剪刀式尾桨配平至相同负载状态。在此基础上，采用Ffowcs Williams-Hawkings（FW-H）方程计算剪刀式尾桨的气动噪声特性。应用该方法对不同前飞速度下剪刀式尾桨的气动力和噪声进行计算分析，着重研究了剪刀角和轴向间距两个重要构型参数的影响。计算结果表明:剪刀角对剪刀式尾桨气动力和气动噪声特性均有重要影响，而轴向间距在合理的变化范围内，对尾桨影响较小。与常规尾桨相比，前飞状态下剪刀式尾桨的噪声指向性变化较小，但噪声幅值变化显著。      

改进型CLOR桨尖旋翼气动特性试验研究及数值分析

 通过风洞试验及数值模拟对具有改进型CLOR(CLOR-Ⅱ)桨尖的旋翼悬停和前飞状态气动特性开展研究。在CLOR桨尖旋翼试验及数值分析的基础上,考虑旋翼非定常流场特点,兼顾旋翼悬停和前飞气动性能,对旋翼桨叶的气动外形进行了改进,主要包括采用多种翼型优化配置以综合改善旋翼前行侧压缩性及后行侧桨叶失速特性,并考虑旋翼前飞状态对其桨叶动力学特性的需求,重新设计了桨尖前后缘的外形。在风洞中分别对3种旋翼进行多种状态条件下的试验研究,为从流动细节上获得不同桨尖旋翼的气动特性差别,采用计算流体力学(CFD)方法对试验状态进行了数值模拟对比。对更高转速状态进行模拟,结果表明相对于其他两种旋翼,CLOR-Ⅱ桨尖旋翼在改善跨声速特性和提高失速迎角等方面具有明显优势,而且综合提高了旋翼悬停和前飞气动性能。     

倾转旋翼典型飞行状态气动特性的CFD分析

针对倾转旋翼特有的悬停、巡航、过渡3种飞行模式,基于CFD方法分别对倾转旋翼的上述3种典型飞行状态进行了数值模拟.首先,基于运动嵌套网格思想建立了适合倾转旋翼运动模拟的多层运动嵌套网格系统.然后,基于RANS(Reynolds averaged Navier-Stokes)方程建立了适用于非定常流场分析的CFD模拟方法.为了提高整个流场的计算效率,时间推进采用高效的隐式LU-SGS(lower upper symmetric Gauss Seidel)方法,并运用SPMD(single program multiple data)并行计算模式.最后,采用所建立的网格生成方法和流场求解器,进行了不同飞行状态下的倾转旋翼气动特性数值分析,验证了建立CFD方法的有效性并得出了多种飞行状态下的气动特性.结果表明:悬停状态下倾转旋翼相对于直升机旋翼有更加明显的下洗流影响区域.巡航状态两旋翼中心间距为2.5R时,倾转旋翼之间存在明显的气动干扰.过渡状态下,倾转旋翼提供的拉力在倾转角为45°前下降较快,且在倾转角较小时,倾转旋翼桨叶表面拉力分布仍能表现出直升机旋翼的特性.      

直升机旋翼非定常气动力及结构动力响应计算

以二维非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程为模型方程,采用Ｊａｍｅｓｏｎ有限体积并引入变系数隐式残值光顺及固接网格计算了直升机旋翼不同展向位置处翼剖面的非定常气动力,进而采用插值的方法模拟直升机在悬停及前飞状态下,旋翼的桨叶所受的非常气动力载荷,并将之与桨叶的结构振动方程相耦合,求解其时域内的结构响应,非定常下气动力计算结果与实验结果做了对比,两者比较吻合,旋翼结构响应计算结果合理。     

基于涡环尾迹模型的共轴刚性旋翼直升机飞行动力学建模

共轴刚性旋翼直升机上下旋翼间距小,旋翼间气动干扰较为复杂,影响飞行动力学特性。针对这一问题,利用涡环单元动态尾迹方法构建了共轴旋翼气动力模型,通过与风洞试验结果比对说明该模型能够准确地计算存在气动干扰时共轴旋翼的气动力特性。以该共轴旋翼气动力模型为基础,建立了共轴刚性旋翼直升机飞行动力学模型,并以XH-59A直升机为研究对象,计算了前进比为0~0.4时的配平特性。通过与飞行试验数据的比对发现：该飞行动力学模型与飞行试验结果比对良好;且模型计算速度较快。通过对配平结果以及旋翼尾迹运动的分析发现：共轴刚性旋翼直升机旋翼间气动干扰会增加悬停和低速前飞时的配平总距和总距差动;低速前飞时的纵向周期变距负梯度现象是由于旋翼间气动干扰与刚性旋翼挥舞运动特性叠加而造成的。     

共轴刚性旋翼气动干扰及操纵特性分析

建立了一套基于高效配平策略和嵌套网格技术的共轴刚性旋翼气动数值分析方法。流场求解采用雷诺平均Navier Stokes方程，空间离散采用Roe monotone upstream centered schemes for conservation laws（MUSCL）格式，湍流模型为Spalart Allmaras模型；配平计算中引入“差量修正”策略。分析了共轴刚性旋翼模型在不同飞行状态下的气动干扰现象和操纵特性。结果表明:所建立的配平方法在同等精度下，比传统方法可显著节省计算时间；共轴刚性旋翼间距较小，其气动干扰明显强于普通共轴旋翼；悬停状态下，由于受上旋翼下洗流影响，下旋翼总距较大但拉力较小；小速度前飞时气动干扰引起下旋翼气流不对称性增大，双旋翼周期变距差也变大；采用适当的升力偏置量，可以提高共轴刚性旋翼气动效率。      

悬停状态剪刀式尾桨气动/噪声特性优化分析

建立了一套基于CFD(computational fluid dynamics)方法和FW-H(Ffowcs Williams-Hawkings)方程的剪刀式尾桨气动噪声预估技术和组合优化算法的尾桨外形参数设计方法。基于嵌套网格方法采用雷诺平均Navier-Stokes(RANS)方程作为尾桨流场求解控制方程,采用了一套适用于剪刀式尾桨悬停气动特性模拟的高效CFD方法。在流场分析的基础上,采用FW-H方程预测剪刀式尾桨在典型观察位置处的气动噪声。分别在控制总距以及控制拉力系数不变的情况下,以提高尾桨悬停效率同时降低气动噪声为目标,对剪刀式尾桨的剪刀角和轴间距两个主要外形参数进行优化设计。将基于拉丁超立方(LHS)方法和径向基函数(RBF)的代理模型方法与遗传算法过程相结合,建立了一种有效的组合优化算法。结果表明:剪刀角和轴间距的不同组合可以通过削弱桨-涡干扰现象从而实现降低旋翼桨-涡干扰噪声的目的。在当前的计算状态下,优化得到的剪刀式尾桨的悬停效率比常规尾桨高16%,其平均声压级比常规尾桨降低了2.3 dB。      

前飞状态旋翼动稳定性及自由度间相互作用研究

为研究直升机在前飞状态下的稳定性,建立了前飞状态下旋翼的结构模型、气动力模型和平衡方程组。对方程组的周期系数处理以及前飞时的相关参数对旋翼动稳定性的影响等进行了研究。在前飞动稳定性分析中引入动力入流自由度。从分析各个自由度间的相互作用入手揭示旋翼动不稳定性的物理本质及影响因素。