基于贝叶斯因子和二阶概率方法的圣地亚热传导模型确认挑战问题

针对美国圣地亚国家实验室为促进模型确认方法的发展所提出的热传导挑战问题,根据现代模型确认思想总结了该挑战问题解答中应包含的核心内容。在确认度量环节采用贝叶斯因子考察了实验数据是否支持所给模型,在此基础上通过二阶概率方法得到了模型预测的分布,以此计算出模型预测结论的置信度。该过程中考虑了模型参数的随机不确定性和认知不确定性,最后通过灵敏度分析辨识了模型参数的不确定性对模型预测的影响。研究表明该挑战问题中的实验数据支持所给模型,模型预测受导热系数的不确定性影响最大,模型预测材料在调控条件下的失效概率不满足调控要求,该结论的置信度为99.97%。     

基于密度权重的可靠性灵敏度分析方法

为了提高可靠性灵敏度求解数字模拟法的效率,提出了一种变量空间确定性低偏差均匀抽样与样本点处联合概率密度函数构造权重相结合的方法,来估计可靠性灵敏度。该方法通过均匀样本点处联合概率密度函数的权重保证了可靠性灵敏度的估计值收敛于真值,而由低偏差抽样代替原问题中的联合概率密度抽样则可以保证更低的误差阶以及在小失效概率条件下抽得的样本有更高的可能性落入失效域,从而保证了所提方法具有更高的收敛速度。另外,所提方法可以采用与独立变量相同的步骤来估计相关变量情况下的可靠性灵敏度,计算简便,适用范围广。算例充分证明了所提方法的优越性。     

基于径向基函数改进的Delaunay图映射动网格方法

为了提高大变形下动网格生成的效率和质量,基于径向基函数插值技术对原始Delaunay图映射动网格方法进行了改进。首先通过带多项式的径向基函数插值方法计算出背景网格远场边界点的位移;然后将背景网格位移插值到计算网格;最后利用衰减函数将计算网格远场位移衰减为零。基于矩形旋转网格变形实例比较了改进方法与原始Delaunay图映射动网格方法之间的差别,并研究了径向基基点数目和衰减函数类型对变形网格质量的影响。矩形旋转网格变形实例说明该方法能够精确恢复出Delaunay背景网格的旋转特性。NACA 0012翼型、NLR 7301两段翼和M6机翼网格变形算例进一步证明,通过添加控制点,该方法能够不重构背景网格实现大变形下高质量动网格的生成。     

区域杂波估计的多目标跟踪方法

高斯粒子概率假设密度（PHD）滤波往往假定杂波密度参数已知,这种做法对于实际应用是不现实的。此外,杂波的参数值通常依赖于环境条件,可能随时间发生变化。因此,多目标跟踪算法中需要实时准确估计杂波密度的参数。基于此,提出了一种多目标跟踪的区域杂波估计方法。首先根据量测信息在线估计出场景中的杂波数目,然后估计落入目标附近感兴趣区域的杂波数,并估计每个目标感兴趣区域杂波强度。仿真结果表明,在复杂场景下算法的跟踪性能明显优于未进行杂波估计的多目标跟踪算法,提高了跟踪的实时性和跟踪精度。     

一种改进的航空发动机结构概率风险评估方法

针对航空发动机适航条款FAR33.75中关于发动机限寿件(ELLP)结构失效概率要求,提出了一种基于Kriging和蒙特卡罗半径外重要抽样(MCROIS)混合的结构概率风险评估方法。该方法针对ELLP高维、小失效概率事件以及极限状态函数为隐式、高度非线性的特点,利用Kriging元模型模拟隐式极限状态函数,然后通过主动学习迭代算法,计算最优点(MPP,最接近设计验算点的样本点),更新实验设计(DOE)并提高Kriging元模型的模拟精度。在此基础上,利用Kriging元模型确定最优抽样半径,构造半径外重要抽样密度函数,在最优抽样半径确定区域进行抽样,通过构造主动学习函数,使样本点更多落在抽样半径确定的球区域附近,加速失效概率计算的收敛,并构建了ELLP风险概率模型,解决了高维、小失效概率事件以及隐式、非线性极限状态函数的发动机结构概率风险评估难题,以某型发动机低压压气机轮盘为应用示例,与传统的蒙特卡罗仿真(MCS)方法进行了对比,验证了该方法的高效率、鲁棒性和仿真精度。     

硫化压力对橡胶／金属热硫化粘接剥离强度的影响

研究 NR1155天然橡胶在不同硫化压力下热硫化后的物理力学性能。基于90°粘合剥离强度测试方法分别测量了不同硫化压力下制备的橡胶／金属复合结构试样的粘合剥离强度,并采用低场 NMR 橡胶交联密度测定仪分析试样中橡胶部分的交联密度,同时采用冷场发射电子扫描显微镜（SEM）研究剥离破坏表面。结果表明：NR1155天然橡胶胶料与 Chemlok205／Chemlok220热硫化胶粘体系具有很好的相容性和协同硫化作用。随着硫化压力的提高,橡胶／金属的粘合剥离强度呈现先降低后增加的规律,且最大剥离力呈现对硫化压力的函数依赖关系,而硫化胶的交联密度则呈现相反的趋势。     

对偶型最大熵概率密度函数模型及其优化解法

针对经典型最大熵概率密度函数模型及其计算目前存在的非线性程度高,优化不收敛,求解效率低等问题,提出了一种对偶型最大熵概率密度函数模型+逐次优化的方法.根据优化过程不稳定,重新推导了拉格朗日系数的线性变换公式.针对几种常见及一种复杂的概率密度函数,采用经典型与对偶型最大熵概率密度函数模型分别计算概率密度及可靠度的对比表明:与经典型最大熵概率密度函数模型相比,对偶型最大熵概率密度函数模型优化函数形式简单,非线性程度低.逐次优化法求解拉格朗日系数不仅克服了初始值敏感性问题,而且计算效率高.对偶型最大熵概率密度函数模型+逐次优化法与其他方法相比,计算精度最高,且能很好的应用于复杂概率分布及可靠性问题.      

PIV技术在复杂二相流场中的应用

光学元器件随飞行器在大气中飞行时,其工作性能越来越多地受到大气悬浮汇聚微粒的影响。大气微粒在复杂流场中呈现何种运动汇聚效应,对于合理准确评估机载光学元器件的工作效能具有十分重要的工程意义,而复杂气动流场中微粒分布状态的预估一直是飞行器外界环境研究中的一个难点。气动问题的复杂性、大气中微粒的多样性一直是制约各种试验手段展开、数值模型建立的主要因素。利用先进的激光粒子图像技术,在风洞中对舵面旋涡主导的复杂流场中的微粒速度及分布特性进行了实验研究。在测量舵面翼梢脱落旋涡特性的基础上,通过激光片光扫描流场全域,同时高帧频CCD相机同步曝光,利用PIV 拍摄到的流场中涡流截面内微粒分布的瞬态图像。结合图像后处理技术,对原始粒子图像进行互相关、二值化处理,通过对图像区域内的灰度值计算,统计相对流场截面内的粒子浓度系数,得到在复杂旋涡结构流场内瞬态粒子的分布特性规律。研究结果表明,利用大气中微粒在激光片光下的米氏散射原理,可以有效地拍摄到复杂流场结构下粒子光学散射及分布的特性图像,解决了传统环境测试设备无法对复杂条件下流场内粒子分布进行实时测量的缺陷;在旋涡为主导的流场中,大气中的微粒由向心力牵引,在涡核周围达到平衡运动状态,微粒环绕涡核形成一条环状带,这一区域中的粒子浓度系数要远大于自由流场中的微粒,涡核中心粒子呈“空洞”状态。     

详细化学反应机理对预混丙烷钝体熄火模拟影响

航空发动机熄火预测是重要关键问题之一，湍流和化学反应的非线性相互作用使预测非常困难。本文采用大涡模拟（LES）对湍流进行高精度模拟，采用概率密度函数输运方程湍流燃烧模型（TPDF）耦合JL4、Z66和H73三种化学反应机理，对预混丙烷钝体熄火现象和规律进行研究。JL4的反应机理最简单，反应释热快，局部放热高，火焰宽度大，火焰两侧温度梯度大，燃烧更加趋于稳定，无法模拟出熄火状态。H73机理绝热火焰温度低，火焰温度低，回流区中部OH含量高；在近熄火状态，大量CO被氧化，释放热量过高导致无法模拟出熄火现象。Z66机理可以模拟出火焰正常状态，在低当量比下也可以模拟出熄火状态。本文算例中，局部Da数大于1的区域超过35%则会发生熄火。     

C/C-SiC复合材料致密度影响因素

利用先驱体转化法制备C/C-SiC复合材料,对试样进行微观结构分析和性能测试,研究渗硅温度、保温时间、真空度和裂解周期对C/C-SiC复合材料致密度的影响。结果表明:随着渗硅温度的升高,材料的致密度呈先加速升高后快速下降趋势;随着保温时间的延长,材料的致密度先快速升高,保持一段时间稳定后再缓慢降低;随着烧结真空度的提高,材料的致密度加速升高;随着裂解周期的增加,材料的致密度不断增大,但增速逐步降低。经过11周期的“浸渍-固化-裂解”过程后,所制备的C/C-SiC复合材料获得最大密度2.09 g/cm3、最小孔隙率7.6%,其综合力学性能最为优异:弯曲强度468 MPa、拉伸强度242 MPa、断裂韧度19.6 MPa?m1/2、维氏硬度17.2 GPa。