化学非平衡高超声速电离流动数值研究

本文对非平衡化学区应高超声速绕流零攻角轴对称钝头体电离空气粘性流动进行了研究，通过带有化学源项的Ｎａｖｉｃｒ－Ｓｔｏｋｅｓ方程对流动作数值模拟，总的连续方程由单个组元守恒方程代替。高温电离空气中存在７个主要的组元，它们依次为Ｎ＿２，Ｏ＿２，ＮＯ，ＮＯ￣＋，Ｎ，Ｏ，和ｅ￣－。化学动力学模型由７个组元和这些组元之间６个可能的反应来描述。发展了一种数值方法以求解这组完全耦合的与时间有关的偏微分方程组。利用一类新的有限差分方法来求解这组以守恒形式表示的控制方程。这是一种以ＴＶＤ格式来离散粘性流通量的方法。在数值实验中采用了ＴＶＤＹｅｅ－Ｈａｒｔｅｎ格式，对全部算例都是用激波捕捉法。关于化学组元，我们试图处理完全催化壁、非催化壁和有限催化壁等三种不同的表面边界层条件。计算是对ＲＡＭＣ－Ⅱ试验条件进行的，得出了７１ｋｍ高空情形的典型流场数值解。对表面压力，Ｓｔａｎｔｏｎ数和热流率给出了化学反应气体结果与完全气体结果的比较。计算的流场电子密度与ＲＡＭＣ－Ⅱ实验数据进行了比较，计算结果与实验数据符合是相当好的。     

航空发动机涡轮气动技术发展研究

本文对航空发动机涡轮气动技术的发展做了概要论述，并对典型的涡轮气动设计方法重点进行了介绍其中包括先进的三维分析和ＣＡＤ技术的应用。     

多向碳／碳复合材料超高温力学性能测试技术研究

本文系统地运用超高温力学性能测试的理论，通过对Ｇｌｅｅｂｌｅ－１５００热模拟实验机的改装，建立了高达２８００℃碳／碳复合材料加热、测控温、加载、数据采集和处理一体化的快速通电加热测试技术，着重介绍了各分系统改装的原理和方法。     

基于SQP方法的常推力月球软着陆轨道优化方法

月球软着陆是未来月球探测中的一项关键技术。针对这项技术，本文给出了一种基于SQP方法的常推力月球软着陆轨道优化方法。该方法通过将常推力月球软着陆轨道离散化，利用离散点处状态连续作为约束条件，把常推力月球软着陆轨道优化问题归结为一个非线性规划问题，对于此问题的求解，其初值均为有物理意义的状态和控制量，从而避免了采用传统优化方法在解决此优化问题时对没有物理意义变量初值的猜测。最后，利用SQP方法求解了此轨道优化问题。仿真计算结果表明这种离散化的方法应用于此轨道优化问题可以避免传统轨道优化方法对初值敏感的问题。     

小样本或单次仿真运行的统计分析方法

提出用Bayes方法对小样本或单次运行的仿真输出进行分析,给出了几种条件下的Bayes方法和伽玛分布置信区间的估计方法,提出了Bayes方法中先验分布的确定方法和准则。并用M/M/1排队系统对几种分析方法进行了实例研究。实验结果表明,Bayes方法有明显的优越性和应用前景。     

混合网格在求解三维N—S方程中的应用

用三棱柱／四面体组成的混合网格求解三维Navier-Stokes方程，在物面附近采用三陵柱网格，而其他地方则采用四面体网格，与单一类型网格相比，混合网格易于调节粘性层内网点数目，并能减少存储量，提高解的效率，Navier-Stokes方程的求解采用有限体积法，对人工耗散项进行了改进，使得求解的结果更接近真实的物理解，本文对上述方法开展了数值实验研究。结果表明，用三棱柱／四面体组成的混合网格求解N－S方程是非常有效的。     

球头铣刀在虚轴机床上的刃磨

对虚轴刃磨机床的工作原理进行了分析.建立球头铣刀刀刃的曲线方程,利用活动标架法,推导球头铣刀前后刀面的数学模型.根据球头铣刀的数学模型,求解球头铣刀的广义刀位轨迹,利用虚轴机床的反解模型,把广义的刀位轨迹换算成连杆坐标系下连杆的长度,控制连杆长度的变化,完成刃磨加工运动.进而论证了利用虚轴机床进行刀具加工的优势.     

基于路径跟踪原理求解6-SPS并联机构位置正解

以6-SPS并联机构为例,基于路径跟踪原理,提出一种求解并联机构位置正解的预估-校正法.该方法通过求解并联机构位姿方程的微分初值问题,经过路径跟踪,可求出指定杆长动平台的正确位姿,并可求出从初值到待求值的中间结果.仿真结果表明,该方法具有计算速度快,收敛范围大,计算结果可靠的特点.     

选择性激光烧结PC粉末的建模研究

建立选择性激光烧结粉末快速成形工艺过程的数值计算模型,并用于估算烧结深度.在该计算模型中,综合考虑了材料的热物性参数随温度变化的情况和激光光强分布不均匀的因素.采用有限差分方法求解烧结过程中的瞬态传热微分方程,对烧结深度进行模拟计算,结果与实验值吻合较好.     

对称矩阵反问题解的稳定性

考虑如下对称矩阵反问题：问题Ⅰ：给定Ｘ，Ｂ∈Ｒ（ｎ×ｎ），求Ａ∈ＳＲ（ｎ×ｎ）使得ＡＸ＝Ｂ，其中ＳＲ（ｎ×ｎ）＝｛｛Ａ｜Ａ∈Ｒ（ｎ×ｎ），ＡＴ＝Ａ｝。问题Ⅱ：给定，求使得其中是矩阵的Ｆｒｏｂｅｎｉｕｓ范数，ＳＡ表示问题Ⅰ的解集合。本文讨论问题Ⅰ解的稳定性，给出问题Ⅱ解的扰动界。