航空发动机非线性模型实时计算的迭代方法研究

针对航空发动机非线性模型将作为未来机载应用而需要严格实时性的问题,研究了求解模型的迭代方法,综述其历史发展情况,并讨论收敛性问题,同时提出一种改进应用方案,即模型在过渡态仿真时可将上一状态点雅可比矩阵直接应用于下一点计算,能够在典型应用方案基础上直接减少气动热力模型计算次数,从而提高整体模型执行效率.通过实际算例对这些应用方案进行仿真比较与分析,结果表明:①Broy-den秩1法执行效率最高,更适宜用于模型实时计算;②改进应用方案简单、有效,为进一步满足非线性模型的严格实时性要求奠定基础.      

Effect of interface debonding on matrix multicracking evolution of fiber-reinfroced ceramic-matrix composites

An analytical methodology was developed to investigate the effect of fiber/matrix interface debonding on matrix multicracking evolution of fiber-reinforced CMCs(ceramic-matrix composites). The Budiansky-Hutchinson-Evans shear-lag model was adopted to analyse the micro-stress field of the damaged composites. The critical matrix strain energy criterion, which presupposes the existence of an ultimate or critical matrix strain energy with matrix, was obtained to simulate the matrix multicracking evolution of CMCs. With the increase of the applied stress, the matrix multicracking and fiber/matrix interface debonding occurred to dissipate the additional energy entered into the composites. The fiber/matrix interface debonded length under matrix multicracking evolution was obtained by treating the interface debonding as a particular crack propagation problem. The conditions for no-debonding and debonding during the evolution of matrix multicracking were discussed in terms of two interfacial properties, i.e., the interface shear stress and interface debonded toughness. When the fiber/matrix interface was bonded, the matrix multicracking evolution was much more intense compared with the interface debonding; when the fiber/matrix interface was debonded, the matrix crack density increased with the increasing of interface shear stress and interface debonded energy. The theoretical results were compared with experimental data of unidirectional SiC/CAS(calcium alumina silicate), SiC/CAS-Ⅱ and SiC/borosilicate composites.      

高声强下多狭缝共振腔的吸声性能

为深入理解高入射声强下多狭缝共振腔吸声机理,采用低频散低耗散的计算气动声学方法对二维多狭缝共振腔开展直接数值模拟研究.首先对标准单狭缝共振腔计算结果进行验证,随后相同的数值模拟方法被应用于相同穿孔率的多狭缝共振腔的数值模拟中.结果显示:高声强下涡脱落对吸声系数的贡献占据了主导地位,各入射频率下均超过68%.多狭缝低频时会导致脱落涡总能量的下降,而高频时升高;而黏性耗散作用随着共振腔狭缝数目的增加而增强.因此综合作用下低频时多狭缝共振腔对吸声效果影响不大,但在高频时多狭缝共振腔有更好的吸声效果.      

基于多项式平方和规划的航空发动机鲁棒LPV/PI控制

针对航空发动机常规(proportion integration,PI)控制器设计过程中难以保证鲁棒性及参数适应性差等问题,提出了一种基于线性变参数(linear parameter varying,LPV)模型及多项式平方和(sum of squares,SOS)规划的控制器设计方法.结合传递函数模型下的鲁棒稳定条件及弱对偶定理给出了多项式描述的LPV模型鲁棒稳定条件,并转化为便于求解的SOS规划问题.根据发动机非线性模型获取不同转速下的传递函数模型,并利用多项式拟合的方法建立发动机LPV模型.根据所提出的定理构造出SOS规划问题,并求解得出LPV/PI控制器.最终以某型双轴涡扇发动机为被控对象,在包线内不同点进行了阶跃仿真,结果表明:高压转子转速控制系统的稳态误差为0,调节时间小于3s.      

X射线脉冲星导航硬件脉冲轮廓累积研究

为了在X射线脉冲星地面实验系统仿真源模拟产生X射线的基础上,能够快速稳定地得到脉冲轮廓,采用硬件历元叠加的方法获得脉冲轮廓。研究了用硬件实现历元叠加及其数据整合的算法,该算法首先在MATLAB现场可编程逻辑阵列(FPGA)中实现,再通过MATLAB硬件描述语言(HDL)代码生成模块把算法转换成HDL,经编译后获得配置硬件的Bit文件,最终在开发板FPGA上实现数据处理的硬件模块。一段时间内的光子到达时间数据通过MATLAB算法得到的脉冲轮廓数据与通过硬件模块处理后得到的数据结果存在误差,在单个时间窗口内误差最大值为2个光子数,误差平均值占光子数统计平均值的0.084%;两组统计的脉冲轮廓数据中不同数据占总数据个数的9.481%,这样的误差不影响后端模拟导航模块的导航。利用硬件实现的历元叠加及其数据整合模块具有处理速度快、设备紧凑、功耗低的特点,为航天器利用X射线脉冲星导航提供了一种可行的硬件数据处理技术上的支持。     

存在观测站位置误差的转发式时差无源定位

无源定位中,由于观测站安放在运动平台等原因造成的观测站位置误差会影响无源定位精度性能。另外到达时间差(简称时差)(TDOA)的转发式测量需要将不同观测站截获到的辐射源信号都转发到同一位置,如主观测站。针对这两个问题,提出了基于约束总体最小二乘(CTLS)的无源定位算法。首先将转发式时差的非线性定位方程转化为不需要中间变量的直接线性方程,再基于CTLS算法依次转化为约束优化问题和无约束优化问题,最后推导给出定位近似闭式解。仿真实验表明在观测站误差较大时,该算法与其他算法相比定位精度性能较好。     

Gust response analysis and wind tunnel test for a high-aspect ratio wing

A theoretical nonlinear aeroelastic response analysis for a flexible high-aspect ratio wing excited by harmonic gust load is presented along with a companion wind tunnel test. A multidisci-plinary coupled numerical calculation is developed to simulate the flexible model wing undergoing gust load in the time domain via discrete nonlinear finite element structural dynamic analysis and nonplanar unsteady vortex lattice aerodynamic computation. A dynamic perturbation analysis about a nonlinear static equilibrium is also used to determine the small perturbation flutter bound-ary. A novel noncontact 3-D camera measurement analysis system is firstly used in the wind tunnel test to obtain the spatial large deformation and responses. The responses of the flexible wing under different static equilibrium states and frequency gust loads are discussed. The fair to good quanti-tative agreements between the theoretical and experimental results demonstrate that the presented analysis method is an acceptable way to predict the geometrically nonlinear gust response for flex-ible wings.     

油膜干涉测量翼型壁面摩阻低速风洞试验技术

油膜干涉试验可以定量测量得到模型壁面的摩擦应力,是开展摩阻特性研究的有力工具。通过对油膜干涉试验技术中的试验原理、试验方法、误差修正、数据结果分析等关键问题的研究,初步建立了试验系统,具备了开展模型壁面摩阻特性研究的试验能力。综合分析表面热膜试验结果与数值计算结果,在以模型弦长为特征长度,试验雷诺数 Re ＝1．3×106的条件下,翼型表面摩擦应力沿模型弦向的分布规律具有较好的一致性。     

正十七烷与正十八烷对航空煤油理化性能及液滴平均直径的影响

讨论了正十七烷与正十八烷(C₁₇,C₁₈)对燃料运动黏度、密度、表面张力理化性能的影响,探究了直链烷烃对SMD(Sauter mean diameter)的影响.将Key's混合规律、Tat混合规律和多项式拟合应用于密度公式的拟合,并将Key's混合规律与对数拟合应用于运动黏度公式的拟合,得出的理化性能拟合公式精度较高;通过测量纯煤油及C₁₇,C₁₈混合油雾化的SMD,实验发现随着直链烷烃的质量分数的增加,混合油的运动黏度增加,SMD也随之增加,雾化效果与纯航空煤油的相比较差,而C₁₈混合油的SMD要比C₁₇混合油要大,并且利用密度、运动黏度的拟合公式对SMD进行预测的偏差小于0.3%.      

低裂解度正癸烷物性快速计算方法

为研究正癸烷轻微裂解时物性计算方法,基于广义对应态法则建立起正癸烷裂解过程中的密度、黏度、导热系数和定压比热容的三种计算方法,包括直接计算、物性库插值计算、物性子库加权计算,在验证直接计算方法的准确度基础上,对比三种方法的计算精度、计算内存和计算时间来对其综合评估,并探究裂解度对高温物性的影响。计算温度变化范围为300~1020K,压力变化范围为1~15MPa,裂解度变化范围为0~0.25。结果表明:不同物性计算方法均能定性预测热物性的特殊变化趋势;以物性直接计算结果为基准,物性子库加权计算的误差与裂解度大小成正比;在不考虑计算内存的情况下,物性库插值计算方法能无损加速计算,低裂解度下快速工程计算可选择物性子库加权方法;裂解度对高温物性的影响显著且不可忽略。