CN112257880B - 一种k/n负载均担系统的维修决策方法及系统 - Google Patents

一种k/n负载均担系统的维修决策方法及系统 Download PDF

Info

Publication number
CN112257880B
CN112257880B CN202011130675.XA CN202011130675A CN112257880B CN 112257880 B CN112257880 B CN 112257880B CN 202011130675 A CN202011130675 A CN 202011130675A CN 112257880 B CN112257880 B CN 112257880B
Authority
CN
China
Prior art keywords
time
probability
calculating
detection
component
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
CN202011130675.XA
Other languages
English (en)
Other versions
CN112257880A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202011130675.XA priority Critical patent/CN112257880B/zh
Publication of CN112257880A publication Critical patent/CN112257880A/zh
Application granted granted Critical
Publication of CN112257880B publication Critical patent/CN112257880B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/20Administration of product repair or maintenance
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • General Physics & Mathematics (AREA)
  • Marketing (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Operations Research (AREA)
  • General Business, Economics & Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Quality & Reliability (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Debugging And Monitoring (AREA)
  • Testing And Monitoring For Control Systems (AREA)

Abstract

本发明公开了一种k/n负载均担系统的维修决策方法,包括:基于一描述部件相关性的负载平均分配规则,以及一构建的描述系统部件退化过程的多阶段随机过程模型,计算部件故障时间的分布函数;根据离线检测时的系统状态和系统可靠度,定义一个检测周期内所有可能发生的状况,并根据所述部件故障时间的分布函数,计算每一个可能发生的状况的概率;计算一个系统生命周期内的系统总运行成本和系统生命周期长度;计算系统单位时间内的运行成本,得到系统的最优检测时间序列和最优可靠度阈值。通过该发明提出一种新的维修决策方案,保证了系统设备以最低维修成本运行在可靠的状态下。

Description

一种k/n负载均担系统的维修决策方法及系统
技术领域
本发明涉及设备维修技术领域,尤其涉及一种k/n负载均担系统的维修决策方法及系统。
背景技术
应用在各种场景下的设备由于温度、湿度和磨损等原因会发生退化,影响了设备的可靠性,甚至导致设备发生重大事故。预防性维护是指设备出现故障前运用统计分析、维护经济型分析等制定合适的维护计划以降低和避免故障所到来的危害和损失,预防性维修最重要的目的是在设备发生故障之前,及时对设备进行维护或替换,使设备故障所带来的损失将到最低。
在传统的针对k/n系统的维修决策方法中,通常忽略系统中各部件的相关性,只考虑各个部件相互独立的情况。然而在k/n负载均担系统中各个部件共同承担系统总负载,某一部件的故障会导致其他存活部件的加速退化,各个部件之间需要考虑相关性,因此传统的k/n系统维修决策方法并不适用于k/n负载均担系统,对k/n负载均担系统来说需要构建新的维修决策技术方案。另一方面,目前很多系统的维修决策的方法基于在线检测的,但是在实际应用中,大多数系统没有办法在线检测各个部件的退化状况,或者即使能够在线检测也导致系统的维修成本很高。
发明内容
基于此,本发明目的在于提供一种针对所述k/n负载均担系统的维修决策方法及系统,基于离线检测的方法,对所述k/n负载均担系统提出一种新的维修决策方案,保证了系统设备以最低维修成本运行在可靠的状态下。
为实现上述目的,本发明提供了一种k/n负载均担系统的维修决策方法,所述方法包括:
S1、基于一描述部件相关性的负载平均分配规则,以及一构建的描述系统部件退化过程的多阶段随机过程模型,计算部件故障时间的分布函数;
S2、根据离线检测时的系统状态和系统可靠度,定义一个检测周期内所有可能发生的状况,并根据所述部件故障时间的分布函数,计算每一个可能发生的状况的概率,其中,所述系统状态包括正常态、加速退化态和故障态;
S3、根据所述每一个可能发生的状况的概率,计算一个系统生命周期内的系统总运行成本和系统生命周期长度;
S4、根据所述系统总运行成本和系统生命周期长度,计算系统单位时间内的运行成本,得到系统的最优检测时间序列和最优可靠度阈值。
优选的,所述步骤S1包括:
根据所述负载平均分布规则,存活部件退化率μα与部件标准负载ωα的关系用式(1)表示;
lnμα=β01ωn-α (1);
其中,α为故障部件数量,并且
Figure GDA0002783381750000021
β0与β1为待拟合参数。
优选的,所述步骤S1包括:
利用多阶段维纳模型进行数据拟合,所述多阶段维纳模型表示为x=xα(t)(α=0,1,...,n-k),其中,α个部件故障的情况下,存活部件的退化水平与时间的变化关系为xα(t);
Figure GDA0002783381750000022
其中,σ为多阶段维纳过程的噪声参数。
优选的,所述步骤S1还包括:
基于所述多阶段维纳模型,计算第α次(α=1,2,...,n-k+1)部件故障时间tα的分布函数fα(t)和部件在任意时间t的退化水平dt的概率密度函数
Figure GDA0002783381750000034
Figure GDA0002783381750000031
Figure GDA0002783381750000032
其中,H是部件退化水平上限。
优选的,所述步骤S2包括:
采用离线检测的方式,根据i-1离线检测得到系统各部件的退化水平{dj|j=1,2,...,n},通过式(5)计算在i检测时的系统可靠度R,
Figure GDA0002783381750000033
定义一个检测周期内所有可能发生的状况,包括:
第一情况:当第i-1次检测时系统处于正常态以及所述系统可靠度大于预设的系统可靠度最低阈值L,并且第i次检测时系统处于正常态以及所述系统可靠度R大于所述系统可靠度最低阈值L;
第二情况:当第i-1次检测时系统处于正常态以及系统可靠度大于所述系统可靠度最低阈值L,并且第i次检测时系统处于正常态以及系统可靠度小于所述系统可靠度最低阈值L;
第三情况:当第i-1次检测时系统处于正常态以及系统可靠度大于所述系统可靠度最低阈值L,并且第i次检测时系统处于加速退化态;
第四情况:当第i-1次检测时系统处于正常态以及系统可靠度大于所述系统可靠度最低阈值L,并且第i次检测时系统处于故障态。
优选的,所述步骤S2还包括:
通过式(6)计算所述第一情况发生的概率
Figure GDA0002783381750000041
PA i=PA i-1×Pr{N(Ti)=0}×Pr{R(Ti|N(Ti)=0)>L} (6);
其中,
Figure GDA0002783381750000042
为1,Ai-1表示在第i-1次检测和第i次检测之间,第一情况发生,
故障数量N为α的概率为:
Figure GDA0002783381750000043
系统可靠度大于所述系统可靠度最低阈值L的概率为:
Figure GDA0002783381750000044
通过式(9)计算所述第二情况发生的概率PB i
PB i=PA i-1×Pr{N(Ti)=0}×Pr{R(Ti|N(Ti)=0)<L} (9);
通过式(10)计算所述第三情况发生的概率PC i
Figure GDA0002783381750000047
其中Pj i=PA i-1×Pr{N(Ti)=j};
通过式(11)计算所述第四情况发生的概率PD i
Figure GDA0002783381750000053
优选的,所述步骤S3包括:
根据所述概率PB i、概率PC i以及概率PD i,计算一个系统生命周期内系统检测次数的期望EIN
Figure GDA0002783381750000054
根据所述概率PB i和概率PC i,计算系统生命周期以预防性维修结束的概率PPM
Figure GDA0002783381750000055
根据所述概率PD i,计算系统生命周期以事后维修结束的概率PCM
Figure GDA0002783381750000056
根据所述部件故障时间的分布函数,计算系统在各个系统状态下运行时间的期望Tj
Figure GDA0002783381750000057
其中,
Figure GDA0002783381750000061
优选的,所述步骤S3还包括:
根据所述一个系统生命周期内系统检测次数的期望EIN、系统生命周期以预防性维修结束的概率PPM、系统生命周期以事后维修结束的概率PCM,计算得到系统一个生命周期的总运行成本期望EC
Figure GDA0002783381750000062
其中,CIN为一次检测的费用,CCM为一次故障性替换的费用,CPM为一次预防性替换的费用,CADDα为系统在α个部件故障的情况下运行单位时间的费用;
计算得到系统生命周期长度ET
Figure GDA0002783381750000063
优选的,所述步骤S4包括:
根据所述系统的总运行成本期望EC和系统生命周期长度ET,计算得到系统单位时间运行成本:
Figure GDA0002783381750000064
对所述系统单位时间运行成本进行优化处理,得到所述最优检测时间序列T*和最优可靠度阈值L*
(L*,T*)=minCPT (20)。
为实现上述目的,本发明提供了一种k/n负载均担系统的维修决策系统,所述系统包括:
模型构建模块,基于一描述部件相关性的负载平均分配规则,以及一构建的描述系统部件退化过程的多阶段随机过程模型,计算部件故障时间的分布函数;
概率计算模块,根据离线检测时的系统状态和系统可靠度,定义一个检测周期内所有可能发生的状况,并根据所述部件故障时间的分布函数,计算每一个可能发生的状况的概率,其中,所述系统状态包括正常态、加速退化态和故障态;
成本计算模块,根据所述每一个可能发生的状况的概率,计算一个系统生命周期内的系统总运行成本和系统生命周期长度;
模型输出模块,根据所述系统总运行成本和系统生命周期长度,计算系统单位时间内的运行成本,得到系统的最优检测时间序列和最优可靠度阈值。
与现有技术相比,本发明一种k/n负载均担系统的维修决策方法及系统,所带来的有益效果为:针对k/n负载均担系统这种复杂的系统,通过构建模型和离线检测的方式,制定维修策略,并给出系统最佳维修时间和最佳可靠度阈值,以对系统检修时间和可靠度阈值进行优化,从而保证包括检测费用、维修费用等在内的系统运行费用达到最低,并且能够保证系统运行在可靠的状态;本发明针对k/n负载均担系统提出针对性的预防性维护策略,对系统的最佳检修时间进行研究,建立了有针对性的周期性维护策略,使系统维护的成本降至最低。
附图说明
图1是根据本发明的一个实施例的k/n负载均担系统的维修决策方法的流程示意图。
图2是根据本发明的一个实施例的k/n负载均担系统的维修决策系统的系统示意图。
具体实施方式
以下将结合附图所示的具体实施方式对本发明进行详细描述,但这些实施方式并不限制本发明,本领域的普通技术人员根据这些实施方式所做出的结构、方法、或功能上的变换均包含在本发明的保护范围内。
如图1所示的本发明的一个实施例,本发明提供一种k/n负载均担系统的维修决策方法,所述方法包括:
S1、基于一描述部件相关性的负载平均分配规则,以及一构建的描述系统部件退化过程的多阶段随机过程模型,计算部件故障时间的分布函数;
S2、根据离线检测时的系统状态和系统可靠度,定义一个检测周期内所有可能发生的状况,并根据所述部件故障时间的分布函数,计算每一个可能发生的状况的概率,其中,所述系统状态包括正常态、加速退化态和故障态;
S3、根据所述每一个可能发生的状况的概率,计算一个系统生命周期内的系统总运行成本和系统生命周期长度;
S4、根据所述系统总运行成本和系统生命周期长度,计算系统单位时间内的运行成本,得到系统的最优检测时间序列和最优可靠度阈值。
本发明针对k/n负载均担系统,基于离线检测的方式制定维修策略,通过构建单位时间运行成本模型,得到系统的最优检测时间和最优可靠度阈值,从而使系统维修运行费用达到最低。该系统具有以下特征:系统由n个相同的不可维修部件组成,其中k个部件正常工作就可以保证系统正常工作;系统高度集成,当某个部件需要替换时,必须对所有部件同时进行替换;系统负载平均分配到所有正常工作的部件上。
采用负载平均分配规则描述部件的相关性,负载平均分配是指系统总负载平均分配到所有工作的负载,根据所述负载平均分布规则,存活部件退化率μα与部件标准负载ωα的关系用式(1)表示;
lnμα=β01ωn-α (1);
其中,α为故障部件数量,并且
Figure GDA0002783381750000091
β0与β1为待拟合参数。该式子建立了部件负载和部件退化率之间的关系,基于该计算式在后续的多阶段随机过程模型中,减少了模型中的参数。通过多阶段随机过程模型描述系统部件的退化过程。所述多阶段随机过程包括多阶段维纳过程、多阶段伽马过程、多阶段逆高斯过程。本发明以多阶段维纳过程为例来进行详细说明。获取系统各部件退化水平的历史数据,利用多阶段维纳模型进行数据拟合,所述多阶段维纳模型表示为x=xα(t)(α=0,1,...,n-k),其中,α个部件故障的情况下,存活部件的退化水平与时间的变化关系为xα(t);
Figure GDA0002783381750000092
其中,σ为多阶段维纳过程的噪声参数。
基于所述多阶段维纳模型,计算第α次(α=1,2,...,n-k+1)部件故障时间tα的分布函数fα(t)和部件在任意时间t的退化水平dt的概率密度函数
Figure GDA0002783381750000093
Figure GDA0002783381750000101
Figure GDA0002783381750000102
其中,H是部件退化水平上限。
根据部件故障数量α定义系统状态包括正常态、加速退化态和故障态,即系统状态sα为:
Figure GDA0002783381750000103
预先设置可靠度最低阈值L,通过离线检测部件退化水平计算系统的可靠度。具体地,采用离线检测的方式,根据i-1离线检测得到系统各部件的退化水平{dj|j=1,2,...,n},通过式(5)计算在i检测时的系统可靠度R,
Figure GDA0002783381750000104
根据系统状态和系统可靠度,定义一个检测周期内所有可能发生的状况,状况包括:
第一情况:当第i-1次检测时系统处于正常态以及所述系统可靠度大于预设的系统可靠度最低阈值L,并且第i次检测时系统处于正常态以及所述系统可靠度R大于所述系统可靠度最低阈值L;
第二情况:当第i-1次检测时系统处于正常态以及系统可靠度大于所述系统可靠度最低阈值L,并且第i次检测时系统处于正常态以及系统可靠度小于所述系统可靠度最低阈值L;
第三情况:当第i-1次检测时系统处于正常态以及系统可靠度大于所述系统可靠度最低阈值L,并且第i次检测时系统处于加速退化态;
第四情况:当第i-1次检测时系统处于正常态以及系统可靠度大于所述系统可靠度最低阈值L,并且第i次检测时系统处于故障态。
根据概率密度函数,分别计算上述四种情况发生的概率。
通过式(6)计算所述第一情况发生的概率
Figure GDA0002783381750000111
PA i=PA i-1×Pr{N(Ti)=0}×Pr{R(Ti|N(Ti)=0)>L} (6);
其中,
Figure GDA0002783381750000112
为1,Ai-1表示在第i-1次检测和第i次检测之间,第一情况发生,
故障数量N为α的概率为:
Figure GDA0002783381750000113
系统可靠度大于所述系统可靠度最低阈值L的概率为:
Figure GDA0002783381750000114
通过式(9)计算所述第二情况发生的概率PB i
PB i=PA i-1×Pr{N(Ti)=0}×Pr{R(Ti|N(Ti)=0)<L} (9);
通过式(10)计算所述第三情况发生的概率PC i
Figure GDA0002783381750000121
其中Pj i=PA i-1×Pr{N(Ti)=j};
通过式(11)计算所述第四情况发生的概率PD i
Figure GDA0002783381750000123
根据所述每一个可能发生的状况的概率,计算一个系统生命周期内的系统总运行成本和系统生命周期长度。
根据所述概率PB i、概率PC i以及概率PD i,计算一个系统生命周期内系统检测次数的期望EIN
Figure GDA0002783381750000124
根据所述概率PB i和概率PC i,计算系统生命周期以预防性维修结束的概率PPM
Figure GDA0002783381750000125
根据所述概率PD i,计算系统生命周期以事后维修结束的概率PCM
Figure GDA0002783381750000126
根据所述部件故障时间的分布函数,计算系统在各个系统状态下运行时间的期望Tj
Figure GDA0002783381750000127
其中,
Figure GDA0002783381750000131
根据所述一个系统生命周期内系统检测次数的期望EIN、系统生命周期以预防性维修结束的概率PPM、系统生命周期以事后维修结束的概率PCM,计算得到系统一个生命周期的总运行成本期望EC
Figure GDA0002783381750000132
其中,CIN为一次检测的费用,CCM为一次故障性替换的费用,CPM为一次预防性替换的费用,CADDα为系统在α个部件故障的情况下运行单位时间的费用;
计算得到系统生命周期长度ET
Figure GDA0002783381750000133
根据所述系统总运行成本和系统生命周期长度,计算系统单位时间内的运行成本,得到系统的最优检测时间序列和最优可靠度阈值。根据所述系统的总运行成本期望EC和系统生命周期长度ET,计算得到系统单位时间运行成本:
Figure GDA0002783381750000134
对所述系统单位时间运行成本进行优化处理,得到所述最优检测时间序列T*和最优可靠度阈值L*
(L*,T*)=minCPT (20)。
最优可靠度阈值L*对系统可靠度最低阈值L进行了优化。基于最优检测时间和最优可靠度阈值,制定相应的系统预防性维修策略。
如图2所示的本发明的一个实施例,本发明提供了一种k/n负载均担系统的维修决策系统,所述系统包括:
模型构建模块20,基于一描述部件相关性的负载平均分配规则,以及一构建的描述系统部件退化过程的多阶段随机过程模型,计算部件故障时间的分布函数;
概率计算模块21,根据离线检测时的系统状态和系统可靠度,定义一个检测周期内所有可能发生的状况,并根据所述部件故障时间的分布函数,计算每一个可能发生的状况的概率,其中,所述系统状态包括正常态、加速退化态和故障态;
成本计算模块22,根据所述每一个可能发生的状况的概率,计算一个系统生命周期内的系统总运行成本和系统生命周期长度;
模型输出模块23,根据所述系统总运行成本和系统生命周期长度,计算系统单位时间内的运行成本,得到系统的最优检测时间序列和最优可靠度阈值。
模型构建模块采用负载平均分配规则描述部件的相关性,获取系统各部件退化水平的历史数据。构建多阶段维纳模型,利用多阶段维纳模型进行数据拟合,根据存活部件的退化水平与时间的变化关系,计算第α次(α=1,2,...,n-k+1)部件故障时间的分布函数和部件在任意时间退化水平的概率密度函数。
概率计算模块中,定义系统状态包括正常态、加速退化态和故障态。预先设置可靠度最低阈值L,通过离线检测部件退化水平计算系统的可靠度,根据系统状态和系统可靠度,定义一个检测周期内所有可能发生的状况,并根据所述部件故障时间的分布函数,计算每一个可能发生的状况的概率。
成本计算模块根据所述每一个可能发生的状况的概率,计算一个系统生命周期内系统检测次数的期望,并根据所述部件故障时间的分布函数,计算系统在各个系统状态下运行时间的期望和系统生命周期长度。
模型输出模块根据所述系统总运行成本和系统生命周期长度,计算系统单位时间内的运行成本,对所述系统单位时间运行成本进行优化处理,得到所述最优检测时间序列和最优可靠度阈值。
尽管为示例目的,已经公开了本发明的优选实施方式,但是本领域的普通技术人员将意识到,在不脱离由所附的权利要求书公开的本发明的范围和精神的情况下,各种改进、增加以及取代是可能的。

Claims (4)

1.一种k/n负载均担系统的维修决策方法,该系统由n个相同的不可维修部件组成,其中k个部件正常工作保证系统正常工作,其特征在于,所述方法包括:
S1、基于一描述部件相关性的负载平均分配规则,以及一构建的描述系统部件退化过程的多阶段随机过程模型,计算部件故障时刻的分布函数;
S2、根据离线检测时的系统状态和系统可靠度,定义一个检测周期内所有可能发生的状况,并根据所述部件故障时刻的分布函数,计算每一个可能发生的状况的概率,其中,所述系统状态包括正常态、加速退化态和故障态;
S3、根据所述每一个可能发生的状况的概率,计算一个系统生命周期内的系统总运行成本和系统生命周期长度;
S4、根据所述系统总运行成本和系统生命周期长度,计算系统单位时间内的运行成本,得到系统的最优检测时间序列和最优可靠度阈值;
所述步骤S1包括:
利用多阶段维纳模型进行数据拟合,所述多阶段维纳模型表示为x=xα(t),α个部件故障的情况下,存活部件的退化水平与时刻t的变化关系为xα(t);
Figure FDA0003862939200000011
其中,α=0,1,…,n-k+1,σ为多阶段维纳过程的噪声参数,tα为α个部件故障时刻,μα为α个部件故障后存活部件的退化率,t0=t-1=0,μn-k+1=μ-1=0,W(t)为标准布朗运动;
所述步骤S2包括:
采用离线检测的方式,根据i-1次离线检测时刻Ti-1得到系统各部件的退化水平
Figure FDA0003862939200000012
通过式(5)计算在i次检测时刻Ti的系统可靠度R(Ti):
Figure FDA0003862939200000013
定义一个检测周期内所有可能发生的状况,包括:
第一情况:当第i-1次检测时系统处于正常态以及所述系统可靠度大于预设的系统可靠度最低阈值L,并且第i次检测时系统处于正常态以及所述系统可靠度R大于所述系统可靠度最低阈值L;
第二情况:当第i-1次检测时系统处于正常态以及系统可靠度大于所述系统可靠度最低阈值L,并且第i次检测时系统处于正常态以及系统可靠度小于所述系统可靠度最低阈值L;
第三情况:当第i-1次检测时系统处于正常态以及系统可靠度大于所述系统可靠度最低阈值L,并且第i次检测时系统处于加速退化态;
第四情况:当第i-1次检测时系统处于正常态以及系统可靠度大于所述系统可靠度最低阈值L,并且第i次检测时系统处于故障态;
通过式(6)计算所述第一情况发生的概率
Figure FDA0003862939200000021
PA i=PA i-1×Pr{N(Ti)=0}×Pr{R(Ti|N(Ti)=0)>L} (6);
其中,
Figure FDA0003862939200000022
Ai-1表示在第i-1次检测和第i次检测之间,第一情况发生,用N(t)表示t时刻故障部件数量,则t时刻N(t)为α的概率为:
Figure FDA0003862939200000023
已知无部件故障,系统可靠度大于所述系统可靠度最低阈值L的概率为:
Figure FDA0003862939200000024
其中,函数1{R(t|(N(t’)=0)>L}表示无部件故障且t’时刻每个部件退化量已知情况下系统可靠度R(t)大于可靠度最低阈值L的概率,t’为离线检测时刻,t为离线检测后需要预测的任一时刻,因此t’<t;
通过式(9)计算所述第二情况发生的概率PB i
PB i=PA i-1×Pr{N(Ti)=0}×Pr{R(Ti|N(Ti)=0)<L} (9);
通过式(10)计算所述第三情况发生的概率PC i
Figure FDA0003862939200000031
其中Pj i=PA i-1×Pr{N(Ti)=j};
通过式(11)计算所述第四情况发生的概率PD i
Figure FDA0003862939200000032
所述步骤S3包括:
根据第二情况发生的概率PB i、第三情况发生的概率PC i以及第四情况发生的概率PD i,计算一个系统生命周期内系统检测次数的期望EIN:
Figure FDA0003862939200000033
其中,i表示第i次检测;
根据概率PB i和概率PC i,计算系统生命周期以预防性维修结束的概率PPM
Figure FDA0003862939200000034
根据概率PD i,计算系统生命周期以事后维修结束的概率PCM
Figure FDA0003862939200000035
根据部件故障时刻的分布函数,计算系统在各个系统状态下运行时间的期望Tj
Figure FDA0003862939200000036
Tj i是第i次检测至第i+1次检测内系统运行在状态sα下运行时间,由式(16)计算得到,Pj i是第i次检测至第i+1次检测内系统发生j次部件故障的可能性;
Figure FDA0003862939200000037
其中,Ti+1是第i+1次离线检测时刻,函数fj(t)是tj的概率密度函数,函数Fj(t)是tj的累积分布函数,可以根据
Figure FDA0003862939200000038
计算得到,根据部件故障数量α定义系统状态包括正常态、加速退化态和故障态,即系统状态sα为:
Figure FDA0003862939200000041
根据一个系统生命周期内系统检测次数的期望EIN、系统生命周期以预防性维修结束的概率PPM、系统生命周期以事后维修结束的概率PCM,计算得到系统一个生命周期的总运行成本期望EC
Figure FDA0003862939200000042
其中,CIN为一次检测的费用,CCM为一次故障性替换的费用,CPM为一次预防性替换的费用,CADDj为系统在j个部件故障的情况下运行单位时间的费用,当n-k+1次部件故障后,系统自动停止运行,CADD n-k+1为0;
计算得到系统生命周期长度的期望ET
Figure FDA0003862939200000043
所述步骤S4包括:
根据系统的总运行成本期望EC和系统生命周期长度的期望ET,计算得到系统单位时间运行成本:
Figure FDA0003862939200000044
对系统单位时间运行成本进行优化处理,得到所述最优检测时间序列T*和最优可靠度阈值L*
(L*,T*)=minCPT (20)。
2.如权利要求1所述的k/n负载均担系统的维修决策方法,其特征在于,所述步骤S1包括:
根据所述负载平均分配规则,存活部件退化率μα与部件标准化负载ωα的关系用式(1)表示;
lnμα=β01ωn-α (1);
其中,α为故障部件数量,α=0,1,…,n-k+1,并且
Figure FDA0003862939200000051
β0与β1为待拟合参数。
3.如权利要求2所述的k/n负载均担系统的维修决策方法,其特征在于,所述步骤S1还包括:
基于多阶段维纳模型,计算α个部件故障时刻tα的概率密度函数fα(t)为;
Figure FDA0003862939200000052
其中,H是部件退化水平上限,
Figure FDA0003862939200000056
为标准正态分布的概率密度函数,Ф(·)为标准正态分布的分布函数,当没有部件故障即α=0时有f0(t)=1;
将系统中的n个部件编号为1,2,···n,无部件故障情况下在任意时刻t编号为m的部件的退化水平dm t的概率密度函数
Figure FDA0003862939200000053
Figure FDA0003862939200000054
其中,
Figure FDA0003862939200000055
是编号为m的部件在t时刻前的t’时刻的退化量。
4.一种k/n负载均担系统的维修决策系统,其特征在于,所述系统执行权利要求1-3任一所述的k/n负载均担系统的维修决策方法,所述系统包括:模型构建模块,基于一描述部件相关性的负载平均分配规则,以及一构建的描述系统部件退化过程的多阶段随机过程模型,计算部件故障时刻的分布函数;
概率计算模块,根据离线检测时的系统状态和系统可靠度,定义一个检测周期内所有可能发生的状况,并根据所述部件故障时刻的分布函数,计算每一个可能发生的状况的概率,其中,所述系统状态包括正常态、加速退化态和故障态;
成本计算模块,根据所述每一个可能发生的状况的概率,计算一个系统生命周期内的系统总运行成本和系统生命周期长度;
模型输出模块,根据所述系统总运行成本和系统生命周期长度,计算系统单位时间内的运行成本,得到系统的最优检测时间序列和最优可靠度阈值。
CN202011130675.XA 2020-10-21 2020-10-21 一种k/n负载均担系统的维修决策方法及系统 Active CN112257880B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011130675.XA CN112257880B (zh) 2020-10-21 2020-10-21 一种k/n负载均担系统的维修决策方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011130675.XA CN112257880B (zh) 2020-10-21 2020-10-21 一种k/n负载均担系统的维修决策方法及系统

Publications (2)

Publication Number Publication Date
CN112257880A CN112257880A (zh) 2021-01-22
CN112257880B true CN112257880B (zh) 2022-11-22

Family

ID=74263732

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011130675.XA Active CN112257880B (zh) 2020-10-21 2020-10-21 一种k/n负载均担系统的维修决策方法及系统

Country Status (1)

Country Link
CN (1) CN112257880B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117875947A (zh) * 2024-03-11 2024-04-12 浙江大学 k/n负载均担系统的可靠性评估和维修决策方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101807265A (zh) * 2008-09-30 2010-08-18 洛克威尔自动控制技术股份有限公司 机器选择、集成和使用的动态多目标优化的系统和方法
CN107229979A (zh) * 2017-04-17 2017-10-03 北京航空航天大学 一种可修退化系统周期性预防维修策略的优化方法
CN109102085A (zh) * 2018-07-09 2018-12-28 华中科技大学 一种数控机床的非周期性预防维修方法
CN110073301A (zh) * 2017-08-02 2019-07-30 强力物联网投资组合2016有限公司 工业物联网中具有大数据集的数据收集环境下的检测方法和系统
CN110909442A (zh) * 2019-10-11 2020-03-24 北京航空航天大学 一种考虑维修负面效果的维修效果建模方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101807265A (zh) * 2008-09-30 2010-08-18 洛克威尔自动控制技术股份有限公司 机器选择、集成和使用的动态多目标优化的系统和方法
CN107229979A (zh) * 2017-04-17 2017-10-03 北京航空航天大学 一种可修退化系统周期性预防维修策略的优化方法
CN110073301A (zh) * 2017-08-02 2019-07-30 强力物联网投资组合2016有限公司 工业物联网中具有大数据集的数据收集环境下的检测方法和系统
CN109102085A (zh) * 2018-07-09 2018-12-28 华中科技大学 一种数控机床的非周期性预防维修方法
CN110909442A (zh) * 2019-10-11 2020-03-24 北京航空航天大学 一种考虑维修负面效果的维修效果建模方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
Opportunistic maintenance based on multi-dependent components of manufacturing system;IUNG B;《CIRP Annals》;20161231;第65卷(第1期);第401-404页 *
基于云计算的设备状态维护关键技术研究;杨明;《中国优秀硕士学位论文全文数据库 (信息科技辑)》;20180115(第1期);全文 *
基于服务性能合同的多部件系统维修策略模型优化;许飞雪;《计算机应用》;20200922;第1184-1191页 *
状态检修决策中的电气设备故障率推算;潘乐真;《电力自动化设备》;20101231;第30卷(第2期);第91-94页 *
非完美特性下的多状态系统检测与维修优化;李志栋;《北京航空航天大学学报》;20171231;第43卷(第5期);第951-960页 *

Also Published As

Publication number Publication date
CN112257880A (zh) 2021-01-22

Similar Documents

Publication Publication Date Title
CN110991666B (zh) 故障检测方法、模型的训练方法、装置、设备及存储介质
CN109033499B (zh) 一种多阶段一致性检验的航空发动机剩余寿命预测方法
US6981182B2 (en) Method and system for analyzing fault and quantized operational data for automated diagnostics of locomotives
JP4792047B2 (ja) システム可用性の計算方法
CN108764551B (zh) 基于系统级寿命信息的成组维修决策方法
JP3968656B2 (ja) プラント機器の保守支援装置
Frostig et al. On the availability of R out of N repairable systems
CN112257880B (zh) 一种k/n负载均担系统的维修决策方法及系统
CN111382494A (zh) 用于检测工业机器的感知数据中的异常的系统及方法
JP6482743B1 (ja) リスク評価装置、リスク評価システム、リスク評価方法、及び、リスク評価プログラム
CN115270078B (zh) 一种机电设备平均修复时间的计算方法和系统
KR20150084613A (ko) 데이터 마이닝을 이용한 타워크레인 고장예측방법
CN113487182B (zh) 设备健康状态评估方法、装置、计算机设备和介质
CN115453356A (zh) 一种动力设备运行状态监测分析方法、系统、终端及介质
CN116820820A (zh) 服务器故障监测方法及系统
JP6574533B2 (ja) リスク評価装置、リスク評価システム、リスク評価方法、及び、リスク評価プログラム
CN107730070B (zh) 一种基于状态监测的智能变电站继电保护系统检修方法
Koutras et al. Semi-Markov availability modeling of a redundant system with partial and full rejuvenation actions
CN114841382A (zh) 设备故障检测方法、装置及存储介质
CN109739210B (zh) 设备部件健康状态的评估方法及装置
JP6482742B1 (ja) リスク評価装置、リスク評価システム、リスク評価方法、及び、リスク評価プログラム
CN112799911A (zh) 一种节点健康状态检测方法、装置、设备及存储介质
CN115163513B (zh) 一种地铁潜污泵的故障诊断方法和装置
CN113835976B (zh) 针对互联网数据中心进行告警处理的方法及装置
CN110232482B (zh) 基于神经网络的设备管理方法及装置

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