CN115310048B - 一种设备在期望时间内完成修复概率的计算方法和系统 - Google Patents

一种设备在期望时间内完成修复概率的计算方法和系统 Download PDF

Info

Publication number
CN115310048B
CN115310048B CN202211223796.8A CN202211223796A CN115310048B CN 115310048 B CN115310048 B CN 115310048B CN 202211223796 A CN202211223796 A CN 202211223796A CN 115310048 B CN115310048 B CN 115310048B
Authority
CN
China
Prior art keywords
time
repair
probability
component
calculating
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202211223796.8A
Other languages
English (en)
Other versions
CN115310048A (zh
Inventor
翟亚利
阮旻智
李松涛
李华
冯杭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Naval University of Engineering PLA
Original Assignee
Naval University of Engineering PLA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Naval University of Engineering PLA filed Critical Naval University of Engineering PLA
Priority to CN202211223796.8A priority Critical patent/CN115310048B/zh
Publication of CN115310048A publication Critical patent/CN115310048A/zh
Application granted granted Critical
Publication of CN115310048B publication Critical patent/CN115310048B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/20—Administration of product repair or maintenance

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Operations Research (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Quality & Reliability (AREA)
  • Probability & Statistics with Applications (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Evolutionary Biology (AREA)
  • Strategic Management (AREA)
  • Marketing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Algebra (AREA)
  • Economics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种设备在期望时间内完成修复概率的计算方法和系统，属于设备故障指标量化领域。包括：在任务时间内，结合各部件的累计工作时间，对其寿命服从的对数正态分布密度函数积分计算，得到任务时间内各部件发生故障的概率；按照检查次序，根据任务时间内各部件发生故障的概率，计算任务时间内各部件的修理权重系数；按照检查次序，根据各部件的状态检查消耗时间服从的正态分布和修理各失效部件消耗时间服从的正态分布，分别计算修复时间分布参数数组；根据任务时间内各部件的修理权重系数和修复时间分布参数数组，计算设备在期望时间内完成修复概率。本发明实现在期望时间内完成修复的概率的预测，能更具体、详尽地描述装备维修性性能。

Description

一种设备在期望时间内完成修复概率的计算方法和系统
技术领域
本发明属于设备故障指标量化领域，更具体地，涉及一种设备在期望时间内完成修复概率的计算方法和系统。
背景技术
当装备出现某种故障现象后，首先对多个可能引起该故障现象的零部件逐一进行检查，直至找到失效的零部件，然后对该失效件采取更换备件等修理方式完成修复。当故障现象和故障原因是一对多的关系时，因失效件的不确定性，会导致每次完成修复的时间也并不相同。目前常采用平均修复时间（mean time to repair，MTTR）来描述装备维修性。
对于海军舰船装备而言，舰员级修理是一种在海上执行任务期间、装备发生故障后在装备现场进行的修理，也是一种在修理设施、修理工具、修理人员数量和水平等方面都极为有限的修理。舰员级MTTR指标对战时恢复装备作战能力至关重要，受到装备生产方和军方的高度重视。生产方采取各种措施来满足军方的MTTR指标，例如采用自动测试技术帮助舰员快速找到故障原因，广泛采用模块化技术来设计装备，使得舰员能快速拆除失效件、更换备件从而修复装备。当前在使用MTTR时，存在着两个较大问题：一是在落实MTTR指标时，装备设计/生产方和军方大多采用“针对双方约定好的某具体故障进行MTTR指标考核”的方式。这种方式的背后原因，在于不能在更一般、更广泛的情况下估计MTTR，只好退而求其次，通过“实现”部分或有代表性故障的平均修复时间来“体现”装备的整体MTTR性能。二是MTTR的数学本质是均值，是一种在宏观、总体层面进行描述的指标，但实际上即便同一种故障现象，由于失效件不同、故障排查时间的不确定性，因此修复时间实际上是分布在某个范围内的。在实际工作中，即便知道类似“修复该故障的平均耗时46分钟”这样的结论，也仍然更希望能得到“在0.5小时内完成修复的概率有多大”这类问题的答案。
发明内容
针对现有技术的缺陷，本发明的目的在于提供一种设备在期望时间内完成修复概率的计算方法和系统，旨在解决现有技术无法预测在期望时间内完成修复的概率的问题。
为实现上述目的，第一方面，本发明提供了一种设备在期望时间内完成修复概率的计算方法，所述设备包括多个部件，所述部件的寿命均服从于对数正态分布，整个任务时间内任意时刻最多一个部件发生故障，故障排查时各部件的状态检查的次序独立不相关，该方法包括：
S1.获取各部件的寿命服从的对数正态分布密度函数、状态检查消耗时间服从的正态分布和累计工作时间，获取修理各失效部件消耗时间服从的正态分布和故障发生后对所有部件的检查次序，并将设备的一段工作时期作为任务时间；
S2.在任务时间内，结合各部件的累计工作时间，对其寿命服从的对数正态分布密度函数积分计算，得到任务时间内各部件发生故障的概率；
S3.按照检查次序，根据任务时间内各部件发生故障的概率，计算任务时间内各部件的修理权重系数；
S4.按照检查次序，根据各部件的状态检查消耗时间服从的正态分布和修理各失效部件消耗时间服从的正态分布，分别计算修复时间分布参数数组；
S5.根据任务时间内各部件的修理权重系数和修复时间分布参数数组，计算设备在期望时间内完成修复概率。
优选地，步骤S2包括：
S21.设置部件编号
Figure 865473DEST_PATH_IMAGE001
；
S22.计算任务时间
Figure 164868DEST_PATH_IMAGE002
内部件
Figure 142051DEST_PATH_IMAGE003
发生故障的概率
Figure 459769DEST_PATH_IMAGE004
：
Figure 864205DEST_PATH_IMAGE005
当
Figure 650896DEST_PATH_IMAGE006
时，
Figure 900611DEST_PATH_IMAGE007
；
当
Figure 151464DEST_PATH_IMAGE008
时，
Figure 664485DEST_PATH_IMAGE009
；
其中，
Figure 789DEST_PATH_IMAGE010
表示部件的数量，
Figure 54195DEST_PATH_IMAGE011
表示部件
Figure 97238DEST_PATH_IMAGE012
的条件概率，
Figure 843477DEST_PATH_IMAGE013
、
Figure 582588DEST_PATH_IMAGE014
分别表示部件
Figure 705265DEST_PATH_IMAGE015
的寿命服从的对数正态分布密度函数中的对数均值和对数根方差，
Figure 602814DEST_PATH_IMAGE016
表示部件
Figure 254375DEST_PATH_IMAGE017
的累计工作时间；
S23.
Figure 565271DEST_PATH_IMAGE018
，若
Figure 163742DEST_PATH_IMAGE019
，进入步骤S22，否则，进入步骤S3。
优选地，步骤S3包括：
S31.设置部件检查序号
Figure 978115DEST_PATH_IMAGE020
；
S32.计算任务时间内检查序号
Figure 269419DEST_PATH_IMAGE021
对应的部件的修理权重系数：
Figure 802031DEST_PATH_IMAGE022
并按照以下方式赋值四个中间变量：
Figure 931DEST_PATH_IMAGE023
，
Figure 856761DEST_PATH_IMAGE024
，
Figure 381283DEST_PATH_IMAGE025
，
Figure 338875DEST_PATH_IMAGE026
；
其中，
Figure 341466DEST_PATH_IMAGE027
表示部件的数量，
Figure 68113DEST_PATH_IMAGE028
，
Figure 763537DEST_PATH_IMAGE029
表示编号为
Figure 942845DEST_PATH_IMAGE030
的部件任务时间内发生故障的概率，
Figure 749127DEST_PATH_IMAGE031
表示故障发生后对所有部件的检查次序，
Figure 845128DEST_PATH_IMAGE032
分别表示编号为
Figure 445874DEST_PATH_IMAGE033
的部件的状态检查消耗时间服从的正态分布中的均值和根方差，
Figure 174795DEST_PATH_IMAGE034
分别表示修理编号为
Figure 456872DEST_PATH_IMAGE033
的失效部件的消耗时间服从的正态分布的均值和根方差；
S33.
Figure 220429DEST_PATH_IMAGE035
，若
Figure 929759DEST_PATH_IMAGE036
，进入步骤S32，否则，进入步骤S4。
优选地，步骤S4包括：
S41.设置部件检查序号
Figure 145977DEST_PATH_IMAGE037
；
S42.分别计算修复时间分布参数数组
Figure 294061DEST_PATH_IMAGE038
和
Figure 849807DEST_PATH_IMAGE039
；
S43.
Figure 57935DEST_PATH_IMAGE040
，若
Figure 948399DEST_PATH_IMAGE041
，进入步骤S42，否则，进入步骤S5。
优选地，步骤S5包括：
S51.设置部件检查序号
Figure 634596DEST_PATH_IMAGE042
，初始化设备在期望时间
Figure 107165DEST_PATH_IMAGE043
内完成修复概率
Figure 158298DEST_PATH_IMAGE044
；
S52.更新设备在期望时间
Figure 349108DEST_PATH_IMAGE045
内完成修复概率：
Figure 776678DEST_PATH_IMAGE046
其中，
Figure 369333DEST_PATH_IMAGE047
表示任务时间内检查序号
Figure 856947DEST_PATH_IMAGE048
对应的部件的修理权重系数；
S53.
Figure 535053DEST_PATH_IMAGE049
，若
Figure 750002DEST_PATH_IMAGE050
，进入步骤S52，否则，终止计算，输出概率
Figure 197164DEST_PATH_IMAGE051
。
优选地，该方法还包括：
S6.选择多个不同的期望时间，重复步骤S5，计算得到修复时间的概率分布。
为实现上述目的，第二方面，本发明提供了一种设备在期望时间内完成修复概率的计算系统，包括：包括处理器和存储器；所述处理器用于存储计算机执行指令；所述处理器用于执行所述计算机执行指令，使得第一方面所述的方法被执行。
总体而言，通过本发明所构思的以上技术方案与现有技术相比，具有以下有益效果：
本发明提供了一种设备在期望时间内完成修复概率的计算方法和系统，按照检测次序，确定任务时间内各部件的修理权重系数，根据各部件的状态检查消耗时间服从的正态分布和修理各失效部件消耗时间服从的正态分布，分别计算修复时间分布参数数组，综合修理权重系数和修复时间分布参数数组，计算设备在期望时间内完成修复概率，实现在期望时间内完成修复的概率的预测，能更具体、详尽地描述装备维修性性能。
附图说明
图1为本发明实施例提供的一种设备在期望时间内完成修复概率的计算方法流程图。
图2为本发明实施例提供的分别采用仿真法和本发明方法得到的修复时间在26~280分钟范围内完成修复的概率分布结果。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
本发明涉及的设备包括多个部件，所述部件的寿命均服从于对数正态分布，整个任务时间内任意时刻最多一个部件发生故障，故障排查时各部件的状态检查的次序独立不相关。图1为本发明实施例提供的一种设备在期望时间内完成修复概率的计算方法流程图。如图1所示，该方法包括：
步骤S1.获取各部件的寿命服从的对数正态分布密度函数、状态检查消耗时间服从的正态分布和累计工作时间，获取修理各失效部件消耗时间服从的正态分布和故障发生后对所有部件的检查次序，并将设备的一段工作时期作为任务时间。
在工程上，对数正态分布用于半导体器件的可靠性分析和某些种类的机械零件的疲劳寿命。对数正态型单元指该单元寿命服从对数正态分布
Figure 855678DEST_PATH_IMAGE052
，其密度函数
Figure 21081DEST_PATH_IMAGE053
，其中，
Figure 790453DEST_PATH_IMAGE054
为对数均值，
Figure 826543DEST_PATH_IMAGE055
为对数根方差。
本发明约定：
（1）某装备由多个对数正态型单元组成，为便于描述，以时间来描述各单元的寿命。
（2）在任意时刻，至多有1个单元发生故障。当某单元发生故障时会影响装备的正常工作，装备会出现某些故障现象，此时需要进行开展修理工作。
（3）在进行故障确认时，对这些单元进行状态检查的次序是独立不相关的，即：不存在“必须先检查单元A、然后再检查单元B”这类对检查次序有特定要求的情况。
（4）已知各单元的寿命分布规律、对每个单元进行（正常与否的）状态检查所消耗时间的分布规律、修理各失效单元消耗时间的分布规律、各单元的累积已工作时间、即将执行任务的时间和某故障现象发生后对所有相关单元的检查次序。
本发明的相关变量约定如下：单元数量记为
Figure 718275DEST_PATH_IMAGE056
；对相关单元的检查次序构成了故障确认方案，按照该次序依次检查单元状态直至找到失效单元为止；为便于论述，按照检查次序对单元进行编号，即：第
Figure 308657DEST_PATH_IMAGE003
个进行检查的单元其编号记为
Figure 678458DEST_PATH_IMAGE003
；单元
Figure 287163DEST_PATH_IMAGE003
的寿命服从对数正态分布
Figure 349797DEST_PATH_IMAGE057
；单元
Figure 161895DEST_PATH_IMAGE003
的累积已工作时间记为
Figure 335387DEST_PATH_IMAGE058
；对单元
Figure 611648DEST_PATH_IMAGE003
的状态检查时间服从正态分布
Figure 782866DEST_PATH_IMAGE059
，其密度函数
Figure 144577DEST_PATH_IMAGE060
；对失效单元
Figure 590602DEST_PATH_IMAGE003
的修复时间服从正态分布
Figure 924632DEST_PATH_IMAGE061
，任务时间记为
Figure 329068DEST_PATH_IMAGE002
。
步骤S2.在任务时间内，结合各部件的累计工作时间，对其寿命服从的对数正态分布密度函数积分计算，得到任务时间内各部件发生故障的概率。
优选地，步骤S2包括：
S21.设置部件编号
Figure 365026DEST_PATH_IMAGE001
。
S22.计算任务时间
Figure 880321DEST_PATH_IMAGE002
内部件
Figure 865595DEST_PATH_IMAGE003
发生故障的概率
Figure 378615DEST_PATH_IMAGE004
：
Figure 714919DEST_PATH_IMAGE005
当
Figure 971588DEST_PATH_IMAGE006
时，
Figure 76947DEST_PATH_IMAGE009
；
当
Figure 760869DEST_PATH_IMAGE008
时，
Figure 318890DEST_PATH_IMAGE062
；
其中，
Figure 175987DEST_PATH_IMAGE010
表示部件的数量，
Figure 322804DEST_PATH_IMAGE011
表示部件
Figure 239944DEST_PATH_IMAGE012
的条件概率，
Figure 488523DEST_PATH_IMAGE013
、
Figure 883732DEST_PATH_IMAGE014
分别表示部件
Figure 698104DEST_PATH_IMAGE015
的寿命服从的对数正态分布密度函数中的对数均值和对数根方差，
Figure 989408DEST_PATH_IMAGE016
表示部件
Figure 522021DEST_PATH_IMAGE017
的累计工作时间。
S23.
Figure 455342DEST_PATH_IMAGE018
，若
Figure 327483DEST_PATH_IMAGE019
，进入步骤S22，否则，进入步骤S3。
步骤S3.按照检查次序，根据任务时间内各部件发生故障的概率，计算任务时间内各部件的修理权重系数。
优选地，步骤S3包括：
S31.设置部件检查序号
Figure 852005DEST_PATH_IMAGE020
。
S32.计算任务时间内检查序号
Figure 799145DEST_PATH_IMAGE021
对应的部件的修理权重系数：
Figure 67315DEST_PATH_IMAGE022
并按照以下方式赋值四个中间变量：
Figure 793962DEST_PATH_IMAGE063
，
Figure 958228DEST_PATH_IMAGE064
，
Figure 465432DEST_PATH_IMAGE065
，
Figure 209397DEST_PATH_IMAGE066
；
其中，
Figure 852868DEST_PATH_IMAGE027
表示部件的数量，
Figure 453614DEST_PATH_IMAGE028
，
Figure 120219DEST_PATH_IMAGE029
表示编号为
Figure 730191DEST_PATH_IMAGE067
的部件任务时间内发生故障的概率，
Figure 680699DEST_PATH_IMAGE031
表示故障发生后对所有部件的检查次序，
Figure 186767DEST_PATH_IMAGE068
分别表示编号为
Figure 668563DEST_PATH_IMAGE067
的部件的状态检查消耗时间服从的正态分布中的均值和根方差，
Figure 754331DEST_PATH_IMAGE069
分别表示修理编号为
Figure 106815DEST_PATH_IMAGE067
的失效部件的消耗时间服从的正态分布的均值和根方差。
S33.
Figure 49363DEST_PATH_IMAGE035
，若
Figure 690560DEST_PATH_IMAGE036
，进入步骤S32，否则，进入步骤S4。
步骤S4.按照检查次序，根据各部件的状态检查消耗时间服从的正态分布和修理各失效部件消耗时间服从的正态分布，分别计算修复时间分布参数数组。
优选地，步骤S4包括：
S41.设置部件检查序号
Figure 642336DEST_PATH_IMAGE037
。
S42.分别计算修复时间分布参数数组
Figure 52588DEST_PATH_IMAGE038
和
Figure 166038DEST_PATH_IMAGE039
。
S43.
Figure 356848DEST_PATH_IMAGE040
，若
Figure 33686DEST_PATH_IMAGE041
，进入步骤S42，否则，进入步骤S5。
步骤S5.根据任务时间内各部件的修理权重系数和修复时间分布参数数组，计算设备在期望时间内完成修复概率。
优选地，步骤S5包括：
S51.设置部件检查序号
Figure 626341DEST_PATH_IMAGE070
，初始化设备在期望时间
Figure 910692DEST_PATH_IMAGE043
内完成修复概率
Figure 260902DEST_PATH_IMAGE044
。
S52.更新设备在期望时间
Figure 554480DEST_PATH_IMAGE045
内完成修复概率：
Figure 939325DEST_PATH_IMAGE046
其中，
Figure 394577DEST_PATH_IMAGE047
表示任务时间内检查序号
Figure 559979DEST_PATH_IMAGE048
对应的部件的修理权重系数。
S53.
Figure 329352DEST_PATH_IMAGE071
，若
Figure 631020DEST_PATH_IMAGE072
，进入步骤S52，否则，终止计算，输出概率
Figure 257174DEST_PATH_IMAGE051
。
优选地，该方法还包括：步骤S6.选择多个不同的期望时间，重复步骤S5，计算得到修复时间的概率分布。
本发明提供了一种设备在期望时间内完成修复概率的计算系统，包括：包括处理器和存储器；所述处理器用于存储计算机执行指令；所述处理器用于执行所述计算机执行指令，使得上述方法被执行。
实施例：已知某部件由10个对数正态型单元组成，各单元的相关信息如表1，即将执行100小时的任务。约定发生故障后，依次按照单元序号1、3、7、8、6、2、4、5、9、10进行状态检查，直至找到失效单元后，对该单元进行修理，完成修复。采用上述方法，计算该故障发生后在180分钟内完成修复的概率。
表1 各单元的相关信息
Figure 831243DEST_PATH_IMAGE073
1）遍历计算各单元发生故障的概率
Figure 466624DEST_PATH_IMAGE074
，单元1至单元10单元发生故障的概率分别为：0.061、0.075、0.047、0.084、0.247、0.065、0.120、0.055、0.100、0.130。
2）按照上述单元检查次序
Figure 826061DEST_PATH_IMAGE075
，遍历计算修理权重系数
Figure 623116DEST_PATH_IMAGE076
为0.062、0.047、0.122、0.056、0.066、0.076、0.085、0.251、0.102、0.132，
Figure 497531DEST_PATH_IMAGE077
、
Figure 608706DEST_PATH_IMAGE078
、
Figure 619388DEST_PATH_IMAGE079
、
Figure 852923DEST_PATH_IMAGE080
见表2。
3）计算修复时间分布参数数组
Figure 152317DEST_PATH_IMAGE081
、
Figure 129501DEST_PATH_IMAGE082
，结果见表2。
4）计算在180分钟时间内完成修复的概率
Figure 712798DEST_PATH_IMAGE083
为0.47。
表2各单元的相关指标
Figure 851655DEST_PATH_IMAGE084
可建立仿真模型验证上述方法的正确性，仿真模型简述如下：
（1）产生
Figure 435083DEST_PATH_IMAGE085
个随机数
Figure 153640DEST_PATH_IMAGE086
，
Figure 404493DEST_PATH_IMAGE087
，
Figure 714252DEST_PATH_IMAGE088
服从单元
Figure 722659DEST_PATH_IMAGE089
的寿命分布规律，且要求所有的
Figure 41645DEST_PATH_IMAGE090
成立，则各单元的剩余寿命
Figure 84687DEST_PATH_IMAGE091
。
（2）在所有
Figure 565347DEST_PATH_IMAGE092
中寻找最小数，对应的序号记为
Figure 388947DEST_PATH_IMAGE093
，即：
Figure 432995DEST_PATH_IMAGE094
，
Figure 658440DEST_PATH_IMAGE095
。
（3）若
Figure 310001DEST_PATH_IMAGE096
成立，则本次仿真有效，按照单元检查次序，检查到单元
Figure 293001DEST_PATH_IMAGE097
，按照相关单元的检查时间分布规律产生对应的随机数，用于模拟其检查时间；按照单元
Figure 953789DEST_PATH_IMAGE098
的修理时间分布规律产生随机数，用于模拟修理该失效件的时间；所有检查时间和修理时间之和，即为本次模拟的修复时间。
在大量多次模拟后，可统计得到修复时间的概率分布。图2为本发明实施例提供的分别采用仿真法和本发明方法得到的修复时间在26~280分钟范围内完成修复的概率分布结果。考虑到仿真的随机性，图2表明二者的结果较为一致。仿真结果表明：该故障的平均修复时间为164.8分钟，修复时间的根方差为60.6分钟。因其修复时间的变化波动较大，以平均修复时间来开展维修管理计划等方面的工作还是过于粗略。
大量仿真验证结果表明：本发明方法能同时考虑装备的可靠性（各单元的寿命分布规律）、装备的健康状态（累积已工作时间）、装备基本组成单元的维修性（各单元的状态检查时间分布规律和修理时间分布规律）和任务时间等因素的影响，准确估计修复时间的概率分布，相比MTTR指标，能更具体、详尽地了解装备的修理性，可用于装备设计阶段的维修性设计方案评估、装备使用阶段的维修方案优化。
本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (5)

1.一种设备在期望时间内完成修复概率的计算方法，其特征在于，所述设备包括多个部件，所述部件的寿命均服从于对数正态分布，整个任务时间内任意时刻最多一个部件发生故障，故障排查时各部件的状态检查的次序独立不相关，该方法包括：
S1.获取各部件的寿命服从的对数正态分布密度函数、状态检查消耗时间服从的正态分布和累计工作时间，获取修理各失效部件消耗时间服从的正态分布和故障发生后对所有部件的检查次序，并将设备的一段工作时期作为任务时间；
S2.在任务时间内，结合各部件的累计工作时间，对其寿命服从的对数正态分布密度函数积分计算，得到任务时间内各部件发生故障的概率；
S3.按照检查次序，根据任务时间内各部件发生故障的概率，计算任务时间内各部件的修理权重系数；
S4.按照检查次序，根据各部件的状态检查消耗时间服从的正态分布和修理各失效部件消耗时间服从的正态分布，分别计算修复时间分布参数数组；
S5.根据任务时间内各部件的修理权重系数和修复时间分布参数数组，计算设备在期望时间内完成修复概率；
步骤S4包括：
S41.设置部件检查序号
Figure 336466DEST_PATH_IMAGE001
；
S42.分别计算修复时间分布参数数组
Figure 402511DEST_PATH_IMAGE002
和
Figure 984802DEST_PATH_IMAGE003
；
S43.
Figure 769087DEST_PATH_IMAGE004
，若
Figure 727816DEST_PATH_IMAGE005
，进入步骤S42，否则，进入步骤S5；
步骤S5包括：
S51.设置部件检查序号
Figure 523733DEST_PATH_IMAGE006
，初始化设备在期望时间
Figure 401559DEST_PATH_IMAGE007
内完成修复概率
Figure 548507DEST_PATH_IMAGE008
；
S52.更新设备在期望时间
Figure 310927DEST_PATH_IMAGE009
内完成修复概率：
Figure 226930DEST_PATH_IMAGE010
其中，
Figure 7149DEST_PATH_IMAGE011
表示任务时间内检查序号
Figure 641392DEST_PATH_IMAGE012
对应的部件的修理权重系数；
S53.
Figure 800978DEST_PATH_IMAGE013
，若
Figure 837067DEST_PATH_IMAGE014
，进入步骤S52，否则，终止计算，输出概率
Figure 932062DEST_PATH_IMAGE015
。
2.如权利要求1所述的方法，其特征在于，步骤S2包括：
S21.设置部件编号
Figure 53602DEST_PATH_IMAGE016
；
S22.计算任务时间
Figure 282458DEST_PATH_IMAGE017
内部件
Figure 907474DEST_PATH_IMAGE018
发生故障的概率
Figure 173371DEST_PATH_IMAGE019
：
Figure 641261DEST_PATH_IMAGE020
当
Figure 549174DEST_PATH_IMAGE021
时，
Figure 28697DEST_PATH_IMAGE022
；
当
Figure 855708DEST_PATH_IMAGE023
时，
Figure 686261DEST_PATH_IMAGE024
；
其中，
Figure 132285DEST_PATH_IMAGE025
表示部件的数量，
Figure 731894DEST_PATH_IMAGE026
表示部件
Figure 464227DEST_PATH_IMAGE027
的条件概率，
Figure 782075DEST_PATH_IMAGE028
、
Figure 766212DEST_PATH_IMAGE029
分别表示部件
Figure 610540DEST_PATH_IMAGE030
的寿命服从的对数正态分布密度函数中的对数均值和对数根方差，
Figure 389140DEST_PATH_IMAGE031
表示部件
Figure 194285DEST_PATH_IMAGE032
的累计工作时间；
S23.
Figure 106746DEST_PATH_IMAGE033
，若
Figure 680947DEST_PATH_IMAGE034
，进入步骤S22，否则，进入步骤S3。
3.如权利要求1所述的方法，其特征在于，步骤S3包括：
S31.设置部件检查序号
Figure 896028DEST_PATH_IMAGE035
；
S32.计算任务时间内检查序号
Figure 784874DEST_PATH_IMAGE036
对应的部件的修理权重系数：
Figure 376392DEST_PATH_IMAGE037
并按照以下方式赋值四个中间变量：
Figure 539520DEST_PATH_IMAGE038
，
Figure 50136DEST_PATH_IMAGE039
，
Figure 829873DEST_PATH_IMAGE040
，
Figure 959503DEST_PATH_IMAGE041
；
其中，
Figure 242717DEST_PATH_IMAGE042
表示部件的数量，
Figure 924234DEST_PATH_IMAGE043
，
Figure 925688DEST_PATH_IMAGE044
表示编号为
Figure 859009DEST_PATH_IMAGE045
的部件任务时间内发生故障的概率，
Figure 855784DEST_PATH_IMAGE046
表示故障发生后对所有部件的检查次序，
Figure 849148DEST_PATH_IMAGE047
分别表示编号为
Figure 337898DEST_PATH_IMAGE049
的部件的状态检查消耗时间服从的正态分布中的均值和根方差，
Figure 933964DEST_PATH_IMAGE050
分别表示修理编号为
Figure 191770DEST_PATH_IMAGE049
的失效部件的消耗时间服从的正态分布的均值和根方差；
S33.
Figure 90456DEST_PATH_IMAGE051
，若
Figure 925557DEST_PATH_IMAGE052
，进入步骤S32，否则，进入步骤S4。
4.如权利要求1所述的方法，其特征在于，该方法还包括：
S6.选择多个不同的期望时间，重复步骤S5，计算得到修复时间的概率分布。
5.一种设备在期望时间内完成修复概率的计算系统，其特征在于，包括：包括处理器和存储器；
所述处理器用于存储计算机执行指令；
所述处理器用于执行所述计算机执行指令，使得权利要求1至4任一项所述的方法被执行。
CN202211223796.8A 2022-10-08 2022-10-08 一种设备在期望时间内完成修复概率的计算方法和系统 Active CN115310048B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211223796.8A CN115310048B (zh) 2022-10-08 2022-10-08 一种设备在期望时间内完成修复概率的计算方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211223796.8A CN115310048B (zh) 2022-10-08 2022-10-08 一种设备在期望时间内完成修复概率的计算方法和系统

Publications (2)

Publication Number Publication Date
CN115310048A CN115310048A (zh) 2022-11-08
CN115310048B true CN115310048B (zh) 2023-01-03

Family

ID=83866846

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211223796.8A Active CN115310048B (zh) 2022-10-08 2022-10-08 一种设备在期望时间内完成修复概率的计算方法和系统

Country Status (1)

Country Link
CN (1) CN115310048B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116579494B (zh) * 2023-05-23 2024-03-19 中国人民解放军海军工程大学 基于维修耗时下机电设备的备件库存预测方法及系统
CN116955910B (zh) * 2023-07-19 2024-04-02 中国人民解放军海军工程大学 一种机械串联部件的保障任务成功率计算方法和系统
CN116843231B (zh) * 2023-07-20 2024-03-19 中国人民解放军海军工程大学 一种考虑维修耗时的机械设备使用可用度量化方法和系统
CN117852844B (zh) * 2024-03-07 2024-05-28 中国人民解放军海军工程大学 设备修理任务分配方法、系统、设备及存储介质

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113326612A (zh) * 2021-05-25 2021-08-31 广西电网有限责任公司电力科学研究院 配电网可靠性评估方法及装置

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR3044787A1 (fr) * 2015-12-03 2017-06-09 Electricite De France Estimation de la fiabilite d'un systeme industriel
CN107016448B (zh) * 2017-04-14 2021-01-05 中国人民解放军海军航空大学 一种故障/维修时间服从任意分布的周期性检查装备可用性评估方法
CN109118108A (zh) * 2018-08-29 2019-01-01 电子科技大学 一种适用于新型复杂机电系统的可靠性预计方法
CN110209515B (zh) * 2019-06-06 2023-05-16 广东电网有限责任公司 一种可靠性评估方法、装置、设备及存储介质
CN112200324A (zh) * 2020-09-11 2021-01-08 中国人民解放军空军工程大学 装备意外损伤平均修复时间测定方法、评定方法及应用

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113326612A (zh) * 2021-05-25 2021-08-31 广西电网有限责任公司电力科学研究院 配电网可靠性评估方法及装置

Also Published As

Publication number Publication date
CN115310048A (zh) 2022-11-08

Similar Documents

Publication Publication Date Title
CN115310048B (zh) 一种设备在期望时间内完成修复概率的计算方法和系统
Liu et al. Dynamic reliability assessment for nonrepairable multistate systems by aggregating multilevel imperfect inspection data
JP4943653B2 (ja) 多事象同時発生時の競合する原因事象の確率及び／又はシステムアベイラビリティを決定するプロセス
RU2757436C9 (ru) Устройство и способ контроля указаний на неисправность от транспортного средства, компьютерочитаемый носитель
US6195624B1 (en) Corrosion analysis system and method
CN115374658B (zh) 一种电子设备最少耗时故障排查次序优化方法和系统
CN115270078B (zh) 一种机电设备平均修复时间的计算方法和系统
Mendoza et al. Optimal life-cycle mitigation of fatigue failure risk for structural systems
Cheng et al. Review and categorization of existing studies on the estimation of probabilistic failure metrics for Reactor Coolant Pressure Boundary piping and steam generator tubes in Nuclear Power Plants
CN115688025A (zh) 一种设备故障修复时间的概率分布的估计方法和系统
Zhao et al. Statistical analysis of time-varying characteristics of testability index based on NHPP
CN116955912A (zh) 一种电子设备的备件保障任务的成功率评估方法及系统
CN115759479A (zh) 一种基于综合值的复杂设备故障定位优化方法和系统
Crosetti Fault tree analysis with probability evaluaticn
Pietruczuk et al. Block inspection policy model with imperfect maintenance for single-unit systems
CN115879720A (zh) 一种采用通用检测工具检测的检查次序优化方法和系统
Coble et al. Identifying suitable degradation parameters for individual-based prognostics
CN110688152B (zh) 一种结合软件开发质量信息的软件可靠性定量评估方法
Bhatti et al. Reliability Analysis of Industrial Model Using Redundancy Technique and Geometric Distribution
Kumar Reliability and probabilistic safety assessment in multi-unit nuclear power plants
CN115293296B (zh) 一种机械设备故障定位优化方法和系统
CN116029692A (zh) 一种针对复杂设备的平均维修时间估计方法和系统
CN116843231B (zh) 一种考虑维修耗时的机械设备使用可用度量化方法和系统
CN115796838A (zh) 一种针对多人合作的故障排查时间估计方法和系统
Hines et al. Application of monitoring and prognostics to small modular reactors

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant