共轴刚性旋翼气动干扰数值计算方法

建立了一个适用于共轴刚性旋翼气动特性分析的数值模拟方法。该方法采用任意拉格朗日欧拉方法(Arbitrary Lagrange Euler,ALE)描述的可压缩Navier-Stokes(N-S)方程求解流场,采用低数值耗散的Roe格式进行空间离散;使用多重嵌套网格方法以模拟双旋翼的运动。针对共轴刚性旋翼配平,引入"差量修正"策略解决了传统配平中雅克比矩阵计算复杂的问题。首先,对Harrington-2共轴双旋翼的悬停气动性能进行了计算,然后,对某2 m直径共轴双旋翼的悬停及前飞状态进行了计算,并与试验值进行了对比。结果表明:在典型状态下拉力系数的计算结果与试验值误差在3%以内,扭矩系数的计算结果与试验值误差基本在5%以内;所采用的数值计算方法对旋翼涡尾迹特征具有较高的捕捉精度,可以有效模拟共轴刚性旋翼悬停和小速度前飞下的复杂流场及其细节特征。     

基于Van Leer＋AUSM混合格式求解三维N-S方程

基于多块粘性结构网格,开展了三维N-S方程数值算法的研究。控制方程的空间离散采用有限体积法,在前人工作的基础上,发展了Van Leer＋AUSM混合格式并应用到对流通量的离散中,粘性项采用中心格式离散并利用格林定理计算粘性通量中的导数项,时间推进采用五步R-K法,湍流模型为S-A一方程模型。最后,以M6机翼和某超声速弹丸的粘性流场作为数值算例,计算表明：发展的数值算法对跨声速、超声速流场均具有较高的分辨率,适用于跨、超声速流场的数值模拟。     

Maxwell方程的高阶间断有限元数值解法

采用高精度方法求解时域Maxwell方程,方程的空间离散采用基于计算流体力学(Computational fluid dynamics,CFD)领域的高阶间断有限元格式,非定常时间迭代采用四步龙格-库塔格式。为了提高计算效率,本文采用了Quadrature-free implementation和网格分区并行技术。数值结果表明,采用高阶格式的情况下,采用稀疏网格便可以得到高精度数值解。另外由于本文的方法基于非结构网格,因此非常适合计算复杂外形的情况。     

双转盘偏心槽球体研磨均匀性定量评价方法

硬质合金球在航空航天、石油矿业等领域有着广泛应用,其研磨方法一直是研究的热点。评价硬质合金球研磨均匀性是研究其加工方法的一种主要手段。为进一步提高研磨均匀性评价方法的准确性,考虑表面原始形貌和研磨压力的影响,本文采用单因素实验法,使用MATLAB数学分析软件推导了球体研磨的材料去除方程,并结合材料去除方程提出了双转盘偏心V形槽球体研磨均匀性的定量评价方法与仿真模型。采用正交实验,对研磨均匀性定量评价方法进行了实验验证。实验表明,双转盘偏心V形槽球体研磨均匀性定量评价方法能够有效预测研磨压力、上盘转速和下盘转速对研磨效果的影响。     

软管-锥套式空中加油动态建模与性能分析

提出了一种基于凯恩方程的自主空中加油软管-锥套动态模型及性能分析方法。该方法将软管视为由有限段以铰链形式连接的刚性 杆组成,锥套为软管末端的一个质点。本文定义了描述系统状态的广义坐标及广义速率,导出了软管段位置多级递推公式和系统动态方程,估算了软管在加油机尾流、定常流和大气扰动下的气动载荷。通过数值仿真分析了平稳大气中锥套在加油机不同飞行条件下的稳态阻力和软管拖拽轨迹,验证了模型的正确性及系统稳定性。最后研究了大气扰动对锥套运动的影响及不同软管段的受扰运动。     

基于迹向控制量的编队星群相位控制

为了得到控制量与控制时机对相位改变的直观公式,基于Hill方程以及迹向控制的约束条件,定性分析共面绕飞编队中控制量、控制时机与编队卫星相位之间的关系,分析结果表明:当需要改变的相位角小于90°时,可实现最省燃料的控制;当需要改变的相位角大于90°时,可进行趋于最小控制量的控制;当在相对运动椭圆上下点实施控制时,相位无变化.     

基于Riccati方程解的高增益未知输入观测器设计

针对具有未知输入和测量噪声的一类Lipschitz非线性系统,研究了状态估计、噪声估计及未知输入重构问题.通过将输出噪声看作扩展状态,把原系统转化为描述系统.针对描述系统,首先基于Riccati方程的解,提出了一种高增益观测器设计方法,实现对系统状态的估计和测量噪声的重构;之后,设计二阶高增益滑模观测器精确估计输出的微分,并利用状态和输出微分的估计,提出了一种未知输入的重构方法.在一Riccati方程有解的前提下,所提出的未知输入和测量噪声的重构,均适用于强时变信号.最后,对一个实际模型仿真,验证所提出方法的有效性.      

三维可压缩 Navier-Stokes 方程的间断Galerkin 有限元方法研究

拓展了二维间断 Galerkin(DG)有限元方法研究,将该数值方法用于三维可压缩欧拉方程和 Navier-Stokes方程的求解。基于六面体网格单元,采用插值方法将物面的四边形面网格单元构造为弯曲面网格单元,更好地表述了真实物面特征;物面边界相邻体网格单元相应构造为高阶体网格单元,其余体网格单元采用八节点六面体单元,以较小的计算代价使网格满足 DG 方法计算需求。通过对三维带 bump 管道内流、圆球绕流以及旋转流线体绕流进行的数值求解,验证了边界弯曲方法的可行性及 DG 方法的高精度特性。此外,由于采用了隐式计算方法,仅需较少的时间步就能迭代收敛。     

旋转对管内流体流动的影响

本文描述了旋转管流的基本运动方程。在笛卡尔坐标系和柱坐标系中,针对两种基本的几何流动,推导出旋转情况下管内的流动方程,定义了影响旋转管流的流动和换热的两个无量纲准则数J和Ra,分析了二次流产生的原因,从而说明旋转对流动的影响以及旋转管流的运动特征。