非牛顿特性对重载点接触热弹流润滑的影响

采用近似完全数值解法对HKD-1型航空润滑油进行了重载点接触热弹流分析,将两种不同流变模型下的计算结果进行比较,并进行了试验验证.结果表明:重载情况下,润滑油的非牛顿特性较显著,非牛顿特性对于弹流润滑中润滑油的温度、剪应力和拖动力影响较大,对压力分布影响较小.采用牛顿模型计算误差非常大,采用Evans-Johnson模型计算精度较高.滑滚比增加,油膜温升增加显著,压力变化不大,只是二次压力峰的峰值有所减小.      

冰雹的数值模拟方法初探

冰雹的撞击对于飞机的整个结构存在着巨大的威胁,偶尔的一次遭遇冰雹袭击也是对飞机结构性能的严峻考验。由于实验测试需要投入大量的人力财力,因此,建立冰雹的数学模型显得尤为重要。目前所建立的数值模型主要有三种：有限元（FE）模型,任意拉格朗日一欧拉（ALE）模型和光滑粒子流体动力学（SPH）模型。首先论证了这些数值模型与实验测试结果吻合的情况,接着模拟了冰雹撞击飞机发动机进气道的过程,最后通过对这三种冰雹模型的比较,得出了SPH模型是分析冰雹撞击飞机问题的最有效、最好的模型的结论。     

基于直接配点法的滑翔轨迹快速优化设计

介绍了基于五次Gauss-Lobatto多项式的直接配点法在再入飞行器三维轨迹最优化问题中的应用。首先给出了再入飞行器轨迹优化问题模型,其中运动方程为三自由度模型,性能指标选为到达指定地点飞行时间最短,控制变量则为无量纲升力系数和倾侧角。再入飞行过程中受到加热率、过载和动压约束,终端状态受到目标位置约束。然后,应用直接配点法将最优控制问题离散化为非线性规划问题,将动态优化问题转化为静态参数最优化问题。选取各节点和配点上的状态量和控制量作为优化参数。最后应用SNOPT软件包对参数最优化问题进行求解。仿真结果表明直接配点法对于再入飞行器轨迹初始参数取值不敏感,且求解过程具有一定的实时性。     

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

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

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

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

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

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

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

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

无网格法耦合RNG k-ε湍流模型在亚、跨声速翼型黏性绕流中的数值模拟

基于移动最小二乘无网格方法,耦合RNG(Re-Normalisation Group)k-ε湍流模型求解雷诺平均Navier-Stokes方程。采用AUSM(Advection Upstream Splitting Method)+-up迎风格式求解数值通量,应用在高度各向异性点云结构中取得良好结果的点云重构技术结合移动最小二乘法拟合空间导数,并用三阶SSP(Strong Stability Preserving)型Runge-Kutta显式时间推进格式求解离散后的控制方程。在此基础之上,实现了对NACA0012、RAE2822翼型亚、跨声速黏性绕流的数值模拟,给出了翼型表面压力系数分布曲线、不同位置处的平均速度剖面、马赫数等值线等计算结果,并与实验值及相关文献数值模拟结果进行比较,结果吻合较好。表明所发展的结合点云重构技术的无网格方法耦合RNGk-ε湍流模型能够成功模拟翼型亚、跨声速黏性绕流,验证了所提算法的有效性,并拓展了无网格方法求解湍流流动的途径。     

Numerical analysis and optimization of boundary layer suction on airfoils

Numerical approach of hybrid laminar flow control(HLFC) is investigated for the suction hole with a width between 0.5 mm and 7 mm. The accuracy of Menter and Langtry’s transition model applied for simulating the flow with boundary layer suction is validated. The experiment data are compared with the computational results. The solutions show that this transition model can predict the transition position with suction control accurately. A well designed laminar airfoil is selected in the present research. For suction control with a single hole, the physical mechanism of suction control, including the impact of suction coefficient and the width and position of the suction hole on control results, is analyzed. The single hole simulation results indicate that it is favorable for transition delay and drag reduction to increase the suction coefficient and set the hole position closer to the trailing edge properly. The modified radial basis function(RBF) neural network and the modified differential evolution algorithm are used to optimize the design for suction control with three holes. The design variables are suction coefficient, hole width, hole position and hole spacing. The optimization target is to obtain the minimum drag coefficient. After optimization,the transition delay can be up to 17% and the aerodynamic drag coefficient can decrease by 12.1%.