航空发动机高压转子“可容模态”设计及实验验证

为降低航空发动机高压转子的振动，提高高压转子对复杂工况的可容度，本文建立带双阻尼器的高压转子动力学模型，开展模态分析和响应分析，研究高压转子关键结构参数对“可容模态”的影响规律。构造了可容度评价函数η，建立了“可容模态”设计方法。设计了高压转子模拟实验器，通过模态校核实验验证了计算模型和计算方法的准确性，通过“可容模态”实验验证了设计方法的可行性。研究发现，利用当量临界转速ωˉ估计转子系统前两阶临界转速的方法可行，高压转子实验器的临界转速满足ωcr1≤ωˉ，ωˉ≤ωcr2≤2ωˉ；高压转子实验器在两阶“可容模态”下连续“共振”时长达到7125s，实验器转子关键部位的振动位移峰值不超过79.42μm，振动速度有效值不超过2.33mm/s，最大波动不超过10.59%。结果表明，所建立的发动机高压转子“可容模态”设计方法是行之有效的。     

离心力场对内燃波转子火焰传播特性的影响

为了揭示离心力对内燃波转子火焰传播特性的影响，建立了内燃波转子燃烧过程的简化模型，数值模拟了不同离心力场强度下波转子通道内的燃烧过程。研究结果表明，离心力只在一定范围内对火焰传播起到加速作用，且这种影响与预混气当量油气比Φ有关，当Φ=1.0时，作用范围为［0～250g］，在Φ=0.6，0.8时，作用范围为［0～700g］；离心力场对内燃波转子火焰传播特性的影响机制是导致火焰锋面变形，进而引起火焰传播速度变化。     

基于改进联合稀疏电容层析成像算法滑油监测研究

由于电容层析成像（Electrical Capacitance Tomography，ECT）的“软场”特性，使其应用于航空发动机滑油管道监测时难以准确获得被测对象的位置和大小信息。针对图像重建过程中的病态问题，在对称交替乘子法（S-ADMM）算法基础上，提出一种采用L2,1-范数作为核范数、(L1-εL2)范数作为正则化器的改进联合稀疏电容层析成像算法。仿真实验采用COMSOL5.3和MATLAB2014a搭建仿真实验平台，仿真结果表明，改进算法与Landweber迭代算法相比，成像误差降低了35.82%，相关系数提高了56.51%；与S-ADMM算法相比成像误差降低了35.02%，相关系数提高了41.82%，图像误差降低至0.1435、相关系数提升至0.9308，且成像时间保持在0.1s。实验结果表明，改进联合稀疏成像算法提高了滑油监测方面的成像质量，具有实效性和适用性。     

放电室长径比对射频离子推力器影响的数值模拟研究

为获得推力器的性能与推力器工况之间的变化规律，利用等离子体密度平衡方程、中性气体密度平衡方程、中性气体功率平衡方程、电子功率平衡方程构成的全域模型对氙工质射频离子推力器开展了数值模拟研究。结果表明，放电室长径比（L/R）的改变引起放电室壁上消耗功率的变化，从而影响推力器的性能。在放电室体积和工质流率不变的情况下，放电室的开口面积越大，越有利于射频功率的转换和推力器性能的提升。此外，改变工质流率导致离子间碰撞平均自由程和碰撞频率发生变化，进而使推力器的射频功率转换效率、推力效率、推功比随工质流率的增加而增加，工质利用率随着工质流率的增加而减小。     

固体姿轨控发动机喉栓运动时序实验设计方法

固体姿轨控发动机通过伺服机构控制喉栓位移实现推力矢量调节。为实现推力精确调控，需拟合大量地面试验数据，得出发动机喉栓位移与推力矢量间的数学模型。提出了一种面向喉栓式固体轨控发动机喉栓运动时序实验设计方法，实现了单次试验中不同燃烧室理论压强及喉栓位移组合最大覆盖，为最低试验成本下精准构建喉栓位移与推力输出模型提供训练样本。以四阀固体轨控发动机作为对象，基于零维内弹道实现燃气发生器压强求解。采用Φp准则评估喉栓位置组合遍历性，并对其进行优化，获得合理的喉栓运动方案。算例结果证明实验设计后方案在整个运动周期内几乎遍历各种组合状态，较实验设计前Φp值提升近40%，验证了本文实验设计方法在时序工程问题中的优良性能。采用实验设计方法可提高固体姿轨控发动机地面试车数据利用效率。     

一种民航发动机滚动轴承故障预警方法

为了提前感知滚动轴承故障，避免民用航空发动机非计划维护带来的损失，提出了一种故障预警方法。对轴承振动信号进行特征工程，提取其时域和频域特征，引入梯度提升决策树（Gradient boosting decision tree，GBDT）算法，量化了特征重要度；在特征相关性分析的基础上，利用核主成分分析（Kernel principal component analysis，KPCA）方法实现特征融合与主元提取，再次结合GBDT构建了故障预警模型，使用交叉验证法实现了模型泛化能力评估。结果表明：KPCA+GBDT模型的泛化性能显著，模型的F1分数高达0.991，对应的受试者工作特性（Receiver operating characteristic，ROC）曲线下面积的值为0.998，体现出该方法用于支撑航空发动机健康管理与维护决策工作的工程应用价值。     

基于叶端定时的转子叶片裂纹监测与诊断

为了实现转子叶片裂纹非接触式监测，提出基于??0范数的叶端定时欠采样信号稀疏重构算法。针对叶端定时信号欠采样问题，建立欠采样信号稀疏表示模型，利用分块正交匹配追踪算法对叶端定时欠采样信号进行分段重构，进而监测在变转速非稳态工况下转子叶片动频的变化情况，据此判断转子叶片是否产生裂纹。针对应变片信号间隔采样的特点，研究基于同步提取变换的时频分析方法，以提高时频图的时频聚集性。开展旋转叶片裂纹扩展试验，同时采用叶端定时和应变片系统测量叶片振动。对比叶端定时和应变片的分析结果，二者均可在裂纹叶片断裂之前至少20min，监测到裂纹叶片和正常叶片的动频发生分离的现象，即裂纹萌生致使叶片动频发生偏移，表明所提出的叶端定时欠采样信号重构方法能够实现转子叶片裂纹的监测与诊断。     

电帆推进在等离子体环境下推力影响因素分析和试验研究

为了分析电帆推力的影响因素，根据电帆推进原理设计了电帆地面模拟装置。首先分析了真空中带电导线和等离子作用机理，得出了电帆推力公式。利用朗缪尔探针测量了电帆模拟试验环境中等离子体参数，计算产生推力的大小，最后通过微推力测量系统测量电帆产生的推力，与计算结果进行了对比。根据推力公式和仿真结果表明，提高电压、增加带电粒子密度和速度、增加电帆长度均可增加推力。在试验中电帆上电压从300V增加到800V时，推力从9μN增加到26μN。仿真和试验结果表明，电压、带电粒子密度和速度、电帆导线长度都对推力有影响。根据试验结果，电帆推力公式的经验推力系数在1～2取值。     

流量比对高位预旋进气涡轮转静盘腔换热影响的数值研究

为提高冷气与涡轮盘的换热效果，针对高位预旋30°进气的转静腔耦合涡轮盘系统进行了数值模拟，研究了不同转静腔入口流量比率和转速对系统换热效果的影响，并对涡轮盘壁面平均努塞尔数Nuˉ、盘体最高温度、盘体温度均匀性进行了分析。结果表明：转速为1×10⁴r/min时，随着转静腔1和2的入口流量比率P从1增加到6，左壁面Nuˉ稍有上升，右壁面Nuˉ先增后减，在P=3时最大；盘体最高温度降低，且P<3时降幅较大，之后较小。盘体温度均匀性受P影响较小，P>2后，左右两侧径向温度均匀性几乎一致，轴向温度均匀性均优于径向，盘缘处轴向温度均匀性最差，越靠近轴心均匀性越好。在转速为6×10³~1.2×10⁴r/min时，随转速增大，系统换热效果增强。当转速<1×10⁴r/min时，随P增大，冷气与涡轮盘之间换热不断改善；转速>1×10⁴r/min后，P在3附近综合换热效果较好。     

基于光滑粒子流体动力学方法的液滴蒸发问题数值模拟研究

基于光滑粒子流体动力学方法（Smoothed Particle Hydrodynamics，SPH），开展了SPH新算法在蒸发燃烧领域的研究。建立了适用于SPH方法的蒸发数值模型，推导了基于傅立叶热传导公式和菲克扩散定律的SPH离散方程；借鉴VOF方法（Volume of Fluid）的思想，提出了SPH粒子的液相质量分数的概念，以有效表征蒸发过程中的相变问题。采用SPH方法对高温环境中单个液滴的蒸发过程进行数值模拟，结果符合D2定律，与理论模型相一致；在强迫对流环境中，液滴的蒸发过程受到对流作用及表面张力的影响，蒸发速率加快；进一步对双液滴在静止、对流环境中的蒸发过程进行数值模拟研究。结果表明，液滴的间距、滴径对多个液滴的蒸发过程影响至关重要，液滴间距至少在两倍的液滴直径以上，相互之间的影响才可以近似忽略。通过本文研究，拓宽了SPH方法在蒸发相变领域的应用范围，研究结果也能够为进一步的燃烧问题研究奠定基础。     

带有闭式布雷顿循环的预冷发动机特性研究

为获取带有闭式布雷顿循环的预冷发动机的飞行包线及性能,同时为提高发动机工程实现可行性,本文基于带有闭式布雷顿循环的预冷发动机基础循环及现有部件技术水平,构建了一种适度预冷发动机方案。对该方案下发动机沿着SABRE3飞行轨迹下的性能和部件匹配规律进行了分析。然后通过对发动机的高度、速度、调节特性进行研究,得到了该方案下发动机的飞行包线及整个包线内的性能。计算结果表明,本文所提出的适度预冷方案与SABRE3方案相比,核心机的比冲基本相当,但单位推力有所降低,工程可实现性提高;通过分别控制氦循环最低、最高温度为目标值,可保证发动机各部件在马赫数0~5的整个飞行过程中均处于稳定工作区间内,发动机比冲在1359 s~2099 s之间,地面点单位推力最大,达到1.9 kN/(kg/s);特性研究发现发动机推力与比冲在高度0~15 km、马赫数1~3之间最高,而单位推力最高的区域主要集中在包线的左侧低马赫数区,随马赫数的增加逐渐降低;发动机对氦压气机前温度的调节十分敏感,而对氦涡轮前温度的调节敏感性较低。综合研究表明,本文所给出的适度预冷方案的预冷发动机具有较好的宽域工作能力。     

旁路放气循环高速单轴涡喷发动机安装性能模拟

针对Ma3.5旁路放气循环单轴涡喷发动机,提出一种压气机旁路放气计算方法,建立基于进/排气系统特性数据库的涡轮发动机安装性能计算模型,分析压气机旁路放气对压气机共同工作线和发动机高空高速推力性能的影响,给出压气机旁路放气量的调节原则,计算发动机安装性能、进/排气安装阻力沿飞行轨迹的变化规律。计算结果表明:当飞行马赫数大于2.3后需打开压气机旁路放气,旁路放气阀门面积和放气流量均随着飞行马赫数的增大基本呈线性增大趋势;通过压气机旁路放气,可显著改善单轴涡喷发动机在高空高速飞行条件下的稳定性和安装推力性能,在飞行马赫数3.0附近,可实现安装推力提高30%以上;在跨声速至飞行马赫数2.0区间内,推力安装损失最大,约为非安装推力的25%～30%。     

飞行状态下脉冲爆震发动机性能估计分析

在理想单循环脉冲爆震发动机分析模型的基础上,通过引入进气道建立了一定飞行状态下的自吸气的脉冲爆震发动机理论性能计算模型,研究了当飞行状态改变时,脉冲爆震发动机的性能变化趋势,计算结果表明,对给定模型,飞行高度和飞行马赫数对脉冲爆震发动机性能的影响是不一样的,决定于进气道和爆震室的推进特性,同时模型也考虑了爆震频率和爆震室直径对性能的影响,这为脉冲爆震发动机的理论设计和分析提供了一定依据。      

旋转爆震冲压发动机总体性能分析

为了快速可靠地评估旋转爆震冲压发动机的总体性能,针对冲压模态下的旋转爆震发动机建立了性能分析模型。模型以飞行条件和冲压发动机关键几何参数作为输入参数,结合气体动力学和C-J爆震理论,获得旋转爆震燃烧室的流场参数分布以及发动机喷管排气参数,输出发动机推力以及燃料比冲,建立了基于连续旋转爆震的冲压发动机性能评估方法。模型参与反应的燃料和氧化剂分别为煤油以及空气,主要研究了燃料温度、喷管喉部面积、燃烧室环面面积、反应物当量比、飞行马赫数以及飞行高度对发动机燃料比冲、推力的影响趋势。研究结果表明,控制其它变量不变,发动机推力与燃料比冲随燃料温度上升而提高;随喷管喉部面积、燃烧室环面面积减小而增大;随飞行高度增加而降低;燃料比冲随当量比、马赫数增大而减小,而推力随当量比、马赫数增大而增大。在高度为25 km、马赫数为4、当量比为0.6的工况下,发动机燃料比冲可达到1 740 s。分析结果表明,模型计算方法可靠,可快速计算出旋转爆震冲压发动机的推力性能,为旋转爆震冲压发动机的设计提供可靠参考。     

脉冲爆震发动机的飞行性能计算与分析

为了获得脉冲爆震发动机的飞行性能参数随飞行高度和飞行速度的变化规律,在应用最小G ibbs自由能理论建立求解爆震燃烧的平衡成分及其浓度、平衡参数的基本关系式,据此对爆震波特性参数进行数值模拟的基础上,通过引入进气道和尾喷管建立了一定飞行状态下的多管自吸气式脉冲爆震发动机理论性能计算模型,对脉冲爆震发动机的性能随飞行高度、速度的变化趋势进行研究。结果表明,在不同飞行高度和飞行马赫数的状态下,脉冲爆震发动机性能有一定的变化规律,同时模型也考虑了加入尾喷管对性能的影响,为脉冲爆震发动机的性能分析与理论设计提供了依据。     

换热预冷发动机预冷特性和发动机性能数值研究

为深入了解换热预冷发动机换热器预冷效果和发动机性能，首先设计了叉排管束式细小通道换热预冷器，并采用非定常数值模拟方法对换热效果进行了仿真，然后结合换热预冷发动机性能计算程序对预冷发动机性能进行了研究。结果表明，飞行马赫数为2.5~4.0且氢气空气质量流比为0.03~0.09条件下，换热预冷器能将来流空气预冷90.6~471.2 K，低温氢气经吸热后温度升幅为266.1~455.3 K，换热效果良好。来流空气经预冷后涡轮发动机的飞行包线最高被拓展至马赫数4.0，达到了与超燃冲压发动机的衔接速域。相比于传统涡轮发动机，氢气空气质量流比为0.03时，本文换热预冷发动机加力状态推力能恢复至设计点推力水平；当氢气空气质量流比增加至0.09时，加力状态推力最高达到设计点推力的两倍左右。马赫数2.6以下换热预冷措施能小幅改善发动机比冲和耗油率（仅计算用于燃烧的氢气量）性能，而飞行速度大于马赫数2.6后换热预冷措施也难以抑制比冲和耗油率迅速恶化的趋势。     

SABRE预冷器计算模型及其在整机模型中的应用

预冷器的性能对SABRE等预冷组合循环发动机具有重要影响,为实现发动机方案设计阶段预冷器的快速设计与评估,建立了预冷器准二维快速评估模型。将SABRE预冷器的几何结构简化为一个扇环形区域,沿径向和周向将该区域划分为二维节点。应用守恒方程及传热关联式完成单个节点计算,再求解节点矩阵的平衡方程组,计算内外流体特定节点上参数的二维分布,得到预冷器出口参数。将模型嵌入发动机总体性能程序中,实现了发动机设计及非设计点的预冷器性能计算的功能。与文献数据对比结果表明,预冷器模型传热计算误差小于5%,摩擦阻力误差小于10%。整机计算结果显示,Ma0~5范围内,预冷器空气侧温降范围为143K至932K,温降随飞行马赫数升高单调上升。预冷器传热有效度范围为0.896-0.945,空气侧总压恢复系数范围为0.852-0.904。     

某民用涡扇发动机飞行包线内吸雨量计算分析

分析吸雨对发动机工作性能的影响,首要的是确定能被发动机吸入的雨量。通过综合考虑适航规章要求的大气雨水分布、环境温度、飞行速度、发动机功率和聚集效应对发动机吸雨的影响,结合完整的推力调节计划,计算分析了某民用大涵道比涡扇发动机在不同推力等级下全飞行包线内风扇进口水气比(WAR)。结果表明:在同一马赫数下,风扇进口水气比在6 100m处达到最大;在该高度以下,水气比随着高度的增加而增大;在该高度以上,水气比随着高度的增加而减小;飞行马赫数越大、环境温度越高,水气比也越大;并且标准天空中慢车推力下的风扇进口水气比最大可达到7.38%。该计算方法可为民用涡扇发动机吸雨适航取证的关键点分析提供参考。     

吸气式激光推进飞行器发动机特性分析

介绍了吸气式脉冲激光爆轰推进的环聚焦构形飞行器概念,建立了进气道流动模型;利用Chapman-Jouguet爆轰理论和Sedov自相似律,建立了冲量发生与推力传递模型,给出了冲量耦合系数及等效比冲与飞行参数、激光脉冲参数之间关系的解析表达式。计算了发动机冲量耦合系数及等效比冲随飞行高度、马赫数、入射激光强度和脉宽等因素变化的曲线,结果可为相关构形飞行器的弹道与激光器工作模式设计提供参考。     

空天机复合循环推进技术新概念研究——液化空气循环在其中的应用

研究了空天机复合循环推进技术中的一个新概念——以液化空气循环为基础的SERJ组合发动机性能。计算了空气液化后分离为液氧和液氮后组合发动机在不同高度、不同飞行马赫数下的性能曲线,并考虑了液化空气量对组合发动机性能的影响。对计算结果的分析表明,液化空气循环虽然使组合发动机推力下降,但能使比冲有较大幅度的上升;液化空气量过大与过小均不利于进一步提高比冲性能,而是存在一个最佳值。