数值模拟喷焰2.7微米红外辐射特性

以固体火箭发动机尾喷焰的红外探测为主要应用背景,利用反向蒙特卡罗法和有限体积法计算了火箭尾喷焰红外辐射特性,获得了尾喷焰在探测波长为2.7μm下的红外辐射的空间、光谱分布特性,两种方法的计算结果也非常接近,具有重要的理论意义和实际应用背景。     

考虑基体开裂的陶瓷基复合材料应力-应变曲线模拟方法及验证

提出了蒙特卡罗方法模拟陶瓷基复合材料基体随机开裂过程,采用剪滞模型描述了复合材料出现损伤时细观应力场,并推导得到了考虑基体开裂时复合材料拉伸应力-应变曲线计算公式.开展了室温环境下C/SiC复合材料的单轴拉伸试验,并将理论预测应力-应变曲线与试验结果进行对比.同时,采用该方法对SiC/CAS,SiC/Si₃N₄复合材料应力-应变曲线进行了模拟,并与国外提供的相关试验数据进行比较,发现两者吻合得较好,从而证实了蒙特卡罗法可有效地模拟考虑基体随机开裂过程的陶瓷基复合材料应力-应变曲线.此外,还分析了Weibull模量、残余热应力和初始开裂应力对应力-应变曲线的影响.研究表明:Weibull模量越大,应力-应变曲线非线性越明显;热残余应力越大,应力-应变曲线偏转越早,非线性越明显;初始开裂应力与Weibull模量对应力-应变曲线影响规律相似.      

混合变量系统基于MCMC的自适应重要抽样法

 系统关键故障的发生,会导致系统处于各种离散性能降级状态.针对传统的基于马尔可夫链蒙特卡罗(Markov chain Monte Carlo,MCMC)的自适应重要抽样法只适用于连续变量系统的不足,提出考虑混合变量的基于MCMC的自适应重要抽样法,以支持系统性能可靠性的高效仿真.该方法首先将由关键故障产生的不同失效域组成失效空间,并通过初始样本点在失效空间中随机游走构造马尔可夫链模拟样本;然后综合考虑连续变量与离散变量,利用核密度估计构建混合核抽样密度函数;再根据该密度函数进行重要抽样仿真并计算系统的性能可靠性;最后对该方法的仿真效率进行理论分析.通过电液舵机(Electro-Hydrostatic Actuator,EHA)案例对方法的正确性和仿真效率进行验证.     

基于加速退化模型的卫星组件寿命与可靠性评估方法

针对卫星产品长寿命、高可靠性的特点,提出了使用加速退化试验的方法来评估卫星用微波接收机的寿命与可靠性指标.首先,在确定微波接收机关键模块与相应失效机理的基础上,选择了适用的加速模型.然后,通过分析产品性能的实际退化过程,选择漂移布朗运动对此进行拟合,并结合漂移布朗运动首达时(first passage time)服从逆高斯分布的特点,建立了产品的可靠性模型.最后,采用极大似然结合最小二乘的评估方法,对该关键模块的恒定应力加速退化试验的蒙特卡罗仿真结果,进行了产品的寿命与可靠性评估并且得到了合理的结果,从而验证了本文方法的正确性.     

地-月低能耗转移轨道中途修正问题研究

采用地-月低能耗转移轨道的探测器从地球停泊轨道转移到极月轨道一般需要3~4个月时间,这类转移轨道对入轨精度有较高的要求。本文对地月转移轨道中途修正问题进行了研究。文中结合地-月低能耗转移轨道的特点,给出一种分段式多目标多次中途修正方案。利用显式制导结合牛顿迭代,分别以地球和月球作为中心天体求解兰伯特问题,在假设探测器各种轨道误差的基础上进行了蒙特卡罗仿真。采用该方法一般需要3~5次中途修正能够满足月球探测器环月轨道入轨精度要求,整个转移过程燃料消耗小于传统地月转移轨道。文中给出的仿真结果验证了该方案的可行性。     

航空发电机健康特征参数与老化模式分析

航空三级式无刷发电机中旋转整流器整流二极管的老化是发电机外特性退化的重要来源.本文选择二极管正向导通压降的变化作为其主要老化模式,对二极管老化和发电机外部特性参数间的关联性进行研究.采用Saber软件搭建航空发电机老化模型,利用蒙特卡罗仿真功能,获取大量发电机外部特性参数数据,从中筛选出与旋转整流器二极管老化状态存在较...     

考虑不确定性的飞机总体参数优化方法

由于在各种设计问题包括飞机概念设计中都存在一定的不确定性,因此在总体参数优化时有必要考虑这种不确定性。以大型客机总体参数优化设计为例,定义了考虑不确定性的飞机总体参数优化问题,该问题与传统飞机总体参数优化的区别是要进行不确定性分析。而不确定性分析的计算量过大,为此提出了一种渐进代理模型方法来解决这一难题。在建立代理模型时,通过连续成批地在设计空间的全局和局部均加入新样本点,不断提高代理模型的全局拟合精度,直至获得满意的代理模型为止。然后在优化过程中使用计算量小的代理模型。大型客机总体参数优化问题中含有5个设计变量,目标函数为起飞重量最轻,并需满足4个性能约束。考虑了不确定性后,不仅使目标值（起飞重量）对总体参数变化的敏感度有所减小,而且满足约束（设计要求）概率显著提高。     

不确定性颤振风险定量分析

 针对小风险颤振定量分析,提出一种概率鲁棒颤振初步分析方法。该方法结合鲁棒颤振分析和概率颤振分析,针对不同层次风险和不确定性水平提供颤振边界特性全面信息,给决策制定提供更好参考。在鲁棒颤振分析中,采用振型迭代的方法,通过结构奇异值( μ )分析和特征值求解得到零风险的鲁棒颤振临界速度。在概率颤振分析中,采用标准蒙特卡罗模拟(MCS)方法得到较大风险时的颤振速度分布和概率临界速度。采用概率比较和二分法估计概率鲁棒颤振不确定度半径和概率鲁棒颤振临界速度,同时采用样本重新利用的随机方法估计颤振风险和不确定性定量关系。通过一个考虑集中质量不确定性的大展弦比双梁式机翼的颤振计算分析表明,确定型μ方法只能用于零风险颤振分析,标准的MCS方法只能用于较大风险颤振速度分布问题。概率鲁棒颤振分析结果表明增加1%颤振风险可以使不确定度半径增大62%或颤振速度边界增大5%,从而放宽了设计要求。     

高温高速稀薄流的DSMC算法与流场传热分析

在非结构网格下,研究了热力学碰撞传能的6种方式和化学反应的8种类型,详细给出了上述碰撞传能和化学反应碰撞类型的子程序计算框图.在可变硬球(VHS)分子模型、Borgnakke-Larsen唯象模型、Bird的化学反应几率模型以及壁面CLL(Cercignani-Lampis-Lord)反射模型的基础上,用Fortran语言编制了能够模拟内能松弛、热力学非平衡和化学非平衡的稀薄气体DSMC(direct si mulation Monte-Carlo)源程序,在地球大气层和火星大气层中完成了Ballute减速装置的8个工况计算(其中Knudsen数从0.05变到30.0,飞行Mach数从26.3变到11.2),并与NASA Langley研究中心2007年发表的计算结果作了比较,本文的结果令人满意.另外,这里用|T-T_v|/T去分析热力学非平衡,用Damk hler数去分析化学反应非平衡,用Stan-ton数去分析飞行器壁面的传热效果,所有这些分析对高温部件的热防护设计十分有益.      

DSMC方法的边界采样与统计技术

研究了适合在三维直接模拟蒙特卡罗(DSMC,Direct Simulation Monte Carlo)方法计算中根据入口边界条件合理地布置进入计算区域内粒子参数的边界采样技术.提出了一种根据穿越入口边界、出口边界以及撞击壁面边界的粒子统计数据获取边界位置处宏观流场参数的新型边界统计技术,弥补了DSMC计算中传统统计技术仅能够获取流场内部参数,不能够获取界面处流场参数的不足.这种界面统计技术具有普遍的适用性,不仅限于DSMC方法.以无穷远来流圆球绕流为例进行计算,验证了DSMC方法的边界采样与统计技术的有效性.