面向对象的通用航空发动机建模技术研究

针对目前发动机性能仿真中燃气热力性质计算方法适用范围和通用性差的问题,加入一种变成分的求解方法,并利用面向对象技术构建了通用的航空发动机性能仿真系统.建模中考虑了部件特性耦合、引气、冷却和功率提取等影响.同时在程序中加入了一些异常处理机制,以便于快速准确地进行调试,增强系统的可靠性.最后利用该仿真系统进行了某型双转子混排涡扇发动机的稳态特性计算,并与试验数据进行了比较.结果表明采用的变成分法进行发动机性能仿真的推力误差一般不超过3%,耗油率误差不超过2.5%,取得了良好的精度.尤其在高空小马赫数时变成分法的精度高于传统的拟合法.      

竖直平板Blasius流层流边界层流动与传热耦合解

对考虑辐射传热、流体黏度随温度变化、有和无滑移边界条件下的可渗透竖直平板Blasius流层流边界层的无量纲速度场与温度场进行了深入研究.经相似变换将描述速度场与温度场耦合的偏微分方程组转换成非线性常微分方程组,用Runge-Kutta法对常微分方程组进行了数值求解.研究了无量纲参数对无量纲速度场及温度场的影响,着重分析了滑移边界条件下速度和温度随无量纲参数的变化规律.结果表明:吸入时边界层变薄,喷注时边界层加厚;对比于无滑移边界条件,滑移边界条件下速度、温度边界层变薄;随着变黏度参数a或喷注与吸入参数的增大,壁面摩擦因数、局部努塞尔数Nu增大,速度和温度边界层变薄;随着普朗特数Pr、热辐射参数R的增加,或毕渥数Bi,布林克曼数Br的减小, 温度边界层变薄.      

采用火焰面/反应进度变量方法模拟湍流燃烧

在自主开发的软件平台上添加了火焰面/反应进度变量方法,选取了两种反应进度变量分别对甲烷/空气的同轴射流燃烧室进行了数值模拟.将计算结果与稳态火焰面模型、文献参考值以及实验值进行了对比研究,结果表明:采用火焰面/反应进度变量方法能捕捉到稳态火焰面模型所不能预测到的火焰抬举等非稳态现象,得到的结果与文献参考值以及实验值出入很小,能准确反应真实的燃烧过程;另外不同反应进度变量的定义对计算结果有较大的影响,采用火焰面/反应进度变量方法对主要中间产物一氧化碳的预测比文献参考值更好,更接近于实验值.      

风扇叶片丢失激励下转子-支承系统结构安全性设计策略

以高涵道比涡扇发动机风扇叶片丢失载荷激励下的转子-支承系统为研究对象,提出了一套结构安全性设计策略,即通过支承方案与载荷分配、变刚度支承结构和支承结构变形控制,结合转子结构动力学特性设计,实现转子结构安全性设计.研究表明:风扇后支点采用变刚度支承结构设计,能够在转子减速停车过程中减小风扇局部振动临界转速与相应振幅.通过设计滚珠轴承支承锥壳锥角,能够使支承具有较高的轴向承载能力并减小转轴变形对滚珠轴承的影响.轴承座底部与转轴间采用鼓形配合面连接设计,能够在大弯矩作用下通过配合面相互滑移避免支承随转轴变形,保证轴承安全.研究结果可为恶劣载荷作用下高涵道比涡扇发动机结构安全性设计提供依据.      

基于自校正Broyden拟牛顿法的航空发动机模型数值计算

以Broyden拟牛顿法为基础结合计算发散判断和校正机制,提出自校正Broyden拟牛顿法.该算法结合牛顿法(Newton-Raphson method)平方收敛和Broyden拟牛顿法超线性收敛特性,通过自适应调整计算步长和校正函数,在非线性系统中具有更好的计算性能.以变循环发动机部件级模型为对象,应用自校正Broyden拟牛顿法进行稳态及动态仿真计算,并与牛顿法和Broyden拟牛顿法作对比.结果表明:自校正Broyden拟牛顿法对恶劣的初始计算条件适应性更高,计算速度更快且收敛能力性更强,动态计算中部件模型计算调用次数为牛顿法的15%,模型动态误差低于Broyden拟牛顿法的15%,同时也低于牛顿法的28%,动态计算最大残差量低于其他两种算法的25%.验证结果表明了自校正Broyden拟牛顿法的优越性.      

考虑来流边界层的超燃燃烧室喷射特性

针对超燃燃烧室中的燃料掺混问题,采用基于雷诺平均Navier-Stokes的数值模拟方法分析了考虑来流边界层条件下的燃料横向射流流场特征及其掺混特性.研究发现:对于确定的来流边界层,燃料喷射存在一个临界动压比.当动压比低于该临界动压比时,增大来流边界层能明显提高燃料的穿透深度和掺混效率.而当动压比大于该临界动压比时,来流边界层厚度对燃料的穿透深度和掺混效率几乎没有影响.对于所研究的流动状态,该临界动压比约为0.900.在相同动压比下,所选厚来流边界层条件下的总压恢复系数仅约为薄来流边界层的0.93倍.其中,来流边界层内的摩擦损失是造成超燃燃烧室低总压恢复的主要因素,而改变来流边界层厚度对喷流及下游流场总压损失造成的影响相对较小.     

对B-P本构模拟材料率相关塑性变形和循环硬化的改进

对Bodner-Partom统一黏塑性本构理论（B—P模型）进行了修正，改善了其对材料加载速率相关的黏塑性和循环硬化的模拟精度。改进模型包含两个各向同性硬化内变量：一个用于描述小非弹性应变范围内快速的各向同性硬化。另一个以缓慢的演化速率来表述大的累积非弹性变形，如循环硬化过程中非弹性应变累积。利用ABAQUS的用户材料子程序将改进模型编成有限元计算代码，通过对镍基合金Udimet 720Li率相关黏塑性特征的本构模拟对其进行了验证，并与B-P模型作了对比。结果表明改进模型能够以更高精度、更加灵活地对材料的非弹性力学行为进行建模。     

结构非概率可靠性灵敏度分析方法

 结构非概率可靠性的灵敏度分析对工程中非概率问题的分析预测和优化具有重要的指导作用。基于结构非概率可靠性指标理论,提出了结构非概率可靠性的灵敏度分析方法。针对极限状态函数为线性的情况,推导出了非概率可靠性灵敏度的解析解。在此基础上,结合函数线性化理论,提出了求解极限状态函数为非线性情况的非概率可靠性灵敏度的近似解析方法。基于求解非概率可靠性指标的优化方法,并结合差分理论,提出了求解非概率可靠性灵敏度的优化方法。给出了算例分析,结果表明了所提方法的可行性。     

变厚度雷达罩电磁性能仿真三种方法的比较

对机载雷达罩提出一种分段变厚度的设计思想,实现变厚度机载雷达罩的建模。分别采用低频的多层快速多极子方法（MLFMM）、高频的物理光学法（PO）及高低频相结合的物理光学混合矩量法（PO/MoM）,利用这三种方法对一个中等尺寸机载雷达天线罩在分段变厚度情况下进行远场辐射方向图的计算,通过对三者所得结果的比较得出：PO/MoM混合方法在电大尺寸雷达罩电磁性能仿真计算中所具有的优越性及机载雷达罩分段变厚度方法的可操作性。     

高超声速二元进气道顶板移动变几何方案数值模拟

针对一种工作于马赫数Ma=4~6范围内的高超声速二元进气道,探索了一种进气道部分压缩顶板可移动的简单变几何方案,利用数值模拟研究了接力点变几何进气道的自起动性能和Ma=4~6的调节方法。结果表明:部分压缩顶板可移动简单变几何进气道方案在Ma=3.7能够实现自起动;变几何所形成的自适应放气槽放气量很小,最大放气量在2%以内,关闭自适应放气槽接力点Ma=4流量系数达到0.77;在整个工作范围内流量系数较高、总体性能较优,该变几何方案的调节方法是切实可用的。