一种求解固体推进剂粘弹性参数的通用数值积分法

根据粘弹性理论,用拉氏逆变换推导出了由拉伸松弛模量E(t)和体积松弛模量K(t)求固体推进剂剪切松弛模量G(t)的一种积分方程式,并用迭代和数值积分法给出了求解G(t)的数值积分算法。算例表明,该法简单实用,便于用ANSYS等软件进行粘弹分析时输入参数的确定。经比较,可将此源于梯形积分公式的数值积分算法作为一种求解固体推进剂粘弹性参数的通用数值积分法。     

剪切模量的近似函数和精确数值积分表达式

以粘弹性理论为基础，采用拉氏逆变换、迭代和数值积分的方法，分别推导出由应力松弛模量E(t)、粘弹泊松比v(t)计算固体推进剂剪切模量G(t)的近似函数和精确数值积分表达式。     

基于大应变黏弹性理论的药柱结构可靠性分析

基于大应变黏弹性理论,用积分随机有限元方法计算结构功能函数的各阶中心矩,用Gram-Charlier级数拟合概率密度函数,再由数值积分得到药柱结构的可靠度。算例表明：基于大应变黏弹性理论的计算结果与实际情况吻合。     

基于关机点状态的航天器落点计算及精度分析

当全程测量数据不完备时,如何根据关机点状态进行事后精确落点计算对于航天器精度评估非常重要。首先建立了航天器自由段和再入段的动力学方程,再入段充分考虑了空气阻力。接着采用复合数值积分法对动力学方程进行求解,使用四阶Runge-Kutta法对积分进行起步,采用Admas-Cowell积分减少误差。通过建立误差系数矩阵,研究了关机点状态估计信息和再入段空气密度偏差对落点计算的影响。实际数据计算和仿真结果表明,采用的复合数值积分法能提高落点计算的精度和速度,建立的误差系数矩阵能分析关机点状态和空气密度偏差对落点位置的影响,具有良好的实际应用价值。     

固体火箭药柱泊松比的表征与诊断

本文提出表征固体火箭药柱泊松比的三个精确表达式和一种数值积分算法。通过实验与计算给出一种固体推进剂药柱黏弹性泊松比的结果，分析了药柱黏弹性泊松比的表征与诊断以及诊断精度的有关问题。该文的结果对固体火箭发动机装药结构完整性分析计算有重要的实用价值。     

基于速度增量修正的卫星入轨段外测数据定轨新方法

在传统的数值积分精密定轨方法基础上,提出基于器箭分离速度增量修正的定轨方法,实现了器箭分离前后两段外测数据的联合定轨。仿真计算结果和任务实战数据验证结果表明,该方法大大提高了地月转移轨道入轨段外测定轨精度,比传统初轨确定方法提高约一个数量级。     

体积松弛模量与体积蠕变柔量的研究

简述了体积松弛模量Ｋ（ｔ）与体积蠕变柔量Ｂ（ｔ）之间的精确与近似转换关系，提出一种根据实测Ｂ（ｔ）诊断Ｋ（ｔ）的数值积分算法。通过体积蠕变实验室出一种丁羟固体推进剂Ｂ（ｔ）和Ｋ（ｔ）的实测与转换计算结果。     

土星探测器飞行动力学与轨道计算

根据牛顿第二定律和万有引力定律，建立了空间探测器的运动微分方程，并以土星探测器为例，利用精确的数值积分方法对其轨道进行了设计和计算，在土星探测器的轨道设计和计算中，充分利用了有关天体的引力机动，有目的地使探测器飞过一系列行星来增加探测器的能量，从而大大降低了探测器的发射条件。     

水星近日点的进动

计算绕太阳运行的航天器轨道的长期进动需要考虑广义相对论。然而,作者发现爱因斯坦
用他的广义相对论解释行星近日点进动的文章有误,他的计算行星近日点进动的公式是建立
一系列错误上的杜撰。本文详细分析了爱因斯坦的失误,尤其是他的积分错误,并用计算机
的数值积分直接验证了他的失误。     

柔性多体航天器单向递推组集建模方法

本文提出一种面向计算机建立柔性多体航天器动力学方程的新方法，其特点是完成运动学正向递推同时实现动力学方程的组集。这种方法不仅可提高航天器、空间机械臂动力学仿真的计算效率，而且可有效地克服数值积分病态的困难。文末一个空间四杆柔性机械臂的算例验证了这种建模方法的上述优点。     

探月返回飞行器跳跃式再入轨迹优化

基于节点自适应稀疏配点法,提出一种高精度求解探月返回飞行器跳跃式再入轨迹优化问题的方法。该方法的基本策略是：首先,应用节点自适应稀疏配点法对完整的跳跃式再入轨迹进行优化;然后,根据优化得到的控制变量对再入动力学方程进行数值积分;当积分至跳跃轨迹的最高点时,以积分得到的状态变量值作为新的初始条件,对二次再入轨迹重新优化。仿真结果表明：1)对二次再入轨迹重新优化能够显著提高跳跃式再入轨迹的优化精度,否则轨迹优化精度低,终端误差较大;2)在跳跃轨迹的最高点进行的二次优化是一种准实时优化,在跳跃式再入轨迹的制导领域具有潜在应用价值。     

小步长低阶Collocation方法轨道积分性能研究

介绍了Collocation方法的原理及误差构成,以Collocation方法在逐历元滤波轨道解算中的应用为切入点对小步长低阶Collocation方法轨道积分的性能进行了研究。分别就二体问题及真实GPS星座卫星轨道为对象进行了一系列验算,得到了一些有益结论：就积分方法本身性能而言,Collocation积分方法较传统数值积分法精度要高且稳定性要好;在一定长度的积分步长范围内,将小步长低阶Collocation方法用于滤波轨道解算是可行且有效的。 但在积分步长很小的情形下,从计算效率方面考虑,不推荐使用Collocation方法;利用小步长低阶Collocation方法取代大步长高阶方法能提高效率且可抑制蚀现象引起的积分结果发散。经试验表明：选用积分步长为900s的4阶方法是合适的。     

空地一体化信息对抗系统构架设想

首先介绍“空天一体化”作战体系概念及信息对抗系统的组成和系统特点,并基于目前信息对抗系统的现状和对未来信息战的推测,对空地一体化信息对抗系统构架做了初步的设想。重点介绍空地一体化信息对抗系统的需求、系统构架、组成、特点及关键技术。     

有限空间轴对称射流流场的数值模拟

本文在к－ε紊流模型下，采用罚有限元法对有限空间轴对称射流流场进行数值模拟。引用罚函数优化迎风格式和退化积分，建立有限元解法的数值矩阵，在不同截面上求得工，压力和紊流参数分布，计算结果与实验数据十分接近。本文研究结果可应用于火箭发动机的喷射设计和射流量论的研究。     

The method of nearby problems for estimation of numerical error in orbital-mechanics simulations

The method of nearby problems has previously been applied to validate the accuracy of numerical solutions of the problems in fluid dynamics. In this paper, the application of this method is extended to the problems in orbital mechanics. Using this method, the error in the numerical simulations of two-body and three-body equations of motion is estimated. A fixed step-size and a variable step-size numerical integration scheme are used for the simulations. The method of nearby problems is found to provide reliable error estimates and has advantages of cost efficiency and ease of use.     

线圈炮——电磁同轴发射器及其系统分析

线圈炮在宇航与国防中有重要应用前景。本文先简述其分析模型与方程。方程为非线性变系数常微分方程的初值问题，有时呈刚性。用Ｔｒｅａｎｏｒ方法求解，所得结果与国外已有装置的试验数据符合良好。并对一个小炮做了初步的原理设计。     

三维药柱固体火箭发动机内弹道数值模拟

介绍了一种适用于工程内弹道计算的数值模拟新方法，用直接积分法计算三维药柱燃面，能得到精确的计算结果，在内弹道计算中，设燃烧产物准定常流动，使用分段解析法进行内弹道的数值模拟。对某翼柱型三维药柱发动机内弹道进行计算，取得了接近实验的结果。     

低重力环境下三维液体非线性晃动的数值模拟

本文讨论低重力环境下三维液体非线性晃动问题。采用ＡＬＥ（任意的拉格朗日－欧拉）运动学描述跟踪自由液面；对于ＡＬＥ描述的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，在时间域上采用分步离散方法中的速度修正格式进行数值离散，在空间域上利用有限元方法进行了数值离散；推导了考虑表面引力效应时有限元边界条件的弱积分形式；给出了表面张力的数值计算公式。模拟了考虑表面张力情况下圆筒型贮腔中液体的非线性晃动，并得到了贮腔壁面处     

Unsteady laminar separation over impulsively moved cylinders

A numerical method is employed to integrate the unsteady laminar boundary layers about circular and elliptic cylinders started impulsively from rest. An upwind differencing scheme permits the integration through thin layers of reversed flow that are embedded in the boundary layer. The results are in agreement with numerical solutions of the full Navier-Stokes equations and with small time inner-outer expansion solutions. At later times, a singularity of the Goldstein type appears; Sears and Telionis defined this singularity as unsteady boundary-layer separation. For the circular cylinder, this singularity travels upstream and for large times, arrives asymptotically at the point of steady-state separation.