航天器交会两点边界值问题

从绝对运动和相对运动两方面讨论近地空间航天器交会中的两点边界值问题。其中，绝对运动涉及多圈Lambert问题，以Lagrange时间方程为研究工具，而相对运动Lambert问题应用C—W线性解。对圆轨道之间的双冲量转移，给定转移角与转移时间，研究最小变轨速度增量所对应的转移圈数与轨道参数的求解方法，提出满足最小变轨速度增量要求的多圈转移的工程图解法，并从工程应用出发，在飞行时间约束下，按最小速度增量要求，阐述航天器交会两点之间飞行轨道（轨迹）设计方法。这种方法将飞行轨迹划分为初始漂移段、轨道转移段与终端停泊段三部分，应用两点边界值问题的解，选择两次冲量机动时刻，使速度增量之和最小。模拟算例表明，这种方法对航天器交会设计是适用的、有效的。     

蜜蜂悬停飞行控制的仿生力学研究

研究蜜蜂悬停飞行的控制问题.用数值求解N-S方程的方法计算拍动翅体的控制导数;用特征模态分析方法分析控制特性.文中δu,δω,δq和δθ分别为水平方向速度,垂直方向速度,俯仰角速度和俯仰角度的扰动量;δ(φ)和δ(φ)分别为拍动幅度和拍动平均角增量;δα1表示上、下拍迎角同时增加的增量;δα2表示下拍迎角增加(或减小)而上拍迎角减小(或增加)的增量.获得以下结果:(1)悬停飞行时,改变Φ和α1主要产生垂直力的变化;改变φ主要产生俯仰力矩的变化;改变α2主要产生水平力及俯仰力矩的变化.(2)蜜蜂悬停飞行的纵向扰动运动由3个特征模态构成:不稳定振荡模态,快衰减模态和慢衰减模态;为实现稳定的悬停飞行,不稳定振荡模态和慢衰减模态需要控制.为控制不稳定振荡模态,以δ(φ)(或δα2)反馈δu,δq和δθ这3个量的某种组合便可;为控制慢衰减模态,以δ(φ)(或δα1)反馈δω便可.这就是说,该昆虫只需用控制变量δ(φ)和δ(φ)(或δ(φ)和δα2,δα1和δ(φ)1,δα1和δα2),便可稳定地悬停(当然也可4个控制变量都用).     

Supersonic Two-Dimensional Minimum Length Nozzle Design at High Temperature. Application for Air

When the stagnation temperature of a perfect gas increases, the specific heat ratio does not remain constant any more, and start to vary with this temperature. The gas remains perfect, its state equation remains always valid, except it will name in more calorically im-perfect gas or gas at High Temperature. The goal of this work is to trace the profiles of the supersonic Minimum Length Nozzle with centered expansion when the stagnation temperature is taken into account, lower than the threshold of dissociation of the molecules and to have for each exit Mach number several nozzles shapes by changing the value of the temperature. The method of characteristics is used with a new form of the Prandtl Meyer function at high temperature. The resolution of the obtained equations is done by the second order of finite differences method by using the predictor corrector algorithm. A study on the error given by the perfect gas model com-pared to our model is presented. The comparison is made with a calorically perfect gas for goal to give a limit of application of this model. The application is for the air.     

FEM Analysis of Spring-backs in Age Forming of Aluminum Alloy Plates

The age forming technology, characterized by huge spring-backs, has been developed to manufacture large integral wing-skin panel parts, which necessitates devising a method of predicting spring-backs. A 7B04-T7451 aluminum alloy creep test in tension is accomplished at 155 ℃, and the creep curves are obtained. The material constants of the mechanism-based creep constitutive equations are determined through experiments. The age forming process and the spring-backs of 7B04 aluminum alloy plates are analyzed using the commercial finite element software ABAQUS. The effects of plate thickness and forming time on spring-backs are researched. The spring-backs decrease with the increase of plate thickness and forming time. The test results verify the reliability of the finite element method （FEM） analysis.     

预处理方法在叶轮机械三维数值模拟中的应用

应用预处理方法求解三维可压缩流动Navier-Stokes方程, 发展了一种适用于叶轮机械的既可以求解跨声速流动也可以求解不可压流动的数值模拟方法, 对NASA Rotor67和北航低速大尺寸压气机算例进行了计算, 计算结果与实验数据相吻合, 表明本文采用的预处理方法具有较大的工程实用性.      

基于Navier-Stokes方程/自由尾迹/全位势方程的旋翼流场模拟混合方法

为充分考虑旋翼尾迹对流场的影响和减少尾迹的数值耗散,建立了一个基于Navier-Stokes方程／自由尾迹分析／全位势方程的旋翼流场求解的新的混合方法。该方法的求解域由三部分组成：一是围绕旋翼桨叶周围的粘性区域．采用可压Navier-Stokes方程来捕捉近场信息,包括激波及尾迹;二是离桨叶较远、粘性可以忽略的等熵流区域,以全位势方程来描述其流动;三是在无粘区域中嵌入自由尾迹模型,模拟桨尖涡从粘性区域进入势流范围的发展变化。为便于流场分区求解和信息传递,采用了嵌套网格方法,并给出了不同区域之间的信息传递方法。以两叶的Caradorma＆Tung模型旋翼和四叶的UH-60A直升机旋翼为算例,计算给出了旋翼桨叶表面的压强分布以及桨尖涡的位置,并与可得到的试验数据及无尾迹模型方法的计算结果进行了对比,表明本文的混合方法能够很大程度地减少旋翼尾迹的数值耗散。     

共轴刚性旋翼气动干扰及操纵特性分析

建立了一套基于高效配平策略和嵌套网格技术的共轴刚性旋翼气动数值分析方法。流场求解采用雷诺平均Navier Stokes方程,空间离散采用Roe monotone upstream centered schemes for conservation laws（MUSCL）格式,湍流模型为Spalart Allmaras模型;配平计算中引入“差量修正”策略。分析了共轴刚性旋翼模型在不同飞行状态下的气动干扰现象和操纵特性。结果表明:所建立的配平方法在同等精度下,比传统方法可显著节省计算时间;共轴刚性旋翼间距较小,其气动干扰明显强于普通共轴旋翼;悬停状态下,由于受上旋翼下洗流影响,下旋翼总距较大但拉力较小;小速度前飞时气动干扰引起下旋翼气流不对称性增大,双旋翼周期变距差也变大;采用适当的升力偏置量,可以提高共轴刚性旋翼气动效率。      

桨叶外形对共轴刚性旋翼悬停性能影响的CFD分析

桨叶外形是影响共轴刚性双旋翼直升机悬停性能的主要因素之一,研究其影响机理对提高共轴刚性旋翼的悬停性能具有重要意义。为了提高外形参数对旋翼气动特性影响的分析精度,首先,基于运动嵌套网格技术和可压RANS方程,建立了一套适用于共轴刚性双旋翼悬停流场模拟的CFD方法;然后,进行了悬停状态共轴旋翼流场气动特性的计算分析,得到上、下旋翼拉力和扭矩沿桨叶展向和周向的分布特征;最后,着重进行了桨叶平面外形参数对共轴刚性旋翼悬停性能影响的分析,包括后掠角、后掠起始位置及尖削非线性弦长分布等,得到不同参数下共轴刚性旋翼悬停效率和压强分布等特征。结果表明,桨尖后掠能够使负压中心外移,延缓气流分离,从而提高悬停效率;尖削及桨叶面积集中在中部段的非线性弦长分布能够有效对桨尖进行卸载,优化气动载荷分布,从而改善悬停性能。     

可降低气动热效应的类凹腔外形优化设计

针对飞行器上由凸起物形成的类凹腔气动加热问题,采用数值方法求解三维N-S方程,研究了类凹腔外形结构的高超声速气动加热规律,获得了三维高超声速流场和局部热流分布,并详细分析了局部的流场结构和气动加热机理。针对凹腔前壁面热流密度过高的问题,提出并验证了一种降低前壁面边缘热流密度的优化外形,将前壁面的热流密度降低至优化前的20%左右。     

有限元线性方程组改进的全稀疏存贮方法

有限元线性方程组的系数矩阵一般具有稀疏性和对称性的特点,全稀疏存贮方法就是利用这些特点,只存贮对称部分的非零元素,采用链表式管理,即节省存贮空间,又便于动态更改;但在实际求解计算时,查找元素花费的时间较多;这里给出了一种改进的方法,是在链表的基础上生成相应的线性表,能够方便快速查询元素,使存贮结构同时具有链表和线性表的优点.数值算例证明该存贮方案在时间和存贮上都较为占优,可靠高效,能够应用于有限元线性方程组的求解.