1、弹用涡扇发动机飞行可靠性评估论文关键词:弹用涡扇发动机 可靠性评估 寿命分布 AMSAA 模型 环境因子论文摘要:对发动机进行可靠性评估是检验发动机是否满足可靠性要求的重要手段.根据某型弹用涡扇发动机在研制过程中的各种地面试验数据和少量的飞行试验数据,采用一种基于 AMSAA 模型的方法建立了环境因子模型,给出了环境因子的求解方法,并利用求得的环境因子获得了等效的飞行数据.在此基础上分别建立了发动机在地面环境和飞行环境下的寿命分布模型,并实现了可靠性评估.弹用涡扇发动机作为导弹的重要组成部分,其可靠性直接影响到导弹的性能、寿命和费用,因此在发动机定型前必须对其可靠性进行定量的评估,以检验是否满
2、足提出的可靠性要求1.确定发动机寿命分布及可靠性评估的依据是发动机在各种试验中得到的时间失效数据.弹用涡扇发动机在研制过程中所进行的试验主要包括各种地面条件下的性能实验(磨合试验、加减速试验、二次起动试验、节流试验)、可靠性摸底试验、可靠性增长试验和飞行试验2.如何利用地面试验数据和少量的飞行试验数据来对发动机的飞行可靠性做出合理评估,成为发动机飞行可靠性评估中一个急需解决的问题3.本文首先给出了基于 AMSAA 模型的环境因子的确定方法,并利用优化方法求得了不同地面试验环境相对于飞行试验环境的环境因子,在此基础上实现了弹用涡扇发动机飞行环境下的寿命分布建模和可靠性评估.1 环境因子的确定在可
3、靠性数据的统计分析中可以使用环境因子将各种试验环境下的产品的试验时间转化为使用环境下的试验时间.11 基于 AM SAA 模型的环境因子设某产品在研制阶段共经历了 m 个(m2)试验项目;产品在研制期间总故障数为 N;产品第 1 个试验项目起止节点为 T0, T1,第 j 个试验项目的起止节点为 Tj-1, Tj(j=1,2,m);产品第 q 个故障的累计试验时间为 tq(q =1,2,N);各试验项目的环境因子为 k1,k2,km.则产品折合后在整个研制试验过程中的第 q 个故障的累计故障时间 tn,q 可由下式计算:式中:Tj-1tqTj;Tn,0 为产品折合后整个研制试验的起始节点;Tn
4、,j-1,Tn,j 分别为产品折合后整个研制试验中的第 j 个试验项目的起止节点.此时环境因子的求解过程可表示为一个多约束极小形式的最优化过程:min(C2N). (2)约束条件为式中:U1-/2 和 C2N, 为产品在显著性水平为 下的增长趋势检验统计量的临界值及 AMSAA 模型拟合优度检验统计量的临界值,可查文献4的增长趋势分析临界值表和克莱默-冯梅赛斯检验临界值表.U 为产品增长趋势检验统计量;C2N 为产品 AMSAA 模型拟合优度检验统计量,可由下式求得:式(4)适用于定时截尾情形,其中 Tn,m 为产品折合后的整个研制试验的截止时间,-b 为产品的增长参数.式(5)适用于定数截尾
5、情形,其中 tn,N 为产品折合后的整个研制试验中最后一个失效时间,式中的产品增长参数-b 为12 环境因子的求解方法1)环境应力类型的划分拟将所进行的试验项目分为:性能试验(包括磨合试验、加减速试验、验证试验、湿启动、干启动试验)、可靠性试验和试飞 3 种环境应力类型.2)环境因子取值范围的选取本文取性能试验的环境因子 k1 的取值范围为 01,可靠性试验的环境因子 k2 的取值范围为 13;试飞的环境因子 k3 的取值为 1.3)环境因子的求解算法式(2)和式(3)可写成如下的一般形式:采用外部罚函数法求解式(8)确定的多约束极小值问题.对可行域作一惩罚项 p(x),有:式中: 为大于等于
6、 1 的常数,取 =2.构造一个增广目标函数: 式中: 为罚因子,0.于是式(8)所示的带约束的极小值问题转化为无约束的极小值问题: 2 寿命分布模型寿命分布建模是发动机可靠性评估的前提和基础.根据某型弹用涡扇发动机的故障性质,选择威布尔分布来描述发动机的寿命分布.因两参数威布尔分布是三参数威布尔分布的特例,本文仅给出三参数威布尔分布参数的点估计方法.三参数威布尔分布的累积分布函数为: 式中:、和 分别为形状参数、尺度参数和位置参数,且 、 和 均大于零.采用基于图估法的极大似然估计法进行威布尔参数的点估计.概率密度函数和失效率函数为: 设数据列为(t1, t2, tn),既含有失效的数据,也
7、含有截尾的数据.令 为待估的模型参数列(,),则对数似然函数为: iCF 和C 分别表示失效数据和截尾数据的子集合.将上式分别对 、 和 求偏导数,并令结果为 0,可得如下非线性方程组: 解方程组(15)可得 , 的估计值.解方程组时的 , 初值由图估法给出.3 飞行可靠性评估结果收集到的原始数据包括 15 个地面试验数据和 5 个飞行数据,如表1 所示.其中地面试验数据包括 8 个故障数据, 7 个截尾数据,飞行数据均为截尾数据.数据后面带+表示截尾数据.由表 1 可知,发动机的累积故障时间为:275, 363, 827, 41998, 45558, 46318,54808,94648 mi
8、n.性能试验起始节点时间为 0,终止节点时间为 41998 min;可靠性试验起始节点时间为41998 min,终止节点时间为 94648 min.利用上述数据求得环境因子k1=095,k2=103.将地面试验中性能试验数据乘以环境因子 k1、地面试验的可靠性试验数据乘以环境因子 k2、试飞的试验数据乘以试飞环境因子 k(k=1).折合后的数据见表 2.带+的为截尾数据,带*的为实际试飞数据.对表 1 和表 2分别进行拟合优度检验,确定了表 1 中的数据符合三参数威布尔分布,表2 中的数据符合两参数威布尔分布.对两种寿命分布分别进行参数估计,结果如表 3 和表 4 所示.发动机地面环境下的可靠
9、性参数和飞行环境下的可靠性参数如表5 和表 6 所示.其中,可靠度是指在时间为 90min 时的可靠度,可靠寿命是指可靠度为 08 时的可靠寿命.发动机地面试验环境下的可靠度函数曲线和飞行环境下的可靠度函数曲线如图 1 所示.4 结 论1)基于 AMSAA 模型所确定的环境因子,是符合工程实际的.2)利用地面试验数据和经环境因子折合后的数据得到的发动机的寿命分布模型是不同的.3)通过引入环境因子,可将地面试验数据折合成飞行数据,从而增大了发动机的飞行样本容量,提高了发动机飞行可靠性评估的可信度.参考文献:1韩庆田,杨兴根,张杰.弹用涡喷发动机可靠性参数的选择和评定方法J.飞航导弹,2004(3):43-46.2肖波平.弹用涡喷发动机性能监视与诊断系统软件研制J.推进技术, 2001, 22(2): 111-113.3傅博,杜振华,赵宇,等.一种基于遗传算法的环境因子的确定方法J.北京航空航天大学学报, 2004, 30(5): 466-468.4陈昭宪.AMSAA 模型的增长分析.电子产品可靠性与环境试验, 19954(1): 66-72.
