航天器四冲量固定时间最优交会

本文研究了航天器的四冲量固定时间最优交会问题及仿真，在两航天器相对运动线性化方程组和最优交会基本方程组基础上，利用特征矢量的对称性。求得了最优点火时间、最优推力和速度增量的计算方法，仿真结果表明，该方法是可行的。     

火箭喷管的二维两相流

本文介绍了编制一组解轴对称喷管气体/粒子流动方程组的计算机程序的进展情况。已经编写了两个程序,是用来解喷管喉部跨声速区的方程组的。第一个程序把两相流体看成是一种具有修正等熵指数和修正分子量的重理想气体,并解跨声速等熵流动方程组,得到的初始流场形状作为第二个程序的输入。第二个程序包含有气体/粒子混合物非平衡效应。这两个程序为进行超声速计算提供准确的初始线数据,这些数据将用作第三个程序的输入。第三个程序是解喷管超声速区两相流动方程组的,不久即将写成。     

基于Reddy简化高阶剪切理论的复合材料对称角铺设矩形板横向弯曲一般解析解

采用复数级数法求解基于Reddy简化高阶剪切理论的复合材料对称角铺设矩形板横向弯曲问题。将待定位移函数展开为复数级数,代入该弯曲问题控制偏微分方程组,确定特征根和挠度待定常数与其他位移函数待定常数之间关系式。首次给出了该弯曲问题实数形式的一般解析解。将该一般解析解代入矩形板弯曲边界条件和角点条件,根据正弦级数的正交性建立关于挠度函数待定常数的线性代数方程组,求解此线性代数方程组可确定挠度函数待定常数。建立了该问题解析求解模式。将Reddy高阶剪切理论解析解与经典理论、一阶剪切理论解析解进行对比计算,验证了一般解析解,并给出数值算例。     

基于时域有限体积法的计算电磁

介绍了基于计算流体力学(CFD)中的时域有限体积法(FVTD)的电磁计算.根据Maxwell方程组,给出了计算电磁(CEM)中FVTD的时间离散、MUSCL迎风型格式插值、空间离散和边界条件的处理,以及雷达散射截面(RCS)的计算公式.对无限长理想导体圆柱和线天线二维完全导电体的表面诱导电流与雷达散射截面计算结果表明:FVTD是一种高效准确的时域有限方法.     

弹道式导弹的落点计算

本文采用了一种计算弹道式导弹落点参数的新方法。在地球为椭球的情况下,综合利用椭圆轨道理论和中间轨道理论,建立弹道式导弹的落点方程组。利用数值计算方法对方程组进行求解,得到弹道式导弹的落点参数。文中导出了一组落点参数相对于轨道根数的偏导数,以及相对于关机点(被动段起始点)的状态参数的偏导数,从而较方便地计算落点散布。分析了地球表面形状对落点参数的影响。在计算机上进行数字仿真结果表明:该方法是正确的、可行的;在满足精度条件下,计算简便,速度快;能够适用各种型号的弹道式导弹的设计。     

柔性接头界面端角点应力奇异性

基于界面理论,建立了柔性接头增强件、弹性件两种结合材料形成的界面端角点附近奇异应力场方程,依据界面连续条件和自由边界条件,得到了与奇异应力场关联的特征值λ的线性齐次方程组,通过该方程组的求解和公式推导,给出了柔性接头界面端角点附近含应力强度系数的奇异应力场和位移场;并根据特征值λ的取值范围,获得了应力奇异性消失的柔性接...     

基于粒子群算法的多波位高度估计方法

为解决雷达斜视模式下对目标区域的高度估计问题,提出了一种基于粒子群算法的多波位高度估计方法。利用粒子群优化算法建立约束条件和目标函数,通过多次迭代使目标函数趋于一个最优值,同时得到高度的最优解。利用粒子群优化算法求解多波位测高方程组可减小方程组近似处理误差。此外,利用粒子群优化算法可随时更改波位数目,增强了该测高方法使用的灵活性,有效提高了高度估计精度。通过理论仿真和实测数据仿真分析,验证了粒子群优化算法在求解多波位测高方程组时的有效性。结果表明:该方法具有较高的高度估计精度。     

基于稳定度设计的空间站姿态TEA稳定LQR控制器

TEA姿态求解和控制器设计是空间站姿态研究的重要问题。将TEA求解转化为非线性方程组求解,推导了准确的TEA附近线性化的空间站姿态动力学方程,采用LQR的稳定度设计进行空间站的TEA稳定控制器设计。研究表明:采用非线性方程组迭代法求解TEA其求解过程更简单、结果更准确;基于稳定度设计的LQR控制器使系统的闭环极点按要求远离虚轴,因此控制器的稳定性更好。     

一种多标校源的高轨伴星时差频差定位算法

针对高轨伴星时差频差无源定位系统中卫星位置、速度、时差和频差等参数的测量系统误差严重影响定位精度的问题,提出了一种基于四个或四个以上已知位置的地面标校源的高斯-牛顿定位算法。该算法首先利用差分法消除星间时差、频差测量的系统误差,再利用标校源的时差频差测量方程组估计出主星和伴星的相对位置和相对速度误差,最后结合时差、频差、地球球面以确定非合作辐射源位置。理论和数字仿真均表明在时差和频差测量的随机误差较小时,本文算法的均方根误差（MSE）接近克拉美-罗下限(CRLB)。
     

基于减小初始扰动的弹管间隙优化设计方法

针对火箭/发射装置系统,用牛顿-欧拉法建立了火箭和发射管的动力学方程组,基于MATLAB平台编写了方程组求解程序,在此基础上分析了弹管间隙对火箭初始扰动的影响,并利用二次回归模型求出了弹管间隙的最优解。计算结果表明,初始扰动角和角速度是弹管间隙的二次函数,而最优弹管间隙为该二次函数构成的齐次方程的解。该方法能为火箭发射系统的设计提供参考。     

液体火箭发动机带有隔板和无隔板的两相冷态解数值模拟

本文对液体火箭发动机燃烧室内带有隔板和无隔板的两种情况进行了两相冷态解数值模拟，气相控制方程用欧拉坐标系下的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组描述，液相控制方程在Ｌａ－ｇｒａｎｇｉａｎ坐标系下进行描述。气、液两相作用通过方程组的源项互相耦合，通过计算比较了有隔板和无隔板时发动机内部流场的参数值，表明了隔板体的数值处理方法对计算两相冷态解的有效性。     

利用真空热试验进行轨道卫星污染研究

利用北京卫星环境工程研究所拥有的各种规模的空间环境模拟器在进行真空热试验时,测量卫星表面返回质量流密度和发射质量流密度,通过加文(GarWin)方程组来确定轨道卫星的污染预估,从而解决了当前型号研制中污染控制的问题.     

航天器损伤评估和安全管理中的多尺度力学问题

文章指出航天器的安全性涉及到材料、元件、器件、子系统、系统和整机6个层次。虽然预测和防止失效的任务是在整机层次,但是损伤却是起源于材料的微损伤这一最底层。损伤从原子层次到整机层次的在多个尺度层次上的演化诱致了整体的破坏。由此引入多尺度力学,即研究各种不同长度尺度和时间尺度相互耦合现象的科学。文章试图将力学的宏观运动方程和微结构转变的动力学方程组成统一的方程组,从而得到宏观标量损伤D与微损伤数密度n的关系。文章还指出对于航天器的抗辐射加固问题,微损伤成核和发展的特征时间,与宏观特征空间尺度,以及外加载荷的特征时间之间的耦合关系,即德博拉(Deborah)数,是判断这类损伤演化诱致失效问题的一个关键。     

嵌金属丝推进剂燃烧起始阶段增速特性研究

考虑了使用嵌金属丝装药的发动机工作过程中的多种因素,建立了金属丝与药柱、燃气间相互作用的控制方程组,并对方程组进行了数值求解,得到了嵌金属丝推进剂在燃烧起始阶段,增速比随时间的详细变化过程。结果表明,金属丝的热扩散率、熔点、直径和黑度以各自不同的方式影响金属丝的增速效果,熔点相对较低的金属丝可能会在内弹道曲线的初始阶段形成一个轻微的压强峰,固有燃速较低、着火温度较低的推进剂对嵌入金属丝较为敏感。文中数值方法计算量小,结果可靠,适于工程应用。     

一种基于现有发动机的数值分析方法在液体火箭发动机设计过程中的应用

本文采用多元回归技术对历史上的和目前正在使用的液体火箭发动机求出拟合方程组。发动机的通用方程组仅能用来估计发动机的结构尺寸、干质量、燃科和氧化剂质量流量,以及涡轮泵需要的扬程功率。而本文得到的方程仅需以所要求的发动机推力和推进剂种类作为输入条件,就可以求得新发动机设计的初步近似结果.以现有的发动机为基础建立的设计方程,估计误差大约为±10～20%.本文不仅讨论了数据拟合技术,而且列举了现有的14种液体火箭发动机的30个设计参数.     

Lawrence-Gurney模型对于激光驱动飞片的适用性分析

Lawrence-Gurney模型给出了求解激光驱动飞片速度的方程组。文章从理论和试验两个方面对Lawrence-Gurney模型的适用性进行了分析,发现该模型不能较好地描述厚度为10μm以下飞片的驱动过程。主要原因是Lawrence-Gurney模型对激光与飞片的具体相互作用过程进行了简化假设,依据能量守恒和动量守恒定律建立数学方程组,而这些假设条件和近似求解条件对于厚度小于10μm的飞片均不能很好地成立,导致对飞片速度、能量耦合效率、冲量耦合系数和飞片烧蚀深度随激光能量变化的理论预估结果与试验结果偏差较大,需要进行修正。     

冲压空气涡轮泵供应系统方案设计及特性分析

为了把冲压空气为动力的涡轮泵供应系统从亚燃冲压发动机拓展应用至超燃冲压发动机,基于煤油燃料的双燃烧室冲压发动机(DCR)提出了一种冲压空气涡轮泵供应系统方案。供应系统的设计方案中,对涡轮泵选型、系统的调控策略及取气/排气方案进行了初步设计。同时,建立了供应系统的静态模型,通过系统压力、流量及功率平衡组成非线性方程组,使用牛顿迭代法对非线性方程组进行数值求解,得到了冲压空气涡轮泵供应系统在不同工况下的静态特性。最后,分析了飞行Ma范围在3.5~5.5下涡轮泵的性能和调节的变化规律。结果表明,涡轮所需的空气流量约占DCR发动机捕获空气总流量的3%,取气方案对发动机气动性能影响不大;离心泵的特性参数相对稳定,可以一直处于高效率工况下工作,但系统对增压后的燃料利用不足,造成涡轮功率利用率较低。     

基于CAD二次开发的电加热片快速设计

对基于CAD二次开发的点加热片快速设计进行了研究.针对方形薄膜电加热片,用建立的阻值计算数学模型,对不同并联数的布线电阻值方程组进行了数值求解.根据布线方程组求解结果,基于CAD二次开发技术,对方形电加热片进行了布线设计及绘图.结果表明:本文基于CAD二次开发构造的方法可用于薄膜电加热片快速布线设计.     

基于相平面控制的固体火箭快速精确入轨策略

对于仅能控制推力方向的固体火箭入轨段,为保证火箭以同时满足六轨道根数约束直接入轨,提出一种基于相平面控制的弹道设计方法。首先,证明经过2次恒定方向推力使火箭同时满足位置和速度矢量(即六轨道根数)约束直接入轨的可行性,并给出此时火箭初始点与入轨点状态需要满足的条件。然后,推导火箭初始点及入轨点状态与推力方向角的关系方程组,采用数值方法求解该方程组以获得保证入轨的推力方向角。仿真结果表明,与传统的弹道设计方法相比,提出的方法对推力方向角的解算更为简便,缩短了弹道设计的耗时,且可以使火箭终端位置和速度矢量6个参数同时满足入轨精度要求。     

卫星返回运动方程欧拉角奇点的处理与轨道计算

用欧拉角描述人造卫星返回运动的再入舱姿态有一定的长处,但在相应的运动方程组中存在欧拉角奇点,给计算机运算带来困难。这种奇点仅是数学上的奇点。本文提出,建立对应两组欧拉角的两套运动方程组及其解的互换关系式,在奇点附近转换运动方程组进行运算,可以克服奇点困难。本文给出具体的计算方法。这种方法不增加计算机运算时间,不减低计算精度,对各种条件下的返回轨道计算普遍适用。