比例式变推力固体姿控发动机非稳态特性研究

为了研究比例式变推力固体姿控发动机的内流场非稳态特性,建立了比例针栓推力器的二维轴对称计算模型,基于动网格技术模拟入口压强随喉部面积变化而变化的推力器工作模式,得到了内流场各性能参数的变化规律。结果表明:在非稳态工作过程中,内流场会出现典型的亚音速回流区、斜激波和流动分离等特征,入口压强、针栓壁面及喷管壁面压强均随针栓靠近喉部而增大,推力器推力逐渐上升,实现了推力连续调节。开关频率会加剧针栓前进过程中头部压强波动。针栓头部收敛角越大,其头部回流区越小。当喉部面积一定时,燃速压强指数越高,发动机压强与推力变化范围越大,为实现预设的推力调节范围,需要选择合适的燃速压强指数。     

非结构网格二维理想MHD绕流逆风格式解法

在非结构网格上发展了针对二维理想磁流体方程组的逆风格式求解方法.控制方程中的对流项采用AUSM格式处理,时间推进采用显式5步龙格-库塔方法.为了消除计算中产生的磁场散度的影响,引入了双曲型散度清除方法.通过对磁流体激波管问题的求解验证了该方法对激波的捕捉能力,对有均匀磁场干扰下的喷管流动情况进行了数值模拟,并与文献中结果进行了对比.计算结果显示了磁场对磁流体流动的干扰效应,该结果与参考文献中的数值模拟结果相吻合.     

星载大型网状天线非线性结构系统有限元分析

论述了大型网状可展开天线的非线性理论模型及其非线性模拟有限单元模型,并以一副4．2m天线为例进行了有限元建模分析,获得了这类天线的动力学振动特性。研究结果表明：拉索几何非线性和间隙接触非线性对大型网状天线系统的动力学特性具有较大影响。     

基于能量等效的直齿锥齿轮时变啮合刚度计算

针对直齿圆锥齿轮啮合刚度的计算问题,从微元思想出发,将变截面齿廓划分为若干等截面微段齿段,基于能量等效建立微段齿段啮合刚度计算模型,并利用积分方法获得单齿啮合刚度。此外,基于力平衡和变形协调条件进一步提出了齿轮时变啮合刚度计算模型,同时根据几何关系导出了相应传动误差计算式。利用有限元分析对解析计算模型进行了验证,并分析了误差来源。结果表明:利用该模型不仅能够将直齿锥齿轮啮合刚度计算精度保证在2%以内,还达到了快速求解的目的。      

动网格技术在增升装置优化中的应用

将基于Denaulay背景图的动网格方法应用于增升装置二维缝道参数优化设计过程中。针对二维增升装置的构型特点和翼型运动的特点,采用了添加辅助点的双层Denaulay背景图动网格方法,并详细论述了辅助点添加、运动技术。为了进一步提高动态网格质量,提出了针对前后缘运动位置中比较极端的两种位置生成两个初始网格,并采用两种初始网格进行动态网格的生成。结果表明新动态网格方法可以为二维增升装置提供高质量的动态网格。     

W型地下井发射环境数值模拟与分析

为研究导弹地下井发射环境,采用计算流体力学方法对W型井导弹发射过程的流场进行了数值模拟。通过将仿真数据与理论计算及试验结果进行对比,验证了计算模型和数值方法的准确性。采用动网格技术处理流场计算区域随导弹运动的变化过程。仿真结果表明,在整个发射过程中燃气流具有反射、引射和充分发展3个阶段,造成发射环境参数有较大变化。分析了压力和温度参数的变化规律。该数值算法和仿真结果可推广应用于发射井的设计、优化及改进。     

基于CFD动网格的内装式空射分离研究

研究了使用飞机发射运载火箭时的机箭分离过程.提出了一种可以解决内装式空射分离问题的整体解决方案,主要基于CAD三维实体建模方法、多体动力学方法和结合六自由度运动的计算流体力学(CFD)动网格方法.以CFD动网格计算为核心,进行了空射分离过程关键问题的研究,得到了火箭在载机气动干扰环境下的运动特性,为轨道卫星和弹道导弹的快速、廉价发射打下了研究基础.     

Adaptive mesh refinement computation of acoustic radiation from an engine intake

A block-structured adaptive mesh refinement (AMR) method was applied to the computational problem of acoustic radiation from an aeroengine intake. The aim is to improve the computational and storage efficiency in aeroengine noise prediction through reduction of computational cells. A parallel implementation of the adaptive mesh refinement algorithm was achieved using message passing interface. It combined a range of 2nd- and 4th-order spatial stencils, a 4th-order low-dissipation and low-dispersion Runge–Kutta scheme for time integration and several different interpolation methods. Both the parallel AMR algorithms and numerical issues were introduced briefly in this work. To solve the problem of acoustic radiation from an aeroengine intake, the code was extended to support body-fitted grid structures. The problem of acoustic radiation was solved with linearised Euler equations. The AMR results were compared with the previous results computed on a uniformly fine mesh to demonstrate the accuracy and the efficiency of the current AMR strategy. As the computational load of the whole adaptively refined mesh has to be balanced between nodes on-line, the parallel performance of the existing code deteriorates along with the increase of processors due to the expensive inter-nodes memory communication costs. The potential solution was suggested in the end.     

Shape control of spacecraft formation using a virtual spring-damper mesh

This paper derives a distance-based formation control method to maintain the desired formation shape for spacecraft in a gravitational potential field. The method is an analogy of a vir-tual spring-damper mesh. Spacecraft are connected virtually by spring-damper pairs. Convergence analysis is performed using the energy method. Approximate expressions for the distance errors and control accelerations at steady state are derived by using algebraic graph representations and results of graph rigidity. Analytical results indicate that if the underlying graph of the mesh is rigid, the convergence to a static shape is assured, and higher formation control precision can be achieved by increasing the elastic coefficient without increasing the control accelerations. A numerical exam-ple of spacecraft formation in low Earth orbit confirms the theoretical analysis and shows that the desired formation shape can be well achieved using the presented method, whereas the orientation of the formation can be kept pointing to the center of the Earth by the gravity gradient. The method is decentralized, and uses only relative measurement information. Constructing a distributed virtual structure in space can be the general application area. The proposed method can serve as an active shape control law for the spacecraft formations using propellantless internal forces.     

基于离散等价方程的非结构网格有限差分法

激波装配法把解析关系式嵌入流场避免了间断引起的理论问题，使全场一致高精度的实现成为可能，有望解决超声速流动转捩研究中的感受性模拟难题。但装配复杂激波结构以后，被分割的流场空间常出现不规则的几何形状，这给常规的基于结构网格的有限差分法应用带来困难。基于离散等价方程理论，提出了一种新的基于非结构网格的有限差分法，在空间二维离散点附近仅用3条网格线就可以构造出一阶迎风格式。数值算例表明收敛过程对网格质量不敏感，解决了激波装配法模拟激波相交出现小夹角后使用结构网格进行计算存在的难题。根据这种方法的特点展望了未来的应用前景。