磁选态单束铯束管束光学参数的模拟计算和设计

为提高铯束管跃迁信号强度,在束光学参数设计中,用蒙特卡罗方法模拟计算铯原子通过率和未跃迁比例。通过改变模拟计算中的束光学参数,分析铯原子通过率、未跃迁比例和倍增器首极电流的变化趋势,得到了主要束光学参数的合理取值范围。在铯束管实际装配和试验中使用束光学参数设计结果,有效地增加了倍增器首极电流。     

考虑轨道摄动的外热流计算分析

为改进传统外热流算法未考虑外热流经多次反射后的吸收和国外商业软件未考虑轨道摄动及太阳矢量变化对外热流的影响等不足,重新定义了外热流角系数和外热流辐射传递因子,综合SGP4轨道模型和Gebhart方法,给出了考虑摄动等因素的辐射传递因子和轨道外热流计算模型,以及相应的计算流程。通过仿真分析论证了考虑多次反射及摄动等因素的重要性。     

多种失效模式相关的机构运动可靠度计算方法

机构运动中必然存在多种相关失效模式,多为影响机构运动的各类随机变量函数.为了在机构运动可靠性分析中有效考虑多种相关失效模式的影响,建立了考虑多种失效模式相关条件下的机构运动可靠度计算模型,并实现了对多种失效模式相关条件下的机构运动可靠度计算模型的求解,为考虑多种相关失效模式下的机构运动可靠性分析提供了一种有效途径.以6408滚动轴承为例,进行100 000抽样计算,考虑失效模式相关条件下运动可靠度为95.998 9%.该方法精度较高,具有一定的理论意义和应用价值.     

Disturbance observer based model predictive control for accurate atmospheric entry of spacecraft

Facing the complex aerodynamic environment of Mars atmosphere, a composite atmospheric entry trajectory tracking strategy is investigated in this paper. External disturbances, initial states uncertainties and aerodynamic parameters uncertainties are the main problems. The composite strategy is designed to solve these problems and improve the accuracy of Mars atmospheric entry. This strategy includes a model predictive control for optimized trajectory tracking performance, as well as a disturbance observer based feedforward compensation for external disturbances and uncertainties attenuation. 500-run Monte Carlo simulations show that the proposed composite control scheme achieves more precise Mars atmospheric entry (3.8?km parachute deployment point distribution error) than the baseline control scheme (8.4?km) and integral control scheme (5.8?km).     

结合跳点引导的无人机随机搜索避撞决策方法

针对无人机在城市空域环境和密集交通流下的避撞决策问题，提出马尔科夫决策过程（MDP）和蒙特卡洛树搜索（MCTS）算法对该问题进行建模求解。蒙特卡洛树搜索算法在求解过程中为保证实时性而使其搜索深度受限，容易陷入局部最优，导致在含有静态障碍的场景中无法实现避撞的同时保证全局航迹最优。因此结合跳点搜索算法在全局规划上的优势，建立离散路径点引导无人机并改进奖励函数来权衡飞行路线，在进行动态避撞的同时实现对静态障碍的全局避撞。经过多个实验场景仿真，其结果表明改进后的算法均能在不同场景中获得更好的性能表现。特别是在凹形限飞区空域仿真模型中，改进后的算法相对于原始的蒙特卡洛树搜索算法，其冲突概率降低了36%并且飞行时间缩短47.8%。     

高能(450 kV)X射线参考辐射场附加过滤研究

利用MCNP程序对450kV X光机建立模型，计算各参考辐射质下的原级谱、平均光子能量。对ISO 4037-1中已有范围的管电压与附加过滤片厚度关系曲线进行二次拟合，确定模拟厚度范围。利用MCNP计算各附加过滤组合下的平均光子能量值，并与外推法得出的平均光子能量值进行对比，将相对偏差最小时的附加过滤作为最终确定值。从部分能量点的模拟与测量结果的对比可看出，模拟得出的平均光子能量在2%以内与PTB及中国计量院测量结果保持一致，附加过滤最大相对偏差为18%。     

Ni,Ti及NiTi合金低能电子散射的Monte Carlo模拟

采用MonteCarlo方法,模拟了低能电子束(能量E0≤5keV)作用下Ni,Ti及NiTi合金固体中的电子散射,分析了3种金属/合金中散射电子的能量与空间分布。研究表明:(1)入射电子能量E0越小,各材料的背散射电子(BSEs)深度、吸收电子(AEs)深度及散射电子能量损失的深度均越小,分辨率越高;(2)相同E0下,各材料的BSEs,AEs及能量损失深度均有NiNiTi>Ti;(3)NiTi合金的BSEs深度、AEs深度及能量损失的深度和分辨率分别近似于Ni,Ti固体相应量的平均值。     

基于蒙特卡洛方法的冰晶撞击特性计算

航空发动机摄入冰晶易发生内部结冰、造成发动机推力损失,严重时引发航空事故。冰晶结冰数值模拟是研究冰晶结冰的重要手段,而撞击收集系数计算是冰晶结冰数值模拟的关键环节。基于NNWICE平台开展冰晶结冰前期的阻力系数、运动轨迹和撞击收集系数的数值模拟与计算研究,实现了基于蒙特卡洛方法的冰晶撞击收集系数的准确计算;在此基础上发展了拉格朗日框架下冰晶粒子运动、传热传质耦合数值计算方法;通过完成二维NACA0012翼型的相关算例,对比分析了不同阻力模型和冰晶相变对撞击收集系数的影响。结果表明阻力模型对冰晶撞击收集系数影响较小,冰晶相变过程会对冰晶运动轨迹及撞击收集系数产生显著影响,为进一步发展冰晶结冰数值模拟计算方法奠定了基础。     

基于蒙特卡罗法的航空发动机空气系统稳态算法优化

针对目前航空发动机空气系统稳态算法中收敛性依赖初值的问题,将蒙特卡罗方法与流体网络法综合应用到空气系统可压缩流体一维网络计算中,提出了一种新的计算方法Monte Carlo-Fluid Network（MC-FN）。该方法将空气系统简化为由节点和元件组成的网络,借助蒙特卡罗方法获得空气系统内各节点压力分配,再根据空气系统中各元件流阻特性和换热特性计算流量、温度。计算中通过将游动次数比较少的蒙特卡罗方法的计算结果作为流量残差法节点压力、温度的初始值,实现快速求得精确收敛解。与流量残差算法相比,MC-FN方法计算精度不变,收敛速度提升了66.5%;与线性求解法相比,MC-FN方法的计算精度提升了25.2%,收敛速度提升了43.8%。