基于动量源方法的纵列式直升机双旋翼/机身干扰流场分析

针对纵列式双旋翼桨叶的几何特征、运动方式及气动特性,建立了基于动量源方法的纵列式直升机双旋翼/机身组合流场的数值模拟方法,该方法中采用的非结构网格能综合考虑双旋翼/机身干扰流场特点且满足动量源计算网格要求。通过算例的计算值与试验值的对比,验证了本文方法的有效性,然后采用该方法分析了悬停及前飞状态下纵列式直升机的双旋翼/机身干扰流场特性。     

直升机旋翼机身流动特性数值模拟

运用CFD数值工具研究旋翼机身发动机空气动力干扰.通过生成多块贴体结构网格,运用有限体积方法直接求解了不可压N-S方程.首先分析ROBIN机型的旋翼机身耦合的流场特性,然后再分析了MI-171V5直升机低速前飞时的流场性质.对两种机型首先进行孤立机身周围粘性流场数值模拟,然后运用作用盘模型分析了旋翼机身空气动力干扰对机身压力分布等方面的影响,与实验数据对比相吻合.     

高速直升机旋翼/螺旋桨/机身干扰特性分析

研究双拉力螺旋桨复合式高速直升机的气动特性可以为高速直升机的设计及气动优化提供参考。基 于动量源方法构建针对双拉力螺旋桨复合式高速直升机旋翼/螺旋桨/机身干扰特性数值计算及分析方法;对 孤立旋翼、旋翼/机身干扰进行算例验证;应用所构建的方法对双拉力螺旋桨高速复合式直升机悬停及前飞状 态的干扰流场进行数值模拟,分析机身对悬停流场影响及不同前飞速度旋翼/螺旋桨/机身干扰特性。结果表 明:悬停时机身对气流的阻塞作用降低了旋翼的升力,螺旋桨对旋翼下洗气流的加速作用使旋翼升力提高;低 速前飞时旋翼/螺旋桨/机身干扰较大,主要体现在旋翼下洗流造成螺旋桨滑流偏折以及机翼上表面压力分布 增大,高速前飞时这种干扰较小。     

直升机旋翼机身发动机耦合流场数值模拟

运用计算流体力学(CFD)数值工具研究旋翼机身发动机空气动力干扰.通过生成多块贴体网格,运用有限体积方法直接求解了不可压N-S方程.分析了某直升机低速前飞时的流场性质.首先对孤立机身周围粘性流场进行数值模拟,然后运用作用盘模型分析了旋翼机身空气动力干扰对机身压力分布等方面的影响;最后运用简单的发动机模型分析了发动机进出口对直升机流场的影响.      

非结构嵌套网格的直升机旋翼/机身前飞流场数值模拟

建立了一个适用于前飞状态直升机单独旋翼、机身以及旋翼/机身组合体流场的数值计算方法及模型.在该方法中,使用非结构运动嵌套网格描述桨叶之间以及桨叶与机身之间的空间位置变化关系,控制方程采用惯性坐标系下的非定常Navier-Stokes(N-S)主控方程,空间方向上使用二阶迎风格式,时间方向上使用Lower-upper symmetric Gauss-Seidel(LU-SGS)格式.对于求解中遇到的旋翼/机身间距离过近所导致的贡献单元生成困难问题提出了一种新的解决方法.应用所建立的方法,分别对单独旋翼、机身及旋翼机身组合体的前飞绕流场进行了数值模拟与分析,计算结果与试验结果的对比表明本文建立的方法对直升机的前飞绕流场的计算和分析是很有效的.      

旋转机翼飞机旋翼/机身干扰流场数值计算分析

基于动量源方法,建立了旋转机翼飞机的旋翼/机身组合流场的数值模拟方法.首先通过Caradonna-Tung旋翼算例计算,验证了本方法的有效性;然后采用该方法对某旋转机翼飞机的悬停及低速前飞状态下的旋翼/机身干扰流场特性进行了详细的计算与分析.结果表明,该方法能够用于计算旋翼与机身的干扰流场,为旋转机翼飞机的设计提供了参考.     

共轴刚性双旋翼/机身干扰流场数值模拟

采用滑移网格技术求解Navier-Stokes (RANS)方程的方法,研究了共轴刚性双旋翼/机身的干扰问题。通过Caradonna-Tung旋翼、Robin直升机、Maryland直升机旋翼/机身干扰和Harrington 2共轴双旋翼等算例,验证了所提出的旋翼流场数值模拟方法的正确性。在此基础上,以Maryland机身为原型,分析了不同桨距的共轴刚性双旋翼与机身之间的干扰特性。结果表明:所提出的数值模拟方法能够很好地模拟共轴刚性双旋翼/机身的气动干扰特性;由于机身对于共轴刚性双旋翼下洗流场的阻滞作用,旋翼的悬停效率增加5%左右,并且随着拉力系数的增大使得悬停效率的增量更加明显;旋翼的悬停效率增加主要来源于下旋翼0°方位角附近的桨叶升力系数的增大,并且拉力系数的增量由桨根向桨尖方向逐渐减小。      

直升机旋翼对机身气动干扰的计算

基于旋翼自由尾迹模型和三维机身源面元模型，建立了一个全耦合的旋翼／机身气动干扰迭代计算方法。为正确模拟旋翼桨尖涡与机身表面间的大定常贴近干扰，采用了一个“分析数值匹配法”的“贴近涡／面干扰模型”。应用该计算方法，分别计算了旋翼干扰状态下的机身表面点的非定常压强分布和机身的非定常气动升力，俯仰力矩随前进比的变化，以及干扰前后的机身定常压强分布。计算表明，机身上由于旋翼的干扰引起的非定常气动载荷呈现出与旋翼相同的周期性。然后，对比了本“贴近涡面干扰模型”和先前“截断涡线模型”对干扰计算的影响，表明了前对计算结果的重要改进。     

直升机旋翼/机身非定常气动干扰数值分析

基于非定常面元/时间步进全展自由尾迹建立了旋翼/机身非定常气动干扰分析方法。方法中耦合了非定常面元法和时间步进自由尾迹,以准确模拟旋翼非定常气动力、旋翼尾迹及桨叶对机身的非定常干扰效应。为模拟前飞状态下具有升力的机身,将机身离散为汇/偶极子面元,并采用涡线镜像法模拟旋翼尾迹靠近机身表面产生的加速效应。通过计算前飞状态的Maryland、ROBIN(Rotor Body INteraction)旋翼/机身干扰下的非定常压力分布,并与可得到的实验值、CFD计算结果对比,验证方法的准确性。随后分析前飞速度、旋翼与机身距离对旋翼/机身非定气动干扰的影响。计算结果表明机身头部和中部非定常压力主要受桨叶的通过性影响,而机身尾梁主要受尾迹/机身干扰影响,机身非定常气动力频率为桨叶片数的倍频。随前飞速度的增加,机身非定常压力幅值增加,尾梁压力幅值先增加后减小;增加旋翼与机身距离将减小机身和尾梁非定常压力幅值。     

旋翼/机身气动干扰的数值模拟

通过求解粘性非定常可压缩N-S方程,发展了旋翼/机身气动干扰的数值模拟方法,主要技术措施包括:(1)通过ADT数据结构,建立了高效的多块对接网格动态重叠技术;(2)在可压缩方程求解中引入低速预处理方法,克服低速流场求解时的刚性问题;(3)加入了多重网格方法,提高收敛效率;(4)湍流模型方程求解需要反复计算壁面最近距离,采用Wigton的优化算法很大程度上提高了计算速度。利用建立的计算方法对佐治亚大学(GIT)旋翼/机身干扰模型进行了数值模拟,分别比较了机身表面的平均压力分布、旋翼的下洗流场以及非定常压力脉动,对干扰现象及其流动机理进行了分析。     

直升机旋翼/机身耦合系统模型导纳的测试研究

直升机旋翼和机身子系统以及其耦合系统的导纳是用动柔度匹配法确定直升机旋翼 /机身耦合系统固有频率必要的数据 ,我们已推导了它们的分析表达式 ,为验证分析表达式的可用性 ,须进行试验测试与研究。由于旋翼的转动 ,给实测带来困难 ,作为起步 ,从模型试验开始。本文介绍旋翼 /机身耦合系统模型及其子系统的导纳的测试方案、试验程序 ,综述以作者自行设计的旋翼 /机身耦合系统模型进行试验所得到的结果及其与理论计算结果的比较。结果表明 ,测试方案与试验方法是可行的 ,理论分析计算公式是正确的     

基于矩阵运算的直升机旋翼/机体耦合动力学模型

针对直升机旋翼/机体耦合动力学建模中,小角度线性简化带来的计算精度下降以及适用范围缩小等问题,在计算桨叶速度及加速度时,基于矩阵运算,保留了运动的非线性。对不同机体运动幅值和频率情况下的小角度线性简化对桨叶速度的影响进行了仿真分析,结果发现,随着机体运动幅值的增加,线性简化影响越来越大,滚转幅值 15°时,简化带来的误差超过 5%;俯仰幅值 15°时,简化带来的误差超过 6%,且随着机体运动频率的减小,误差越来越大。对不同桨距的情况进行仿真,结果表明,所建模型可以很好地模拟直升机的动态响应,采用该模型可进行直升机旋翼/机体耦合动态响应仿真分析。     

基于非结构嵌套网格方法的旋翼地面效应数值模拟

 建立了一个基于非结构嵌套网格的流场求解器,用来精确模拟复杂的旋翼近地流场,为更好地分析地面效应(IGE)对旋翼气动特性的影响提供一套计算方法。在该求解器中,控制方程采用惯性坐标系下的非定常N-S方程,空间方向上采用二阶迎风格式,并用背景网格的一个面模拟地面作用,以方便地面边界条件的处理。应用所建立的模型,首先针对有实验结果可供对比的旋翼无地面效应（OGE）流场进行了数值模拟,以验证计算方法。然后着重计算了IGE下流场中的桨尖涡空间位置和旋翼拉力增益,并对旋翼在小速度前飞状态下的近地流场及地面涡形成过程进行了模拟和分析。在此基础上,得出了一些有意义的结论。     

基于动量源方法的涵道尾桨CFD分析

本文采用CFD方法分析了悬停状态下涵道尾桨的流场特性。该方法以动量源代替桨叶对流场的作用来分析涵道尾桨的气动特性。即通过编写UDF程序把动量源加入到N-S方程中,建立了一套适合于直升机涵道尾桨流场和气动特性的CFD数值模拟方法。然后利用所建立的数值模拟方法对某个涵道尾桨进行分析,并与试验结果进行对比以验证数值方法的有效性。     

直升机旋翼/机身耦合系统气动/机械稳定性分析

考虑了铰接式旋翼桨叶绕挥舞、摆振和变距铰的整体刚性运动与桨叶中等弹性变形之间的动力学耦合作用,将直升机机身和传动轴作为弹性体。通过构造一种特殊的24自由度刚柔混合单元得到旋翼/机身耦合系统的周期时变动力学方程,根据Floquet理论对稳态周期解的稳定性进行研究。计算结果表明,旋翼轴的截面刚度对旋翼/机身耦合系统的气动/机械稳定性有一定的影响。      

直升机旋翼悬停流场的欧拉方程计算

描述了三维Ｅｕｌｅｒ方程数值模拟旋翼悬停流场的过程,欧拉方程的求解采用了风格中心有限体积法、五步Ｒｕｎｇｅ－Ｋｕｔｔａ时间推进格式。     

基于黏性涡模型的旋翼流场数值方法

 建立了一种适用于旋翼非定常流场特性分析的黏性涡数值方法。在该方法中:流场中的大尺度涡被离散为若干微小的涡元,通过求解涡量-速度形式的Navier-Stokes方程模拟涡元的输运等过程;黏性扩散效应采用高精度的粒子强度交换法进行计算,而桨叶附着涡以及新生涡环量采用了Weissinger-L升力面理论进行求解;为显著提高计算效率,在诱导速度及其梯度的计算中还引入了快速多极子算法(FMM)。应用上述方法,对悬停和前飞状态下的多个旋翼流场算例进行了计算,通过对比旋翼尾迹涡量特征和诱导速度分布等,验证了该方法的有效性。此外,还将本方法与旋翼计算流体力学(CFD)方法及传统的自由尾迹方法进行了比较,结果表明黏性涡方法在兼顾效率的同时,还能够更好地捕捉旋翼尾迹运动。