基于故障率的测试性验证试验故障样本分配方案

现有的分层抽样算法中故障率数据不准确是导致测试性验证试验结果不可信的一个原因。为解决这个问题,提出了利用验前信息和Monte Carlo仿真相结合的方法计算故障率。首先考虑故障率的随机特性和量级范围,用Gamma分布来拟合故障率的概率分布,根据上下限分位点确定了Gamma分布中的超参数。然后利用舍取抽样得到了大量的服从Gamma分布的随机数,求平均值作为实际系统的故障率值。以此故障率数据进行分层抽样,得到的故障样本分配方案更合理,试验结果更可信。     

基于模式重要度的航空电源系统可靠性估计

针对航空电源系统可靠性分析的需要,研究现有典型系统的结构特点,提出基于模式重要度的系统可靠性估计方法。该方法首先根据航空电源各子系统的组成,采用蒙特卡罗方法对各子系统可靠性进行分析。根据子系统中各部件的模式重要度,建立子系统可靠性估计模型,进而获取并联系统和总系统的可靠性估计模型。该方法结合了蒙特卡罗方法的优势和系统本身的结构特点,既解决了复杂航空电源系统可靠性分析的难题,又克服了单纯采用蒙特卡罗方法对其进行分析时,仿真效率较低的问题。最后,采用所提出的方法对两种类型的电源系统进行仿真分析,并与单纯蒙特卡罗方法仿真分析结果进行比较,验证了该方法的可行性和有效性。     

多阶段含延缓纠正可靠性增长的MS-NHPP模型及其Bayes分析

针对复杂航空机电系统的多阶段含延缓纠正可靠性增长试验,提出了多阶段非齐次泊松过程(MS-NHPP)模型及其Bayes分析方法。该模型将多个阶段试验看做多个非齐次泊松过程(NHPP),它们具有相同的时间起点和不同的强度函数参数,并满足序化约束关系。对该模型采用Bayes方法进行分析,以Dirichlet分布作为验前分布,并在验后计算中采用基于Metropolis-Hastings准则的马尔可夫链蒙特卡罗(MCMC)方法获得各阶段可靠性参数的验后期望、方差与置信限。实例分析表明,该方法可以有效融合多阶段可靠性增长试验数据以及专家经验和历史信息,尤其适用于小样本情况,参数估计值满足序化约束关系,可靠性增长评估结果更为可信。     

各向异性陶瓷基复合材料涡轮叶片概率性热分析方法

考虑陶瓷基复合材料等纤维增韧复合材料导热系数的各向异性及分散性,建立了基于概率统计的陶瓷基复合材料涡轮叶片热分析方法。研究中以Mark Ⅱ涡轮叶片冷却结构为例,综合利用有限元方法和蒙特卡洛方法,分析了应用陶瓷基复合材料后的温度场均值和波动特性。计算中将导热系数作为随机输入参数,分析了导热系数各向异性及其分散度对叶片前缘滞止点温度、尾缘温度以及高温区域（T>900K）面积的影响。计算中发现在本文的计算工况下,考虑导热系数存在正态波动情况时,叶片前缘滞止点、尾缘温度波动也满足正态分布。前缘滞止点温度在导热系数变异系数为01,导热系数比为2时其温度波动最大,相比12731K的均温,有16%的概率超温913K。尾缘温度在导热系数变异系数为01,导热系数比为10时波动最大,有16%的概率超过均值11529K达527K。计算结果表明:导热系数分散度所带来的波动,会导致叶片内部高温关注区域（T>900K）的面积增大,并且高温关注区域相对增加量ΔShot随导热系数变异系数α的增加而增加。计算结果表明,高温关注区域相对增加量最大发生在导热系数比为2,变异系数为0.1时,此时ΔShot=4.8%。      

两种基于方差的全局灵敏度分析W指标改进算法

在全局灵敏度分析(SA)中,基于方差的灵敏度分析指标(包括Sobol指标和W指标)应用广泛。其中Sobol指标是将输入随机变量固定于特定点时,求得其对输出响应量的平均影响;而W指标求解当输入随机变量在各自分布区间上缩减变化时,输入变量对输出响应量的影响程度。相比Sobol指标,W指标所反映的信息更加全面。但目前对W指标的求解方法还比较欠缺,双层重复抽样蒙特卡罗(DLRS MC)方法和双层一次抽样蒙特卡罗(DLSS MC)方法是两种传统的求解方法。针对W指标的求解问题,提出了两种新算法:改进的蒙特卡罗模拟(AMCS)和基于稀疏网格积分(SGI)的方法。AMCS只需抽取一组样本便可计算出所有变量的各阶W指标,由于该方法是通过筛选策略来计算条件区间上的方差,避免了DLSS MC法中由于小数取整带来的计数误差,从而提高了计算W指标的精度。基于SGI的方法则利用稀疏网格积分来计算三重矩进而得到W指标,由于该方法继承了稀疏网格积分的高效性,因而进一步提高了W指标的计算效率。最后,给出了两个数值算例和一个工程算例,用于验证所提方法求解W指标的准确性和高效性。     

结冰条件下人-机-环系统的飞行风险概率

以结冰条件下的飞行风险量化概率为研究对象,基于蒙特卡罗飞行仿真实验对结冰条件下人-机-环系统的耦合特性进行了分析,并获取了飞行参数极值样本。构建了飞行风险发生的判定条件;对飞行参数极值样本进行了统计特性分析,验证了其厚尾分布特征。一维分布类型辨识结果表明广义极值分布对相对速度和迎角极值的描述精度最高。为描述二维变量对相关性的各自影响程度,提出了一种新的双参数变权重Copula模型;辨识结果表明该Copula模型能以较高的精度通过假设检验。相关性分析的结果表明相对速度和迎角同时出现极大值和极小值的概率较大。基于二维极值样本的Copula分布模型求出了不同结冰程度下的飞行风险概率值,探讨了飞行风险的非线性增长趋势。     

竞争性故障模型可靠性评估的非参数估计方法

针对航空发动机现场数据的可靠性评估面临着小子样与随机删失问题,对可靠性评估中的非参数估计方法进行了对比分析.应用蒙特卡罗方法,在竞争性故障模型下,设计了一种对比可靠性评估方法优劣性的仿真方法.该方法以三重威布尔分布为例,通过模拟产生具有随机删失特点的航空发动机现场故障数据,分别对比4种可靠性评估非参数估计方法的优劣性.最后通过一个工程应用案例验证了模型的有效性.仿真结果表明:在三重威布尔竞争性故障模型下,非参数估计方法的优劣顺序依次是Kaplan-Meier估计、平均秩次法、Herd-Johnson估计、Nelson-Aalen估计.      

二维空间中目标轨迹预测算法研究与分析

本文主要研究移动目标在二维空间内的轨迹预测算法,首先介绍了马尔可夫蒙特卡罗预测算法,包括影响状态阀值的因素,以及转移概率和积累概率的计算方法;接着介绍了线性神经网络的预测算法,如神经元模型和神经网络模型的组成,之后重点介绍了线性神经网络的学习规则,以及学习率对预测算法的影响.最后通过仿真试验,对两种预测方法的精准度和自适应时间进行了对比和分析     

战术导弹多杀伤概率指标检验方法

针对战术导弹多杀伤概率指标大样本检验的实际问题,提出小样本条件下多杀伤概率指标的联合检验、平均检验和蒙特卡洛随机模拟检验等方法,并在实际试验工程中得到应用。     

叶盘结构盘片耦合振动特性的参数敏感性

在参数化建模、求解的基础上重点研究了叶盘结构盘片耦合振动频率转向特性、由结构参数随机失谐引起的振动局部化特性对结构设计参数变化的敏感性.对盘片耦合振动频率转向特性的研究发现,随着轮盘扇区间耦合刚度或叶片刚度的变化,频率转向间隙曲线的改变具有碗状形式;对振动局部化特性的研究表明,以轮盘扇区间耦合刚度和结构参数失谐强度为自变量的失谐响应放大因子曲面呈类马鞍形;响应放大阈值现象存在与否,与轮盘扇区间耦合刚度的大小直接相关.若选择位于鞍部的轮盘扇区间耦合刚度-失谐强度组合参数设计,失谐响应放大因子最多能够被减小37.5%.