1、引用格式:罗赓，张畔，唐平建，等 某型加速度计可靠性寿命试验抽样方法J 电光与控制，2023,30(3):101-106 LUO G,ZHANG P,TANG P J,etal Sampling of a certain accelerometer in reliability life testing J Electronics Optics Control,2023,30(3):101-106某型加速度计可靠性寿命试验抽样方法罗赓1，张畔1，唐平建1，刘恂2(1 中国华阴兵器试验中心，陕西 华阴714000;2 陆军工程大学，南京210000)摘要:为了解决某型加速度计进行可靠性寿命试验最
2、小试验样本量难以确定的难题，基于可靠性统计推断理论和机电产品寿命分布特点，提出一种结合原始数据修正、分布模型建立和抽样方法设计的适合于工程应用的最小样本量抽样方法。通过该方法得到了该型加速度计进行可靠性寿命试验的样本量为 11，并进行了试验验证。结果表明，该抽样方法计算得到的先验参数值与原值之间的相对误差不超过 1 5%，可以满足工程应用。关键词:可靠性;寿命试验;抽样方法;最小样本量;加速度计中图分类号:TP114 3文献标志码:Adoi:10 3969/j issn 1671 637X 2023 03 018Sampling of a Certain Accelerometer in Re
3、liabilityLife TestingLUO Geng1,ZHANG Pan1,TANG Pingjian1,LIU Xun2(1 Ordnance Test Center，Huayin 714000,China;2 Army Engineering University of PLA,Nanjing 210000,China)Abstract:In order to solve the problem that it is difficult to determine the minimum sample size for thereliability life test of an a
4、ccelerometer,based on the reliability statistical inference theory and the lifedistribution characteristics of mechanical and electrical products,a minimum sample size sampling methodsuitable for engineering application is proposed,which combines the original data modification,distributionmodel esta
5、blishment with sampling method design It is obtained through this method that the sample size ofthe accelerometer for reliability life test is 11,which is experimentally verified The results show that therelative error between the prior parameter value calculated by the sampling method and the origi
6、nal value isless than 1 5%,which can meet the engineering applicationKey words:reliability;life testing;sampling method;minimum sample size;accelerometer0引言某型加速度计构造复杂，既为某型弹种控制平台距离修正系统中重要组成部分1，亦为其自然环境贮存条件下考核的薄弱失效环节，因此，需进行可靠性寿命试验2 4 以评估其可靠性寿命。试验中样本量直接决定着评估精度，样本量越多，可靠性评估精度越高。但是对于加速度计这类单体价值高的高新被试品，需采用最少
7、样本量以满足可靠性寿命试验5 的精度与效率。如何确定试验中最小试验样本量是当前可靠性寿命试验领域的研究难点。文献 6 9 基于可靠性相关统计理论，研究了产品服从特定分布时的抽样方法;文献 10 基于机械材料力学性能参数，研究了产品贮存收稿日期:2022-02-22修回日期:2022-03-26基金项目:国家自然科学基金(61471385)作者简介:罗赓(1990)，男，陕西咸阳人，硕士，工程师。环境试验的抽样方法。上述方法均不能有效解决该型加速度计在分布未知情况下进行可靠性寿命试验最小试验样本量。1原始数据修正在贮存时间 ti后的该型加速度计总体中，抽取 ni个样本检测，有 Yi个未失效。针对
8、该次检测，记 D=(ti，ni，Yi)1ik，k1，k 为时点个数。则其失效分布函数为F(t)=P(T t)=1 R(t)(1)式中:T 为该型加速度计可靠性寿命;R(t)为可靠度函数。则pi=Yi/ni。(2)此时偶尔会出现被试品可靠度为关于时间的增函数这一与自然规律相悖的“倒挂”现象，即 pi pj(i j)。因此，采用 PAVA 方法11 修正“倒挂”的原始数据第 30 卷第 3 期2023 年 3 月电光与控制Electronics Optics ControlVol 30No 3Mar 2023以满足相应约束条件。设 ti时刻检测样本的失效频率为fi=1 Yini。(3)若 fifi
9、+1，则f*i=fi。(4)若 fi fi+1，则B=i，i+1(5)B=i+i+1(6)f*i=fB=AV(B)=iBfiiiBi(7)式中，=(1，k)，i0，为给定的权函数。如此重复上述方法，直到 fi fi+1。2分布模型建立设 H0:FP0，P0=F(;)，其中，P0为分布族，为未知参数，为参数空间。相关研究表明机电产品常见的寿命分布族有以下4 种10。1)指数分布族，F(t)=1 expt()t0。(8)2)威布尔分布族，F(t)=1 expt()mt0。(9)3)极值分布族，F(t)=1 exp expt ()t0。(10)4)对数正态分布族，F(t)=t012xe(ln x)2
10、22dxt 0。(11)针对不同分布族之间不好比较的问题，采用基于极小卡方估计和拟合优度检验相结合的贮存可靠性评估方法建立模型。令2()=ki=1ni(pi pi()2pi()(1 pi()。(12)对式(12)极小化求解，即2(n)=inf2()(13)可得极小 2估计 n以及极小 2统计量 2(n)。根据可靠性统计理论易知，同自由度下 2统计量越小，分布族拟合程度越好，然而不同分布族之间自由度也不尽相同。基于此，引入检验的拟合优度 p 来衡量假定分布族与数据的拟合程度。令p=P(2(k s)2(n)(14)式中，s 为参数个数。求解式(14)可得不同分布族对应 p 值，其值越大，分布族拟合
11、程度越好。3抽样方法设计取来自修正后总体 Xt1(贮存时间 t1后)及 Xt2(贮存时间 t2后，t1t2)的两组截尾样本。1)设a1 a2 ar1at是来自Xt1的按照时间序列进行排序的一组样本，定时截尾时间为 ai，失效的样本数量为 r1，样本容量为 n1，则 ai=ainr1atr1 nn1。2)设 q1 q2 qr1qt是来自 Xt2一组按照时间序列进行排序的的样本，定时截尾时间为 qt，失效的样本数量为 r2，样本容量为 n2，则 qi=qinr2qtr2 nn2。31对数正态分布311模型建立设贮存时间 t1后其寿命随机变量为 X1，服从对数正态分布，即 X1 LN(，2)。利用
12、t 统计量，构造枢轴量t1=n1(a)s1。(15)对 t1进行温索化，得到 t1 t(r11)，即t1=r11n11t1(16)式中:a=1n1n1i=1ln ai;s21=1n1 1n1i=1(ln ai a)2。设贮存时间 t2后其寿命随机变量为 X2，服从对数正态分布，然而由于已经贮存一段时间，在环境载荷的冲击作用下其剩余寿命也随之减少。因此需对剩余寿命进行修正，根据产品可靠性寿命试验失效机理不变原理，引入一个时间调节系数 h(0 h1)，即 X2LN(h，2)。利用 t 统计量，构造枢轴量t2=n2(q h)s2(17)对 t2进行温索化，得到 t2 t(r21)，即t2=r21n2
13、1t2(18)其中:q=1n2n2i=1ln qi;s22=1n2 1n2i=1(ln qi q)2=n2n2 12。当 h 有先验分布时，有P t2t1(r21)|h=1 (19)即201第 30 卷电光与控制罗赓等:某型加速度计可靠性寿命试验抽样方法Pr21n21n2(q h)s2t1(r21)h=1 (20)变换形式后，得P 1hqs2(n21)(r21)n2t1 (r21)h=1 。(21)记 L为 h 给定时，置信度为 1 的平均置信下限，则L=qs2(n21)(r21)n2t1 (r21)E(h1)。(22)因为 t2 t(r21)，故其数学期望为E(r21)n2(q h)(n21
14、)s()2=0(23)即=qE(h1)。(24)精度可靠性评估采用 n=EeL ee()表示，则n=(eEE(h1)E(n21)(r21)t1 (r21)1)(eE(h1)(r21)(r21)t1 (r21)1)，(25)r2即为满足一定精度要求的所需最小失效样本数量。设某个精度 0，精度要求 n 满足条件的最小样本容量 r2可由式(26)求出，(eE(h1)(r21)(r21)t1 (r21)1)。(26)3 1 2h 的估计构造 t 统计量t=w y(h)sw1/r1+1/r2 t(r1+r22)(27)式中，S2w=(r11)S21+(r21)S22r1+r22。对于给定的置信度 1 ，
15、有Pw y(h)sw1/r1+1/r2 t/2(r1+r22()=1 。(28)由于 a是 的无偏估计，可以认为 =E()=w。故 h 在置信度 1 的置信区间为ysw1/r1+1/r2t/2(r1+r22)w，y+sw1/r1+1/r2t/2(r1+r22)w()。(29)3 1 3E(h1)的取值设 g(H)是在连续场合下随机变量 H 的函数，若其数学期望存在，则E g(H)=bag(h)pH(h)dh(30)式中:0 a b1;pH(h)为随机变量 H 的密度函数。故E(h1)=bah1pH(h)dh。(31)3 2指数分布3 2 1模型建立设贮存时间 t1后其寿命随机变量为 X1服从指
16、数分布，即 X1 exp(1)。设贮存时间 t2后其寿命随机变量为 X2服从指数分布，同理引入一个时间调节系数 h(0 h1)，即 X2exp(1h)。定时截尾试验场合下，易知:2T1 2(2r1+1)，2T2h2(2r2+1)，且相互独立，其中，T1=r1i=1ai+(n1 r1)at，T2=r2i=1qi+(n2 r2)qt。当 h 有先验分布时，有P2T2h21 (2r2+1)|h=1 。(32)变换形式后，得P 2T221 (2r2+1)h|h=1 。(33)记 L为 h 给定时，置信度为 1 的平均置信下限，则L=2T221 (2r2+1)E(h1)。(34)因为2T2h 2(2r2+1)，故其数学期望为E2T2h()h=2 r2+12()+1 r2+()12(35)即E(2T2h)=2h r2+12()+1 r2+()12。(36)结合重期望公式 E(Z)=EE(Z K)，则E(L)=221 (2r2+1)r2+12()+1 r2+()12E(h1)E(h)。(37)精度可靠性评估采用n=EL()(38)301第 3 期表示，则n=221 (2r2+1)r2+12()+1
