浙江大学学报(工学版), 2020, 54(7): 1298-1307 doi: 10.3785/j.issn.1008-973X.2020.07.007

自动化技术、计算机技术

基于延迟时间模型的不完全检修计划优化模型

刘葛辉,, 陈绍宽,, 金华, 刘爽, 彭宏勤

Optimum imperfect inspection and maintenance scheduling model considering delay time theory

LIU Ge-hui,, CHEN Shao-kuan,, JIN Hua, LIU Shuang, PENG Hong-qin

通讯作者: 陈绍宽,男,教授,博导. orcid.org/0000-0002-2097-9936. E-mail: shkchen@bjtu.edu.cn

收稿日期: 2019-07-2  

Received: 2019-07-2  

作者简介 About authors

刘葛辉(1995—),男,博士生,从事系统可靠度分析和维修计划研究.orcid.org/0000-0003-1270-6011.E-mail:16114221@bjtu.edu.cn , E-mail:16114221@bjtu.edu.cn

摘要

针对现有基于检测的状态维修方法仅适用于短期维修计划,且通常单一地考虑不完全维修或不完全检测的问题,构建适合设备维修现场特征的长期维修计划优化模型,考虑检测与维修活动均不完全情况下的退化过程. 利用延迟时间模型描述退化过程,围绕双重不确定的检修活动,构建基于递推关系的可靠度模型,该模型可以有效表示系统退化速度和故障发生概率. 以单位时间系统费用最小为决策目标,引入可靠度、可用度等约束,建立检修计划优化模型,求解得到最佳检测周期和更换周期. 案例研究表明,利用提出的模型能够有效优化系统检修计划,节省维修成本. 通过与固定周期和固定可靠度阈值的维修策略进行对比,说明检测精度对于检测模型的优化效果具有显著影响,当检测精度属于较高水平时,检修模型明显优于其他维修模型.

关键词: 延迟时间模型 ; 预防性检修 ; 不完全检修 ; 检测周期 ; 可靠度 ; 维修费用

Abstract

The existing condition-based maintenance methods based on inspection activities were only suitable for short-term schedules and involved only imperfect inspection or maintenance. A long-term optimization model was proposed based on on-site features of system maintenance scheduling by incorporating deterioration process with both imperfect maintenance and inspection activities. The delay time theory was introduced to describe deterioration process of system with double contingencies from imperfect maintenance and inspection. The proposed recursive reliability model can accurately evaluate the deterioration rate and failure rate of system. A maintenance scheduling optimization model was applied to minimize the average cost by searching appropriate inspection cycles and replacement strategy. The reliability and availability constraints were considered to meet the requirement of maintenance. Case studies show that the proposed model can attain optimal inspection strategy with minimum system cost. The accuracy of inspection was crucial for the scheduling optimization model compared with age-based and reliability-based maintenance models. The proposed model with a high accuracy of inspection reached the minimum cost compared to models without inspections.

Keywords: delay time model ; preventive maintenance ; imperfect maintenance ; inspection cycle ; reliability ; maintenance cost

PDF (1296KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

刘葛辉, 陈绍宽, 金华, 刘爽, 彭宏勤. 基于延迟时间模型的不完全检修计划优化模型. 浙江大学学报(工学版)[J], 2020, 54(7): 1298-1307 doi:10.3785/j.issn.1008-973X.2020.07.007

LIU Ge-hui, CHEN Shao-kuan, JIN Hua, LIU Shuang, PENG Hong-qin. Optimum imperfect inspection and maintenance scheduling model considering delay time theory. Journal of Zhejiang University(Engineering Science)[J], 2020, 54(7): 1298-1307 doi:10.3785/j.issn.1008-973X.2020.07.007

预防性维修(preventive maintenance)是保证系统状态、提高运行质量的有效方法,合理的预防性维修计划对于减少系统故障和提高可用度有重要意义[1]. 状态维修(condition-based maintenance)方法根据检测的状态调整维修计划,以精确地实施预防性维修[2]. 检测是状态维修方法的重要内容,已有研究表明预防性维修活动中约80%的任务都是由检测识别并解决的,检测间隔过短或过长会引起维修过度或维修不足等问题,所以合理的检测计划对系统的维护和运行有重要作用[3].

在状态维修模型中,随机过程方法可以直接描述退化过程[4-5]. 在实际的维修工作中,常常难以明确具体的退化状态,需要对退化状态进行离散处理[6]. 延迟时间模型(delay time model)在多阶段退化方面应用广泛,该模型在1984年首次提出并应用于生产线检修问题[7]. 目前,该模型被广泛应用于检修计划优化,如优化检测间隔[8]、检测和更换策略[9],并开始应用于多部件系统[10]. 这一类研究通常假设维修效果是完全的,即检出缺陷或发生故障时进行更换. 在实际情况中维修通常是不完全的(imperfect maintenance)[11],即维修后系统状态不能恢复至全新. 有效役龄方法是表示不完全维修的常用方法[12],并已在延迟时间模型中得到应用[6, 13-14]. 此外,实际维修工作中检测不一定是准确的,不完全检测(imperfect inspection)概念最早由Christer等[15]提出,并逐渐发展为多等级检测[16]、时变检测误差[17]、多阶段延迟时间等模型[18]. 目前,综合考虑不完全检测和维修的研究较少,且均假设故障时进行更换,无法有效表达系统的故障概率情况,难以实现长期检修计划的优化[19].

在既有研究的基础上,本文提出同时考虑不完全检测和不完全维修的可靠度模型,以有效描述系统退化过程和故障发生概率. 以最小化系统费用率为目标,构建检修计划优化模型,通过分析检测结果和故障概率,考虑可靠度和可用度等实际约束,同时优化系统检修和更换策略. 通过与既有研究中不考虑检测的维修策略进行对比,获得应用检修策略的必要条件.

1. 问题描述和维修策略

根据延迟时间模型[6],设正常阶段时间为u,延迟时间阶段时间为v,此时系统寿命为u+v. 定义随机变量uv的密度函数和累计分布分别为f1u)、F1u)和f2v)、F2v),并且2个阶段对应的发生率分别为λu)和hu),称为缺陷率和缺陷发生后的故障率. 系统维修计划取决于检测的状态,当检测需要成本时,有必要首先确定最优的检测策略. 此外,维修不完全性导致系统经多次维修后须对系统进行更换,因此检修计划的另一个重要部分为更换策略.

设系统以固定周期T进行检测,检测只能以概率r(0<r≤1)识别已发生的缺陷,但不会将正常情况识别为缺陷. 当检测到缺陷时,系统进行预防性维修,维修后缺陷将被修复. 若检测不到缺陷,则不进行维修. 维修是不完全的,所以维修后系统不能恢复到全新状态,采用有效役龄方法描述维修后的状态[20],记役龄回退因子为a(0≤a≤1),表示维修后系统有效役龄的缩减比例. 若系统累计运行时间达到最大可用时间TC或累计检测次数到达次数上限τ,则进行预防性更换(preventive replacement),更换后系统的状态为全新. 最大可用时间为系统的固有属性,也称为技术寿命,若在之前某时刻进行更换可降低费用,则这一时刻称为经济寿命[21]. 系统在运行中发生故障时,通过最小维修(minimal repair)恢复到可运行状态,最小维修不改变系统退化情况,缺陷须在下一次检测时通过预防性维修排除.

以上过程中,检测、预防性维修、预防性更换、最小维修(即恢复性维修)的费用分别记为cincpcrcc,相应停机时间分别记为tintptrtc. 最小维修费用包括了故障产生的期望损失费用,因此高于预防性维修和更换的费用. 在TC已知的情况下,检修策略定义为检测周期T和检测上限次数τ. 当检测周期过长时,缺陷不能被有效识别,系统易发生故障;若周期过短,则增加额外的维修成本. 需要通过确定合理的检测间隔和更换策略,以最小化系统平均维修费用.

为了方便模型构建,提出以下假设.

假设1 系统的初始状态为全新.

假设2 系统在一个更新周期内存在最低可靠度约束,表示实际中最大故障概率约束[1-3].

假设3 检修或故障引起的系统停机时间产生停机损失费用,为了保证系统的可使用时间,系统存在最低可用度约束[19, 22].

假设4 系统检修时间和故障处理时间较短,因此不考虑停机时间对系统退化的影响.

2. 不完全检修可靠度模型

维修优化模型通常包括可靠度模型(即退化模型)和维修计划优化模型[22]. 首先构建考虑不完全检修的可靠度模型. 检测的不完全性使预防性维修成为随机事件,且在不同时间维修后系统的退化速率不同. 另外,不完全检测使已经发生的缺陷可能延长至下一个检修间隔,系统状态不满足马尔可夫性,因此需要构建基于递推关系的可靠度模型.

2.1. 可靠度模型

2.1.1. 役龄回退模型

通过役龄回退因子描述不完全维修的效果[20],设系统在第k次检修(时间为tk)时累计寿命为Tk,则不完全维修后系统的累计寿命(即有效役龄)恢复为aTk. 假设役龄回退因子同时作用于正常阶段和延迟时间阶段,则系统在tk进行维修后的缺陷率λk+1u)和独立故障率hk+1u)变化分别如下.

${\lambda _{k + 1}}(u) = \left\{\!\!\! {\begin{array}{*{20}{l}} {{\lambda _k}(u + a{T_k}) = \lambda (u + a k T)} {\text{,}}&{{\text{在}}{t_k}{\text{时刻维修}}};\\ {{\lambda _k}(u + {T_k}){\text{,}}}&{{\text{不维修.}}} \end{array}} \right.$

${h_{k + 1}}(v) = \left\{\!\!\! {\begin{array}{*{20}{l}} {{h_k}(v + a{T_k}) = h(v + a k T)} {\text{,}}&{{\text{在}}{t_k}{\text{时刻维修}}};\\ {{h_k}(v + {T_k}) {\text{,}}}&{{\text{不维修.}}} \end{array}} \right.$

式中:Tk为上次维修后系统的累计使用时间.

2.1.2. 缺陷检出概率

考虑在tii>0)时刻对系统进行检测,检测到缺陷并进行预防性维修的概率为Pdti). 设上一次发生预防性维修的时间为tk(1<ki−1),如图1所示. 在(tl−1tl)某一时刻发生缺陷(tktl−1<tlti),且在检测时刻ti前没有发生故障的概率如下.

图 1

图 1   检测缺陷过程示意图

Fig.1   Sketch map of defect detection


$\begin{array}{*{20}{l}} {{P_{\rm{d}}}({t_i}|{t_k}) = }\\ {r\sum\limits_{l = k + 1}^i {\left\{ {{{(1 - r)}^{i - l}}\int_{{t_{l - 1}} - {t_k}}^{{t_l} - {t_k}} {{g_{k + 1}}(u)[1 - {F_{k + 1}}({t_i} - {t_k} - u)]{\rm{d}}u} } \right\}} .} \end{array}$

若检出概率r=1,则式(3)改写为

${P_{\rm{d}}}({t_i}|{t_k}) = \int_{{t_{i - 1}} - {t_k}}^{{t_i} - {t_k}} {{g_{k + 1}}(u)[1 - {F_{k + 1}}({t_i} - {t_k} - u)]{\rm{d}}u} ,$

即被检出的缺陷只能发生在时段(ti−1ti).

2.1.3. 故障发生概率

考虑上一次预防性维修在tk时刻发生,系统在检测间隔(ti−1ti)中的td时刻发生故障的情况,如图2所示. 缺陷可能发生在(tl−1tl)中的某一时刻,其中tktl−1<tlti. 该情形发生的概率如下:

图 2

图 2   发生故障过程示意图

Fig.2   Sketch map of failure occurrence


$\begin{array}{l} {P_{\rm{f}}}({t_i}|{t_k}) = \sum\limits_{l = k + 1}^i {\left\{ {{{(1 - r)}^{i - l}}\int_{{t_{l - 1}} - {t_k}}^{{t_l} - {t_k}} {{g_{k + 1}}(u)} \times } \right.} \\ \left. {[{F_{k + 1}}({t_i} - {t_k} - u) - {F_{k + 1}}({t_{i - 1}} - {t_k} - u)]{\rm{d}}u} \right\}. \end{array}$

系统在(ti−1ti)时段中发生故障后,会在ti时刻通过预防性维修排除缺陷. 在ti时刻进行预防性维修的条件概率为

${P_{\rm{m}}}({t_i}|{t_k}) = {P_{\rm{d}}}({t_i}|{t_k}) + {P_{\rm{f}}}({t_i}|{t_k}).$

ti时刻进行预防性维修的概率为

$\begin{split} {P_{\rm{m}}}({t_i}) = & \sum\limits_{k = 0}^{i - 1} {{P_{\rm{m}}}({t_k}){P_{\rm{m}}}({t_i}|{t_k})} = \\ & \sum\limits_{k = 0}^{i - 1} {{P_{\rm{m}}}({t_k})[{P_{\rm{d}}}({t_i}|{t_k}) + {P_{\rm{f}}}({t_i}|{t_k})]} . \end{split}$

在系统开始使用时(k=0),有Pmt0)=1.

2.1.4. 可靠度函数

可靠度Rt)表示从初始时刻至t时刻系统不发生故障的概率. 为了计算任意时刻的可靠度,考虑ti−1t时刻中(ti−1<t<ti)不发生故障的概率. 在上一次维修发生在tk时刻的情况下,系统在时段(ti−1t)发生故障的概率如下:

${P_{\rm{f}}}(t|{t_{i{\rm{ - }}1}}) = \sum\limits_{k = 1}^{i - 1} {{P_{\rm{m}}}({t_k})} {P_{\rm{f}}}(t|{t_k};{t_{i{\rm{ - }}1}}).$

式中:

ti−1t)不发生故障的概率为

${P_{\rm{R}}}(t|{t_{i - 1}}) = 1 - {P_{\rm{f}}}(t|{t_{i - 1}}).$

因此,系统在tt>ti−1)时刻的可靠度为

$R(t) = R({t_{i - 1}}){P_{\rm{R}}}(t|{t_{i - 1}}).$

2.2. 参数估计方法

在不完全检修模型中,除了缺陷分布和故障分布参数,还需确定检出概率和役龄回退因子. 基于文献[14],提出维修和检测间隔均不固定的情况下的参数估计方法.

在实际检测中存在维修和不维修2种情况,记系统的实际维修时间序列为 $(t_1^0,t_2^0, \cdots ,t_n^0)$,设相邻维修 $t_k^0{\text{、}}t_{k + 1}^0$之间共有mk)次检测未检测到缺陷,即 $t_{k,m(k)}^0$k<n)的下一次检测时刻为 $t_{k + 1}^0$. 在检测时刻 $t_{k,i}^0,i = 1,\cdots,m(k),k = 1,\cdots,n - 1$检测到缺陷的概率为

$\begin{split} & {P_{\rm{d}}}(t_{k,i}^0|t_k^0) = r\sum\limits_{l = 1}^i {\left\{ {{{(1 - r)}^{i - l}}\int_{t_{k,l - 1}^0 - t_k^0}^{t_{k,l}^0 - t_k^0} {{g_{k + 1}}(u) \times } } \right.} \\ & {[1 - {F_{k + 1}}(t_{k,i}^0 - t_k^0 - u)]{\rm{d}}u} \Bigg\};\;i = 1,\cdots,m(k). \end{split}$

对于检出概率的似然函数为

${L_1} = \sum\limits_{k = 1}^{n - 1} {\left[\sum\limits_{i = 1}^{m(k)} {\ln\; (1 - {P_{\rm{d}}}(t_{k,i}^0|t_k^0))} + \ln\; ({P_{\rm{d}}}(t_{k + 1}^0|t_k^0))\right]}. $

记系统的实际故障序列为 $(x_1^0,x_2^0, \cdots ,x_m^0)$,实际中故障发生次数较少,为了便于计算,将故障发生时间转换为在检测间隔中是否发生故障,即 $x_{k,i}^0$=1表示在相邻2次检测 $t_{k,i}^0{\text{、}}t_{k,i + 1}^0$中发生故障,这一情形的发生概率为

$\begin{split} & {P_{\rm{f}}}(x_{k,i}^0|t_k^0) = \\ & \sum\limits_{l = 1}^i {\left\{ {{{(1 - r)}^{i - l}}\int_{t_{k,l - 1}^0 - t_k^0}^{t_{k,l}^0 - t_k^0} {{g_{k + 1}}(u)[{F_{k + 1}}(t_{k,i + 1}^0 - t_k^0 - u) - } } \right.} \\ & \left. {{F_{k + 1}}(t_{k,i}^0 - t_k^0 - u)]{\rm{d}}u} \right\};\; i = 0,\cdots ,m(k). \end{split}$

对于故障发生的似然函数为

$\begin{split} {L_2} = & \sum\limits_{k = 1}^{n - 1} {\left\{ {\left( {\sum\limits_{i = 0}^{m(k)} {\ln\; (1 - x_{k,i}^0)} (1 - {P_{\rm{f}}}(x_{k,i}^0|t_k^0))} \right) + } \right.} \\ & \left. {x_{k,i}^0\ln\; ({P_{\rm{f}}}(x_{k,i}^0|t_k^0))} \right\}. \end{split}$

综合以上,构建似然函数:

$L = {L_1} + {L_2}.$

3. 检修计划优化模型

3.1. 检修策略模型

在检修策略下,系统的检测间隔为T,检测上限次数为τ. 系统第i次检测的时间为

${t_i} = i T;\quad i = 1,2,\cdots,\tau .$

系统的更新周期受到τ和TC的约束,定义示性变量σσ=1表示系统的更新周期为TC.

$\sigma = \left\{ {\begin{array}{*{20}{c}} 1,&{T \tau \geqslant {\rm {TC}},\; R({\rm {TC}}) \geqslant {R_{{\rm{min}}}}}; \\ 0,&{\text{其他}}. \end{array}} \right.$

在更新周期内,系统的总停机时间TDown包括检测时间、预防性维修时间、最小维修时间和更换时间. 系统一共进行τ−1次检测,因为在第τ次检测或TC时刻将直接执行更换.

$\begin{split} T_{\rm{Down}} = & (\tau - 1){t_{\rm{in}}} + \sum\nolimits_{i = 1}^{\tau - 1} {{P_{\rm{d}}}({t_i}) {t_{\rm{p}}} - } \\& {\left[ {\sigma \ln R({\rm{TC}}) + (1 - \sigma )\ln R({t_\tau })} \right] {t_{\rm{p}}} + {t_{\rm{r}}}.} \end{split}$

可用度A为可运行时间与更新周期的比值:

$A(T,\tau ) = 1 - \frac{{{T_{{\rm{Down}}}}(T,\tau )}}{{\sigma \times {\rm{TC}} + (1 - \sigma ) T \tau }}.$

根据上述分析和更新收益理论[9],建立以系统平均费用率最小为目标、考虑可靠度和可用度约束的检修计划优化模型(P1),目标函数与约束条件如下.

$ \begin{split} \min\; \bar C(T,\tau ) = & E({C_{{\rm{sum}}}})/E(L) = \\ & \left\{ {(\tau - 1){c_{\rm{in}}} + \sum\limits_{i = 1}^{\tau - 1} {{P_{\rm{d}}}({t_i}) {c_{\rm{p}}}} - [\sigma \ln R({\rm{TC}}) + } \right.\\ & {(1 - \sigma )\ln R({t_\tau })] {c_{\rm{c}}} + {c_{\rm{r}}} + {T_{{\rm{Down}}}}{c_{\rm{d}}}} \Bigg\}\Bigg/\\ & [\sigma \times {\rm{TC}} + (1 - \sigma )T\tau] . \end{split} $

$\begin{split} & {\rm{s.t.}}\\& \;\;\;{R_{{\rm{min}}}} \leqslant \left\{ {\begin{array}{*{20}{c}} {R({t_\tau })},&{{\rm{ }}T \tau < \rm TC}; \\ {R(\rm TC)},&{{\rm{ }}T \tau \geqslant \rm TC}; \end{array}} \right. \end{split}$

$A(T,\tau ) \geqslant {A_{\min }};$

$1 \leqslant T \leqslant {T_{\max }};$

$1 \leqslant \tau \leqslant \left\lceil {{\rm{TC}}/T} \right\rceil ;$

$T,\tau \in \bf N.$

式(20)为目标函数,表示总期望费用Csum与更新周期L的比值,费用包括检测费用、预防性维修费用、最小维修费用、更换费用和停机损失费用. 式(21)表示最低可靠度约束,系统在任意时刻的可靠度均不低于Rmin. 式(22)表示最低可用度约束. 式(23)~(25)为决策变量约束,检测间隔和检测上限次数均为正整数,且存在上限约束. 其中,Tmax为系统无维修退化至最低可靠度的时间.

3.2. 非检测策略模型

模型(P1)同时考虑了检测和维修活动. 为了验证检修模型的效果,基于文献[2324]构建2种不考虑检测的预防性维修模型,与模型(P1)的结果进行比较.

3.2.1. 固定周期维修模型

定周期维修策略按照固定时间间隔直接对系统进行预防性维修,其余内容均与模型(P1)相同. 维修策略定义为维修周期T1和最大维修次数τ1.

系统的检测时刻、更新周期和可用度的计算分别类似于式(16)、(17)、(19),其中停机时间的计算公式为

$\begin{split} & {T_{{\rm{Down}}}} = ({\tau _1} - 1) {t_{\rm{p}}} + {t_{\rm{r}}} + {t_{\rm{c}}} \times \\ & \left[ {\sum\limits_{i = 1}^{{\tau _{_1}} - 1} { - \ln R({t_i})} - \sigma \ln R({\rm TC}) - (1 - \sigma )\ln R ({t_{{\tau _1}}})} \right]. \end{split}$

优化模型(P2)构建如下:

$ \begin{split} & \min {{\bar C}_1}({T_1},{\tau _1}) = {{E({C_{{\rm{sum}}}})}/{E({L_1})}} =\\ & \left\{ {({\tau _1} - 1){c_{\rm{p}}} + \Bigg[\sum\limits_{i = 1}^{{\tau _1} - 1} { - \ln R({t_i})} -\sigma \ln R({\rm{TC}}) - } \right.\\ & \left. {\ln R({t_{{\tau _1}}})(1 - \sigma )\Bigg]{c_{\rm{c}}} + {c_{\rm{r}}} + {T_{{\rm{Down}}}}{c_{\rm{d}}}} \right\}\Bigg/\\ & [\sigma \times {\rm{TC}} + (1 - \sigma ){T_{\rm{1}}}{\tau _{\rm{1}}}]. \end{split} $

约束条件与式(21)~(25)相同.

3.2.2. 固定可靠度阈值维修模型

在固定可靠度阈值模型中,当上一次维修后的可靠度下降至R2时,系统需要进行预防性维修或更换,其余内容均与模型(P1)相同. 维修策略定义为可靠度阈值R2和最大维修次数τ2.

ρit)为系统第i次预防性维修后的故障率函数,则第i个维修间隔的表达式为

$\int_0^{{T_i}^*} {{\rho _{i - 1}}(t){\rm{d}}t} = - \ln {R_2}; \quad i = 1,2,\cdots,{\tau _2}.$

实际维修间隔为可靠度对应间隔的下界整数值:

${T_i} = \left\lfloor {{T_i}^*} \right\rfloor ;\quad i = 1,2,\cdots,{\tau _2}.$

由此可得系统第i次维修的时间:

${t_i} = \sum\limits_{j = 1}^i {{T_j}} .$

用于计算更新周期的示性变量σ表达式为

$\sigma = \left\{ {\begin{array}{*{20}{l}} 1,&{\displaystyle\sum\limits_{i = 1}^{{\tau _2}} {{T_i}} \geqslant {\rm{TC}},\;R({\rm{TC}}) \geqslant {R_{{\rm{min}}}}}; \\ 0,&{\text{其他}} . \end{array}} \right.$

在更新周期内,系统的可用度计算类似于式(19),其中停机时间的计算公式为

$\begin{split} & {{T_{{\rm{Down}}}} = ({\tau _2} - 1) {t_{\rm{p}}} + {t_{\rm{r}}} + {t_{\rm{c}}} \times }\\& {[{\tau _2}( - \ln {R_2})(1 - \sigma ) -\sigma \ln R({\rm{TC}}) ].} \end{split}$

优化模型(P3)构建如下:

$ \begin{split} \min \;{{\bar C}_2}({R_2},{\tau _2}) = & {{E({C_{{\rm{sum}}}})}/{E({L_2})}} = \\ & \{ {({\tau _2} - 1){c_{\rm{p}}} + [{\tau _2}( - \ln {R_2})(1 - \sigma ) - } \\ & \left.\sigma {\ln R({\rm{TC}}) ]{c_{\rm{c}}} + {c_{\rm{r}}} + {T_{{\rm{Down}}}}{c_{\rm{d}}}} \right\}/\\ & \left[ {\sigma \times {\rm{TC}} + (1 - \sigma )\sum\limits_{i = 1}^{{\tau _2}} {{T_i}} } \right]. \end{split} $

$\begin{array}{*{20}{l}} {{\rm{s.t.}}}\\ \;\;\;\;\;\;{{R_{\min }} \leqslant {R_2} \leqslant 1;}\\[-10pt] \end{array}\qquad\qquad\qquad\qquad\qquad\qquad\;\;\;$

$1 \leqslant {\tau _2} \leqslant {\tau _{\max }},\quad {\tau _2} \in {\bf{N}}.\qquad\qquad\qquad\qquad$

其余约束条件与式(22)相同.

式(34)为决策变量约束,可靠度阈值不能小于最小可靠度. 式(35)类似于式(24),为维修次数上限约束.

3.3. 求解算法

为了求解模型(P1)的最优检测周期T*检测上限次数τ*,考虑可靠度单调递减的特点,设计如图3所示的枚举优化算法. 在由式(20)计算目标函数值的过程中,需要由式(7)计算每次检测时预防性维修的概率Pmti),并由式(10)计算可靠度.

图 3

图 3   枚举优化算法流程图

Fig.3   Flow diagram of enumeration optimization algorithm


4. 案例分析

4.1. 案例背景

以某铁路专用线配属的SS4B机车中空气管路系统为研究对象,选取其中5个子系统进行独立的检修计划优化. 假设各系统的正常阶段和延迟时间阶段均服从威布尔分布,参数估计的结果如表1所示. 表中,m1m2分别为第1阶段和第2阶段分布的形状参数,l1l2分别为第1阶段和第2阶段分布的尺度参数.退化参数均通过参数检验. 各系统的检修成本和最低可靠度参数如表2所示. 设各系统的停机损失费用为300 元/h,最大可用时间均为730 d,可用度约束为0.98.

表 1   检修和退化参数估计结果

Tab.1  Parameters of inspection and maintenance and deterioration process

系统序号 检修参数 退化参数 参数验证
a r m1 l1−1 m2 l2 统计量实际值 统计量参照值 通过检验
1 0.05 0.68 1 0.003 5.347 6 126.344 0 18.364 18.548
2 0.02 1.00 1 0.011 1.857 1 124.111 0 16.465 16.748
3 0.05 0.78 1 0.020 5.678 7 112.360 3 20.936 21.954
4 0.05 0.57 1 0.006 3.063 3 180.541 0 15.475 16.748
5 0.05 0.87 1 0.011 5.639 2 99.360 3 18.325 18.548

新窗口打开| 下载CSV


表 2   系统检修成本和可靠度参数

Tab.2  Parameters of maintenance costs and reliability

系统序号 ci/元 cp/元 cr/元 cc/元 ti/h tp/h tr/h tc/h Rmin
1 100 280 1 800 4 000 1.5 3.0 6.0 20.0 0.94
2 80 400 1 600 3 500 0.5 1.0 2.5 25.0 0.94
3 50 150 840 2 750 0.5 1.5 4.5 8.0 0.93
4 100 320 2 150 7 200 2.0 4.0 10.0 15.0 0.92
5 120 500 1 670 9 000 1.5 3.5 6.5 10.0 0.94

新窗口打开| 下载CSV


4.2. 最优求解和策略比较

通过枚举优化算法求解,得到各系统的最优检测间隔和更换时间,如表3所示. 表中,Ga为求解结果与最优可用度的差异.如系统1最优检测间隔为41 d,在第11次检测时进行预防性更换,更新周期为451 d,平均费用率为24.27 元/d. 各系统的可用度均处于较高水平,并且与可实现的最大可用度差异很小.

表 3   检修计划优化模型求解结果

Tab.3  Optimal solution of inspection model

系统 Tmax/d $\bar C$/(元·d−1 T/d τ L/d A Ga/%
1 134 24.27 41 11 451 0.997 66 0.02
2 66 18.84 24 30 720 0.998 61 0
3 93 14.55 27 27 729 0.998 56 0
4 144 40.44 42 7 294 0.996 12 0.06
5 88 34.38 30 23 690 0.996 65 0.01

新窗口打开| 下载CSV


以系统1为例,分析目标函数与检测间隔和次数的关系,如图4所示. 在相同的检测间隔下,目标值随着检测次数的增大而减小. 各检测间隔对应的最优值(最优解曲线)随着检测间隔的增加先增后减,并在41天时达到最优目标值. 该曲线在检测间隔增大过程中存在波动,且最优检测次数下降. 因为当检测间隔较大时,系统的可靠度下降速度较快,最大检测次数受到最低可靠度的限制. 每个检测间隔下受最低可靠度和最大可用时间约束的最大维修次数分别如图4的点线和虚线所示. 当检测间隔超过32 d时,可靠度下降较快,最优解曲线受到最低可靠度的限制. 当维修间隔较小时,可靠度下降缓慢,较多次数的检测能够保证较高的可靠度,此时最优解曲线受到最大使用时间的约束. 此外,可靠度约束曲线表明,只有当检测间隔和次数均较小时,可用度较低. 在考虑停机损失费用时,可用度与费用的优化方向接近. 在以费用率为优化目标时,全局最优解均满足可用度约束,这与表3的结果一致.

图 4

图 4   枚举优化算法求解过程分析

Fig.4   Analysis of solution procedure


不同维修策略的最优结果如表4所示. 表中,Gc 为求解结果与检修策略费用的差异. 系统1和4的最优策略为可靠度阈值策略,其余系统的最优策略均为检修策略. 由系统参数可知,系统1和4的检出概率较低且检测成本较高,所以会产生大量无效检测,导致系统成本增加. 在2个非检测策略中,固定可靠度阈值策略优于固定周期策略. 因为最低可靠度约束使得固定可靠度阈值策略平均增加一个检测间隔的可用时间,从而延长系统的更新周期.

表 4   不同维修策略求解结果与比较

Tab.4  Comparisons of different maintenance strategies

系统 检修策略(P1) 对比策略1:固定周期维修(P2) 对比策略2:固定可靠度阈值维修(P3)
$\bar C$/(元·d−1 T/d L/d ${\bar C_1}$/(元·d−1 Gc/% T1/d L1/d ${\bar C_2}$/(元·d−1 Gc/% R2 L2/d
1 24.27 41 451 21.81 −10.14 90 450 19.73 −18.71 0.990 512
2 18.84 24 720 39.89 +111.67 42 126 33.95 +80.14 0.984 152
3 14.55 27 729 16.81 +15.59 61 366 15.04 +3.39 0.988 370
4 40.44 42 294 32.24 −20.28 92 368 28.80 −28.77 0.986 477
5 34.38 30 690 37.58 +9.31 65 260 34.39 +0.01 0.984 263

新窗口打开| 下载CSV


4.3. 检出概率的灵敏度分析

检出概率是检修模型的重要参数. 由4.2节的策略比较可知,当检出概率较低时,检修策略不是最优策略,因此需要分析检出概率对检修计划的影响. 如图5所示为不同检出概率r下模型最优解和目标值的变化趋势. 随着检出概率的增加,目标值均呈现先稳定后下降的特点;相应地,检测间隔首先不变,之后迅速减小,最后波动增加;更新周期先稳定后逐渐增加,最后在小范围内波动. 最优检测间隔和更新周期的波动变化是因为维修间隔和检测次数均为离散变量,非线性变化的可靠度和线性变化的更新周期均需要满足约束条件. 随着检出概率的增加,会出现检测间隔减小且检测次数增大、检测间隔增大且检测次数减小的情况,分别导致更新周期的增大和减小.

图 5

图 5   检出概率对目标函数和维修计划的影响

Fig.5   Influence of inspection accuracy on objective values and maintenance schedules


图5可知,各系统均存在一个检出概率的变化阈值,如系统1的阈值为0.47. 当检出概率小于该值时,进行检测反而增加成本. 当检出概率超过该值时,检测可以减缓可靠度下降速度,从而可以在较长更新周期下满足最低可靠度约束,所以系统通过频繁检测而延长更新周期. 随着检出概率的进一步增大,检测间隔逐渐增大,因为检测更容易识别系统运行中的缺陷,同时更新周期开始受到最大使用时间约束. 综上所述,检出概率要大于相应阈值时才可以执行检修策略. 检出概率越大,不仅能够降低系统的维修成本,而且能够延长系统更新周期,避免系统的频繁更换.

以系统1、2为例,分析检出概率对不同策略的影响,如图6所示. 2个非检测策略的结果不受检出概率变化的影响. 当检出概率分别超过0.83、0.65之后,检修策略成为最优策略,说明只有当检出概率保持在较高水平时,通过检修确定维修计划才具有应用价值. 基于表1的检出概率,对比理想情况(r=1)下的优化结果,如图7所示. 可知,当检出概率为1时各系统的费用率均为全部策略的最优结果,相比可靠度阈值策略(P3)平均可以降低费用率36.78%.

图 6

图 6   检出概率对不同维修策略的影响

Fig.6   Influence of inspection accuracy on different maintenance strategies


图 7

图 7   理想情况下检测模型的目标值及其比较结果

Fig.7   Comparisons of objective values among different maintenance strategies


4.4. 检修概率和故障概率的变化

随着检测次数的增加,检测费用的增长是线性的,累计预防性维修概率、故障期望次数的变化是非线性的. 不同检出概率对执行预防性维修和出现故障有较大的影响. 以系统1为例,设检测周期为20 d,不同检出概率下预防性维修概率、故障期望次数Ef随检测次数的变化如图8所示.

图 8

图 8   预防性维修和故障期望随着检测次数的变化

Fig.8   Variation of probability of preventive maintenance and expected number of failures


随着检测次数的增加,预防性维修的发生概率(包括检出缺陷概率和上一间隔的故障概率)先增后减,当检出概率为1时不存在增加的过程. 故障期望次数基本呈递增趋势,当检出概率较小(r<0.4)时存在一个先减后增的过程,该过程是不完全检测产生的累计效应造成的. 当检出概率较低时,未被检出的缺陷会累计到检测之后的间隔并引发故障,故障期望呈现递增趋势. 相对密集的故障增大了使检修后系统状态处于较新状态的概率,之后发生缺陷和故障的概率均会下降. 预防性维修的发生概率和故障次数期望之差可以近似看成检出缺陷的概率,呈现先增后减的趋势,且下降速度逐渐增大,因为不完全维修的累计效应使得缺陷更容易以故障的形式出现. 在同一检出概率下,不完全维修使故障费用的比例逐渐增加,预防性维修比例逐渐减小,为了实现最优费用率,系统必然在某次检测时进行更换.

随着检出概率的增大,同一次检测对应的预防性维修概率逐渐增大,故障次数期望逐渐减少,不完全检测的累计效应逐渐不明显. 增大检出概率能够降低目标函数值,较高的检出概率使得缺陷更容易在检测时被发现,减少了故障的发生.

5. 结 语

本文基于延迟时间模型,建立同时考虑不完全检测和不完全维修的可靠度模型和检修计划优化模型,应用实际数据进行参数标定和案例分析. 结果表明:利用可靠度模型能够有效地描述系统退化和故障情况,且检测间隔对于可靠度变化有明显的影响,较短的间隔能够有效地降低可靠度下降速度. 利用检修计划优化模型可以有效地确定最优检测间隔和更换策略,降低平均维修成本. 与既有研究中基于固定周期和固定可靠度阈值的维修模型相比,只有当检出概率处于较高水平时,检修模型能够有效地降低维修成本,检测无误差时可以平均节省成本36.78%. 本文研究为单一系统维修计划的制定与优化提供了有效方法,今后进一步研究多部件系统的协同检修计划,有助于提高检修计划模型的适用性.

参考文献

CHEN S, HO T, MAO B, et al

A bi-objective maintenance scheduling for power feeding substations in electrified railways

[J]. Transportation Research Part C: Emerging Technologies, 2014, 44: 350- 362

DOI:10.1016/j.trc.2014.04.016      [本文引用: 2]

ALASWAD S, XIANG Y

A review on condition-based maintenance optimization models for stochastically deteriorating system

[J]. Reliability Engineering and System Safety, 2017, 157: 54- 63

DOI:10.1016/j.ress.2016.08.009      [本文引用: 1]

吕文元, 郑睿

基于时间延迟的维修类型优化组合模型及案例分析

[J]. 系统工程理论与实践, 2013, 33 (7): 1654- 1660

DOI:10.3969/j.issn.1000-6788.2013.07.004      [本文引用: 2]

LV Wen-yuan, ZHENG Rui

Delay-time-based model of optimization grouping maintenance type and case study

[J]. System Engineering: Theory and Practice, 2013, 33 (7): 1654- 1660

DOI:10.3969/j.issn.1000-6788.2013.07.004      [本文引用: 2]

周炳海, 陶红玉, 綦法群

带随机突变的两阶段退化系统视情维修建模

[J]. 哈尔滨工业大学学报, 2016, 48 (1): 30- 35

[本文引用: 1]

ZHOU Bing-hai, TAO Hong-yu, QI Fa-qun

Condition-based maintenance modelling for two-stage deteriorating systems with random change

[J]. Journal of Harbin Institute of Technology, 2016, 48 (1): 30- 35

[本文引用: 1]

郭驰名, 郭波, 王文彬, 等

非周期不完全检测下的维修优化

[J]. 国防科技大学学报, 2013, 35 (4): 176- 181

DOI:10.3969/j.issn.1001-2486.2013.04.031      [本文引用: 1]

GUO Chi-ming, GUO Bo, WANG Wen-bin, et al

Maintenance optimization under non-periodic imperfect inspections

[J]. Journal of National University of Defense Technology, 2013, 35 (4): 176- 181

DOI:10.3969/j.issn.1001-2486.2013.04.031      [本文引用: 1]

赵斐, 刘学娟

考虑不完美维修的定期检测与备件策略联合优化

[J]. 系统工程理论与实践, 2017, 37 (12): 3201- 3214

DOI:10.12011/1000-6788(2017)12-3201-14      [本文引用: 3]

ZHAO Fei, LIU Xue-juan

Joint optimization of periodical inspection and spare parts policies considering imperfect maintenance

[J]. Systems Engineering: Theory and Practice, 2017, 37 (12): 3201- 3214

DOI:10.12011/1000-6788(2017)12-3201-14      [本文引用: 3]

CHRISTER A H, WALLER W M

Delay time models of industrial inspection maintenance problems

[J]. Journal of the Operational Research Society, 1984, 35 (5): 401- 406

DOI:10.1057/jors.1984.80      [本文引用: 1]

LI Y, MA X, ZHAI Q, et al

A delay time model for a mission-based system subject to periodic and random inspection and postponed replacement

[J]. Reliability Engineering and System Safety, 2016, 150: 96- 104

DOI:10.1016/j.ress.2016.01.016      [本文引用: 1]

ARTS J, BASTEN R

Design of multi-component periodic maintenance programs with single-component models

[J]. IISE Transactions, 2018, 50 (7): 606- 615

DOI:10.1080/24725854.2018.1437301      [本文引用: 2]

刘学娟, 赵斐. 基于延迟时间理论的n中取k系统检测区间模型[J/OL]. 控制与决策, 2020, 35(6): 1469-1475 [2019-03-07]. https://doi.org/10.13195/j.kzyjc.2018.1395.

[本文引用: 1]

LIU Xue-juan, ZHAO Fei. Delay-time-based inspection model for k-out-of-n systems [J/OL]. Control and Decision, 2020, 35(6): 1469-1475 [2019-03-07]. https://doi.org/10.13195/j.kzyjc.2018.1395.

[本文引用: 1]

STAPELBERG R F. Handbook of reliability, availability, maintainability and safety in engineering design [M]. London: Springer, 2009.

[本文引用: 1]

WANG W

An overview of the recent advances in delay-time-based maintenance modeling

[J]. Reliability Engineering and System Safety, 2012, 106: 165- 178

DOI:10.1016/j.ress.2012.04.004      [本文引用: 1]

胡海军, 程光旭, 段权, 等

一种包含非理想维修的延迟时间模型

[J]. 西安交通大学学报, 2009, 43 (6): 103- 107

DOI:10.3321/j.issn:0253-987X.2009.06.022      [本文引用: 1]

HU Hai-jun, CHENG Guang-xu, DUAN Quan, et al

Delay time model based on imperfect maintenance

[J]. Journal of Xi’an Jiaotong University, 2009, 43 (6): 103- 107

DOI:10.3321/j.issn:0253-987X.2009.06.022      [本文引用: 1]

WANG L, HU H, WANG Y, et al

The availability model and parameters estimation method for the delay time model with imperfect maintenance at inspection

[J]. Applied Mathematical Modelling, 2011, 35 (6): 2855- 2863

DOI:10.1016/j.apm.2010.11.070      [本文引用: 2]

CHRISTER A H, WANG W, BAKER R D, et al

Modelling maintenance practice of production plant using the delay-time concept

[J]. IMA Journal of Mathematics Applied in Business and Industry, 1995, 6 (1): 67- 83

[本文引用: 1]

WANG H, WANG W, PENG R

A two-phase inspection model for a single component system with three-stage degradation

[J]. Reliability Engineering and System Safety, 2017, 158: 31- 40

DOI:10.1016/j.ress.2016.10.005      [本文引用: 1]

DRIESSEN J P C, PENG H, VAN HOUTUM G J

Maintenance optimization under non-constant probabilities of imperfect inspections

[J]. Reliability Engineering and System Safety, 2017, 165: 115- 123

DOI:10.1016/j.ress.2017.03.020      [本文引用: 1]

刘勤明, 李永朋, 叶春明. 基于三阶段时间延迟模型的设备预防维修策略研究[J/OL]. 控制与决策, 2020, 35(7): 1780-1786 [2019-03-08]. https://doi.org/10.13195/j.kzyjc.2018.1692.

[本文引用: 1]

LIU Qin-ming, LI Yong-peng, YE Chun-ming. Research on preventive plan of equipment based on three-stage time delay model [J/OL]. Control and Decision, 2020, 35(7): 1780-1786 [2019-03-08]. https://doi.org/10.13195/j.kzyjc.2018.1692.

[本文引用: 1]

YANG L, YE Z S, LEE C G, et al

A two-phase preventive maintenance policy considering imperfect repair and postponed replacement

[J]. European Journal of Operational Research, 2019, 274 (3): 966- 977

DOI:10.1016/j.ejor.2018.10.049      [本文引用: 2]

MALIK M A K

Reliable preventive maintenance scheduling

[J]. AIIE Transactions, 1979, 11 (3): 221- 228

DOI:10.1080/05695557908974463      [本文引用: 2]

PARGAR F, KAUPPILA O, KUJALA J

Integrated scheduling of preventive maintenance and renewal projects for multi-unit systems with grouping and balancing

[J]. Computers and Industrial Engineering, 2017, 110: 43- 58

DOI:10.1016/j.cie.2017.05.024      [本文引用: 1]

BOUVARD K, ARTUS S, BERENGUER C, et al

Condition-based dynamic maintenance operations planning and grouping: application to commercial heavy vehicles

[J]. Reliability Engineering and System Safety, 2011, 96 (6): 601- 610

DOI:10.1016/j.ress.2010.11.009      [本文引用: 2]

杨春节, 童晟, 孙长生, 等

基于可靠度约束的混合预防性维修模型

[J]. 浙江大学学报:工学版, 2008, 42 (8): 1376- 1379

[本文引用: 1]

YANG Chun-jie, TONG Sheng, SUN Chang-sheng, et al.

Hybrid preventive maintenance model based on reliability constraint

[J]. Journal of Zhejiang University: Engineering Science, 2008, 42 (8): 1376- 1379

[本文引用: 1]

LIU G, LONG X, TONG S, et al

Optimum consecutive preventive maintenance scheduling model considering reliability

[J]. Journal of Shanghai JiaoTong University: Science, 2019, 24 (4): 490- 495

DOI:10.1007/s12204-019-2089-z      [本文引用: 1]

/