多体系统振动的传递矩阵法

本文是文献工作的继续。导得了空间运动的任意刚体的场传递矩阵，建立了多体系统振动的传递矩阵法，为多体系统振动特性的研究提供了新的理论和方法，从而有效地简化多体系统振动特性的研究，大大减小了计算工作量，并且便于解析讨论。     

关于薄膜反射镜成形机理的研究

聚酰亚胺薄膜反射镜是近几年发展起来的新概念。以Karman方程为基础讨论了薄膜反射镜的成型原理，分析了外载荷对薄膜反射镜面型的影响；给出了非线性耦合Karman方程的求解方法；以一种聚酰亚胺圆薄膜为例，用有限元法对Karman方程进行了数值求解，用zemike多项式拟合方法得到了圆薄膜在外载荷作用下的近似解析解，证明了把Karman方程解析解表示为多项式的合理性；用有限元法得到的不同外载荷作用下的薄膜的变形分析证明了通过控制外载荷的加载方式和载荷强度可以控制薄膜反射镜的面形，改善反射镜像差。     

飞机最快上升转弯轨迹的计算

70年代初,Miele A等人提出SGRA方法用于优化飞行轨迹计算。但方法存在时间加权系数难于确定、长时间的飞行轨迹计算积分精度不易控制及要计算数字微分等不足。对此提出一改进梯度方法,并计算了某型战斗机从亚音速三维转弯升高到超音速的最短时间飞行轨迹。结果表明,本方法稳定有效,为解决飞机最快机动轨迹的优,乜提供了一种行之有效的手段。     

基于变域变分有限元的翼型反设计

反问题作为一种自由(未知)边界问题近年来得到了广泛的研究。变域变分通过把自由边界结合在变分泛函中，使其与求解流场的控制方程结合起来，从而自由边界求解和流场分析可以完全耦合进行。本文采用变域变分和可变形有限元来设计可压缩和不可压缩翼型，计算显示即使在初始翼型的一侧为三角形时仍能得到满意的结果。文中还采用了两种计算网格(H型和O型)来设计翼型，比较了它们的精度，计算表明本方法完全能用于翼型的反设计。     

六类高超声速激波-激波干扰的数值模拟研究

基于粘性全N-S方程，利用有限体积法进行数值离散，对Ⅰ到Ⅵ六种类型的激波一激波干扰现象进行了数值模拟，其中采用了二阶Harten-TVD格式，LU隐式算法以及Baldwin-Lomax代数湍流模型。每种类型都得到了清晰的干扰流场结构，其中Ⅳ型干扰的模拟结果与实验结果进行了比较，符合较好。数值模拟结果的对比指出，由入射位置差异引起的不同类型激波-激波干扰结构有很大不同，对壁面压强分布的影响也存在较大差异。     

顺序函数法求解二维非稳态热传导逆问题

基于有限控制体积法对二维非稳态热传导问题的数值模拟，导出了处理二维非稳态热传导逆问题的顺序函数法，并用该方法来对一典型的圆环域边界条件反演问题进行了反演计算，结果表明，顺序函数法是求解二维非稳态热传导逆问题的有效方法，当测量噪声比较小时，顺序函数法能得出较高精度的反演结果；当测量噪声增大时，反演结果仍能较好地再现出精确解在空间和时间方向上的变化趋势，具有一定精度。     

DEVELOPMENT AND APPLICATIONS OF WENO SCHEMES IN CONTINUUM PHYSICS

INTRODUCTIONHigh order accurate weighted essentiallynon- oscillatory( WENO) schemes have been de-veloped to solve a hyperbolic conservation lawut+ .f( u) =0 ( 1 )The first WENO scheme was constructed in Ref.[1 ]for a third order finite volume version in onespace dimension.In Ref.[2 ],third and fifth or-der finite difference WENO schemes in multiplespace dimensions are constructed,with a generalframework for the design of the smoothness indi-cators and nonlinear weights.L ater,second…     

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

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

后摆心潜入式喷管柔性接头力学特性数值仿真

针对后摆心潜入式喷管柔性接头在燃气压力与摆动联合载荷下,弹性件与增强件的力学特性计算问题,建立了对称刚体-柔性体结合的几何模型;采用了有限元方法,获取了柔性接头在燃气压力、摆动载荷及联合载荷下应力、弹性比力矩分布规律。研究表明:随着燃气压力增大,弹性比力矩降低;随摆角增大,弹性比力矩增大。文中的数值仿真方法,可以应用于柔性接头力学特性的计算,为柔性接头强度的安全性评估提供指导。     

非开普勒轨道下卫星交会变轨决策研究

针对卫星在非开普勒轨道下的运动特征,提出了研究卫星交会变轨决策的方法。首先,在Hill坐标下,采用比例导引方法求解卫星进行变轨机动时所需要的速度增量。然后,再通过Lagrange插值算法解出卫星在任意特征点上变轨机动所需要的速度增量,根据数字仿真得出大量的仿真数据并画出满足变轨决策条件的决策曲面。最后,对决策曲面的涵义做出分析,说明了该方法对于研究非开普勒轨道下卫星变轨决策问题具有一定的意义。