基于EFFD方法的自然层流短舱优化设计

采用extended free-form deformation(EFFD)方法研究了自然层流(natural laminar flow,NLF)短舱的气动外形优化设计方法.使用基于Bernstein基函数的EFFD方法完成了NLF短舱剖面的参数化,利用基于k-εSST(shear stress transport)两方程湍流模型的γ-θ转捩模型进行自然转捩预测,结合EFFD、一种混合动网格方法、Kriging代理模型和改进的粒子群算法(particle swarm optimization,PSO)建立了针对NLF短舱气动外形的优化设计框架.采用该框架分别对通气NLF短舱和带动力NLF短舱进行优化设计.单独通气NLF短舱优化结果的外表面实现48%的层流,阻力系数比初始通气NLF短舱减小了0.0003.带动力NLF短舱的优化结果外表面保持了41%的层流.这些结果表明采用相关技术建立的优化设计框架在NLF短舱设计中具有一定应用价值.     

交会对接任务碰撞规避策略制定研究简

针对交会对接任务目标飞行器与追踪器轨道运行特性,综合考虑规避策略计算方法与工程实际相结合的问题,提出高度规避、时间规避以及与正常轨控相结合的碰撞规避策略计算方法等三种空间目标碰撞规避策略计算方法.高度规避计算方法采用了Lambert飞行原理,用简化二体开普勒模型取代高精度轨道预报方法,迭代求解规避机动速度增量,实现了通过约束过交点与目标径向距离差得到速度增量的最优解;时间规避计算方法通过轨道周期与速度增量的关系,实现了通过约束过交点与目标的时间差得到速度增量的最优解;与正常轨控相结合的碰撞规避策略计算方法,在正常控制考虑冗余控制量的基础上,对控制策略的控制开始时间或沿迹方向的速度增量进行较小的修正,使两者通过碰撞点的时刻或径向距离错开,达到碰撞规避的目的,该方法不仅可以节省燃料、而且对任务的影响较小.通过对三种空间目标碰撞规避策略计算方法仿真分析结果表明,完全适用于交会对接任务,可为我国载人航天任务飞行安全提供技术保障.     

基于人工免疫算法的载人登月任务地月转移双脉冲中止策略研究彭 坤,果琳丽,向开恒,王 平,杨 雷

载人登月任务中,任务中止策略设计是确保航天员安全返回的重要基础。首先结合"星座"计划飞行方案分析了载人登月任务各飞行阶段的中止策略;其次针对地月转移巡航段进行了双脉冲中止策略设计,以速度增量数值、方位角以及变轨时间间隔为控制变量,加入轨道同向、近地点高度、偏心率以及飞行时间约束,提出双脉冲变轨计算流程;最后采用人工免疫算法对该问题进行了求解和优化。仿真算例表明,双脉冲中止策略存在多组解,其全局分布特性为:飞行时间越短速度增量需求越大;飞行时间相近时,大偏心率中止轨道对应的速度增量小;故障点离地月加速点越近,所需速度增量越小。同时也验证了人工免疫算法求解双脉冲中止策略问题的有效性。     

无控武器对地攻击精度分析

针对目前无控武器对地攻击精度影响因素来源单一,分析影响精度的多种误差源。建立直升机无控武器命中精度分析模型,针对平飞方式下的无控武器命中精度问题,采用蒙特卡洛法多次仿真。针对命中精度与各个误差源的关系方程,提出采用多元线性回归法对仿真结果进行分析,给出误差源对命中精度的影响程度。对关系方程和误差影响程度进行显著性检验,使其更加准确、可信。     

k/N（G）结构系统初始备件配置方法

针对现有备件冗余度建模中无限呼叫总体假设不合理、供应可用度计算困难的问题,提出一种k/N（G）结构系统初始备件配置方法。首先,采用有限源排队论获取备件短缺数概率分布。然后,给出计算k/N（G）结构系统供应可用度的方案,并证明了方案的正确性。该方案在计算随机挑选的一套系统未因备件短缺而停机的概率基础上,利用k（或更多）套系统良好的二项分布概率估计系统供应可用度。最后,应用边际分析算法确定各项备件的数量。实例分析表明：该方法能优化k/N（G）结构系统各组成部件的库存数量,且与单项法相比具有明显的优势,能以较少的备件费用达到相同的供应可用度目标。     

定长补偿电火花铣削加工中锥形电极的形成及稳定性

电火花铣削可以利用简单电极加工复杂工件,对铣削过程中电极损耗进行补偿是保证工件加工精度的重要措施。使用圆柱电极进行定长补偿铣削加工,加工过程中电极端面会形成一个圆锥形。为保证补偿精度,对电极圆锥形端面的形成及稳定性进行了研究。通过对定长补偿下工件被加工情况的研究,阐述了圆锥形电极形成的具体过程;并且论证了在加工进入稳定阶段以后,圆锥端面角度会基本保持在一个恒定值,研究了初始加工深度和补偿长度对锥形过渡阶段的影响;在电极直径、电参数一定情况下,验证了加工模型中锥形角度与目标加工深度的函数关系,实验结果与模型计算结果相差2.1%;最后给出了加工实例,并获得较好的加工效果。     

Determination of thermal expansion coefficients for unidirectional fiber-reinforced composites

In the present work, the coefficients of thermal expansion（CTEs） of unidirectional（UD）fiber-reinforced composites are studied. First, an attempt is made to propose a model to predict both longitudinal and transverse CTEs of UD composites by means of thermo-elastic mechanics analysis. The proposed model is supposed to be a concentric cylinder with a transversely isotropic fiber embedded in an isotropic matrix, and it is subjected to a uniform temperature change. Then a concise and explicit formula is offered for each CTE. Finally, some finite element（FE） models are created by a finite element program MSC. Patran according to different material systems and fiber volume fractions. In addition, the available experimental data and results of other analytical solutions of CTEs are presented. Comparisons are made among the results of the cylinder model,the finite element method（FEM）, experiments, and other solutions, which show that the predicted CTEs by the new model are in good agreement with the experimental data. In particular, transverse CTEs generally offer better agreements than those predicted by most of other solutions.     

旋转爆震发动机爆震波建立过程实验研究

为了深入研究旋转爆震发动机爆震波建立过程及形成机理,采用小能量火花单次点火的方式进行了一系列旋转爆震发动机起爆实验。发动机采用环缝-喷孔对撞式掺混方式,燃料为H2,氧化剂为空气,实验成功起爆旋转爆震波,并连续旋转稳定传播,爆震波传播频率为5.09~6.45k Hz,传播速度为1286~1644.8m/s。在发动机稳定工作过程中,集气腔与燃烧室相互影响,二者处于平稳的动态平衡。其次,通过对旋转爆震波起爆过程详细分析发现,点火形成的初始火焰在环形燃烧室经历一个类似DDT的火焰发展过程,成功转变为爆震波,且从点火到爆震波建立之间的火焰发展传播过程和时间间隔均表现出很强的随机性。此外,为验证小能量火花点火的可靠性,还进行了小能量点火重复性实验,发现在稳定工况条件下采用小能量点火成功率最高可达100%,各组旋转爆震波传播速度在1440m/s附近波动。     

一种基于历史标校数据的移相值预测方法

地面测控设备在执行测控任务前,需要对和差接收通道的相位差进行标校(简称校相),以满足测角分系统的目标自跟踪要求.针对因标校设备故障而无法进行校相的问题,通过使用数理统计的多元回归分析方法,对积累的历史标校数据进行分析,确定了设备工作的环境温度对跟踪接收机移相值有显著影响.在此基础上,利用滑动窗口多项式拟合方法得到拟合移相值,并与最近一次的标校移相值进行时变加权,可以得到基于历史标校数据的移相值预测方法.仿真试验表明该方法具有较好的预测精度,可将预测移相值作为跟踪接收机参数,提高测控设备应急测控能力.     

一种星敏感器-陀螺组合定姿的实时在轨标定方法

研究了一种星敏感器一陀螺组合定姿方式中的姿态敏感器误差的实时在轨标定方法。首先,选择直观的欧拉角作为姿态描述参数,根据星敏感器和陀螺的测量原理建立星敏感器一陀螺在轨标定的测量方程和状态方程,并以此建立数学模型。其次,采用简单高效的EKF（ExtendedKalmanFilter,扩展卡尔曼滤波）作为估值算法,进行了在轨标定数值仿真。对于航天器姿态定向中出现的姿态角和星敏感器安装角之间的耦合问题,通过在特定姿态通道上施加简单姿态机动实现了解耦。数值结果表明,该实时在轨标定方法,尤其是所提出的姿态角和星敏感器安装角解耦策略,可以实现对航天器姿态的实时精确估计以及对星敏感器安装误差、陀螺常值漂移和相关漂移等误差的实时在轨标定。该方法可用于航天器姿态测量设备的实时在轨标定和航天器姿态的高精度实时确定。