[发明专利]一种威布尔型有寿件备件的保障概率计算方法

摘要

本发明公开了一种威布尔型有寿件备件保障概率的计算方法，本计算方法主要包括如下步骤：获得有寿件备件的正态等效参数，在获得正态等效参数后，利用正态等效参数来计算保障概率。按照本发明实现的计算方法，主要采用了对备件的工作寿命进行正态等效，在现有技术中并未对有寿件的寿命和工作寿命进行区分，更加不会意识到寿命分布和工作寿命分布是可以不相同的，也就不容易想到对工作寿命进行正态等效，而本发明想到了采用正态等效的方法计算备件保障概率，保证了计算结果的准确性。

主权项

一种威布尔型有寿件备件保障概率的计算方法，其特征在于，该方法主要包括如下步骤：步骤一：获得正态等效参数，具体包括如下步骤：1.1)在(0,7a)范围内选择包含m个数的等差数组[x1,…,xm]，1≤i≤m，其中其中a为威布尔分布参数中的尺度参数；1.2)遍历计算f1(xi)，其中b为威布尔分布参数中的形状参数；1.3)遍历计算对100×f2(xi)进行四舍五入取整后得到N2(xi)；1.4)在数组中找到最接近更换周期Trep的数xr，即xr≤Trep且xr+1＞Trep；1.5)由N2(xi)得到数组[N3(x1),…N3(xr)]，1.6)记Trep和a之间的较小值为min A,在范围(0,min A)内生成等差数组[μ1,…,μn]；在范围(0,0.4min A)内生成等差数组[σ1,…,σq]；1.7)以如下遍历的方式，生成包含n×q行的矩阵其中1.8)对权重系数组[P11,…,P1n·q]进行初始化，令1.9)对所述权重系数组[P11,…,P1n·q]进行循环更新，具体如下：1.9.1)令i1＝11.9.2)令i2＝11.9.3)计算其中π为圆周率参数；1.9.4)更新i2，令i2＝i2+1；如果i2≤n·q，则执行1.9.3)，否则执行1.9.5)；1.9.5)更新权重系数组[P11,…,P1n·q]，令1.9.6)更新i1，令i1＝i1+1；如果i1≤r，则执行1.9.2)，否则执行1.10)；1.10)根据权重系数组[P11,…,P1n·q]，计算μ、σ：μ=Σv=1n·qP1vy(v,1),σ=Σv=1n·qP1vy(v,2);]]>步骤二：计算备件保障概率利用所述步骤一获得的正态分布参数，按照下式计算备件保障概率P：P=1-1σ2π(1+c)∫-∞Twe-(x-(1+c)μ)22(1+c)σ2dx]]>其中，c为所述备件数量变量。