基于二级声辐射模型的地面声场高效预测方法

建立了可用于直升机运动状态下噪声传播特性分析的高效地面声场预测方法。该方法包括计算旋翼气动、噪声的自由尾迹方法和时域FW-H方程;计入大气、地面等环境因素的声传播模型及地面声场计算模型;为提高计算效率,在声源计算与噪声传播之间引入基于“紧致球声源”的二级声辐射模型;在此基础上,还提出了通过建立“特征参数声辐射球库”以实现直升机运动状态下噪声实时预测的方法。以AH-1旋翼为算例,通过与地面声场直接计算法对比,说明了方法的有效性及高效性;此外,文中还分析了大气、地面等环境因素对噪声传播以及地面声场特性的影响。      

电控旋翼气动特性建模与风洞试验验证

首先建立了带襟翼翼型的非定常气动力模型,继而基于Peters-He广义动态尾迹理论,考虑襟翼偏转对电控旋翼叶素环境的影响,建立了电控旋翼有限状态尾迹模型;进一步基于Theodorsen理论推导出电控旋翼桨叶挥舞响应与桨叶变距和襟翼操纵量的关系,综合以上建立了电控旋翼气动特性分析模型.以改进型电控旋翼试验系统为平台进行了风洞试验,测量了不同风速、不同襟翼操纵条件下的电控旋翼气动力、桨距、襟翼偏角及旋翼挥舞角的变化情况.理论计算结果与试验数据符合情况良好,验证了所建立的分析模型的正确性,并得出以下结论:旋翼转速一定时,桨叶变距与襟翼操纵基本呈线性关系;旋翼拉力随襟翼总距的增加而逐渐减小,襟翼总距较大时,其实际气动效率略有下降;前飞状态时,襟翼总距操纵会引起桨叶的纵向周期变距.      

地面效应中前飞旋翼的自由尾迹计算

发展了一种自由尾迹模型,用于对地面效应状态下的旋翼前飞流场进行数值模拟,并求解该状态下的旋翼气动特性.为了扩大模型的适用范围,使用了面元法模拟地面对流场的影响.为提高尾迹迭代的收敛性和精度,在计算模型中引入了地下尾迹等比例修正、环式地面面元分布等改进措施.结果显示:前进比增加,地面附近的旋翼流场由环流模式转变为地面涡模式.在前飞速度很小时出现的环流,对旋翼前缘的轴向下洗速度影响十分强烈,导致此时出现随前进比增加,旋翼拉力减小、需用功率几乎不变的特殊现象.将得到的旋翼气动特性参数随前进比变化情况同实验值进行了对比,证实了新构建模型的准确性.      

共轴双旋翼悬停地面效应气动特性分析

基于时间步进算法建立了适用于共轴双旋翼气动特性分析的旋翼自由尾迹模型,并采用面元法对水平和倾斜地面效应影响下共轴双旋翼系统气动特性进行了研究。在建模的过程中,采用Weissinger-L一阶升力面模型模拟了桨叶的三维效应,并充分考虑了旋翼/旋翼和旋翼/地面之间的气动干扰。通过与尾迹几何和诱导速度试验数据的对比,验证了计算模型的可行性。在此基础上,对共轴双旋翼在地面效应影响下的旋翼尾迹几何形状、流场诱导速度矢量分布和上下两旋翼桨叶的拉力系数分布进行了计算及分析。结果表明,悬停状态下,共轴双旋翼上下旋翼间存在强烈的气动干扰,且地面的影响使旋翼尾迹涡线径向扩展且向上卷起,对共轴双旋翼下旋翼拉力产生明显的影响。     

围绕旋翼飞行器的三维结构化运动嵌套网格生成方法

考虑到旋翼飞行器的旋翼各种运动和机身复杂外形特点,建立了一种鲁棒性高的三维结构化嵌套网格生成方法。采用三维Poisson方程对初始网格进行迭代,获得了高质量的旋翼桨叶三维结构化网格。以此方法为基础,以复杂外形几何体表面网格为初始网格,建立了一种抛物型法向外推方法。推进过程中,使用经抛物化处理的新Poisson方程进行光顺迭代,源项由Hilgenstock法确定,有效克服了椭圆型方程生成网格时边界点较难调整的困难。通过网格最小正交度的检测,表明了该抛物型方法的先进性。进一步将高效的挖洞法("Top map"法)和贡献单元搜寻法("Inverse map"法)引入上述三维结构化网格生成中,建立了围绕旋翼飞行器的运动嵌套网格生成方法,采用RANS方程分别对Caradonna-Tung旋翼和V-22倾转旋翼机流场进行计算流体力学(Computational fluid dynamics,CFD)计算,表明该网格生成方法能够满足CFD计算的需求。     

基于目标压力分布的旋翼先进气动外形反设计分析方法

基于雷诺平均Navier-Stokes（RANS）方程、运动嵌套网格、目标压力分布及余量修正方法,构建了一套直升机旋翼桨叶先进气动外形反设计方法。为避免由桨叶气动外形变化导致的网格畸变,发展了一套基于Poisson方程求解的旋翼桨叶结构化贴体正交网格的快速、自动化生成方法,为提高运动嵌套网格的生成质量和通用性,采用剖面间网格插值与桨尖翻折相结合的方法,同时建立了基于“Top Map”和“Inverse Map”相结合的洞边界划定和贡献单元搜寻的新方法。基于Navier-Stokes方程和双时间法建立了旋翼非定常流场模拟方法,通量求解采用Roe-MUSCL格式,并使用低速预处理法来克服前飞旋翼流场收敛中遇到的刚性问题。在计算流体力学（CFD）方法基础上,基于旋翼翼型压力系数余量联立各方位角处的反设计MGM（Modified Garabedia-McFadden）超定方程组,并依据激波分离、失速等约束设置了各方位角处的反设计权重系数,创建了基于MGM超定方程组最小二乘解的旋翼气动外形（翼型）设计方法。应用所建立的方法,分别针对多目标、多状态和前飞时的旋翼（翼型）气动外形进行反设计分析,验证了本文方法的有效性。最后,将该方法拓展应用到旋翼桨尖气动外形设计中,设计得到与UH-60A直升机旋翼气动特性相似的矩形桨叶外形。     

基于非定常面元/黏性涡粒子混合法的旋翼/平尾非定常气动干扰

旋翼/平尾非定常气动干扰是导致直升机纵向“抬头（Pitch-up）”现象的主要原因。为在直升机设计阶段准确分析旋翼/平尾非定常气动干扰以及由此引起的低速纵向操纵特性变化,通过涡量等效原则和Neumann物面边界条件建立了适用于旋翼/平尾气动干扰分析的非定常面元/黏性涡粒子混合法。该方法耦合了考虑尾迹时变效应的非定常面元法、黏性涡粒子法及涡量镜面法,以准确模拟旋翼和平尾的非定常气动载荷、旋翼尾迹的非定常特性以及旋翼尾迹对平尾的气动干扰效应。首先通过计算NASA ROBIN（Rotor Body Interaction）旋翼尾迹几何和诱导速度分布,并与实验测量值、时间精确自由尾迹及CFD计算结果对比验证方法的准确性。相比于时间精确自由尾迹,本文方法计算精度更高。随后分析了旋翼/平尾非定常气动干扰对平尾向下气动载荷和气动导数的影响,并分析了平尾构型对旋翼/平尾非定常气动干扰的影响规律。分析表明：旋翼尾迹与平尾干扰导致低速状态的平尾载荷突增,气动导数反号;低平尾气动载荷突增较大,高平尾较小,但高速气动导数反号;前置平尾载荷突增量减小,但对应速度范围较宽;右旋直升机右平尾载荷突增量较小,但气动导数特性基本不变。     

计入离心力影响的直升机旋翼翼型结冰数值模拟

建立了一套计入离心力影响的直升机旋翼翼型结冰的数值模拟方法.首先生成围绕翼型的贴体正交网格,然后用Navier-Stokes（N-S）方程求解黏性绕流流场.在此基础上,利用拉格朗日法建立水滴运动方程.其中,为提高计算效率,提出了结合位移矢量的水滴所处单元寻觅方法.最后,结合桨叶工作特点,发展了一种计入离心力影响的三维结冰模型.通过与桨叶结冰实验的对比,验证了本文结冰预测方法的可靠性.对比常规结冰模型,桨叶结冰量减少22.3%;若考虑桨叶的挥舞运动影响,桨叶结冰量进一步减少,表明了离心力及桨叶运动在结冰数值模拟中的重要性.通过不同剖面间的结冰量和冰形对比,分析并获得了桨叶结冰特征.结果表明离心力的影响程度随径向位移的增加而增加,下翼面结冰量随挥舞角的增加而减少.      

基于CFD方法的直升机旋翼/尾桨非定常气动干扰计算

采用CFD方法针对直升机旋翼/尾桨气动干扰问题进行了数值计算研究.建立了一个适用于旋翼/尾桨干扰气动力计算的方法.在该方法中,控制方程采用三维非定常雷诺平均Navier-Stokes(RANS)方程,网格系统使用嵌套网格技术以计及旋翼和尾桨的相对运动.基于所建立的方法,分别计算研究了悬停状态、低速前飞状态以及巡航平飞状态时干扰状态的旋翼和尾桨气动力,并深入分析了尾桨旋转方向对干扰特性的影响.计算结果表明:旋翼对尾桨的气动干扰在各计算状态都存在,且干扰影响不总是负面的,这取决于飞行状态;而尾桨对旋翼的干扰作用则在悬停状态时较为明显,但在前飞状态时可忽略不计;尾桨旋转方向对尾桨气动力的影响较大,合理地选择尾桨的旋转方向可以改善其气动特性.      

基于高阶WENO格式的旋翼非定常涡流场数值模拟

为提高直升机旋翼黏性涡流场计算流体力学（CFD）的计算精度及降低数值耗散,发展了一套基于五阶加权本质无振荡（WENO）格式的旋翼非定常涡流场数值计算方法。采用运动嵌套网格方法生成围绕前飞状态旋翼的网格系统,以Navier-Stokes方程作为主控方程,湍流模型采用一方程Spalart-Allmaras模型。为了提高旋翼流场中对涡的形成、演化等发展过程的模拟精度,采用高间断分辨率的Roe-WENO格式计算对流通量。对状态变量的插值,选取适当的加权因子,通过多个重构模板的凸组合来构造五阶WENO格式,在交界面附近获得具有五阶精度的状态变量。为了提高计算效率,采用高效的隐式LU-SGS（Lower-Upper Symmetric Gauss-Seidel）推进方法和并行计算策略。最后,运用该方法分别对悬停状态下C-T（Caradonna-Tung）旋翼和UH-60A旋翼以及前飞状态下C-T旋翼和SA349/2旋翼的涡流场及气动特性进行了数值模拟,并将数值结果与传统二阶精度格式的计算值进行了对比。结果表明：在计算时间增长10%～20%的前提下,五阶WENO格式能够捕捉到更精确的涡流场特性,反映了五阶WENO格式在计算旋翼非定常涡流场时具有更高的计算精度与更低的数值耗散特性。