任务准备期内的军用飞机瞬时可用度

军用飞机是一种复杂的可修系统,其瞬时可用度计算是装备综合保障领域研究的难点。本文构建了基于作战任务和飞机技术状态驱动的随机维修网络(SMTN),SMTN中各项活动的维修时间服从不同分布,定义了SMTN的矩母函数并求解系统维修度函数的数字特征,将得到的期望和方差分别代入不同分布形式函数进行验证,从而选定符合实际情况的SMTN的维修度函数,并对任务准备期内的军用飞机可用度进行仿真。在计算过程中基于蒙特卡罗技术设计了一种新方法,在提高计算效率的同时,可以较好地反映维修保障过程的逻辑性、拓扑性和随机性,仿真结果可以真实反映装备在任务准备期内瞬时可用度的波动规律。该方法可以广泛应用到各种复杂的可修系统,为装备的使用决策提供依据,为保障系统的保障效能评估提供量化指标,对于综合保障领域"动态"指标的研究有着一定探索价值。     

探测器触地关机软着陆稳定性分析

深空探测器软着陆过程中发动机未能正常关闭将极大降低其软着陆稳定性,因此引入了触地关机着陆方案。针对触地关机方案,建立了腿式探测器动力学分析模型及发动机推力控制模型。考虑着陆环境与探测器着陆状态的不确定性,采用蒙特卡罗模拟分别对探测器在主发动机关闭失败时、触地关机着陆方案下、带有姿态控制的触地关机着陆模式下的软着陆稳定性进行了分析。采用描述性采样方法抽取样本点,基于均值估计相对误差建立模拟终止准则,在保证模拟精度的前提下提高蒙特卡罗模拟效率。计算并对比了3种着陆方式下探测器软着陆可靠度,触地关机方案可在主发动机关闭失败后将探测器稳定着陆可靠度提升11.6%,姿态控制的引入可进一步将安全着陆的可靠度提升9.7%。      

关于蒙特卡罗模拟抽样方法的研究

提出重要性抽样与拉丁超立方抽样相结合的MonteCarlo模拟计算方法。算例证实该复合法的计算效率高。     

Investigations on temperature distribution of satellite surfaces affected by solar absorptivity

Based on the Monte Carlo ray-tracing method, the network coefficients of thermal network model describing the radiation heat transfer among satellite surfaces is solved by considering the surface material optical characters. It is superiority to the conventional Gebhart’s method in view of the grey body and the diffuse reflection assumptions. The zone leveling method is used to discrete the governing equations and the solar absorpivity is separated and considered to be an important correction parameter. Effects of the solar incidence round angle, the zenith angle and the ratio of absorpivity to emissivity (RAE) on temperature distribution are numerically simulated and discussed in detail. The higher or the lower the RAE may be lead to the alternative heating and cooling tend with a larger heating or cooling velocity of main body surfaces than the solar array surfaces. Furthermore, maximum temperature of main body is almost larger than solar arrays. Under the same RAE, solar incidence angle make a great effect on the uniform character of temperature distribution.     

液体火箭发动机喷管壁面辐射热流的数值计算

为火箭发动机喷管设计及其冷却系统的设计提供有效参考数据作为应用背景,计算喷管壁所受的红外辐射热流。分别推导了有限体积法和反向蒙特卡罗法适用于不含粒子条件下喷管壁面所受红外辐射热流的计算公式,并针对不同的边界条件设定,重点分析了有限体积法的特性,并提出了有限体积法的局限性。最后,以某三维双圆弧喷管模型为例,用有限体积法和反向蒙特卡罗法自主编写程序计算在不同边界设定下壁面受到的红外辐射热流密度。通过计算结果对有限体积法的局限性进行了论证,同时肯定了反向蒙特卡罗法在壁面辐射热流计算上的优势。     

低能质子辐照下F46/Ag光学性能退化及损伤模型

利用空间综合辐照环境模拟器对F46/Ag二次表面镜进行低能质子辐照试验,发现30 keV质子辐照会导致F46/Ag的光学性能发生显著退化,波长为350 nm处的光谱反射率变化值随累积辐照通量呈指数关系增加,太阳吸收比变化值随累积辐照通量线性增加。通过Monte Carlo模拟计算得到不同深度电离能损和位移能损的分布,发现质子辐照F46/Ag主要通过电离进行能量的传递。当质子能量为30 keV时,F46表面获得质子传递能量的最大值。量子化学计算结果表明30 keV质子辐照F46/Ag时材料表面C—F键和C—C键断裂所需吸收的能量远小于质子传递给F46表面的能量。XPS分析表明质子辐照后材料表层发生了一系列断键和重组反应,生成了自由基和分子片段,表面发生脱氟和碳富集,与量子化学计算结果相吻合。     

DSMC量子动理学化学反应模型的数值模拟

为评估直接模拟蒙特卡罗（DSMC）方法中新近发展的一种量子动理学（QK）化学反应模型,在自研高超声速稀薄气体流动数值模拟软件RariHV中实现该模型,并进行了单个绝热单元内的热泳平衡/非平衡反应测试。结果表明,数值解与理论解吻合良好,QK模型具有很好的预测化学反应流动的能力。为进一步验证QK模型在真实流动中的表现,还计算了高超声速绕圆柱的化学反应流动,QK模型计算的驻点线上温度和组分分布都与总碰撞能（TCE）模型结果吻合良好。相较于TCE模型,QK模型的优势是不再依赖反应速率系数的实验数据,在深空探测等化学反应数据缺乏的领域具有一定的应用前景。     

SiC/SiC mini复合材料拉伸性能分散性的数值仿真方法

建立了一种能够高效计算连续纤维增强SiC/SiC陶瓷基复合材料力学性能统计分布函数的数值仿真方法。建立SiC/SiC mini复合材料简化二维胞元模型,输入组分材料力学性能参数的概率分布函数,得到胞元的力学性能。将胞元作为mini复合材料二维模型的单元对模型进行网格划分,并赋予其材料性能参数,从而得到mini复合材料的拉伸应力-应变曲线。经多次重复计算得到了mini复合材料拉伸应力-应变曲线的概率分布函数。从结果中发现,SiC/SiC mini复合材料拉伸强度满足威布尔分布,且相比于纤维和基体强度的模数,复合材料拉伸强度的模数增大,即分散性减小,同时尺度参数降低,表明mini复合材料出现概率最大的拉伸强度值变小;随着纤维和基体强度概率分布威布尔模数的增大,mini复合材料拉伸强度的尺度参数有增大的趋势。      

雷暴电场对LHAASO探测面宇宙线次级粒子能量的影响

采用Monte Carlo方法,通过模拟研究了不同强度雷暴电场对LHAASO探测面宇宙线次级电子能量的影响.在电场作用下,电子的能量分布发生了变化.在低能段总电子数目增加明显,而在高能段电场的影响不明显.当电场强度为1700V·cm-1（大于逃逸电场）时,能量<120MeV的电子被加速,能量<60MeV的总电子数目呈指数增长（增幅高达约2252%）,雷暴电场对次级粒子的加速机制与相对论电子逃逸雪崩机制（RREA）相符.当电场强度为1000V·cm-1（小于逃逸电场）时,能量<70MeV的电子被加速,其数目明显增加,但是增幅（约86%）远小于逃逸电场时的幅度.对电场强度小于逃逸电场时的雷暴电场加速宇宙线次级粒子的物理机制进行了讨论.研究结果可为理解LHAASO实验数据特点以及研究雷暴期间宇宙线强度的变化等提供参考.      

直升机高速四点接触球轴承接触疲劳可靠性评估方法

为了准确高效地分析高速四点接触球轴承的接触疲劳可靠性,将蒙特卡罗随机有限元(MC SFEM)与应力-强度干涉理论相结合提出了一种有效的可靠性分析方法.在高速、高温和弹流润滑（EHL）的复杂工况下,利用有限元法(FEM)建立了基于热 结构耦合分析的有限元模型.考虑到轴承材料属性、润滑油参数及接触疲劳强度的随机性对可靠性的影响,在借助MC SFEM和K S (Kolmogorov Smirnov)检验确定了高速四点接触球轴承疲劳接触应力分布类型的条件下,根据应力-强度干涉理论建立了高速四点接触球轴承接触疲劳可靠性评估方法.结果表明:与传统的MC SFEM方法相比,该方法的耗时仅是传统MC SFEM的0.243％,两者之间的可靠性误差仅为1.48％,进一步验证了该方法的有效性.