双节套筒式延伸喷管运动分析

用质心运动定理对双节套筒式延伸喷管的展开过程进行分析,得到一个二阶常微分方程;然后用四阶龙格－库塔法进行数值求解,可得到延伸出口锥的轴向位移、速度和加速度及展开时间,并以单节套筒式延伸喷管为例;计算结果与试验基本吻合．     

大气波动传播问题的一种无反射数值边界格式

根据无反射边界条件的基本原则,考虑到高阶的数值边界格式可以减小边界误差和虚假反射,利用曲线拟合中最小二乘法的思想,提出了一种高阶的光滑拟合外推边界格式(SFEBS).因为大气运动的控制方程简化后可以归结为一个对流方程,所以作为边界格式之间的比较和检验,在一维情况下用对流方程和数值模拟中具有代表性的波包和激波作为算例,将其与传统的基于Taylor展开思想构造的边界条件(TEBS)进行了比较.计算结果表明,在高阶情况下,用SFEBS计算波动传播问题的虚假反射约为用同阶TEBS计算的1/6,说明高阶的SFEBS可以作为一种很好的无返射数值边界格式.为今后模拟大气波动提供了一种很好的数值边界格式.      

圆柱形静电分析器的数值模拟

对美国DMSP系列卫星载60°圆柱形分析器的性能进行了数值模拟:以同时过松弛迭代法解Laplace方程,得到分析器内的势分布,由荷电粒子在分析器势场中的运动轨道以及传输概率的MontCarlo积分得出分析器的各种响应和几何因子．变动与分析器几何相联系的输入量,计算机程序可移用于其它圆柱形分析器.      

微波干涉仪测量羽流电子密度反演方法

利用微波干涉仪进行推进器羽流密度测量的反演方法研究.重点针对羽流等离子体密度反演的Abel逆变换数值求解,采用离散化算法、三次样条插值算法、Hankel-Fourier算法以及Nestor-Olsen算法等四种方法进行仿真分析和比较.通过对这四种方法在无噪声和有噪声两种情况下反演误差的比较,找出适用于实际干涉仪数据处理的方法;通过对采样点数和反演之间的关系分析,给出解决欠采样问题的方法和数值参考;在分析比较的基础上,总结数值处理的必要步骤.利用模拟实验测试数据进行验证并对反演结果进行分析.      

空间站大型伸展机构柔性多体动力学数值计算

对用ｋａｎｅ方法建立的空间站大型伸展机构柔性多体动力学数学模型,从数值计算方法和软件实现方面进行了研究。文章采用局部附着浮动坐标系描述各构件的刚体大范围伸展运动;用截断的子结构正则模态描述各构件的柔性变形;用求解刚性方程的ＴＲＥＡＮＯＲ方法进行数值积分;用混合法对系统的微分－代数动力学方程组进行求解。还对空间站伸展机构多体系统的约束方程建立、独立坐标选取和动力学方程的ｓｔｉｆｆ病态问题进行了讨论。并针对空间站伸展机构在驱动展开过程中以及展开完成后的动力学响应进行了算例计算。     

数值修约的通用数学模型及其一个VB6．0实现过程

按“四舍五入法”规则修约数值时,其舍入误差不是等概率分布的,它仅仅是GB/T 8170-1987《数值修约规则》的一个“有缺陷的”特例。根据GB/T 8170-1987《数值修约规则》相关内容,建立了数值修约的通用数学模型,用VB6.0编写了其实现过程。经测试,证明该过程正确可靠,它对测控系统中的数值修约有一定的参考价值。     

三角域上的二元2n次插值多项式

本文用调配函数法构造了三角域上的二元2n次插值多项式,讨论了逼近阶,并给出了一种函数展开。     

航天器太阳帆板展开过程的最优控制

本文讨论航天器太阳帆板展开过程中主体姿态的最优控制问题. 利用多体动力学与最优控制理论建立数学模型, 考虑系统的非完整约束特性, 提出太阳帆板展开过程的最优控制算法. 通过数值仿真, 表明该方法对太阳帆板展开姿态控制的有效性.      

大型复杂航天器的柔性附件展开的动力学分析

分析了复杂航天器在柔性附件展开过程中的几个关键问题,并提出了相应的解决方案,编制了应用软件。通过数值模拟详细分析了其附件的展开过程,讨论的内容包括卫星姿态角和姿态角速度的变化规律、柔性附件的展开运动及其弹性振动等。     

太阳电池阵展开运动的数学模拟

分析了太阳电池阵的展开运动过程,并根据力学原理建立了模拟太阳电池阵展开运动规律的数学模型.根据某型卫星的数据,进行了地面模拟展开过程计算,并给出了模拟计算和地面展开运动试验的比较.结果表明,这一简化的模型和数值方法对于刚性太阳电池阵展开运动模拟的可行性.     

Numerical study of traveling ionosphere disturbances with vertical incidence data

High frequency ionosphere vertical sounding, as an important and representative application for detecting the ionosphere and studying the characteristics of radio propagation, can be utilized to monitor the ionosphere continuously variation and to acquire the ionosphere asymmetrical features of diverse scale above the ionosphere vertical sounding stations. This is a first article on real time application of numerical methods to obtain the parameters of traveling ionosphere disturbances (TIDs) using vertical incident ionograms. In this paper, the distribution of ionosphere electron density with TIDs is constructed using a background ionosphere model superimposed a perturbation theory model. The background ionosphere electron density is modelled by the inversion of vertical incident ionograms which are observed before the appearance of the disturbance. Based on the fourth order Adams-Bashforth-Moulton (the so-called ABM) predictor corrector method, instead of Runge-Kutta method, the fast digital ray tracing method is established. According to process of the disturbed trace simulation and parameters inversion, the characteristic parameters of ionosphere disturbance at different detection time can be obtained in real time. The numerical analysis of TIDs is then captured completely.     

基于CLYC闪烁体的中子能谱测量及反演方法

空间中子是影响航天器和航天员安全的重要辐射要素之一。优化中子探测器,提高测量效率,提升反演精度是中子测量的难点。中国空间站将搭载一种基于新型中子探测材料Cs2LiYCl6:Ce（CLYC）闪烁体的中子探测器,该探测器具有同时测量热中子和快中子,以及探测效率高等特点。针对该新型探测器的中子能谱反演,分析了不同能量中子在该探测器中的响应特点,分析了中子反演常用的概率迭代法和非负最小二乘（NNLS）法的优缺点,考虑到这2种方法在CLYC探测器反演应用中的不足,提出了基于增广矩阵的非负最小二乘（AM-NNLS）法。数值实验结果表明:AM-NNLS法具有反演运算效率高和反演相对误差小的特点,验证了所提方法的有效性。      

Durbin法在阻尼梁动响应求解中的运用分析

运用铁木辛柯梁理论和K-V阻尼理论，研究了非比例阻尼梁在冲击载荷作用下的频域振动求解方法。推导采用了传统拉普拉斯正变换和基于Durbin公式的拉普拉斯反变换策略（统称拉普拉斯法），发展了阻尼梁系统的动力学方程解法。拉普拉斯法的推演同时涵盖了3种典型的梁边界条件，具有广泛的适用性。数值法的验证采用了特殊构造的比例阻尼点条件，并与基于模态叠加法的求解结果进行了对比分析，且数值算例充分考虑了数值参数和系统参数的影响。计算结果表明:在不同边界条件和受载状态下，拉普拉斯法与模态叠加法均能合理地计算出基本阻尼梁系统的动响应曲线，且两者的求解精度保持在同一量级；同时，捕捉到拉普拉斯法的求解精度会受到系统长细比等参数的影响。拉普拉斯法具有比传统实、复模态叠加法更易操作的特性，但其精度受到了算法固有参数和阶跃外载型式的影响，稳定性仍需进一步提高。      

高超声速飞行器BTT非线性控制器设计与仿真

高超声速飞行器的气动特性比一般的飞行器更为复杂,选用BTT(Bank-to-Turn)技术,即倾斜转弯技术可以满足其对于气动外形的要求,但随之给动力学系统带来了快时变、严重非线性及强烈耦合的特点.针对高超声速飞行器倾斜转弯非线性控制器的这一特点,采用了一种更为有效的非线性系统的控制方法,即非线性动态逆技术.首先建立了高超声速飞行器的非线性数学模型,然后根据奇异摄动理论将动力学系统的受控状态变量分为快变量和慢变量2部分,应用非线性动态逆理论分别对快逆回路和慢逆回路进行设计,其中慢逆回路控制器的输出作为快逆回路控制器的输入指令,最后对于所设计的系统在高超声速下的倾斜转弯运动进行了仿真验证.仿真结果表明该控制系统可以实现倾斜转弯,并可以满足高超声速飞行器稳定飞行的要求.      

马尔可夫更新方程的迭代算法研究

冷储备系统的动态结构变化使得系统可用度分析必须基于瞬态求解,针对频域两部件冷储备系统马尔可夫更新方程解算困难的问题,提出基于时域伏尔特拉方程的迭代算法,证明了该迭代方程的唯一性和收敛性,并给出迭代算法的误差分析.理论和实例分析表明,当冷储备系统部件的工作寿命和维修时间分布是有界函数时,利用迭代算法可以容易实现时域马尔可夫更新方程的求解.      

A high order method for orbital conjunctions analysis: Sensitivity to initial uncertainties

A high order method to quickly assess the effect that uncertainties produce on orbital conjunctions through a numerical high-fidelity propagator is presented. In particular, the dependency of time and distance of closest approach to initial uncertainties on position and velocity of both objects involved in a conjunction is studied. The approach relies on a numerical integration based on differential algebraic techniques and a high-order algorithm that expands the time and distance of closest approach in Taylor series with respect to relevant uncertainties. The modeled perturbations are atmospheric drag, using NRLMSISE-00 air density model, solar radiation pressure with shadow, third body perturbation using JPL’s DE405 ephemeris, and EGM2008 gravity model. The polynomial approximation of the final position is used as an input to compute analytically the expansion of time and distance of closest approach. As a result, the analysis of a close encounter can be performed through fast, multiple evaluations of Taylor polynomials. Test cases with objects ranging from LEO to GEO regimes are considered to assess the performances and the accuracy of the proposed method.     

Gravity field models from kinematic orbits of CHAMP,GRACE and GOCE satellites

The aim of our work is to generate Earth’s gravity field models from GPS positions of low Earth orbiters. Our inversion method is based on Newton’s second law, which relates the observed acceleration of the satellite with forces acting on it. The observed acceleration is obtained as numerical second derivative of kinematic positions. Observation equations are formulated using the gradient of the spherical harmonic expansion of the geopotential. Other forces are either modelled (lunisolar perturbations, tides) or provided by onboard measurements (nongravitational perturbations). From this linear regression model the geopotential harmonic coefficients are obtained.     

蒸发液滴空间实验研究的图像反馈控制系统

提出了一种利用图像反馈控制系统测量液滴蒸发速率的方法.该系统主要由图像采集、图像处理、反馈控制三个部分组成.其工作原理是对CCD采集到的液滴图像进行分析,得到液滴的物性参数,利用控制注液器动态注入来维持液滴的大小,由此得到单位时间内的注入量,即液滴蒸发速率.或通过对液滴图像几何尺寸的计算,得到单位时间内液滴的变化量,进而得到液滴的蒸发速率.为了从动态变化的图像中准确找到液滴的轮廓,简单介绍了基于拉普拉斯方程,运用牛顿法和龙格库塔法等数值方法,对图像轮廓进行拟合的算法理论.以拟合得到的液滴轮廓为基础,利用数值积分计算液滴的表面积和体积,从而由体积的变化量确定蒸发速率.介绍了图像反馈控制系统的软件结构和硬件结构.并给出了利用其进行液滴蒸发测量得到的实验结果.该系统是为我国SJ-10返回式卫星上蒸发与流体界面效应空间实验研究项目专门开发的.      

用整数小波变换进行无失真SAR图像压缩

介绍了利用提升方案(lifting scheme)构造小波变换的3个主要优点,给出了几个整数到整数可逆小波变换对,讨论了SAR图像的特点,提出了一种利用整数到整数小波变换进行无失真SAR图像压缩的方法.该方法具有速度快、可进行同址运算和能进行分级传输等特点.比较了对ERS-1单视SAR图像和IEEE标准测试图的压缩结果,表明性能要略优于DPCM(Differential Pulse Code Modulation)方法,且只须对单视SAR图像作一级小波分解,分解级数增加会使压缩效果变差.     

一种跨声速定常流场求解加速方法

跨声速定常流场的隐式求解相当于使用牛顿迭代法求解一个非线性方程组。为满足牛顿迭代收敛性的要求,通常需要对所求解问题进行全局化处理。在同伦延拓的框架内,提出了一种基于拉普拉斯算子的方程延拓方法,提高了定常流场隐式求解收敛速度。针对定常流场通常初始化为均匀来流的特点,一方面利用拉普拉斯算子的椭圆性加快边界条件信息向流场内部的传播,另一方面利用拉普拉斯算子的线性和正定性改善延拓问题的正则性,综合两者增加拟牛顿算法的稳定性,提高可用CFL数,最终达到提高流场求解效率的目的。由于流场问题的复杂性和非线性,难以通过理论分析得出先验的最优非线性求解策略。因此,通过无黏NACA0012翼型、湍流RAE2822翼型和三维ONERA M6机翼等算例的数值实验,研究了拉普拉斯项参数对收敛效率的影响,给出了效率较优的参数组合,验证了本文方法在跨声速情况下相对于经典伪时间推进法可以节约20%以上的CPU计算时间。