突加不平衡下GTF发动机转子系统瞬态响应分析

针对GTF发动机低压转子系统突加不平衡瞬态响应问题,基于有限元法建立了考虑啮合单元、行星架单元、膜盘联轴器单元的多体接触、多转子耦合动力学模型,介绍了考虑非线性因素下突加不平衡瞬态响应计算方法,求解了低压转子系统的突加不平衡瞬态响应,并分析了关键部件刚度对突加不平衡的影响规律。结果表明：当风扇转子发生突加不平衡时,各部件均表现为振动幅值突增,后迅速趋于稳定,低压涡轮转子发生明显拍振,突加不平衡载荷主要由行星架结构和第1.5支点刚性支承承担;行星架刚度主要影响风扇转子突加不平衡响应、第1.5支点和行星架处的外传力;膜盘联轴器刚度主要影响增压级突加不平衡响应、转子在突加不平衡之后稳定运转情况。     

守恒型可压缩一维湍流方法及其在超声速标量混合层中的应用

一维湍流（ODT）方法是一种能在一维计算域上遵循湍流基本物理规律的湍流建模方法。通过结合确定性和随机性求解方法,能够在一维计算域上准确捕捉到湍流统计规律,且降维建模可显著减小计算量。ODT方法主要被广泛用于不可压湍流和湍流燃烧研究,若要将其拓展用于模拟高速可压缩湍流,需对建模方法进行深度改进。相比于不可压ODT方法,本文基于欧拉参考框架,针对可压缩湍流的特性,将因变量由原始变量改为有利于减小可压缩湍流模拟误差的守恒通量,并加入了组分求解模块。对确定性和随机性求解模块均进行了相应的深度改进,开发出具有标量混合模拟功能的守恒型可压缩ODT方法。在确定性模块中改为求解以守恒通量为变量的一维截断控制方程,在随机性模块中构造一维涡时,将三联映射的作用对象也相应地由原始变量改为守恒通量,并选用了可保证变密度情况下动量守恒的双核变换。通过模拟空间发展超声速平面湍流混合层并将自相似阶段结果与实验结果比对,验证该方法对可压缩剪切湍流场中标量混合的捕捉精度。守恒型可压缩ODT方法模拟得到的速度场和组分场的平均剖面和脉动强度分布与实验结果准确吻合,精度明显优于传统的耦合梯度扩散亚格子模型的大涡模拟方法（LES-GRAD.DIFF.）以及耦合线性涡（LEM）亚格子模型的大涡模拟方法（LES-LEM）,且该方法的降维处理使其在降低计算成本方面具有显著优势。     

“汽电双驱”系统电机瞬态性能仿真研究

以上海电机厂有限公司“汽电双驱”项目中三相异步电动机/发电机为依据,基于以Ansoft软件,利用场路耦合法对系统中异步电机的瞬态性能进行分析。该电机为6 800 kW,750 r/min。针对其3种不同的工况,即空载起动、风机负载起动、大负载起动,分别进行分析计算,模拟系统中电机由电动机状态切换成发电机状态的过程,获得该三相异步电机在不同状态下的瞬时起动电流、电磁转矩、气隙磁密等电气参数,进而分析该异步电机的瞬态性能和相应的电磁场分布。     

融合柔性基础传递函数的转子系统建模方法及验证

航空发动机、火箭涡轮泵等现代高端涡轮机械对性能有着极致的需求,其轻柔化的支承结构导致转子与柔性基础振动耦合特性明显,必须将其柔性支承结构纳入整体动力学分析之中。然而,柔性基础的高精度建模过程需要耗费大量时间,同时简化的基础模型也难以表达其真实的动力学特性。为此,提出了融合柔性基础传递函数的转子系统动力学建模仿真方法：将实测得到的柔性基础传递函数通过状态子空间法拟合获得低维时域和频域表征的数学模型,通过支承力模型将柔性基础模型耦合到转子有限元模型之中,最终形成基于实物传递函数的混合转子动力学模型。同时,提出了柔性基础从稳态到瞬态的模型转换策略,并采用线性-非线性节点（显-隐）分离的快速瞬态数值积分方法进行求解。基于所提模型和计算方法,开展了柔性基础-转子系统的稳态、瞬态动力学数值仿真,以及稳定转速扫频激励和降速不平衡激励的振动测试试验,结果表明：实测的柔性基础振动特性,能够有效地体现在整机动力学分析之中,且考虑柔性基础特性的转子动力学模型所预测的瞬态动力学响应与试验结果更为吻合。所提方法为含复杂柔性基础的转子系统提供了行之有效的整机建模及瞬态仿真方法。     

本征正交分解法对正十二烷喷雾射流的解析

应用大涡模拟方法对航空煤油单组分替代燃料正十二烷冷态喷雾现象进行模拟计算,探究喷雾射流中的多物理场耦合现象。从宏观层面分析燃油破碎蒸发及油气混合过程,发现流场涡团在其中起着重要的作用。运用本征正交分解方法对正十二烷射流质量分数场和涡量场数据库进行加工,从“微观”层面深入探究两者的联系。结果表明:正十二烷质量分数场喷雾结构与射流涡量场涡团结构存在显著的空间分布相关性,反应出两者结构演变具有重要关联特征。湍流脉动致使涡团尺度发生变化,从而影响燃油破碎蒸发,最终使得喷雾不同尺度结构分布形态发生改变。      

基于放大因子与Spalart-Allmaras湍流模型的转捩预测

为了验证放大因子输运方程与Spalart-Allmaras(S-A)湍流模型耦合对转捩现象的模拟精度,选取Schubauer and Klebanoff(S-K)平板、S809低速翼型、30p30n多段翼型以及复杂的三维HiLiftPW-1构型进行自由转捩计算,并将计算结果与实验进行比较分析,其中针对S809算例,还与Langtry-Menter(L-M)转捩模型进行了比较.算例结果表明:放大因子输运方程与S-A湍流模型的耦合能够较好的捕捉转捩位置以及转捩发展过程,对分离泡诱导的转捩的模拟相比L-M转捩模型更精确,转捩位置的捕捉精度提升了10%;对比实验,多段翼转捩位置的捕捉误差最大为6.5%;针对三维高升力增升构型,以实验作为参考,全湍流计算与考虑边界层转捩的对比显示考虑边界层转捩能够更加精确的模拟气动力系数,升力和表面摩擦阻力系数的模拟精度精度提升1%.      

基于雷诺平均Navier-Stokes方程的表面传热系数计算

采用有限体积法数值求解控制二维绕流的雷诺平均Navier-Stokes(RANS)方程组,计算了光滑和粗糙NACA0012翼型以及圆柱表面的局部表面传热系数.分析了近壁面网格间距、湍流模式和表面粗糙度模型对数值计算结果的影响.结果表明:切应力输运(SST)湍流模型能够区分层流和湍流边界层的对流传热特性,并能预测转捩的发生;采用Spalart-Allmaras(S-A) 扩展模型能够计算粗糙壁面的对流传热系数,但采用忽略转捩函数的S-A模型不能有效计算层流边界层的传热系数.当近壁面网格间距接近10-5量级的黏性子层时,在光滑和粗糙壁面都能得到准确的传热系数分布.结合合适的近壁面网格间距,湍流模式和表面粗糙度模型可以得到与实验数据十分接近的表面传热系数曲线.通过与求解不可压缩RANS方程得到的结果比较后发现,不可压缩RANS方程主要忽略了压缩和黏性耗散效应,这种效应可以通过绝热升温项的形式并入总体热分析.      

基于LPV状态空间模型的小型航空煤油二冲程发动机喷油控制策略

为了简化基于平均值模型的喷油控制模型的计算,根据线性变参数(LPV)状态空间模型原理,在平均值模型的基础上建立了LPV状态空间喷油控制模型,利用Matlab/Simulink工具对稳态工况和瞬态工况的喷油脉宽控制模型进行仿真,最后在发动机试验台架上进行试验验证.仿真和试验结果表明:①基于LPV状态空间模型的喷油控制模型能够满足瞬态工况和稳态工况控制的要求,计算的喷油脉宽精度略低于平均值模型,由于LPV状态空间模型计算简单,能够直接应用于控制算法的设计,因此在工程中容易实现.②与插值算法相比,基于LPV状态空间模型的控制策略使发动机的动力输出有所下降,主要是由于模型简化造成进气空气流量减小引起的.      

基于涡方法生成大涡模拟进口条件的数值计算

探索改进涡方法来生成大涡模拟的非定常进口条件.改进方法中,为了避免局部漩涡数量过多或者过少,采用密度分布方式放置漩涡场.并用Langevin方程控制漩涡运动,模拟实验方法中的蜂窝器,使改进后的涡方法生成的脉动速度场更加符合湍流的特征.在已知雷诺平均的流场结果下,利用涡方法产生漩涡场,进而生成能满足大涡模拟所需要的非定常进口流场.为了检验改进后的涡方法在生成脉动速度过程中的效果,在槽道中进行了对比数值试验,并借助于直接数值模拟数据做对比,对比分析槽道进出口的平均速度、涡量以及雷诺应力统计,证明改进后的涡方法在生成大涡模拟进口条件下是非常有效的.