风扇/增压级带间隙三维粘性流场数值模拟

使用三维数值模拟软件NUMECA(Numerical Mechanics Application)对双涵道风扇/增压级带转子叶尖间隙的粘性流场进行了三维定常流动数值模拟,获得了其设计转速下的流场特性.结果显示,随着外涵背压的升高,风扇顶部激波后阻塞区周向和轴向尺寸增大,激波变形加剧,继续提高外涵压比,则外涵效率和进出口流量下降.内涵背压的升高,使内涵各转子叶片顶部的泄漏流动影响范围向上游扩散.内涵背压过高,效率和进出口流量同样会有下降的趋势.同时通过和实验结果的比较表明,数值模拟软件模拟的风扇/增压级的三维定常粘性流场能够较真实的反映出转子叶尖间隙泄漏流场的结构特征,以及间隙对流场参数的影响.      

北京航空航天大学学报2007年第33卷总目录（总第167期～总第178期）

北京航空航天大学学报2007年第33卷第1期(总第167期)3477后掠翼对细长体中段大迎角绕流的影响……………………………………………………………李国辉邓学蓥3478合成射流作用机理及其对共轴射流掺混的影响………………………………刘艳明伍耐明董金钟周盛3479基于实体单元的转子动力特性计算方法………………………………………陈萌洪杰朱彬刘书国3480跨声风扇转子的BVF气动优化方法……………………………………………吴宏李秋实郭明周盛3481蜂窝夹芯板超高速碰撞仿真………………………………………………………………徐小刚黄海贾光辉3…     

喷管分离流动及其侧向载荷

利用商业软件CFX对某液体火箭大面积比喷管地面条件下的分离流动进行了三维数值模拟.计算获得了喷管入口总压从8MPa减小到1MPa时的流场参数分布和侧向力载荷情况.结果表明,随着入口总压的降低,喷管内流场会依次经历自由激波分离和受限激波分离两种分离激波模态.受限激波分离模态下喷管壁面压强具有较大波动,再附着点压强甚至高于环境压强.流动分离情况下,喷管将受到一定侧向载荷作用,载荷方向随机分布.入口总压为4MPa时计算得到的侧向载荷最大,实际侧向载荷峰值可能出现在自由激波分离与受限激波分离转换瞬间.     

扩张型双喉道喷管的流动特性和起动方法

利用数值模拟方法,对二元扩张型双喉道喷管的流动特性和起动方法进行了研究.结果表明:扩张型双喉道喷管内会出现正激波系,产生了很大的总压损失,使第2喉道壅塞,喷管不能起动.在低落压比条件下,喉道注气可以形成大的分离区,使激波强度减弱、喷管可以起动;在大落压比条件下,喉道注气不能形成大的分离区,喷管不能起动.扩张段注气可以在喷管内形成大的分离区,使正激波转变成斜激波系,减小了总压损失,使第2喉道流通能力增强、喷管起动.      

基于IFA-ELM的航空发动机自适应PID控制新方法

针对大涵道比涡扇发动机强非线性、变参数的特点,提出了一种基于优化极端学习机(ELM)对发动机参数进行预测的自适应PID控制方法.为提高ELM的预测精度和实时性,采用适用于多峰值寻优的改进萤火虫算法(IFA)优化ELM网络参数,形成优化的ELM训练方法IFA-ELM.该算法在保证预测精度的前提下,有效简化了网络规模,并提高了其泛化能力.利用该算法建立发动机风扇转速预测模型,基于该模型,采用梯度下降法在线调整PID参数,提升发动机动态性能.数字仿真验证表明,与常规PID控制相比,基于IFA-ELM的自适应PID法调节时间减少了0.2~1.4s,超调量降低了0.2%~1.5%,验证了该控制方法的有效性.      

用GAO-YONG湍流模型计算激波/边界层干扰

采用GAO-YONG可压缩湍流方程组,模拟了平板激波/湍流边界层干扰现象.运用SIMPLE算法求解方程组,并分别采用三阶精度的QUICK格式和中心格式离散对流项和扩散项.计算结果较好预测了入射斜激波在平直壁面引起湍流附面层分离的流动特征: 分离点的反射激波、分离包引起的膨胀扇以及再附点的反射激波.对流场的时均参数与实验值进行了比较,计算得到的壁面压力分布、摩阻系数分布和速度型与实验值比较吻合很好.结果表明GAO-YONG可压缩湍流方程组能够高精度模拟平板激波/湍流边界层干扰流动.      

激波碰撞干扰流动非定常效应的数值研究

将AUSM+上风格式应用于数值模拟半球绕流激波碰撞干扰的流动中,揭示了该流动在定常来流条件下的非定常本质.分析弓形激波后超音速"喷流"的流动结构及其非定常效应的形成机制,提出一种新的激波碰撞干扰流动的非定常形成机理.      

航空发动机转子结构布局优化设计方法

结构是航空发动机功能、性能及可靠性设计水平的综合体现,一切技术要求、性能指标、强度指标或者结构的安全性和可靠性都应建立在合理的结构布局设计上。提出了航空发动机转子结构布局并对其进行优化设计的观点,同时以典型高推重比涡扇发动机高压转子的结构构型为例,基于试验设计（DOE）的响应面法,应用有限元计算并通过多目标遗传优化算法,分别从抗变形能力、力学环境适应能力以及转子结构效率综合进行相关优化计算,论证了合理的结构布局形式可以大幅度提升转子的力学特性。研究方法对于航空发动机转子系统的初始结构布局设计具有指导意义,可以显著减少结构设计的迭代次数,缩短设计周期。      

软件流水中的循环展开优化

在软件流水中应用循环展开可以实现分数值的启动间距,提高资源的利用率,同时基于展开的优化技术可以降低程序的资源需求和关键路径的长度.提出了基于程序特性的展开因子算法UTBPC(Unrolling Times Based Program Characteristics),解决了循环展开的核心问题——展开因子的确定,同时提出了基于展开的软件数据预取优化技术,提高了软件数据预取的效率.所有这些都在ORC(Open Research Compiler)中实现,并对SPEC 2000中的程序进行了测试,平均性能提高了2.6%.实验结果表明新提出的算法和基于展开的数据预取优化技术提高了编译器的整体性能.      

柔片式密封数值计算及性能分析

针对高速转子柔片式密封系统,建立了密封区域流场的数学模型,采用一种修正的压力耦合方程组的半隐式方法(Semi-Implicit Method for Pressure-Linked Equations Revised, SIMPLER)算法和有限差分法对数学模型进行数值求解,构建了柔片式密封性能分析程序,得到密封系统的压力场分布、速度场分布及泄漏量等,研究了不同结构参数和工况参数对密封性能的影响.分析结果表明:泄漏量随密封压差的增大呈线性增长趋势,而随转子转速的增大略有减小;增加柔片宽度、减小柔片长度、减小前/背板-转子间距均对降低密封系统的泄漏量有益,而柔片安装角度的变化则对密封系统的泄漏量影响不明显.基于流场的结果,柔片和转子的动力学响应将被研究.      

旋翼结冰模型与纵列式直升机平衡特性分析

为了研究纵列式直升机结冰后的飞行动力学问题,基于CFD(ComputationalFluid Dynam ics)技术提出了旋翼结冰的理论模型;基于旋翼冰风洞试验数据,考虑桨叶表面附着冰脱落及桨叶表面局部温度对结冰的影响,提出了旋翼结冰的工程模型.建立了纵列式直升机结冰后的飞行动力学模型,研究了不同飞行速度条件下结冰对样机CH-47B平衡特性的影响,分析了环境温度、结冰时间、液态水含量以及平均水滴直径对旋翼扭矩的影响.计算结果的对比分析显示:两种结冰模型在计算精度上基本一致,表明结冰模型合理,能够用来进行结冰后的飞行动力学分析,为纵列式直升机结冰后的飞行动力学研究奠定了基础.     

共轴刚性旋翼悬停状态桨叶表面压力测量试验与计算研究

针对共轴刚性模型旋翼悬停状态,开展了桨叶表面压力测量试验与数值模拟研究。试验采用微型压力传感器进行桨叶表面压力测量,不仅获得了桨叶表面压力的试验数据,同时为CFD计算方法计算桨叶表面压力提供了验证数据。计算与试验结果对比吻合度良好,验证了CFD计算方法的有效性。研究获得了共轴刚性旋翼上下旋翼桨叶表面的流动情况和压力特性,结果表明:对于上下各4片桨叶的共轴刚性旋翼,桨叶表面压力随着桨叶旋转呈周期性变化,旋转一周出现8个小周期;在上下旋翼扭矩配平的悬停状态,下旋翼桨叶大部分区域受下洗流影响,下旋翼剖面拉力低于上旋翼;在桨尖区域,下旋翼的桨距角大于上旋翼,受各自上洗流的影响,下旋翼剖面拉力高于上旋翼。      

航空发动机转子动力优化设计软件工具研究

通过对航空发动机多转子系统的转子动力学优化模型的深入分析与研究,采用模糊数学中的正态分布的隶属函数的加权之和,对多阶临界转速相对于多个常用工作转速的分布状态进行了描述,并由此构造了目标函数;按照设计工程需要选定了对临界转速有显著影响的多个因素作为设计变量;按照设计规范的要求选定性能约束; 从而成功地建立了航空发动机转子动力学优化数学模型.在Windows98/NT平台上开发了航空发动机转子动力学优化设计软件工具.该工具可以实现多转子系统的转子动力学优化设计和整体转子的有限元变形与应力分析.以航空发动机转子为实例,验证了转子动力学优化设计软件工具在航空发动机优化设计中的应用,取得了满意的优化结果.     

直升机前飞状态旋翼结冰风洞试验研究

为研究直升机旋翼在不同前飞状态的结冰情况,研制了直升机旋翼模型结冰试验台,在中国空气动力研究与发展中心结冰风洞4.8 m×3.2 m试验段中进行了旋翼结冰风洞试验。以某型直升机2 m直径旋翼缩比模型为试验对象,分析了旋翼转速、初始拉力系数对旋翼结冰的影响。结冰试验过程中保持旋翼操纵恒定,利用天平测量了旋翼拉力和扭矩性能的动态变化,并采用二维冰形切割及三维冰形扫描的方式分别测量了桨叶展向典型剖面的翼型及桨叶的整体结冰形态。试验得到了旋翼结冰关键数据,分析结果表明:旋翼桨叶结冰主要集中在桨叶前缘和下表面,结冰会在降低旋翼升力的同时增大旋翼扭矩和功率;小拉杆杆端轴承的积冰可能造成卡塞,变距拉杆上的积冰可能造成杆端轴承卡塞,从而使旋翼操纵失效。      

悬停状态旋翼固定尾迹分析中三角级数的应用

为了在固定尾迹分析中更好地描述桨叶附着涡环量分布,结合叶素理论,采用一个三角级数表达附着涡环量的连续分布,并应用于悬停状态下单旋翼和共轴双旋翼的气动特性分析,通过小量计算即获得连续光滑的附着涡环量分布,将其应用到推导的公式中计算出连续分布的诱导速度,进而分析了共轴双旋翼上下旋翼气动干扰及扭矩平衡条件.计算结果与实验结果进行了比较分析,验证了该方法的有效性,并得出结论:下旋翼总距变化对上旋翼诱导作用影响微小,但对自身诱导作用明显;当上下旋翼扭矩平衡时,下旋翼总距比上旋翼总距约大1°.     

基于小扰动理论的桨叶叶素气动载荷计算方法

针对直升机飞行仿真应用,采用小扰动理论建立了桨叶叶素气动载荷计算方法。该计算方法集成了弹性桨叶挥舞-摆振-扭转模态,耦合了旋翼动态入流、非定常动态失速气动模型和旋翼配平模型。整个算法采用标准的状态空间格式。为验证算法,以UH-60A直升机为例,在低速和高速2种定常飞行情况下,研究了偏航流效应对旋翼气功性能的影响,结果表明高速飞行时偏航流效应对旋翼气功性能影响较大。模拟了桨盘平面诱导速度分布以及桨叶叶素气动载荷,将仿真结果与飞行试验结果进行对比验证,本文建立的方法可以准确地求解定常飞行时桨叶的非定常气动载荷,捕捉其沿方位角的变化特征,对于大速度前飞仍能适用。      

垂直下降状态下的旋翼三维流场数值模拟

应用激励盘模型,采用计算流体力学(CFD, Computational Fluid Dynamics)方法结合动量叶素理论,计算了旋翼在垂直下降时的流场特征与气动性能.计算过程中,旋翼载荷被描绘成沿桨盘径向分布、与当地流动参数相关的压力源项.在直角坐标系中,运用有限体积方法直接求解了定常不可压N-S(Navier-Stokes)方程,湍流模型为S-A(Spalart-Allmaras)一方程模型,重点使用诱导速度经验公式计算出了压力源项.旋翼流场模拟结果和性能预计结果同实验测量数据吻合较好,表明这种CFD方法可以有效地分析旋翼下降飞行时的空气流动特性,为进一步的涡环状态研究提供了参考.     

表面粗糙或带凸起转盘风阻扭矩实验

为研究壁面不同粗糙程度以及带有螺栓等表面凸起的转盘风阻扭矩特性,建立了一个转盘扭矩测量实验台.自由盘测量结果与经典经验关系式较好的符合度验证了测量方法的准确性.分析研究了在不同粗糙度以及附有6种凸起结构的转盘风阻扭矩规律.结果表明,壁面的粗糙会增大转盘受到的扭矩,在所测量旋转雷诺数范围内,从气动光滑一直到粗糙高度0.2 mm,扭矩系数增大1倍.基于气动光滑自由盘经验关系式进一步总结出相对粗糙度对于风阻扭矩的影响系数.对比6种凸起引起的风阻扭矩特性发现,叶型凸起会大幅降低风阻扭矩.凸起安装迎风角度的改变会带来10%的风阻扭矩差别.结合自由光滑盘扭矩系数经验关系式,总结出不同凸起形状对于扭矩系数的影响因数.