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

提出重要性抽样与拉丁超立方抽样相结合的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.     

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

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

基于多体仿真模型的运动机构可靠性仿真试验系统研究

以蒙特卡罗法为理论核心,利用多体运动学和动力学分析软件--LMS Virtual.Lab作为运动机构可靠性分析平台,通过Visual Basic.NET程序设计语言建立用户界面并实现对LMS Virtual.Lab的调用,以进行运动机构的可靠性仿真试验,采用得到的随机变量数据库和结果数据库对运动机构进行典型失效模式的可靠性分析,从而形成了一套可以实现运动机构的可靠性仿真与分析功能的可靠性仿真试验系统.以对心曲柄滑块机构为例,进行了几何尺寸误差和运动副间隙误差影响下的机构运动精度可靠性分析;以某锁机构为例,研究多失效模式下动力学特性对于机构定位可靠性的影响并进行了灵敏性分析.通过这两个实例验证该可靠性仿真试验系统进行机构可靠性仿真试验和可靠性分析的可行性和有效性.     

低能质子辐照下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％,进一步验证了该方法的有效性.      

姿控发动机喷管内流场的DSMC并行算法

在分析了直接仿真Monte-carlo（DSMC）方法内在并行度的基础上，针对工作站群机并行系统（NOWS）分布主存的特点，构造了适合于DSMC仿的并行算法模型，在PVM环境下将其应用于姿控发动机喷管内流场的DSMC仿真，结合并行算法效率评价准则对如何提高并行计算效率进行研究。仿结果表明了实施DSMC并行计算的必要性与有效性。