纵列式双旋翼直升机气动力和流场的数值模拟

纵列式双旋翼直升机的旋翼尾迹受前飞自由流、机身和旋翼之间的干扰影响,旋翼尾迹不仅影响机身的空气动力,也影响旋翼自身的空气动力特性.应用欧拉方程与计算流体力学(CFD)方法,采用自适应非结构直角网格技术生成计算域网格,数值模拟纵列式双旋翼直升机全机流场,给出了考虑旋翼影响的机身表面压力分布以及全机流场特性.     

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

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

直升机传热与红外特征建模与分析(英文)

直升机蒙皮表面及排气温度分布受到发动机排气系统、旋翼下洗气流和太阳辐射的直接影响。为了能够精确地模拟直升机蒙皮表面及排气温度分布,在进行处于悬停状态的直升机三维流场和温度场建模与分析时,综合考虑了旋翼下洗对喷流的影响和太阳辐射对机身加热的影响。基于温度场和流场数据,通过正反射线踪迹法对直升机红外波段辐射强度的分布进行计算。针对一假想的排气系统与后机身融合的直升机模型,本文对其热特性和红外辐射特性进行了数值计算和分析。     

基于动量源方法的直升机旋翼/机身流场数值模拟

建立了一套基于动量源方法的直升机旋翼/机身流场数值模拟方法.将旋翼对流场的作用以动量源项的形式代表,建立了一个包含动量源项的CFD方法.针对直升机旋翼/机身流场的特点以及动量源方法对网格系统的要求,采用了一种混合网格生成方法.通过结合旋翼桨叶的运动方式、几何特征及气动特性,建立一个包含动量源项的N-S方程的旋翼流场计算方法和迭代流程,并采用该方法进行直升机旋翼、机身和旋翼/机身下洗流场的数值模拟,得到了关于旋翼/机身干扰流场的一些有意义的结论.     

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

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

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

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

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

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

基于配平的共轴刚性双旋翼直升机飞行性能计算与分析

以共轴刚性双旋翼直升机为研究对象,对其飞行性能展开研究。首先,建立了考虑双旋翼间干扰、旋翼以及机身对尾推的干扰、旋翼对机身及尾翼干扰的直升机飞行力学模型。其次,建立了共轴刚性双旋翼直升机需用功率计算模型,以X2直升机为算例,计算了其配平状态下的需用功率,并用文献中的试验数据对需用功率计算结果进行了验证。然后,建立了发动机功率与耗油率模型以及直升机飞行性能计算模型。最后,计算并对比了X2直升机带推力螺旋桨模式和纯直升机模式下的飞行性能。研究结果表明,带推力螺旋桨模式的直升机飞行性能优于纯直升机模式的飞行性能。     

直升机机身、平尾、垂尾的气动干扰

运用面元法模拟直升机机身流场,计算前飞状态的ROBIN(rotor body interaction)模型孤立机身顶部中线的压力系数分布,并与参考数据、CFD计算结果对比,验证了本方法的准确性.采用离散涡系模拟平尾、垂尾和短翼等升力面,以带短翼的UH-60直升机为例,研究了升力面参数变化对机身、平尾、垂尾气动干扰的影响.结果表明:改变升力面的安装角显著改变了气流对其周围的绕流情况,在参数中影响最大;平尾参数尤其是安装角会对垂尾和短翼的压力系数产生较大影响;减小垂尾展长和增大短翼安装角会提升各自的压力系数.     

悬停时的旋翼/机身气动干扰研究

直升机在悬停时旋翼/机身的气动干扰,主要引起机身垂向阻力的增加。本文根据参考文献,给出估算垂向阻力的理论方法,并以某机为例进行了计算,计算结果与该机的旋翼/机身组合风洞试验结果相比,获得了很好的一致。     

Bump进气道中鼓包诱导的激波/边界层干扰特性

为了探索Bump进气道中鼓包诱导的锥形激波和机身发展而来的湍流边界层干扰问题,分析其气动优势,首先选取了半锥和半棱锥这两种与鼓包的流场结构具有一定相似性的构型作为参照,采用数值仿真方法,分别对这三类典型的三维激波/湍流边界层干扰问题进行了流场分析。在此基础之上,设计了三个不同马赫数的鼓包,并研究了设计马赫数对鼓包流场特性的影响。结果表明:当三类构型的无黏激波强度相等时,半锥诱导产生的旋涡强度最强,鼓包次之,半棱锥最弱。尽管鼓包诱导的流场非常复杂,其干扰流场却呈现出准锥形相似的特性。虽然半锥对边界层的排移能力最强,但是综合考虑边界层排移能力及进气道出口流场畸变下,鼓包最具优势,这也是其被选为超声速进气道前缘压缩面的重要原因之一。此外,在设计状态下,适当增加设计马赫数能改善鼓包排移边界层的能力,但设计马赫数太高,边界层排移能力基本不变,反而使得进气道总压损失急剧增加。      

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

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

非线性叶间黏弹减摆器对直升机空中共振的影响分析

 建立带非线性叶间黏弹减摆器的直升机旋翼/机体耦合动稳定性分析模型。与全机飞行力学平衡计算相结合，旋翼/机体耦合动稳定性分析模型考虑前飞状态桨叶变距操纵、机体姿态角和桨毂纵向安装角。针对具有非线性特性的叶间黏弹减摆器，采用基于复模量的非线性VKS改进模型、Simulink时域仿真和多桨叶坐标变换等效阻尼识别法分析直升机悬停、前飞状态下旋翼/机体耦合动稳定性及减摆器双频动幅值，并就减摆器布局、全机总重以及前飞速度对桨叶摆振后退型模态阻尼的影响进行分析。结果表明：由悬停到前飞直升机动稳定性一般均下降，一定速度后又上升；加上减摆器能消除前飞不稳定区；叶间黏弹减摆器抬头连接能提高模态阻尼。     

直升机粒子分离器三维两相流场的数值模拟

本文以某型直升机粒子分离器为研究对象, 在三维非正交曲线坐标系下计算湍流分叉流动和砂粒运动, 结合边界条件特殊处理, 采用分块计算来考虑叶片对流场的影响, 首次对有叶片粒子分离器三维两相流动进行了数值模拟。      

倾转旋翼机多部件对机翼气动干扰的分析及优化

为研究倾转旋翼机各部件对机翼的气动干扰效应，建立了一套适用于倾转旋翼机流场CFD求解方法，提出并生成了一套由三棱柱/四面体组成的适用于倾转旋翼机多种飞行状态的非结构混合网格系统.以三维Navier-Stokes方程为主控方程，采用动量源方法进行倾转旋翼的模拟，并引入了高效的OpenMP并行加速技术等，提出了一套新型的动量源项添加及搜索方法.计算分析了倾转旋翼机旋翼/机身/短舱对于机翼气动特性的干扰影响，得出了一些对设计有指导意义的结论.采用基于径向基函数（RBF）的代理模型方法，根据机翼不同展向位置的干扰程度对各段翼型配置进行气动优化.优化结果表明:考虑气动干扰作用，在巡航速度下优化后的全机升阻比增长达到了36.78%.      

用运动嵌套网格方法数值模拟旋翼前飞非定常流场

通过求解Ｅｕｌｅｒ方程数值模拟了直升机旋翼前飞非定常流场。为了模拟包括旋转、周期性变距和周期性挥舞的非定常运动,本文采用了一种能够快速完成重叠网格间流场信息交换的运动嵌套网格方法。空间采用中心平均的有限体积法进行离散,时间方向采用双时间法隐式推进求解。为了验证方法的正确性,数值模拟了悬停流场,数值结果和实验值吻合很好。尽管前飞状态的计算没有实验验证,但计算结果与其它文献计算结果吻合很好。