删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

基于故障行为的惯导产品贮存寿命试验设计*

本站小编 Free考研考试/2021-12-25

惯导产品(如加速度计、陀螺等)广泛应用于导弹、卫星、火箭等飞行器中,其精度高低和性能优劣对于保障飞行器姿态是否稳定和目标是否准确起着非常重要的作用。由于材料、工艺、制造等因素,惯导产品在受到时间和外界环境的作用影响下,发生内应力释放、蠕变、退磁等变化,导致其关键参数(如加速度计零偏、标度因数)产生漂移,造成在使用过程中普遍存在长期重复性差、使用维护复杂等问题,难以满足整个惯导系统的精度及稳定性要求。惯导产品受使用情况决定,贮存寿命是一个重要指标,其长短会影响武器系统的使用维护性及效能的发挥。因此,如何在短时间内获得惯性产品贮存寿命,为惯导系统的维护标定提供依据是亟待解决的问题。
加速退化试验(ADT)由于不需观测到产品实际失效发生,只需要获取产品在高应力水平下的性能退化数据,通过向正常应力水平下的外推进行寿命评估,使得加速试验技术的应用和高可靠长寿命产品寿命评估得到了长足的发展和广泛应用。加速试验方案优化设计的研究可追溯到20世纪60年代。现有加速退化试验设计方法主要有基于解析的加速退化试验设计方法和基于仿真的加速退化试验设计方法。基于解析的加速退化试验设计方法发展较早,也更为成熟。Yu研究了在一定试验费用约束下,以产品p分位寿命的均方差最小为目标的加速退化试验方案优化设计问题[1]。在此基础上,Yu进一步考虑了产品的退化速率服从逆Weibull分布情况下的优化设计问题[2]。Weaver和Meeker采用大样本近似方法基于随机效应线性退化模型研究了单应力加速退化试验设计和评估问题[3]。Hu等研究了基于Wiener过程的步进加速退化试验方案设计,并证明了多水平步进加速退化试验在一定的准则和约束下会退化为一个恒定应力加速试验[4]。Wang等研究了退化模型为逆高斯过程的试验方案设计问题,得到了在费用约束下的最小p分位寿命渐进方差的优化设计方案,并验证该方法的稳健性[5]。Tsai和Ling等考虑了退化过程中的单调性,采用Gamma过程对产品的退化行为进行建模,同样以费用为约束,以平均故障前时间(MTTF)的最小渐进方差为目标进行优化设计,并进行了敏感性分析[6-7]。而Shi和Meeker则先后用非贝叶斯和贝叶斯方法研究了破坏性加速退化试验方案的优化设计问题[8-9]
基于解析的加速退化试验设计方法虽然能够准确得到试验的最优设计方案,然而对于某些复杂的试验情形则难以得到解析解。随着计算机的计算能力不断提高,基于仿真的加速退化试验设计方法越来越被科学家和工程人员所接受和广泛运用。Li和Kececioglu以LED为对象,考虑应力水平设置因素,给出了基于仿真的加速退化试验设计方法[10]。Nasir和Pan利用仿真方法基于Hellinger距离和DIC准则考虑加速试验设计中的模型鉴别问题[11]。Wang等利用基于仿真的试验设计方法优化了多应力加载条件下的各应力水平的样本数,测量频率及试验终止时间[12]。上述加速退化试验设计方法的研究主要基于统计和优化方法,如果能够利用产品的故障行为模型(FBM),融合物理和其他试验信息,将能快速给出一种试验应力水平设计方法,并且根据文献调研目前还鲜有针对惯导产品的加速退化试验设计方法相关研究。
本文以典型惯导产品加速度计为例,在通过对产品特性和主机理分析基础上,确定表征产品寿命特征参数及试验应力。根据仿真计算、强化试验分析结果以及工程经验,确定各试验应力极限,从而确定试验应力的最高水平和最低水平。根据考虑分散性的故障行为模型,给定置信度下,高应力水平下的置信下限与低应力水平下的置信上限不交叉的原则,在最高应力水平和最低应力水平之间,确定中间应力水平,从而给出了基于故障行为模型的加速度计加速退化试验设计方法。
1 加速度计故障行为模型 石英挠性加速度计是一种机械摆式加速度计,由表头本体和伺服放大电路两部分组成,表头主要由软磁体、磁钢、摆片、力矩器线圈等构成,表头结构见图 1
图 1 石英挠性加速度计表头结构 Fig. 1 Structure of quartz flexible accelerometer header
图选项




故障行为即产品或单元在不同的环境条件和输入条件下,随时间表现出来的,可从外部探知的非正常状态。故障行为模型是描述多种机理作用下产品故障规律的模型,可以用状态变化模型或失效时间模型来表示[13]
在贮存过程中,由于石英挠性加速度计内的石英摆片、环氧树脂胶等材料均耐潮,并且考虑实际贮存剖面为不通电的密封恒温环境(-5~40℃的室温条件),因此,仅考虑贮存过程中的温度因素对加速度计参数的影响,而不再考虑湿度、磁场和电应力等。标度因数是加速度计的重要参数。通过主机理分析可以得知标度因数变化主要由于温度作用下磁钢退磁[14]。加速度计标度因数主要呈现退化状态,由此描述标度因数随时间和应力变化规律的故障行为模型如下[15]
(1)


(2)

式中:ΔBmBm分别为最大磁感应强度变化量和最大磁感应强度;ΔBr为剩余磁感应强度变化量;K1为标度因数;ΔK1t时刻标度因数与t=0时刻标度因数之差;t为时间,h;T为绝对温度,K。
考虑模型参数的分散性,根据式(2),建立标度因数加速退化的随机模型为
(3)

式中:μA=2.31×105σAσB可根据磁钢试验曲线模型参数的分散性确定,这里假设参数AB的变异系数分别为0.05和0.01,即σA=0.05μAσB=0.01μB;根据加速度计标度因数测试误差的随机散布特性(可通过实测数据统计得到),确定σε=5×10-6
2 基于故障行为模型的加速退化试验设计方法 基于故障行为模型的加速退化试验设计方法针对缺乏预试验数据的情况下粗略给出试验应力水平问题,利用故障行为模型可描述参数加速时变规律的特性,根据考虑分散性的故障行为模型,给定置信度下,以高应力水平下的置信下限与低应力水平下的置信上限不交叉为原则,设计给出试验剖面应力水平。
2.1 含置信度的故障行为模型包络 首先确定应力水平极限条件。根据加速度计强化试验分析结果以及工程经验,确定了温度极限,高温环境下加速度计加速退化机理不变的温度上限为90℃。
在确定工作极限条件后,根据考虑分散性的标度因数加速退化模型,即式(3),采用一阶泰勒展开可以近似得到t时刻的方差:
(4)

假设服从正态分布,则t时刻置信度为γ=1-α的单侧置信上限和下限分别为
(5)

(6)

t时刻置信度为(1+γ)/2的的置信区间上下限。uα为标准正态分布的α分位数;σK1t时刻的标准差。
根据考虑分散性的故障行为模型式(5),图 2给出了25、60和90 ℃下10年内的标度因数退化置信曲线(γ=90%),可以看出3个应力水平下的置信曲线没有交叉。
图 2 含置信度的标度因数退化曲线(γ=90%) Fig. 2 Degradation curves of scale factor with confidence coefficient (γ=90%)
图选项




2.2 试验应力水平设计方法 本试验采用基于仿真的加速退化试验方法。在试验样本量、检测时间间隔以及最高应力水平给定的情况下,主要对低应力和中间应力2个水平进行优化。优化时考虑的约束为根据考虑分散性的故障行为模型,给定置信度下,高应力水平下退化量的置信下限与低应力水平下退化量的置信上限不交叉的原则,在最高应力水平和最低应力水平之间,确定中间应力水平。其具体过程如图 3所示。
图 3 T2可选范围确认流程示意图 Fig. 3 Schematic diagram of optional range validation procedure for T2
图选项




步骤1??根据式(3)生成进行N次仿真,按退化量由小到大对N次仿真进行排序,得到,则为其在90%置信度下的双侧置信下限。
步骤2??由图 2可知,其置信下限曲线为一条光滑曲线,将其近似看成某一温度下的退化轨迹,由式(3)反推得到近似等效温度T′1
步骤3??令T2=T′1-1,在T2温度下进行N次仿真,得到退化量在90%置信度下的双侧置信上限。
步骤4??若在任一时刻T′1的置信下限均大于T2的置信上限,则T2此时的值为T2可选温度范围的上限值;若不成立,则令T2=T2-1,重复步骤4,直至满足条件跳出循环为止。
同理可以获得T2可选温度的下限值和T3可选温度的上限值。为了保证一定的加速效果,T3可选温度的下限值已经确定为30℃。
2.3 确定试验最优方案 确定T2T3的可选温度范围并进一步确定备选方案后,需要对备选方案进行优选,其流程如图 4所示,Ti为第i应力水平的温度。
图 4 加速度计标度因数加速退化试验方案优选流程 Fig. 4 Optimization procedure of accelerometer scale factor accelerated degradation test scheme
图选项




步骤1??在进行优化设计前,首先确定备选方案集、优化目标函数、故障行为模型(即加速退化模型)和仿真规模Ωd
步骤2??根据故障行为模型和备选方案dl,仿真得到该方案中的各加速应力下Ωd条标度因数的退化曲线。
步骤3??利用加速度计标度因数的加速退化模型对仿真得到的退化曲线进行拟合,得到备选方案dl下各加速应力下的伪寿命分布及其参数估计。
步骤4??从步骤3中的各加速应力下的伪寿命分布中分别抽取一组伪寿命,构成备选方案dl下一组加速应力下伪寿命仿真数据。
步骤5??利用加速模型对步骤4中的伪寿命仿真数据进行拟合,并外推得到正常温度下可靠寿命估计。
步骤6??重复步骤4和步骤5直至完成Ωd次仿真。
步骤7??对备选方案dl下的Ωd个正常温度下的寿命估计值计算目标函数值。
步骤8??重复步骤2至步骤7直至完成所有仿真计算。
步骤9??对所有备选方案的目标函数值进行排序,确定其中的最优试验方案。
3 案例分析 根据图 3所阐述的方法和流程进行仿真,得到T2可选温度的上限值为75.85℃、下限值为43.85℃,T3可选温度的上限值为61.85℃。其中,T2为可选温度的上限值和T3为可选温度的上限值时,其退化轨迹如图 5所示。
图 5 T2, T3温度可选范围上限时的标度因数退化曲线 Fig. 5 Degradation curves of scale factor under optional range upper limit of temperature T2, T3
图选项




根据上述仿真结果得到的温度选择范围,结合工程实践中所需考虑的温度梯度问题,确定备选方案的应力设置如表 1所示, 1个T2T3构成一组备选方案。
表 1 备选方案T2, T3温度设置 Table 1 Optional plan of temperature setting of T2, T3
T2 T3
75 60, 55, 50, 45, 40, 35, 30
70 55, 50, 45, 40, 35, 30
65 50, 45, 40, 35, 30
60 45, 40, 35, 30
55 40, 35, 30
50 35, 30
45 30


表选项






以备选方案1“T2=75℃, T3=60℃”为例,分别生成3个应力水平下的20 000条退化轨迹,通过非线性拟合得到每条退化轨迹的参数估计值,再分别外推得到每个样本的伪寿命值。在得到的试验温度Ti下的Ωy个伪寿命估计的基础上,进行伪寿命分布的估计,确定伪寿命的分布类型和分布参数。根据图 6显示的概率分布可以看出,对数正态分布可以有效拟合伪寿命数据。
图 6 90℃伪寿命对数正态分布概率图 Fig. 6 Probability graph of lognormal distribution of pseudo life at 90℃
图选项




鉴于伪寿命可用对数正态分布近似描述,那么正常贮存温度下的对数寿命服从正态分布,此时的寿命点估计值同时是样本的中位数,能有效反应群体的平均寿命情况。因此选择正常贮存温度下寿命的对数方差作为优化目标函数比较合适。根据式(3),伪寿命加速模型为
(7)

式中:的失效阈值的绝对值。对多应力水平伪寿命数据(Ti,ln tij),ln tij为该应力水平下的伪寿命数据,采用非线性拟合可以得到参数A′B的估计?′,进一步可获得对数寿命。重复2.3节步骤1、步骤2,Ωd次,则可得到Ωd个对数可靠寿命估计。试验方案dl下的目标函数,即对数方法的计算公式为

式中:为初始寿命估计值;Ul的方差。
优化结果如图 7所示,从图 7中可以看出当T2=75℃, T3=55℃(图中黑点处)时,目标函数达到最小值0.064 2。
图 7 备选方案优化目标函数值 Fig. 7 Optimization objective function value of optional plan
图选项




4 结论 1) 本文提出了一种基于故障行为模型的加速退化试验设计方法。以高应力水平下的置信下限与低应力水平下的置信上限不交叉为原则,采用基于仿真的加速退化试验方法,优化确定了试验方案,从而解决了缺乏预试验数据情况下试验设计问题。
2) 以典型惯导产品加速度计为对象,进行了方法可行性和有效性验证。通过主机理分析得到加速度计主要参数标度因数变化的主退化机理为温度作用下的磁钢退磁,并得到描述标度因数应力作用下时变规律的故障行为模型。在加速度计最高温度应力水平为90℃条件下,基于该模型设计加速度计加速退化试验的优选应力为T2=75℃, T3=55℃。该设计方法为提高基于加速退化试验的产品可靠寿命评估精度提供了前提保障。

参考文献
[1] YU H. Designing an accelerated degradation experiment by optimizing the estimation of the percentile[J].Quality and Reliability Engineering International, 2003, 19(3): 197–214.DOI:10.1002/(ISSN)1099-1638
[2] YU H. Designing an accelerated degradation experiment with a reciprocal Weibull degradation rate[J].Journal of Statistical Planning and Inference, 2006, 136(1): 282–297.DOI:10.1016/j.jspi.2004.06.030
[3] WEAVER B P, MEEKER W Q. Methods for planning repeated measures accelerated degradation tests[J].Applied Stochastic Models in Business and Industry, 2014, 30(6): 658–671.DOI:10.1002/asmb.2061
[4] HU C, LEE M, TANG J. Optimum step-stress accelerated degradation test for Wiener degradation process under constraints[J].European Journal of Operational Research, 2015, 241(2): 412–421.DOI:10.1016/j.ejor.2014.09.003
[5] WANG H, WANG G, DUAN F. Planning of step-stress accelerated degradation test based on the inverse Gaussian process[J].Reliability Engineering & System Safety, 2016, 154: 97–105.
[6] TSAI T, SUNG W, LIO Y L, et al. Optimal two-variable accelerated degradation test plan for Gamma degradation processes[J].IEEE Transactions on Reliability, 2016, 65(1): 459–468.DOI:10.1109/TR.2015.2435774
[7] LING M H, TSUI K L, BALAKRISHNAN N. Accelerated degradation analysis for the quality of a system based on the Gamma process[J].IEEE Transactions on Reliability, 2015, 64(1): 463–472.DOI:10.1109/TR.2014.2337071
[8] SHI Y, MEEKER W Q. Bayesian methods for accelerated destructive degradation test planning[J].IEEE Transactions on Reliability, 2012, 61(1): 245–253.DOI:10.1109/TR.2011.2170115
[9] SHI Y, MEEKER W Q. Accelerated destructive degradation test planning[J].Technometrics, 2009, 51(1): 1–13.DOI:10.1198/TECH.2009.0001
[10] LI Q, KECECIOGLU D B. Design of an optimal plan for an accelerated degradation test:A case study[J].International Journal of Quality & Reliability Management, 2013, 23(4): 426–440.
[11] NASIR E A, PAN R. Simulation-based Bayesian optimal ALT designs for model discrimination[J].Reliability Engineering & System Safety, 2015, 134(2): 1–9.
[12] WANG Y, ZHANG C, ZHANG S, et al. Optimal design of constant stress accelerated degradation test plan with multiple stresses and multiple degradation measures[J].Journal of Risk and Reliability, 2015, 229(1): 83–93.
[13] ZENG Z G, KANG R, WEN M, et al. Measuring reliability during product development considering aleatory and epistemic uncertainty[C]//Proceedings of the Annual Reliability and Maintainability Symposium. Piscataway, NJ: IEEE Press, 2015: 1-6.
[14] XU D, CHEN Y X, KANG R. Quartz flexible accelerometer stability duration prediction on storage condition based on accelerated degradation test[J].Chemical Engineering Transactions, 2013, 33(1): 541–546.
[15] 许丹, 廖寻, 陈云霞, 等. 加速度计贮存稳定期分析[J].北京航空航天大学学报, 2013, 39(2): 173–177.
XU D, LIAO X, CHEN Y X, et al. Analysis of the storage stability of accelerometers[J].Journal of Beijing University of Aeronautics and Astronautics, 2013, 39(2): 173–177.(in Chinese)


相关话题/方案 设计 优化 数据 过程

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 产生匀强磁场的圆柱形线圈组设计方法*
    匀强磁场在精密测量、航空航天等领域都有广泛的应用,在近年来飞速发展的量子传感技术中,高均匀度的匀强磁场更是扮演着重要的角色。在量子传感技术中,原子陀螺仪作为一种新型陀螺仪,有体积小和精度高的优点,已经成为下一代超高精度陀螺仪的发展方向[1-2];原子磁强计是一种利用原子自旋效应来测量磁场的新型磁强计 ...
    本站小编 Free考研考试 2021-12-25
  • 基于Kriging模型的浮空器氦气昼夜温差最优化*
    平流层浮空器具有广泛的应用前景,可以用于对地观测、战区侦查以及远程通信。日本ISAS组织在2002年放飞了总体积达到600000m3的高空气球[1]。2005年,美国洛克希德·马丁公司设计了长152.4m、直径45.7m、总体积105000m3的高空飞艇[2],并且于2011年7月进行了缩比验证艇的 ...
    本站小编 Free考研考试 2021-12-25
  • 高超末段机动突防/精确打击弹道建模与优化*
    高超声速飞行器是一种全球进攻性武器,能进行远距离、防区外作战,最终在末段实现高精度、大威力的定点打击。然而,多层次反导防御系统尤其是末段爱国者-3(PAC-3)低层防御系统的发展,对高超声速飞行器的生存造成很大的威胁。因此,有必要研究末段机动突防和精确打击问题。再入弹头一般利用反导武器拦截时间短、过 ...
    本站小编 Free考研考试 2021-12-25
  • 基于分布式动力的翼身融合飞机整流罩气动设计*
    近年来,翼身融合(Blended-Wing-Body,BWB)布局飞行器由于具有提高气动效率、降低油耗、排放和噪声等方面的潜力[1-3]而逐渐成为研究热点,是一种有希望取代常规布局的未来民机布局形式。为了能进一步发挥其优势,研究者提出在BWB布局上应用一种被称为分布式动力(distributedpr ...
    本站小编 Free考研考试 2021-12-25
  • 随机优化的改进交叉熵方法*
    工程优化设计是以产品结构为基础,根据设计要求、产品功能等,结合现代优化技术,在可行域中搜索最优设计方案,以达到缩短设计周期、提高产品质量和降低生产成本的目的。目前的结构优化设计方法大致可以分为两大类:基于梯度的优化设计方法和随机优化设计方法。前者利用梯度信息从初值点逐渐向最优值点搜索[1];后者则是 ...
    本站小编 Free考研考试 2021-12-25
  • 机载时间触发光网络的设计与实时性分析*
    随着航空电子系统综合化程度的逐渐提高,航空电子网络朝着大吞吐量、强扩展性、强实时性的方向发展。波分复用光(WavelengthDivisionMultiplexing,WDM)网络是目前较好的航空电子网络解决方案之一[1-2],它具有高带宽,强灵活性,对协议和比特率透明,强扩展性等优点,但在强实时性 ...
    本站小编 Free考研考试 2021-12-25
  • 小型飞机自动着舰系统设计准则适用性分析*
    航空母舰(简称航母)搭载有不同种类、功能各异的舰载机以实现多种作战任务。随着无人机技术的发展,航母上应用了越来越多的小尺寸无人机进行战场监控、数据中继等任务。这类飞机通常被划分为Ⅰ类飞机,采用撞网、钩拦阻索等方式回收。这些回收方式都需要飞机在自动着舰系统(AutomaticCarrierLandin ...
    本站小编 Free考研考试 2021-12-25
  • 考虑物料配送的飞机移动生产线调度问题优化*
    飞机移动生产线是丰田生产方式/精益生产理论在飞机制造业的应用,对于降低飞机的装配成本、提升装配效率、缩短装配工期具有重要意义。在飞机整个装配过程中,移动生产线上的飞机以平稳的速度从一个工作站移动到下一个并缓慢通过整个生产线,从而完成飞机所需的所有装配作业,而移动生产线速度的调整,可以设定生产线的不同 ...
    本站小编 Free考研考试 2021-12-25
  • 基于Lamb波气体传感器的数据采集仪*
    超声波作为一种信息载体,已广泛应用于化工工业、航空航天、无损检测、海洋探查及生物医学等领域[1]。Lamb波气体传感器是一种新型的超声波传感器,具有高灵敏度、低损耗、多模式特性,可实现流速矢量、压强、温度、水分等多参数的测量,在气体传感、空气动力学、多参数解耦等方面展示出十分广阔的工业应用前景[2- ...
    本站小编 Free考研考试 2021-12-25
  • 基于飞机客舱空气品质的桥载空调送风优化*
    夏季民航客机在机场过站、短停时,仍需给客舱制冷,为机组人员、中转旅客登机后提供凉爽、舒适的客舱环境。飞机落地后,关闭飞机机载空调,启动悬挂于廊桥下方的桥载空调为客舱制冷,因为飞机机载空调消耗航空煤油,运行成本高且污染环境,而使用桥载空调既经济又不会污染机场空气[1],然而目前的桥载空调控制模式单一, ...
    本站小编 Free考研考试 2021-12-25