航空发动机多变量系统状态反馈解耦

利用状态反馈实现了航空发动机的动态解耦。在解耦时不需要增加新的状态反馈即可实现解耦系统的任意极点配置。最后的算例表明了该方法的有效性。     

高速发动机油膜惯性对活塞裙润滑影响

利用一种新的迭代算法 ,分析了高速发动机油膜惯性对活塞裙润滑特性的影响 .该迭代方法是利用有限元法和差分法交替求解Navier Stocks方程和雷诺方程 ,据此导出了适用于高速发动机活塞裙的润滑计算的混和润滑雷诺模型 .新的模型借助惯性系数 ,引入了油膜惯性项 ;同时给出求解含有油膜惯性项的迭代步骤和有限元表达式。计算结果表明 :随着惯性系数和活塞裙的长径比的不同 ,油膜惯性会对油膜的摩擦力、压力和承载力产生不同的影响 ,这种影响对承载力尤为明显 .该润滑模型也可用于中、低速发动机的活塞裙润滑计算以及不计入惯性项 (惯性系数置为零 )的某些润滑问题求解 .     

基于控制系统工作原理的视线角速度获取方法

针对红外成像导引头弹目视线角速度获取精度较差,无法满足导弹制导指标要求的问题,根据红外导引头角跟踪原理,比较研究了直接微分法、基于马尔科夫模型的方法和基于控制系统工作原理的方法3种弹目视线角速度提取方案。利用试验数据进行了仿真分析。仿真结果表明:直接微分法和基于马尔科夫模型的方法均不适用于红外成像导引头弹目视线角速度提取。将控制系统设计为二阶无静差系统,并以此建立了基于控制系统工作原理的弹目视线角速度提取模型。仿真分析表明:该方法适用于红外导引头弹目视线角速度提取。     

基于DNA弹性细杆模型的Euler-Lagrange方程

采用了分析力学的方法,基于最小势能原理,DNA弹性曲杆的Euler-Lagrange方程组,运用弹性细杆模型拟合了A-,B-,Z-DNA的r-h曲线,讨论了弹性细杆模型描述3种DNA几何构型的可行性。与实验数据对比发现,通过计算,给定合适的曲率κ、挠率τ、扭转角χ,圆截面弹性细杆模型可以用来描述A-,B-DNA的几何构型;当椭圆截面的长宽比k=0.141时,采用椭圆截面弹性细杆模型拟合Z-DNA的r_0-h曲线与实验结果相一致。此研究内容有望为DNA分子构型的研究与探索提供参考。     

基于切齿高阶运动的螺旋锥齿轮齿面修缘方法

提出了一种利用机床柔性使用常规直刃刀具进行螺旋锥齿轮修缘的方法。根据不同刀具建立了螺旋锥齿轮修缘齿面的求解方程,并以此建立了螺旋锥齿轮啮合分析有限元模型分析对比了不同刀具修形对应修缘齿面的传动性能。根据对比结果选取圆弧刃刀具作为切齿高阶运动齿面修缘的优化目标,在齿面方程引入高阶运动系数并建立了基于高阶切齿运动的修缘齿面逼近优化计算模型,求解获得新的切齿调整参数得到近似修缘齿面。采用有限元分析和滚动检查实验对优化获得的修缘齿面和未修缘齿面进行了传动性能验证,其结果表明:切齿高阶运动逼近的修缘齿面啮合时的1阶振动幅值减小约16.6%,可以有效提高齿轮副的传动性能。      

一种关于静压气体轴承节流孔系数的计算方法

基于层流边界层方程的分离变量算法和雷诺方程的解析算法,提出了一种关于单节流孔静压气体止推轴承的节流孔系数的计算方法。该方法通过比较层流边界层方程计算获得的气体轴承的质量流量和雷诺方程计算获得的质量流量计算获得了节流孔系数。将计算获得的节流孔系数和节流孔系数为常数0.8代入单节流孔气体止推轴承的雷诺方程中,计算获得的承载力与分离变量算法求解层流边界层方程获得的承载力进行对比,可以发现,相对于采用节流孔系数为0.8来说, 采用该计算的节流孔系数求解雷诺方程的承载力与分离变量算法求解获得的承载力结果精度最大提高了8%。从而验证了该计算节流孔系数方法的正确性。      

Thermal state calculation of chamber in small thrust liquid rocket engine for steady state pulsed mode

This paper presents a method of thermal state calculation of combustion chamber in small thrust liquid rocket engine. The goal is to predict the thermal state of chamber wall by using basic parameters of engine: thrust level, propellants, chamber pressure, injection pattern, film cooling parameters, material of wall and their coating, etc. The difficulties in modeling the startup and shutdown processes of thrusters lie in the fact that there are the conjugated physical processes occurring at various parameters for non-design conditions. A mathematical model to predict the thermal state of the combustion chamber for different engine operation modes is developed. To simulate the startup and shutdown processes, a quasi-steady approach is applied by replacing the transient process with time-variant operating parameters of steady-state processes. The mathematical model is based on several principles and data commonly used for heat transfer modeling: geometry of flow part, gas dynamics of flow, thermodynamics of propellants and combustion spices, convective and radiation heat flows, conjugated heat transfer between hot gas and wall, and transient approach for calculation of thermal state of construction. Calculations of the thermal state of the combustion chamber in single-turn-on mode show good convergence with the experimental results. The results of pulsed modes indicate a large temperature gradient on the internal wall surface of the chamber between pulses and the thermal state of the wall strongly depends on the pulse duration and the interval.     

一种模拟不可压流的节点间断有限元格子波尔兹曼方法

基于节点间断有限元方法,本文提出了一种求解格子波尔兹曼方程(Lattice Boltzmann equation,LBE)的新方法,即节点间断有限元-格子波尔兹曼方法(Nodal discontinuous Galerkin-Lattice Boltzmann method,NDGLBM)。在该方法中,LBE的碰撞过程和迁移过程被拆分成了两步:碰撞过程用LBM多松弛时间(Multiple relaxation time,MRT)模型进行求解,迁移过程则写成对流方程并采用节点间断有限元方法进行求解。其中,空间离散采用了非结构网格,时间离散采用了四阶、五步龙格-库塔格式。为了验证NDG-LBM的可行性,文中模拟了顶盖驱动方腔流、静止单个圆柱绕流、旋转-静止双圆柱绕流以及高雷诺数NACA0012翼型绕流。计算所得的数值结果与其他文献中的结果吻合度很好。     

喉道对压气机超声叶栅流态及性能的影响

为更深入认识超声叶栅流动机理,以ARL-SL19、CM-1.2和SM-1.5叶栅为研究对象,采用数值模拟和理论分析相结合的方式开展喉道对超声叶栅激波结构和性能影响的研究。研究结果表明:超声叶栅存在两种稳定工作状态,起动状态和溢流状态;在来流马赫数较高时,叶栅只工作于起动状态;在来流马赫数较低时,叶栅只工作于溢流状态;存在一个马赫数区间,叶栅的工作状态由前一个状态决定;对于低马赫数C形超声叶栅,高压比下气动喉道起决定因素;对于高马赫数S形超声叶栅,真实喉道起决定因素;若为气动喉道导致溢流,溢流实现更大的裕度和更低的损失,进口马赫数和气流角会受压比影响;若为真实喉道引进的溢流,溢流会降低裕度并增加损失,叶栅保持唯一进气角流动,但进口气流角和马赫数与起动状态不同。     

基于后缘小翼的旋翼翼型动态失速控制分析

针对后缘小翼(TEF)的典型运动参数对旋翼翼型动态失速特性的控制进行了研究。发展了一套适用于带有后缘小翼控制的旋翼翼型非定常流动特性模拟的高效、高精度CFD方法。通过求解Poisson方程生成围绕旋翼翼型的黏性贴体和正交网格,为保证后缘小翼附近的网格生成质量,建立了基于翼型点重构的方法来描述后缘小翼的偏转运动;为克服变形网格方法可能导致网格畸变的不足,发展了一套适用于带有后缘小翼控制的旋翼翼型运动嵌套网格方法。基于非定常雷诺平均Navier-Stokes(URANS)方程、双时间法、Spalart-Allmaras(S-A)湍流模型和Roe-Monotone Upwind-centered Scheme for Conservation Laws(Roe-MUSCL)插值格式,发展了旋翼翼型非定常气动特性分析的高精度数值方法,并采用Lower-Upper Symmetric Gauss-Seidel(LU-SGS)隐式时间推进方法及并行技术提高计算效率。以有试验结果验证的HH-02翼型和SC1095翼型为算例,精确捕捉了动态失速状态下的气动力迟滞效应,验证了本文方法的有效性。着重针对SC1095旋翼翼型的动态失速状态开展后缘小翼的控制分析,提出了可以体现翼型升力、阻力及力矩综合特性的关系式Po和Pc,揭示了后缘小翼振荡频率、相位差和偏转幅值对动态失速特性影响的规律。研究结果表明:当后缘小翼偏转的相对运动频率为1.0,且小翼运动规律与翼型振荡规律之间的相位差为0°时,后缘小翼能够更好地抑制翼型动态失速现象;在此状态下,当偏转幅值为10°时,SC1095翼型最大阻力系数和最大力矩系数可以分别降低19%和27%。