表决部件的备件量计算方法、仿真方法、终端及存储介质技术

本发明专利技术公开了一种表决部件的备件量计算方法、仿真方法、终端及存储介质，属于表决部件备件，要解决的技术问题为以换件维修的方式对故障单元进行修理时，如何确定机电类表决部件的备件需求量。其方法包括基于机电单元寿命服从威布尔分布，构建备件数量与保障概率之间的关系式，并以机电类表决部件在保障任务期间计划的工作时间、保障概率不低于保障概率阈值为条件，通过上述关系式得到备件数量的预测值。仿真方法用于验证上述计算方法的有效性。终端中处理器被配置用于调用程序指令执行上述计算方法和/或仿真方法。存储介质中程序指令当被处理器执行时处理器执行上述计算方法和/或上述仿真方法。

Calculation method, simulation method, terminal and storage medium of spare parts quantity of voting parts

【技术实现步骤摘要】
表决部件的备件量计算方法、仿真方法、终端及存储介质
本专利技术涉及表决部件备件领域，具体地说是一种表决部件的备件量计算方法、仿真方法、终端及存储介质。
技术介绍
表决部件k/n(G)是指该部件由n个同型单元组成。当部件开始工作后，这n个单元一起开始工作，随着部件中故障单元的逐渐增多，部件的输出性能会下降(此时部件相当于降额使用)，但只要部件中完好单元的数量不小于k时，则该部件仍可视为在正常工作。一旦部件中的完好单元数量下降到小于k时，则认为此时该部件已不能满足工作要求，视为部件出现故障而停止工作。采用表决结构的部件，使得该部件在维修性上可获得的好处之一是能在部件不停机的情况下完成对故障单元的修理，即：当部件中的部分单元出现故障时，只要单元故障数量不大于n-k，就能在部件中其他单元继续工作的同时，对故障单元开展维修工作，从而为修理工作赢得时间。在本文中，以换件维修的方式对故障单元进行修理。机电类表决部件是指组成部件的单元是机电类元器件。威布尔分布常用于描述机电件的寿命分布情况，例如滚珠轴承、继电器、蓄电池、液压泵、齿轮、材料疲劳件等。威布尔型单元指寿命服从威布尔分布的单元，单元寿命X服从威布尔分布W(A,B)，X的密度函数为：某机电类表决部件记为k/n(G)，由n个同型机电单元组成，单元寿命服从威布尔分布W(A,B)。该表决部件在保障任务期间计划的工作时间记为Tw。在保障任务期间，采用即时修理策略，即：一旦出现单元故障立刻开展换件维修工作；如果备件消耗完毕，则部件将继续工作到完好单元数量为k-1时刻或保障任务的终点时刻。若任务期间该部件能始终工作则视为保障任务成功。备件数量记为S，要求在这些备件支持下保障成功的概率不得低于PS。基于上述，以换件维修的方式对故障单元进行修理时，如何确定机电类表决部件的备件需求量，是需要解决的技术问题。
技术实现思路
本专利技术的技术任务是针对以上不足，提供一种表决部件的备件量计算方法、仿真方法、终端及存储介质，来解决如何确定机电类表决部件的备件需求量的问题。第一方面，本专利技术提供一种表决部件的备件量计算方法，用于计算机电类表决部件的备件需求量，所述机电类表决部件记为k/n(G)，由n个同型机电单元组成，机电单元寿命X服从威布尔分布W(A,B)，其中，尺度参数A＞0，形状参数B≥1，所述表决部件正常工作时，至少有k个机电单元正常工作；所述计算方法包括：基于机电单元寿命X服从威布尔分布W(A,B)，构建备件数量S与保障概率Pt之间的关系式，并以机电类表决部件在保障任务期间计划的工作时间Tw、保障概率不低于保障概率阈值Ps为条件，通过上述关系式得到备件数量S的预测值。作为优选，所述计算方法包括如下步骤：S100、初始化备件数量S，S＝0；S200、通过如下步骤计算Ps1：S210、设定数组p1，数组p1中共有S+1个数，数组p1中第j个元素的表达式为：其中，变量x表示机电单元的寿命，1≤j≤S+1，S220、初始化i＝1；S230、如果i＝1，设定数组p2＝p1，否则，计算p1与p2的卷积，并将卷积结果赋值给p2，即：p2＝p1*P2；S240、设定i＝i+1，如果i≤n，执行步骤S230，否则，通过如下公式计算Ps1：S300、遍历生成S个备件被n个机电单元消耗的所有情况，遍历结果为Sall，所述Sall为ml×n的矩阵，Sall中行向量形式为[i1,i2,......,in]，in表示第n号机电单元消耗的备件数量；S400、通过如下步骤计算Ps2：S410、初始化Ps2＝0，i＝1；S420、设定n×2的矩阵P2，矩阵P2中第一列元素P2(r,1)的表达式为：其中，sn表示Sall中第i行向量，snr表示sn中第r个元素，1≤r≤n；矩阵P2中第二列元素P2(r,2)的表达式为：p2(r,2)＝1-p2(r,1)S430、初始化r＝1；S440、如果r＝1，设定行向量pt2＝P2(r,:)，否则，计算P2(r.:)与pt2的卷积，并将卷积结果赋值给pt2，即：pt2＝P2(r,:)*Pt2；其中，P2(r,:)为矩阵P2中的第r行向量；S450、设定r＝r+1，如果r≤n，执行步骤S440，否则，通过如下公式计算Ps2，并执行步骤S460：S460、设定i＝i+1，如果i≤ml，执行步骤S420，否则，执行步骤S500；S500、设定保障概率Pt＝Ps1+Ps2，如果Pt＜Ps，令S＝S+1后，执行步骤S200，否则，执行步骤S600；S600、终止计算，S为满足保障概率要求且最为经济的备件需求量。作为优选，步骤S300基于n-1重for循环实现遍历生成S个备件被n个机电单元消耗的所有情况，其中，第1重循环中的变量i1取值范围为0～S之间的整数；第2重循环中的变量i2取值范围为0～(S-i1)之间的整数；第j重循环中的变量ij取值范围为之间的整数；其中，1≤j≤n-1，变量第二方面，本专利技术提供一种仿真方法，所述仿真方法用于模拟机电类表决部件的备件保障过程，以验证如第一方面任一项所述的表决部件的备件量计算方法的有效性；所述仿真方法包括：按照如下仿真步骤多次进行仿真，统计保障结果标志F并计算其均值，将其均值作为保障概率的数学期望值；仿真步骤为：L100、设定当前机电类表决部件的备件数量N1＝S，按照机电单元的寿命分布规律，产生n个随机数simTi，1≤i≤n；L200、选取上述n个随机数中的最小值，将所述作为最小值的随机数的序号记为g；L300、如果simTg＞Tw，仿真结束，并令保障结果标志F＝1，否则，执行步骤L400；L400、如果N1＞0，令N1＝N1-1，按照机电单元的寿命分布规律产生1个随机数simTt，令simTg＝simTg+simTt，并执行步骤L200，否则，执行步骤L500L500、通过如下方法更新保障结果标志F：对n个随机数simTi，进行升序排序，并选取第n-k+1个数值，如果simTn-k+1≥Tw，则保障成功，令F＝1，否则令F＝0；其中，1≤i≤n。第三方面，本专利技术提供一种终端，包括处理器、输入设备、输出设备和存储器，处理器、输入设备、输出设备和存储器相互连接，存储器用于存储计算机程序，所述计算机程序包括程序指令，所述处理器被配置用于调用所述程序指令执行如第一方面任一项所述的表决部件的备件量计算方法和/或如第二方面任一项所述的仿真方法。第四方面，本专利技术提供一种存储介质，所述存储介质为计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序包括程序指令，程序指令当被处理器执行时所述处理器执行如第一方面任一项所述的表决部件的备件量计算方法本文档来自技高网...

【技术保护点】
1.表决部件的备件量计算方法，其特征在于用于计算机电类表决部件的备件需求量，所述机电类表决部件记为k/n(G)，由n个同型机电单元组成，机电单元寿命X服从威布尔分布W(A,B)，其中，尺度参数A＞0，形状参数B≥1，所述表决部件正常工作时，至少有k个机电单元正常工作；/n所述计算方法包括：基于机电单元寿命X服从威布尔分布W(A,B)，构建备件数量S与保障概率Pt之间的关系式，并以机电类表决部件在保障任务期间计划的工作时间T

【技术特征摘要】
1.表决部件的备件量计算方法，其特征在于用于计算机电类表决部件的备件需求量，所述机电类表决部件记为k/n(G)，由n个同型机电单元组成，机电单元寿命X服从威布尔分布W(A,B)，其中，尺度参数A＞0，形状参数B≥1，所述表决部件正常工作时，至少有k个机电单元正常工作；
所述计算方法包括：基于机电单元寿命X服从威布尔分布W(A,B)，构建备件数量S与保障概率Pt之间的关系式，并以机电类表决部件在保障任务期间计划的工作时间Tw、保障概率不低于保障概率阈值Ps为条件，通过上述关系式得到备件数量S的预测值。


2.根据权利要求1所述的表决部件的备件量计算方法，其特征在于所述计算方法包括如下步骤：
S100、初始化备件数量S，S＝0；
S200、通过如下步骤计算Ps1：
S210、设定数组p1，数组p1中共有S+1个数，数组p1中第j个元素的表达式为：



其中，变量x表示机电单元的寿命，1≤j≤S+1，






S220、初始化i＝1；
S230、如果i＝1，设定数组p2＝p1，否则，计算p1与p2的卷积，并将卷积结果赋值给p2，即：
p2＝p1*P2；
S240、设定i＝i+1，如果i≤n，执行步骤S230，否则，通过如下公式计算Ps1：



S300、遍历生成S个备件被n个机电单元消耗的所有情况，遍历结果为Sall，所述Sall为ml×n的矩阵，Sall中行向量形式为[i1,i2,......,in]，in表示第n号机电单元消耗的备件数量；
S400、通过如下步骤计算Ps2：
S410、初始化Ps2＝0，i＝1；
S420、设定n×2的矩阵P2，矩阵P2中第一列元素P2(r,1)的表达式为：



其中，sn表示Sall中第i行向量，snr表示sn中第r个元素，1≤r≤n；
矩阵P2中第二列元素P2(r,2)的表达式为：
p2(r,2)＝1-p2(r,1)
S430、初始化r＝1；
S440、如果r＝1，设定行向量pt2＝P2(r,:)，否则，计算P2(r.:)与pt2的卷积，并将卷积结果赋值给pt2，即：
pt2＝P2(r,:)*Pt2；
其中，P2(r,:)为矩阵P2中的第r行向量；
S450、设定r＝r+1，如果r≤n，执行步骤S440，否则，通过如下公式计算Ps2，并执行步骤S460：



S460、设定i＝i+1，如果i≤ml...

【专利技术属性】
技术研发人员：孔庆福，赵中元，彭英武，杨明明，刘君，李庆民，
申请(专利权)人：青岛航讯网络技术服务有限公司，
类型：发明
国别省市：山东;37