现代复杂系统通常装备有状态监测单元如传感器、仪表、计算设备[1 ] ,这为基于状态的维修(CBM)实施提供了系统状态数据.相对于基于时间的维修方式,CBM能够有效控制“维修不足”和“维修过度”的问题[2 ] .通常来说,获取系统状态信息的方式可分为定期检测和连续监测两类,相比于定期检测,连续监测能够实时获取系统状态信息,从而更为及时有效地进行CBM决策,降低系统故障风险和维修费用.随着传感器技术的发展,连续监测的应用场景越来越多,特别是对于一些关键部件或系统,如发电机组,其本身价值高、停机损失大,对系统状态进行连续监测更为适用.同时,复杂系统通常存在多个退化过程,且由于共同的运行环境、应力等因素,退化过程间往往是相关的,在CBM决策中需要充分考虑这些因素[3 ] ,这也增加了多元退化系统CBM决策的难度.
对于有多个相关退化过程的系统而言,CBM决策的基础是对不同退化过程间的相关性进行建模描述.目前,常用的方法包括多维退化模型[4 ,5 ] 、退化率相关模型[6 ,7 ] 和基于Copula函数的退化相关模型[8 ,9 ] .文献[4 ,5 ]分别采用多元正态分布和二元Birnbaum-Saunders(B-S)分布来描述不同退化过程间的相关关系,而分布形式也限制了此类方法的应用.文献[6 ,7 ]将不同退化过程间的相关关系表示为一个部件退化水平对其余部件退化速率的影响.在工程实际中,这种影响关系及程度很难采用定量化方法描述和确定.相对而言,Copula函数在定量描述随机变量相关性方面具有很高的灵活性,并且提供了分析相关结构对系统可靠性影响的参数化方法.基于不同种类的Copula相关结构,文献[10 ]利用仿真方法对定期检测策略下的系统维修费用进行分析,其中预防性维修阈值为固定值.文献[11 ]建立了多部件系统检测维修策略优化模型,其中部件故障间的相关性采用Gaussian Copula函数来描述.文献[12 ]采用Levy Copula函数描述系统二元退化间的相关关系,并在此基础上提出系统定期检测与机会维修策略,利用仿真方法获得了最优预防性维修阈值和机会维修概率阈值.在现有的基于Copula函数建立多元退化模型的文献中,很少研究在连续监测条件下的CBM决策优化问题,且没有考虑系统备件约束,而对一些价值较高的重要系统而言,很难保证时刻有可用备件,通常需要提前预定.在连续监测条件下,文献[13 ]研究了单一退化过程的预防性更换阈值优化问题.在此基础上,文献[14 ]考虑二元相关Levy退化过程,以系统长期运行条件下的期望维修费用率和可用度为准则,计算获得在不同系统结构下的最优预防性更换阈值.上述研究均是在维修延迟(维修资源的到达需要一定时间,此处的维修资源也可指备件)条件下,对预防性更换阈值进行决策,没有综合考虑备件订购阈值以及预防性维修阈值的优化问题.
针对上述问题,本文考虑系统不同退化过程间具有的相关性,利用Copula函数建立系统多元退化模型.在此基础上,提出系统维修及备件订购策略,分析直接维修费用、备件库存费用和故障停机损失,建立系统维修费用模型.以系统期望维修费用率准则,引入人工蜂群(ABC)算法搜索最优维修及备件订购策略,从而提高维修工作的经济性,为科学制定系统维修及备件订购计划提供方法支撑.
1 系统描述
系统性能会随着工作时间的增加不断退化,当退化超过阈值时会导致系统出现故障.假设系统存在n 个相关性能退化过程,任意退化过程超过特定阈值时系统即会出现故障.
现有的退化过程模型主要包括退化轨迹模型、退化量分布模型和基于随机过程的退化模型等.相对于前两种模型,基于随机过程的退化模型能够描述性能退化在时间轴上的不确定性,更符合工程实际.在基于随机过程的退化模型中,Gamma过程由于其良好的性质,广泛应用于连续非减的退化过程建模,包括材料的磨损、侵蚀、裂纹、老化等.本文采用Gamma过程来建立系统退化过程模型.
令Xi (t ),i =1,2,…,n 表示系统第i 个性能退化过程,根据Gamma过程平稳独立增量性质可知,对于任意t >u >0,Xi (t )-Xi (u )=Xi (t -u )~ Γ(αi (t -u ),βi ),其中:Γ(·,·)为Gamma分布;αi >0和βi >0分别为形状参数和尺度参数.假设系统初始退化量均为0,那么Xi (t )的概率密度函数和概率密度分布函数可表示为
(1) f X i ( t ) ( x | α i , β i ) = 1 Γ ( α i t ) β i α i t x α i t - 1 e - x β i I x ≥ 0 F X i ( t ) ( x | α i , β i ) = ∫ 0 x f X i ( t ) ( s | α i , β i ) d s = 1 - Γ α i t , x β i Γ ( α i t )
式中:Γ(α )=∫ 0 ∞ rα -1 e- r dr 为Gamma函数;Γ(α ,x )=∫ x ∞ rα -1 e- r dr 为不完全Gamma函数;I { x ≥0} 为示性函数.系统第i 个退化过程的均值为αi βi t ,方差为αi β i 2 t.
如上所述,本文采用Copula函数描述退化过程的相关性.为便于实际应用,这里假设:在不同时间间隔内,不同退化过程的退化过程增量间的相关性可忽略.在文献[8 ]中也采用了类似的假设分析多元Wiener退化过程的相关性.令t >u >0,根据Sklar定理,系统n 个退化过程在时间区间[u ,t ]上退化增量Xi (t -u )的联合概率分布函数可表示为
(2) Ht - u (x )=C (F X 1 ( t - u ) ( x 1 ),F X 2 ( t - u ) ( x 2 ),…,F X n ( t - u ) ( xn );θ )
(3) ht - u (x )=∂ H t - u ( x ) ∂ x = c (F X 1 ( t - u ) ( x 1 ),F X 2 ( t - u ) ( x 2 ),…,F X n ( t - u ) ( xn );θ )×∏ i = 1 n f X n ( t - u ) ( xn )
(4) c (F X 1 ( t - u ) ( x 1 ),F X 2 ( t - u ) ( x 2 ),…,F X n ( t - u ) ( xn );θ )=∂ C ( F X 1 ( t - u ) ( x 1 ) , F X 2 ( t - u ) ( x 2 ) , … , F X n ( t - u ) ( x n ) ; θ ) ∂ F X 1 ( t - u ) ( x 1 ) ∂ F X 2 ( t - u ) ( x 2 ) … ∂ F X n ( t - u ) ( x n )
式中:C 为Copula函数;c 为Copula密度函数;x =[x 1 x 2 … xn ],xi ≥0;F X i ( t - u ) ( xi )为第i 个退化过程在[u ,t ]内增量的分布函数;θ =[θ 1 θ 2 … θm ]为Copula函数中的参数向量,其取值直接影响着随机变量间相关关系强弱.对于Copula函数参数估计问题,可基于退化数据通过最大似然估计或边际推断(IFM)法分两步获得.在此基础上,可采用赤池信息准则(AIC)选择合适的Copula函数描述系统退化相关性.关于Copula函数的参数估计方法参考文献[15 ].
在此基础上,假设系统退化故障阈值向量为 Q L =[ Q L 1 Q L 2 …Q L n ],可获得系统可靠度函数为
(5) R sys (t )=P (X 1 (t )≤Q L 1 ,X 2 (t )≤Q L 2 ,…,Xn (t )≤Q L n ) = Ht (Q L )=C (F X 1 ( t ) ( Q L 1 ) ,F X 2 ( t ) ( Q L 2 ) ,…,F X n ( t ) ( Q L n ) ;θ )
由于Gamma函数和不完全Gamma函数的影响,F X i ( t ) ( Q L i ) 的计算较为复杂,可采用B-S分布对其进行近似计算[16 ] .F X i ( t ) ( Q L i ) 可近似表示为
(6) F X i ( t ) ( Q L i ) ≈1-F BS (t ;Q L i ) = Φ 1 γ i λ i t - t λ i
式中:F BS (t ;Q L i ) 为B-S分布函数;γi =β i / Q L i ;λi =Q L i /αi βi ;Φ (·)为标准正态分布.
令Ui (t )=1 γ i λ i t - t λ i ,由此,系统可靠度函数可近似表示为
(7) R sys (t )=C F X 1 ( t ) ( Q L 1 ) , F X 2 ( t ) ( Q L 2 ) , … , F X n ( t ) ( Q L n ) ; θ ≈C Φ ( U 1 ( t ) ) , Φ ( U 2 ( t ) ) , … , Φ ( U n ( t ) ) ; θ
2 系统维修与备件订购策略
本文的系统维修和备件订购采用控制限策略(CLP),即根据系统退化状态确定相应的备件订购和维修活动.备件订购阈值为Q A =[Q A 1 Q A 2 …Q A n ],预防性更换阈值为Q M =[Q M 1 Q M 2 …Q M n ].为建立系统维修与备件订购决策优化模型,对相关策略做出如下假设.
(1) 对系统在运行过程中的退化状态进行连续监测,系统各性能退化指标水平可以通过连续监测得到;单位时间状态监测费用记为C I .
(2) 初始备件库存为0,当系统任意退化过程达到备件订购阈值时,系统发出订货信号,备件经过时间τ 后到达.此时,若系统不进行更换,会产生相应的备件库存费用;备件订购费用记为C O ,单位时间库存费用记为C S .
(3) 当系统任意退化过程达到预防性更换阈值后,系统发出维修信号,如果此时库存中有备件,直接对系统进行更换;否则,在备件到达后对系统进行更换.更换后,系统修复如新,且之后系统退化过程与更换前的系统无关.如果备件是在系统故障后到达,即会产生停机损失,系统单位时间内的停机损失记为C LO .一般来说,系统维修更换费用与其退化量相关,退化量高的系统通常所需维修费用也较多.系统退化量是由多个退化过程共同决定的,由于各退化过程的故障阈值和量纲可能会不同,进而将影响系统退化量的计算.因此,首先定义各退化过程的相对退化量为
(8) χi (t )=X i ( t ) Q L i , i =1,2,…,n
由式(8)可以看出,当任意χi (t )≥1时,系统处于故障状态.多退化过程引发的系统故障为竞争失效过程,系统状态由多退化过程中的最大相对退化量决定.基于此,系统在t 时刻的维修更换费用可表示为
(9) CMAN =CR +ρEsup (χt )
式中:E sup (χt )=max{E (χ 1 (t )),E (χ 2 (t )),…,E (χn (t ))},E (χi (t ))为χi (t )的期望;C R 为系统更换固定费用;ρ 为更换费用随系统退化量变化系数.
(4) 维修时间相对于系统寿命周期来说非常短,在建模过程中忽略不计.
(5) 为使备件订购阈值和系统预防性更换阈值有意义,令Q A i ≤Q M i ≤Q L i ,i =1,2,…,n.
以单退化过程为例,系统退化过程和更新过程如图1 所示.其中:t A 、t M 和t L 分别为退化过程到达备件订购阈值、预防性更换阈值及故障阈值的时间,且t L ≥t M ≥t A >0.
图1
图1
系统退化和更新过程
Fig.1
System degradation and renew process
对于Gamma过程而言,其退化路径为跳跃过程,因此有P (t M -t A =0)>0,P (t L -t A =0)>0,如图1 中系统第2个更新周期内的退化过程所示.事实上,对于很多左极右连的Levy过程(包含Gamma过程)来说,均存在上述性质.
系统不同维修时机如图2 所示.根据备件到达的时间不同,系统维修情况可划分为以下3种:
(1) 当t A +τ <t M 时,在系统退化状态到达预防性更换阈值时对系统进行更换;
(2) 当t M ≤t A +τ <t L 时,在备件到达时对系统进行预防性更换;
(3) 当t A +τ ≥t L 时,在备件到达时对系统进行故障更换.
图2
图2
系统的不同维修时机
Fig.2
Different maintenance opportunities of the system
3 费用模型
本文以系统长期运行情况下单位时间平均费用(即期望费用率)为指标,对系统预防性维修和备件订购策略进行评价.由前文分析可以看出,备件订购阈值和预防性更换阈值的高低会影响不同维修方式发生概率、维修费用、库存费用、更新周期长度,从而影响系统期望维修费用率.
3.1 期望维修费用率计算
考虑上述3种不同的系统维修情况,在不同情况下系统在一个更新周期内的维修费用C (T )及其相应的更新周期长度T 可表示为
(10) C ( T ) = C S ( t M - t A - τ ) + C I ( t M ) + C R + ρ E sup ( χ t M ) + C O T = t M
(11) C ( T ) = C R + ρ E sup ( χ t A + τ ) + C I ( t A + τ ) + C O T = t A + τ
(12) C ( T ) = C LO ( t A + τ - t L ) + C I ( t A + τ ) + C R + ρ E sup ( χ t A + τ ) + C O T = t A + τ
在此基础上,根据更新报酬原理可得系统的期望维修费用率为
(13) γ ∞ =[C S E (max{t M -t A ,τ })-C LO E (min{t L -t A ,τ })+(C LO -C S )τ +ρ (E sup (χ t A + τ )- E sup (χ t M )) P (t M -t A ≤τ )+ρE sup (χ t M ) +C R +C O ]/[E (max{t M -t A ,τ })+E (tA )]+C I
3.2 具体表达式推导
为获得系统期望维修费用率的解析式,假设G M (u )和G L (u )分别为t M -t A 和t L -t A 的分布函数,相应的生存函数G - M ( u )=1-G M (u ),G - L ( u )=1-G L (u ),可得:
(14) E (max{t M -t A ,τ })=τG M (τ )+∫ τ + ∞ s dG M (s )=τG M (τ )+∫ τ + ∞ ∫ 0 s du dG M (s )=τG M (τ )+ ∫ 0 τ ∫ τ + ∞ dG M (s )du +∫ τ + ∞ ∫ u + ∞ dG M (s )du =τG M (τ )+τ -τG M (τ )+∫ τ + ∞ (1- G M (u ))du =τ +∫ τ + ∞ G - M ( u )du
(15) E (min{t L -t A ,τ })=∫ 0 τ sG L (s )ds +τ G - L ( τ )=G L (τ )τ -∫ 0 τ G L (s )ds +τ G - L ( τ )=∫ 0 τ G - L ( s )ds
(16) P (t M -t A ≤τ )=1-P (t M -t A >τ )=1-G - M ( τ )
(17) P (t A >t )=P (X 1 (t )<Q A 1 ,X 2 (t )<Q A 2 ,…,Xn (t )<Q A n ) =Ht (Q A )
式中:Ht (Q A )可由式(2)获得.在此基础上可得期望E (t A )的表达式为
(18) E (t A )=∫ 0 + ∞ P (t A >u )du =∫ 0 + ∞ Hu (Q A )du
在上述推导结果基础上,当τ为常数时,可得系统期望维修费用率表达式为
(19) γ∞ =[CS ∫ τ + ∞ G - M ( s)d s-CL ∫ 0 τ G - L ( s)d s+CL τ+ρ(Esup (χ t A + τ ) -Esup (χ t M )) [1-G - M ( τ)]+ρEsup (χ t M ) +CR +CO ]/[τ+∫ τ + ∞ G - M ( s)d s+∫ 0 + ∞ Hu (Q A )du ]+C I
通过式(19)可以看出,计算系统期望维修费用率的关键在于推导G - M ( s ),∫ τ + ∞ G - M ( s )ds ,∫ 0 τ G - L ( s )ds ,E sup (χ t M ) 和E sup (χ t A + τ ) 的表达式.
根据定义,G - M ( s )和G - L ( s )分别为t M -t A 和t L -t A 的生存函数,即
(20) G - M (s) = P ( t M - t A > s ) G - L (s) = P ( t L - t A > s )
基于前文中的系统相关退化过程模型,可以通过推导获得G - M ( s ),∫ τ + ∞ G - M ( s )ds 和∫ 0 τ G - L ( s )ds 的表达式为
(21) G - M (s) = ∫ 0 + ∞ ∫ x ∈ [ 0 , Q A ] ∂ H s ( Q M - x ) ∂ s h v ( x ) d x d v + H s ( Q M ) ∫ τ + ∞ G - M (s) d s = ∫ τ + ∞ H s ( Q M ) d s - ∫ 0 + ∞ ∫ x ∈ [ 0 , Q A ] H τ ( Q M - x ) h v (x) d x d v ∫ 0 τ G - L (s) d s = ∫ 0 + ∞ ∫ x ∈ [ 0 , Q A ] H τ ( Q L - x ) h v (x) d x d v - ∫ 0 + ∞ H v ( Q A ) d v + ∫ 0 τ H s ( Q L ) d s
式中:Q L -x =[Q L 1 - x 1 Q L 2 - x 2 …Q L n - xn ],Q M -x 也是类似的表示;∫ x ∈ [ 0 , Q A ] ( ·)dx 为关于[x 1 x 2 … xn ]的n 重积分;Hs (·)和hv (·)可分别由式(2)和(3)获得.
对于E sup (χ t M ) 和E sup (χ t A + τ ) ,由定义可知:
E sup (χ t M ) =max{E ( X 1 ( t M ) ) Q L 1 ,E ( X 2 ( t M ) ) Q L 2 ,…,E ( X n ( t M ) ) Q L n },系统各性能退化均为Gamma过程,在任意时刻退化量联合分布的边缘分布为Gamma分布,根据鞅停时定理[13 ] ,容易得到:
(22) E (Xi (t M ))=αi βi E (t M )
式中:E (t M )与式(18)有类似表达式,此处不再赘述.在此基础上,可得:
(23) E sup (χ t M ) =maxα 1 β 1 E ( t M ) Q L 1 , α 2 β 2 E ( t M ) Q L 2 , … , α n β n E ( t M ) Q L n
(24) E sup (χ t A + τ ) =Max{α 1 β 1 ( E ( t A ) + τ ) Q L 1 ,α 2 β 2 ( E ( t A ) + τ ) Q L 2 ,…,α n β n ( E ( t A ) + τ ) Q L n }
将式(21)、(23)和(24)代入式(19)即可得到系统期望维修费用率的解析表达式.
3.3 期望维修费用率估算
通过前文获得的系统期望维修费用率模型,可以得到相关的数值计算结果.但由于模型中涉及到x =[x 1 x 2 … xn ]的高维积分,虽然在MATLAB等科学计算软件中集成了高维积分的数值算法,可以直接使用,但算法的计算复杂度较高.为简化计算,下面给出模型的近似表达式.
假设系统各退化过程在订货时间点t A 的退化量已知,记为
X t A =[X 1 (t A ) X 2 (t A ) … Xn (t A )],那么容易得到:
(25) G - M ( s )=P (t M -t A >s )=P (t Q M - X t A >s )= P (Xs <Q M -X t A )
由于退化过程均为随机过程,所以在实际中X t A 为随机变量.因此,这里用X t A 的期望值E (X t A ) =[E (X 1 (t A )) E (X 2 (t A )) … E (Xn (t A ))]代替X t A . 由此可得:
(26) G - M ( s )≈P (Xs <Q M -E (X t A )) = Hs (Q M -E (X t A ))
(27) G - L ( s )≈Hs (Q L -E (X t A ))
因此,∫ τ + ∞ G - M ( s )ds 和∫ 0 τ G - L ( s )ds的表达式可近似为
(28) ∫ τ + ∞ G - M (s) d s ≈ ∫ τ + ∞ H s ( Q M - E ( X t A ) ) d s ∫ 0 τ G - L (s) d s ≈ ∫ 0 τ H s ( Q L - E ( X t A ) ) d s
通过这种近似方法,可以避免x =[x 1 x 2 … xn ]的高维积分运算.将式(26)和(28)代入式(19),即可得到系统期望维修费用率的近似表达式,记为γ ~ ∞ . 可以看出,γ ~ ∞ 和γ ∞ 之间的差值主要是由E (X t A ) 和X t A 引起的,所以当X t A 的方差较小时,估计结果也相对较为精确;反之,误差会相应增大.
4 优化模型及算法
在上述费用模型基础上,为制定最优系统维修与备件订购策略,以系统长期运行期望维修费用率最小为优化目标,对预防性维修阈值和备件订购阈值进行优化.基于此,建立如下优化模型
(29) (Q A * ,Q M * ) =argmin γ ∞ (Q A ,Q M )
在上述决策模型中,可以看出目标函数TC∞ 具有非线性、不可微的特点,且模型决策变量较多,当系统存在n个退化过程时,模型有2n个决策变量,因此,难以得到模型的解析解.
为解决这一问题,采用ABC 算法搜索系统最优维修与备件订购策略.ABC 算法是一种基于蜂群搜索蜜源行为的群体智能优化算法,在智能优化算法中属于比较新的算法,该算法控制参数少、易于实现、全局收敛性能好,并且其优化性能优于其他一些传统智能优化算法[17 ] ,如遗传算法(GA )、粒子群优化(PSO )算法、差分进化(DE )算法等,特别是在非线性优化函数求解方面具有良好的性能.ABC 算法已广泛应用于解决各类优化问题,关于ABC 算法的具体机理和优化过程可参考文献[18 ],这里不再赘述.在案例分析中,ABC 算法采用MATLAB 实现.
5 案例分析
为验证上述决策优化模型的有效性,以两部件系统为例,即n=2,对系统维修及备件订购策略进行优化,并分析相关参数对模型优化结果的影响.考虑两种不同的情形,在案例1中假设系统两个退化过程参数相同,单位时间退化量均值为2,方差为4,故障阈值也相同;在案例2中假设系统两个退化过程参数存在差异,单位时间内退化量均值分别为2和1.5,方差分别为4和1,相应的故障阈值也不同.不同案例的退化过程及故障阈值参数如表1 所示.
对于模型中的费用参数,为方便说明,假定单位时间系统状态监测费用为1,其余各项费用值均为与状态监测费用的比值.系统维修费用参数值如表2 所示.
在Copula函数选择方面,以Gaussian Copula函数为例,描述系统两个退化过程间的相关关系,二元Gaussian Copula函数的分布函数为
(30) C (u ,v ;θ )=∫ - ∞ Φ - 1 ( u ) ∫ - ∞ Φ - 1 ( v ) 1 2 π 1 - θ 2 exp2 θxy - x 2 - y 2 2 ( 1 - θ 2 ) dx dy
式中:Φ -1 (·)为标准正态分布逆函数.当表示两个变量间为正相关关系时,θ ∈[0,1],在本案例中令 θ =0.7.当然,这里也可以采用其他类型的Copula函数描述退化间相关关系,其分析过程相同.
5.1 系统维修及备件订购策略优化
在以上参数设置基础上,假设备件订购交付时间为τ =1.为验证期望维修费用率式(13)的准确性,在不同的维修及备件订购策略下采用Monte-Carlo(MC)方法对维修费用进行仿真,具体过程为:首先对系统退化过程进行仿真,得到备件订购时间、系统更换、系统故障时间,进而计算系统库存费用、更换费用及停机损失,最后得到费用具体值.在这里重复上述过程104 次,得到系统维修费用的104 个实现值.在此基础上,求得平均值并与解析结果进行比较,如表3 所示.
由表3 可以看出,解析方法得到的期望费用率与仿真结果很接近,并且在MC仿真结果95%置信区间内,由此说明了上述系统期望维修费用率表达式的准确性.
下面分析当备件订购阈值(Q A 1 ,Q A 2 ) 为常数时,系统期望维修费用率γ ∞ 随预防性更换阈值(Q M 1 ,Q M 2 ) 的变化情况,结果如图3 所示.
图3
图3
γ ∞ 随(Q M 1 ,Q M 2 ) 变化情况
Fig.3
γ ∞ versus (Q M 1 ,Q M 2 )
由图3 可以看出,无论对于案例1还是案例2,当(Q A 1 ,Q A 2 ) 为常数时,随着(Q M 1 ,Q M 2 ) 的增大,γ ∞ 总体上来说呈现先减小后增大的趋势,这是由库存费用、维修费用增加以及系统更新周期增长共同作用的结果.在给定的备件订购阈值下,存在最优的(Q M 1 ,Q M 2 ) 使得γ ∞ 达到最小.同理,下面固定预防性更换阈值(Q M 1 ,Q M 2 ) ,分析系统期望维修费用率γ ∞ 随备件订购阈值(Q A 1 ,Q A 2 ) 的变化情况,结果如图4 所示.
图4
图4
γ ∞ 随(Q A 1 ,Q A 2 ) 的变化情况
Fig.4
γ ∞ versus (Q A 1 ,Q A 2 )
由图4 可以看出,无论对于案例1还是案例2,当(Q M 1 ,Q M 2 ) 为常数时,随着(Q A 1 ,Q A 2 ) 的增大,γ ∞ 总体上来说呈现先减小后增大的趋势,这主要是受预防性更换概率减小、故障更换概率增大以及库存费用减少的影响.同样,在给定的预防性更换阈值下,存在最优的(Q A 1 ,Q A 2 ) 使得γ ∞ 达到最小.
下面基于ABC算法对模型进行求解,得到最优维修及备件订购策略.在ABC算法中,设置种群数量为10,适应度函数为1/γ ∞ (Q A ,Q M ).以案例1为例,考虑到算法中随机因素的影响,运行ABC算法5次,获得的系统维修及备件订购策略优化过程如图5 所示,其中:K 为迭代次数;上标*表示最优.
图5
图5
ABC算法优化过程
Fig.5
Optimization process of ABC algorithm
由图5 可以看出,当迭代次数达到100次后,最优期望费用率值不再变化,且各优化过程均收敛于相同的最优值,对案例2也有类似的结果.当退化过程增多时,可适当增加算法中种群数量的设置,以提高算法收敛速度及精度.在此基础上,得到的系统最优维修及备件订购策略如表4 所示.此外,假如在这里采取文献[14 ]中的策略,即不单独考虑系统备件订购阈值,而是在预防性更换阈值(Q M 1 ,Q M 2 ) 处订购备件,即Q A i =Q M i . 在此策略下,所获得的系统最优维修策略见表4 .
通过对比表4 中的数据可知,相对于文献[14 ]中的维修策略,综合优化备件订购阈值和预防性更换阈值能够有效降低系统期望维修费用率.在案例1和案例2下系统最优维修及备件订购策略对应的系统期望维修费用率γ ∞ * 分别降低了6.9%和5.9%,这是由于同时优化调整Q A i 和Q M i 值可以在降低停机费用的同时,增加更新周期长度,使得系统寿命得到更为充分的利用.
为验证系统期望维修费用率的近似表达式的精确性,下面对期望维修费用率近似值γ ~ ∞ 和精确值γ ∞ 进行对比.由于在推导γ ~ ∞ 时将X t A 用E (X t A ) 代替,所以这里考察不同备件订购阈值(Q A 1 ,Q A 2 ) 下,γ ~ ∞ 和γ ∞ 的差异.如图6 所示.
图6
图6
不同(Q A 1 ,Q A 2 ) 下γ ~ ∞ 和γ ∞ 对比
Fig.6
Comparison of γ ~ ∞ and γ ∞ at different (Q A 1 ,Q A 2 ) values
由图6 可以看出,当(Q A 1 ,Q A 2 ) 较小时,γ ~ ∞ 和γ ∞ 的值较为接近;随着(Q A 1 ,Q A 2 ) 的增大,γ ~ ∞ 和γ ∞ 的差值也逐渐增大,但最大差值仅为0.12,在精度要求不是很严格时,可以用γ ~ ∞ 代替γ ∞ 以简化计算.γ ~ ∞ 和γ ∞ 差值增大的原因主要是因为(Q A 1 ,Q A 2 ) 变大导致退化达到备件订购阈值的期望时间增大,从而退化量在该时间点的方差变大,所以导致估计结果的误差也相应变大.
5.2 灵敏度分析
首先分析备件交付时间对系统最优维修及备件订购策略的影响,当备件交付时间τ 发生变化时,对于案例1和案例2,(Q A 1 * ,Q A 2 * ,Q M 1 * ,Q M 2 * ) 和γ ∞ * 的变化情况如图7 所示.
图7
图7
备件交付时间对优化结果的影响
Fig.7
Influence of spare parts delivery time on optimization results
由图7 可以看出,无论是在案例1还是案例2的情形下,随着备件交付时间τ 的增长,最优备件订购阈值(Q A 1 * ,Q A 2 * ) 和预防性更换阈值(Q M 1 * ,Q M 2 * ) 均呈下降趋势,最优期望维修费用率γ ∞ * 则是逐渐上升的.这是由于总体上来说,随着备件交付时间的增长,导致不能对退化严重系统及时进行更换,维修费用和停机损失会增加,因此,需要降低备件订购阈值.这将会导致库存费用的增加,预防性更换阈值降低是为了平衡库存费用和预防性更换费用.此外,由图7 还可以看出τ 的变化对(Q A 1 * ,Q A 2 * ) 影响较大,而对(Q M 1 * ,Q M 2 * ) 影响相对较小.
下面分析θ 的变化对系统最优维修及备件订购策略的影响,在案例1和案例2规定的情形下,计算结果如表5 所示.
由表5 可以看出,当Q 发生变化时,系统的最优维修及备件订购策略也发生了一定变化.当相关系数增大时,预防性更换阈值(Q M 1 * ,Q M 2 * ) 的值逐渐减小,而备件订购阈值(Q A 1 * ,Q A 2 * ) 的变化趋势不明显,最优期望维修费用率γ ∞ * 值逐渐减小.γ ∞ * 值的减小是由于随着相关系数的增大,系统可靠性也随之提高的原因.
6 结语
多元退化系统CBM决策通常需要考虑不同退化过程间的相关性.在维修策略中同时考虑预防性更换阈值和备件订购阈值,在推导系统退化状态到达各阈值时间联合分布基础上,得到了连续监测条件下的系统期望维修费用率精确和近似表达式.在案例分析中,通过ABC算法得到了最优维修和备件订购策略,验证了模型和算法的有效性.通过和已有策略的对比表明了所提系统维修和备件订购策略的优势.同时,在案例中也检验了系统期望维修费用近似表达式的精度,分析了备件交付时间和退化相关程度参数对决策结果的影响.值得注意的是,本文采用Gamma过程建立退化过程模型,事实上,对于包含Gamma过程的Levy过程来说,上述模型也同样适用.此外,对于一些非Copula相关的多元退化系统在连续监测下的维修阈值优化问题,可参考上述方法建立优化模型.
参考文献
View Option
[1]
LEE J , WU F J , ZHAO W , et al . Prognostics and health management design for rotary machinery systems—Reviews, methodology and applications
[J]. Mechanical Systems and Signal Processing , 2014 , 42 (1/2 ): 314 -334 .
DOI:10.1016/j.ymssp.2013.06.004
URL
[本文引用: 1]
[2]
DE JONGE B , TEUNTER R , TINGA T . The influence of practical factors on the benefits of condition-based maintenance over time-based maintenance
[J]. Reliability Engineering & System Safety , 2017 , 158 :21 -30 .
DOI:10.1016/j.ress.2016.10.002
URL
[本文引用: 1]
[3]
OLDE KEIZER M C A , FLAPPER S D P , TEUNTER R H . Condition-based maintenance policies for systems with multiple dependent components: A review
[J]. European Journal of Operational Research , 2017 , 261 (2 ): 405 -420 .
DOI:10.1016/j.ejor.2017.02.044
URL
[本文引用: 1]
[4]
PAN Z Q , BALAKRISHNAN N . Reliability modeling of degradation of products with multiple performance characteristics based on gamma processes
[J]. Reliability Engineering & System Safety , 2011 , 96 (8 ): 949 -957 .
DOI:10.1016/j.ress.2011.03.014
URL
[本文引用: 2]
[5]
WANG X L , BALAKRISHNAN N , GUO B . Residual life estimation based on nonlinear-multivariate Wiener processes
[J]. Journal of Statistical Computation and Simulation , 2015 , 85 (9 ): 1742 -1764 .
DOI:10.1080/00949655.2014.898765
URL
[本文引用: 2]
[6]
RASMEKOMEN N , PARLIKAD A K . Condition-based maintenance of multi-component systems with degradation state-rate interactions
[J]. Reliability Engineering & System Safety , 2016 , 148 :1 -10 .
DOI:10.1016/j.ress.2015.11.010
URL
[本文引用: 2]
[7]
杨志远 , 赵建民 , 程中华 . 退化相关多部件系统预防性维修决策模型
[J]. 系统工程与电子技术 , 2018 , 40 (4 ): 823 -832 .
[本文引用: 2]
YANG Zhiyuan , ZHAO Jianmin , CHENG Zhonghua . Preventive maintenance decision model of multi-component system with degradation interaction
[J]. Systems Engineering and Electronics , 2018 , 40 (4 ): 823 -832 .
[本文引用: 2]
[8]
WANG X L , BALAKRISHNAN N , GUO B , et al . Residual life estimation based on bivariate non-stationary gamma degradation process
[J]. Journal of Statistical Computation and Simulation , 2015 , 85 (2 ): 405 -421 .
DOI:10.1080/00949655.2013.824448
URL
[本文引用: 2]
[9]
WANG X L , GUO B , CHENG Z J . Residual life estimation based on bivariate Wiener degradation process with time-scale transformations
[J]. Journal of Statistical Computation and Simulation , 2014 , 84 (3 ): 545 -563 .
DOI:10.1080/00949655.2012.719026
URL
[本文引用: 1]
[10]
HONG H P , ZHOU W , ZHANG S , et al . Optimal condition-based maintenance decisions for systems with dependent stochastic degradation of components
[J]. Reliability Engineering & System Safety , 2014 , 121 :276 -288 .
DOI:10.1016/j.ress.2013.09.004
URL
[本文引用: 1]
[11]
ZHANG N L , YANG Q Y . Optimal maintenance planning for repairable multi-component systems subject to dependent competing risks
[J]. IIE Transactions , 2015 , 47 (5 ): 521 -532 .
DOI:10.1080/0740817X.2014.974115
URL
[本文引用: 1]
[12]
LI H P , DELOUX E , DIEULLE L . A condition-based maintenance policy for multi-component systems with Lévy copulas dependence
[J]. Reliability Engineering & System Safety , 2016 , 149 :44 -55 .
DOI:10.1016/j.ress.2015.12.011
URL
[本文引用: 1]
[13]
BÉRENGUER C , GRALL A , DIEULLE L , et al . Maintenance policy for a continuously monitored deteriorating system
[J]. Probability in the Engineering and Informational Sciences , 2003 , 17 (2 ): 235 -250 .
DOI:10.1017/S0269964803172063
URL
[本文引用: 2]
[14]
MERCIER S , PHAM H H . A preventive maintenance policy for a continuously monitored system with correlated wear indicators
[J]. European Journal of Operational Research , 2012 , 222 (2 ): 263 -272 .
DOI:10.1016/j.ejor.2012.05.011
URL
[本文引用: 4]
[15]
CHERUBINI U , LUCIANO E , VECCHIATO W . Copula methods in finance [M]. Oxford, UK: John Wiley & Sons Ltd. , 2004 : 154 -160 .
[本文引用: 1]
[16]
王浩伟 , 徐廷学 , 刘勇 . 基于随机参数Gamma过程的剩余寿命预测方法
[J]. 浙江大学学报(工学版) , 2015 , 49 (4 ): 699 -704 .
[本文引用: 1]
WANG Haowei , XU Tingxue , LIU Yong . Remaining useful life prediction method based on Gamma processes with random parameters
[J]. Journal of Zhejiang University (Engineering Science) , 2015 , 49 (4 ): 699 -704 .
[本文引用: 1]
[18]
KARABOGA D , BASTURK B . A powerful and efficient algorithm for numerical function optimization: Artificial bee colony (ABC) algorithm
[J]. Journal of Global Optimization , 2007 , 39 (3 ): 459 -471 .
DOI:10.1007/s10898-007-9149-x
URL
[本文引用: 1]
Prognostics and health management design for rotary machinery systems—Reviews, methodology and applications
1
2014
... 现代复杂系统通常装备有状态监测单元如传感器、仪表、计算设备[1 ] ,这为基于状态的维修(CBM)实施提供了系统状态数据.相对于基于时间的维修方式,CBM能够有效控制“维修不足”和“维修过度”的问题[2 ] .通常来说,获取系统状态信息的方式可分为定期检测和连续监测两类,相比于定期检测,连续监测能够实时获取系统状态信息,从而更为及时有效地进行CBM决策,降低系统故障风险和维修费用.随着传感器技术的发展,连续监测的应用场景越来越多,特别是对于一些关键部件或系统,如发电机组,其本身价值高、停机损失大,对系统状态进行连续监测更为适用.同时,复杂系统通常存在多个退化过程,且由于共同的运行环境、应力等因素,退化过程间往往是相关的,在CBM决策中需要充分考虑这些因素[3 ] ,这也增加了多元退化系统CBM决策的难度. ...
The influence of practical factors on the benefits of condition-based maintenance over time-based maintenance
1
2017
... 现代复杂系统通常装备有状态监测单元如传感器、仪表、计算设备[1 ] ,这为基于状态的维修(CBM)实施提供了系统状态数据.相对于基于时间的维修方式,CBM能够有效控制“维修不足”和“维修过度”的问题[2 ] .通常来说,获取系统状态信息的方式可分为定期检测和连续监测两类,相比于定期检测,连续监测能够实时获取系统状态信息,从而更为及时有效地进行CBM决策,降低系统故障风险和维修费用.随着传感器技术的发展,连续监测的应用场景越来越多,特别是对于一些关键部件或系统,如发电机组,其本身价值高、停机损失大,对系统状态进行连续监测更为适用.同时,复杂系统通常存在多个退化过程,且由于共同的运行环境、应力等因素,退化过程间往往是相关的,在CBM决策中需要充分考虑这些因素[3 ] ,这也增加了多元退化系统CBM决策的难度. ...
Condition-based maintenance policies for systems with multiple dependent components: A review
1
2017
... 现代复杂系统通常装备有状态监测单元如传感器、仪表、计算设备[1 ] ,这为基于状态的维修(CBM)实施提供了系统状态数据.相对于基于时间的维修方式,CBM能够有效控制“维修不足”和“维修过度”的问题[2 ] .通常来说,获取系统状态信息的方式可分为定期检测和连续监测两类,相比于定期检测,连续监测能够实时获取系统状态信息,从而更为及时有效地进行CBM决策,降低系统故障风险和维修费用.随着传感器技术的发展,连续监测的应用场景越来越多,特别是对于一些关键部件或系统,如发电机组,其本身价值高、停机损失大,对系统状态进行连续监测更为适用.同时,复杂系统通常存在多个退化过程,且由于共同的运行环境、应力等因素,退化过程间往往是相关的,在CBM决策中需要充分考虑这些因素[3 ] ,这也增加了多元退化系统CBM决策的难度. ...
Reliability modeling of degradation of products with multiple performance characteristics based on gamma processes
2
2011
... 对于有多个相关退化过程的系统而言,CBM决策的基础是对不同退化过程间的相关性进行建模描述.目前,常用的方法包括多维退化模型[4 ,5 ] 、退化率相关模型[6 ,7 ] 和基于Copula函数的退化相关模型[8 ,9 ] .文献[4 ,5 ]分别采用多元正态分布和二元Birnbaum-Saunders(B-S)分布来描述不同退化过程间的相关关系,而分布形式也限制了此类方法的应用.文献[6 ,7 ]将不同退化过程间的相关关系表示为一个部件退化水平对其余部件退化速率的影响.在工程实际中,这种影响关系及程度很难采用定量化方法描述和确定.相对而言,Copula函数在定量描述随机变量相关性方面具有很高的灵活性,并且提供了分析相关结构对系统可靠性影响的参数化方法.基于不同种类的Copula相关结构,文献[10 ]利用仿真方法对定期检测策略下的系统维修费用进行分析,其中预防性维修阈值为固定值.文献[11 ]建立了多部件系统检测维修策略优化模型,其中部件故障间的相关性采用Gaussian Copula函数来描述.文献[12 ]采用Levy Copula函数描述系统二元退化间的相关关系,并在此基础上提出系统定期检测与机会维修策略,利用仿真方法获得了最优预防性维修阈值和机会维修概率阈值.在现有的基于Copula函数建立多元退化模型的文献中,很少研究在连续监测条件下的CBM决策优化问题,且没有考虑系统备件约束,而对一些价值较高的重要系统而言,很难保证时刻有可用备件,通常需要提前预定.在连续监测条件下,文献[13 ]研究了单一退化过程的预防性更换阈值优化问题.在此基础上,文献[14 ]考虑二元相关Levy退化过程,以系统长期运行条件下的期望维修费用率和可用度为准则,计算获得在不同系统结构下的最优预防性更换阈值.上述研究均是在维修延迟(维修资源的到达需要一定时间,此处的维修资源也可指备件)条件下,对预防性更换阈值进行决策,没有综合考虑备件订购阈值以及预防性维修阈值的优化问题. ...
... .文献[4 ,5 ]分别采用多元正态分布和二元Birnbaum-Saunders(B-S)分布来描述不同退化过程间的相关关系,而分布形式也限制了此类方法的应用.文献[6 ,7 ]将不同退化过程间的相关关系表示为一个部件退化水平对其余部件退化速率的影响.在工程实际中,这种影响关系及程度很难采用定量化方法描述和确定.相对而言,Copula函数在定量描述随机变量相关性方面具有很高的灵活性,并且提供了分析相关结构对系统可靠性影响的参数化方法.基于不同种类的Copula相关结构,文献[10 ]利用仿真方法对定期检测策略下的系统维修费用进行分析,其中预防性维修阈值为固定值.文献[11 ]建立了多部件系统检测维修策略优化模型,其中部件故障间的相关性采用Gaussian Copula函数来描述.文献[12 ]采用Levy Copula函数描述系统二元退化间的相关关系,并在此基础上提出系统定期检测与机会维修策略,利用仿真方法获得了最优预防性维修阈值和机会维修概率阈值.在现有的基于Copula函数建立多元退化模型的文献中,很少研究在连续监测条件下的CBM决策优化问题,且没有考虑系统备件约束,而对一些价值较高的重要系统而言,很难保证时刻有可用备件,通常需要提前预定.在连续监测条件下,文献[13 ]研究了单一退化过程的预防性更换阈值优化问题.在此基础上,文献[14 ]考虑二元相关Levy退化过程,以系统长期运行条件下的期望维修费用率和可用度为准则,计算获得在不同系统结构下的最优预防性更换阈值.上述研究均是在维修延迟(维修资源的到达需要一定时间,此处的维修资源也可指备件)条件下,对预防性更换阈值进行决策,没有综合考虑备件订购阈值以及预防性维修阈值的优化问题. ...
Residual life estimation based on nonlinear-multivariate Wiener processes
2
2015
... 对于有多个相关退化过程的系统而言,CBM决策的基础是对不同退化过程间的相关性进行建模描述.目前,常用的方法包括多维退化模型[4 ,5 ] 、退化率相关模型[6 ,7 ] 和基于Copula函数的退化相关模型[8 ,9 ] .文献[4 ,5 ]分别采用多元正态分布和二元Birnbaum-Saunders(B-S)分布来描述不同退化过程间的相关关系,而分布形式也限制了此类方法的应用.文献[6 ,7 ]将不同退化过程间的相关关系表示为一个部件退化水平对其余部件退化速率的影响.在工程实际中,这种影响关系及程度很难采用定量化方法描述和确定.相对而言,Copula函数在定量描述随机变量相关性方面具有很高的灵活性,并且提供了分析相关结构对系统可靠性影响的参数化方法.基于不同种类的Copula相关结构,文献[10 ]利用仿真方法对定期检测策略下的系统维修费用进行分析,其中预防性维修阈值为固定值.文献[11 ]建立了多部件系统检测维修策略优化模型,其中部件故障间的相关性采用Gaussian Copula函数来描述.文献[12 ]采用Levy Copula函数描述系统二元退化间的相关关系,并在此基础上提出系统定期检测与机会维修策略,利用仿真方法获得了最优预防性维修阈值和机会维修概率阈值.在现有的基于Copula函数建立多元退化模型的文献中,很少研究在连续监测条件下的CBM决策优化问题,且没有考虑系统备件约束,而对一些价值较高的重要系统而言,很难保证时刻有可用备件,通常需要提前预定.在连续监测条件下,文献[13 ]研究了单一退化过程的预防性更换阈值优化问题.在此基础上,文献[14 ]考虑二元相关Levy退化过程,以系统长期运行条件下的期望维修费用率和可用度为准则,计算获得在不同系统结构下的最优预防性更换阈值.上述研究均是在维修延迟(维修资源的到达需要一定时间,此处的维修资源也可指备件)条件下,对预防性更换阈值进行决策,没有综合考虑备件订购阈值以及预防性维修阈值的优化问题. ...
... ,5 ]分别采用多元正态分布和二元Birnbaum-Saunders(B-S)分布来描述不同退化过程间的相关关系,而分布形式也限制了此类方法的应用.文献[6 ,7 ]将不同退化过程间的相关关系表示为一个部件退化水平对其余部件退化速率的影响.在工程实际中,这种影响关系及程度很难采用定量化方法描述和确定.相对而言,Copula函数在定量描述随机变量相关性方面具有很高的灵活性,并且提供了分析相关结构对系统可靠性影响的参数化方法.基于不同种类的Copula相关结构,文献[10 ]利用仿真方法对定期检测策略下的系统维修费用进行分析,其中预防性维修阈值为固定值.文献[11 ]建立了多部件系统检测维修策略优化模型,其中部件故障间的相关性采用Gaussian Copula函数来描述.文献[12 ]采用Levy Copula函数描述系统二元退化间的相关关系,并在此基础上提出系统定期检测与机会维修策略,利用仿真方法获得了最优预防性维修阈值和机会维修概率阈值.在现有的基于Copula函数建立多元退化模型的文献中,很少研究在连续监测条件下的CBM决策优化问题,且没有考虑系统备件约束,而对一些价值较高的重要系统而言,很难保证时刻有可用备件,通常需要提前预定.在连续监测条件下,文献[13 ]研究了单一退化过程的预防性更换阈值优化问题.在此基础上,文献[14 ]考虑二元相关Levy退化过程,以系统长期运行条件下的期望维修费用率和可用度为准则,计算获得在不同系统结构下的最优预防性更换阈值.上述研究均是在维修延迟(维修资源的到达需要一定时间,此处的维修资源也可指备件)条件下,对预防性更换阈值进行决策,没有综合考虑备件订购阈值以及预防性维修阈值的优化问题. ...
Condition-based maintenance of multi-component systems with degradation state-rate interactions
2
2016
... 对于有多个相关退化过程的系统而言,CBM决策的基础是对不同退化过程间的相关性进行建模描述.目前,常用的方法包括多维退化模型[4 ,5 ] 、退化率相关模型[6 ,7 ] 和基于Copula函数的退化相关模型[8 ,9 ] .文献[4 ,5 ]分别采用多元正态分布和二元Birnbaum-Saunders(B-S)分布来描述不同退化过程间的相关关系,而分布形式也限制了此类方法的应用.文献[6 ,7 ]将不同退化过程间的相关关系表示为一个部件退化水平对其余部件退化速率的影响.在工程实际中,这种影响关系及程度很难采用定量化方法描述和确定.相对而言,Copula函数在定量描述随机变量相关性方面具有很高的灵活性,并且提供了分析相关结构对系统可靠性影响的参数化方法.基于不同种类的Copula相关结构,文献[10 ]利用仿真方法对定期检测策略下的系统维修费用进行分析,其中预防性维修阈值为固定值.文献[11 ]建立了多部件系统检测维修策略优化模型,其中部件故障间的相关性采用Gaussian Copula函数来描述.文献[12 ]采用Levy Copula函数描述系统二元退化间的相关关系,并在此基础上提出系统定期检测与机会维修策略,利用仿真方法获得了最优预防性维修阈值和机会维修概率阈值.在现有的基于Copula函数建立多元退化模型的文献中,很少研究在连续监测条件下的CBM决策优化问题,且没有考虑系统备件约束,而对一些价值较高的重要系统而言,很难保证时刻有可用备件,通常需要提前预定.在连续监测条件下,文献[13 ]研究了单一退化过程的预防性更换阈值优化问题.在此基础上,文献[14 ]考虑二元相关Levy退化过程,以系统长期运行条件下的期望维修费用率和可用度为准则,计算获得在不同系统结构下的最优预防性更换阈值.上述研究均是在维修延迟(维修资源的到达需要一定时间,此处的维修资源也可指备件)条件下,对预防性更换阈值进行决策,没有综合考虑备件订购阈值以及预防性维修阈值的优化问题. ...
... ]分别采用多元正态分布和二元Birnbaum-Saunders(B-S)分布来描述不同退化过程间的相关关系,而分布形式也限制了此类方法的应用.文献[6 ,7 ]将不同退化过程间的相关关系表示为一个部件退化水平对其余部件退化速率的影响.在工程实际中,这种影响关系及程度很难采用定量化方法描述和确定.相对而言,Copula函数在定量描述随机变量相关性方面具有很高的灵活性,并且提供了分析相关结构对系统可靠性影响的参数化方法.基于不同种类的Copula相关结构,文献[10 ]利用仿真方法对定期检测策略下的系统维修费用进行分析,其中预防性维修阈值为固定值.文献[11 ]建立了多部件系统检测维修策略优化模型,其中部件故障间的相关性采用Gaussian Copula函数来描述.文献[12 ]采用Levy Copula函数描述系统二元退化间的相关关系,并在此基础上提出系统定期检测与机会维修策略,利用仿真方法获得了最优预防性维修阈值和机会维修概率阈值.在现有的基于Copula函数建立多元退化模型的文献中,很少研究在连续监测条件下的CBM决策优化问题,且没有考虑系统备件约束,而对一些价值较高的重要系统而言,很难保证时刻有可用备件,通常需要提前预定.在连续监测条件下,文献[13 ]研究了单一退化过程的预防性更换阈值优化问题.在此基础上,文献[14 ]考虑二元相关Levy退化过程,以系统长期运行条件下的期望维修费用率和可用度为准则,计算获得在不同系统结构下的最优预防性更换阈值.上述研究均是在维修延迟(维修资源的到达需要一定时间,此处的维修资源也可指备件)条件下,对预防性更换阈值进行决策,没有综合考虑备件订购阈值以及预防性维修阈值的优化问题. ...
退化相关多部件系统预防性维修决策模型
2
2018
... 对于有多个相关退化过程的系统而言,CBM决策的基础是对不同退化过程间的相关性进行建模描述.目前,常用的方法包括多维退化模型[4 ,5 ] 、退化率相关模型[6 ,7 ] 和基于Copula函数的退化相关模型[8 ,9 ] .文献[4 ,5 ]分别采用多元正态分布和二元Birnbaum-Saunders(B-S)分布来描述不同退化过程间的相关关系,而分布形式也限制了此类方法的应用.文献[6 ,7 ]将不同退化过程间的相关关系表示为一个部件退化水平对其余部件退化速率的影响.在工程实际中,这种影响关系及程度很难采用定量化方法描述和确定.相对而言,Copula函数在定量描述随机变量相关性方面具有很高的灵活性,并且提供了分析相关结构对系统可靠性影响的参数化方法.基于不同种类的Copula相关结构,文献[10 ]利用仿真方法对定期检测策略下的系统维修费用进行分析,其中预防性维修阈值为固定值.文献[11 ]建立了多部件系统检测维修策略优化模型,其中部件故障间的相关性采用Gaussian Copula函数来描述.文献[12 ]采用Levy Copula函数描述系统二元退化间的相关关系,并在此基础上提出系统定期检测与机会维修策略,利用仿真方法获得了最优预防性维修阈值和机会维修概率阈值.在现有的基于Copula函数建立多元退化模型的文献中,很少研究在连续监测条件下的CBM决策优化问题,且没有考虑系统备件约束,而对一些价值较高的重要系统而言,很难保证时刻有可用备件,通常需要提前预定.在连续监测条件下,文献[13 ]研究了单一退化过程的预防性更换阈值优化问题.在此基础上,文献[14 ]考虑二元相关Levy退化过程,以系统长期运行条件下的期望维修费用率和可用度为准则,计算获得在不同系统结构下的最优预防性更换阈值.上述研究均是在维修延迟(维修资源的到达需要一定时间,此处的维修资源也可指备件)条件下,对预防性更换阈值进行决策,没有综合考虑备件订购阈值以及预防性维修阈值的优化问题. ...
... ,7 ]将不同退化过程间的相关关系表示为一个部件退化水平对其余部件退化速率的影响.在工程实际中,这种影响关系及程度很难采用定量化方法描述和确定.相对而言,Copula函数在定量描述随机变量相关性方面具有很高的灵活性,并且提供了分析相关结构对系统可靠性影响的参数化方法.基于不同种类的Copula相关结构,文献[10 ]利用仿真方法对定期检测策略下的系统维修费用进行分析,其中预防性维修阈值为固定值.文献[11 ]建立了多部件系统检测维修策略优化模型,其中部件故障间的相关性采用Gaussian Copula函数来描述.文献[12 ]采用Levy Copula函数描述系统二元退化间的相关关系,并在此基础上提出系统定期检测与机会维修策略,利用仿真方法获得了最优预防性维修阈值和机会维修概率阈值.在现有的基于Copula函数建立多元退化模型的文献中,很少研究在连续监测条件下的CBM决策优化问题,且没有考虑系统备件约束,而对一些价值较高的重要系统而言,很难保证时刻有可用备件,通常需要提前预定.在连续监测条件下,文献[13 ]研究了单一退化过程的预防性更换阈值优化问题.在此基础上,文献[14 ]考虑二元相关Levy退化过程,以系统长期运行条件下的期望维修费用率和可用度为准则,计算获得在不同系统结构下的最优预防性更换阈值.上述研究均是在维修延迟(维修资源的到达需要一定时间,此处的维修资源也可指备件)条件下,对预防性更换阈值进行决策,没有综合考虑备件订购阈值以及预防性维修阈值的优化问题. ...
退化相关多部件系统预防性维修决策模型
2
2018
... 对于有多个相关退化过程的系统而言,CBM决策的基础是对不同退化过程间的相关性进行建模描述.目前,常用的方法包括多维退化模型[4 ,5 ] 、退化率相关模型[6 ,7 ] 和基于Copula函数的退化相关模型[8 ,9 ] .文献[4 ,5 ]分别采用多元正态分布和二元Birnbaum-Saunders(B-S)分布来描述不同退化过程间的相关关系,而分布形式也限制了此类方法的应用.文献[6 ,7 ]将不同退化过程间的相关关系表示为一个部件退化水平对其余部件退化速率的影响.在工程实际中,这种影响关系及程度很难采用定量化方法描述和确定.相对而言,Copula函数在定量描述随机变量相关性方面具有很高的灵活性,并且提供了分析相关结构对系统可靠性影响的参数化方法.基于不同种类的Copula相关结构,文献[10 ]利用仿真方法对定期检测策略下的系统维修费用进行分析,其中预防性维修阈值为固定值.文献[11 ]建立了多部件系统检测维修策略优化模型,其中部件故障间的相关性采用Gaussian Copula函数来描述.文献[12 ]采用Levy Copula函数描述系统二元退化间的相关关系,并在此基础上提出系统定期检测与机会维修策略,利用仿真方法获得了最优预防性维修阈值和机会维修概率阈值.在现有的基于Copula函数建立多元退化模型的文献中,很少研究在连续监测条件下的CBM决策优化问题,且没有考虑系统备件约束,而对一些价值较高的重要系统而言,很难保证时刻有可用备件,通常需要提前预定.在连续监测条件下,文献[13 ]研究了单一退化过程的预防性更换阈值优化问题.在此基础上,文献[14 ]考虑二元相关Levy退化过程,以系统长期运行条件下的期望维修费用率和可用度为准则,计算获得在不同系统结构下的最优预防性更换阈值.上述研究均是在维修延迟(维修资源的到达需要一定时间,此处的维修资源也可指备件)条件下,对预防性更换阈值进行决策,没有综合考虑备件订购阈值以及预防性维修阈值的优化问题. ...
... ,7 ]将不同退化过程间的相关关系表示为一个部件退化水平对其余部件退化速率的影响.在工程实际中,这种影响关系及程度很难采用定量化方法描述和确定.相对而言,Copula函数在定量描述随机变量相关性方面具有很高的灵活性,并且提供了分析相关结构对系统可靠性影响的参数化方法.基于不同种类的Copula相关结构,文献[10 ]利用仿真方法对定期检测策略下的系统维修费用进行分析,其中预防性维修阈值为固定值.文献[11 ]建立了多部件系统检测维修策略优化模型,其中部件故障间的相关性采用Gaussian Copula函数来描述.文献[12 ]采用Levy Copula函数描述系统二元退化间的相关关系,并在此基础上提出系统定期检测与机会维修策略,利用仿真方法获得了最优预防性维修阈值和机会维修概率阈值.在现有的基于Copula函数建立多元退化模型的文献中,很少研究在连续监测条件下的CBM决策优化问题,且没有考虑系统备件约束,而对一些价值较高的重要系统而言,很难保证时刻有可用备件,通常需要提前预定.在连续监测条件下,文献[13 ]研究了单一退化过程的预防性更换阈值优化问题.在此基础上,文献[14 ]考虑二元相关Levy退化过程,以系统长期运行条件下的期望维修费用率和可用度为准则,计算获得在不同系统结构下的最优预防性更换阈值.上述研究均是在维修延迟(维修资源的到达需要一定时间,此处的维修资源也可指备件)条件下,对预防性更换阈值进行决策,没有综合考虑备件订购阈值以及预防性维修阈值的优化问题. ...
Residual life estimation based on bivariate non-stationary gamma degradation process
2
2015
... 对于有多个相关退化过程的系统而言,CBM决策的基础是对不同退化过程间的相关性进行建模描述.目前,常用的方法包括多维退化模型[4 ,5 ] 、退化率相关模型[6 ,7 ] 和基于Copula函数的退化相关模型[8 ,9 ] .文献[4 ,5 ]分别采用多元正态分布和二元Birnbaum-Saunders(B-S)分布来描述不同退化过程间的相关关系,而分布形式也限制了此类方法的应用.文献[6 ,7 ]将不同退化过程间的相关关系表示为一个部件退化水平对其余部件退化速率的影响.在工程实际中,这种影响关系及程度很难采用定量化方法描述和确定.相对而言,Copula函数在定量描述随机变量相关性方面具有很高的灵活性,并且提供了分析相关结构对系统可靠性影响的参数化方法.基于不同种类的Copula相关结构,文献[10 ]利用仿真方法对定期检测策略下的系统维修费用进行分析,其中预防性维修阈值为固定值.文献[11 ]建立了多部件系统检测维修策略优化模型,其中部件故障间的相关性采用Gaussian Copula函数来描述.文献[12 ]采用Levy Copula函数描述系统二元退化间的相关关系,并在此基础上提出系统定期检测与机会维修策略,利用仿真方法获得了最优预防性维修阈值和机会维修概率阈值.在现有的基于Copula函数建立多元退化模型的文献中,很少研究在连续监测条件下的CBM决策优化问题,且没有考虑系统备件约束,而对一些价值较高的重要系统而言,很难保证时刻有可用备件,通常需要提前预定.在连续监测条件下,文献[13 ]研究了单一退化过程的预防性更换阈值优化问题.在此基础上,文献[14 ]考虑二元相关Levy退化过程,以系统长期运行条件下的期望维修费用率和可用度为准则,计算获得在不同系统结构下的最优预防性更换阈值.上述研究均是在维修延迟(维修资源的到达需要一定时间,此处的维修资源也可指备件)条件下,对预防性更换阈值进行决策,没有综合考虑备件订购阈值以及预防性维修阈值的优化问题. ...
... 如上所述,本文采用Copula函数描述退化过程的相关性.为便于实际应用,这里假设:在不同时间间隔内,不同退化过程的退化过程增量间的相关性可忽略.在文献[8 ]中也采用了类似的假设分析多元Wiener退化过程的相关性.令t >u >0,根据Sklar定理,系统n 个退化过程在时间区间[u ,t ]上退化增量Xi (t -u )的联合概率分布函数可表示为 ...
Residual life estimation based on bivariate Wiener degradation process with time-scale transformations
1
2014
... 对于有多个相关退化过程的系统而言,CBM决策的基础是对不同退化过程间的相关性进行建模描述.目前,常用的方法包括多维退化模型[4 ,5 ] 、退化率相关模型[6 ,7 ] 和基于Copula函数的退化相关模型[8 ,9 ] .文献[4 ,5 ]分别采用多元正态分布和二元Birnbaum-Saunders(B-S)分布来描述不同退化过程间的相关关系,而分布形式也限制了此类方法的应用.文献[6 ,7 ]将不同退化过程间的相关关系表示为一个部件退化水平对其余部件退化速率的影响.在工程实际中,这种影响关系及程度很难采用定量化方法描述和确定.相对而言,Copula函数在定量描述随机变量相关性方面具有很高的灵活性,并且提供了分析相关结构对系统可靠性影响的参数化方法.基于不同种类的Copula相关结构,文献[10 ]利用仿真方法对定期检测策略下的系统维修费用进行分析,其中预防性维修阈值为固定值.文献[11 ]建立了多部件系统检测维修策略优化模型,其中部件故障间的相关性采用Gaussian Copula函数来描述.文献[12 ]采用Levy Copula函数描述系统二元退化间的相关关系,并在此基础上提出系统定期检测与机会维修策略,利用仿真方法获得了最优预防性维修阈值和机会维修概率阈值.在现有的基于Copula函数建立多元退化模型的文献中,很少研究在连续监测条件下的CBM决策优化问题,且没有考虑系统备件约束,而对一些价值较高的重要系统而言,很难保证时刻有可用备件,通常需要提前预定.在连续监测条件下,文献[13 ]研究了单一退化过程的预防性更换阈值优化问题.在此基础上,文献[14 ]考虑二元相关Levy退化过程,以系统长期运行条件下的期望维修费用率和可用度为准则,计算获得在不同系统结构下的最优预防性更换阈值.上述研究均是在维修延迟(维修资源的到达需要一定时间,此处的维修资源也可指备件)条件下,对预防性更换阈值进行决策,没有综合考虑备件订购阈值以及预防性维修阈值的优化问题. ...
Optimal condition-based maintenance decisions for systems with dependent stochastic degradation of components
1
2014
... 对于有多个相关退化过程的系统而言,CBM决策的基础是对不同退化过程间的相关性进行建模描述.目前,常用的方法包括多维退化模型[4 ,5 ] 、退化率相关模型[6 ,7 ] 和基于Copula函数的退化相关模型[8 ,9 ] .文献[4 ,5 ]分别采用多元正态分布和二元Birnbaum-Saunders(B-S)分布来描述不同退化过程间的相关关系,而分布形式也限制了此类方法的应用.文献[6 ,7 ]将不同退化过程间的相关关系表示为一个部件退化水平对其余部件退化速率的影响.在工程实际中,这种影响关系及程度很难采用定量化方法描述和确定.相对而言,Copula函数在定量描述随机变量相关性方面具有很高的灵活性,并且提供了分析相关结构对系统可靠性影响的参数化方法.基于不同种类的Copula相关结构,文献[10 ]利用仿真方法对定期检测策略下的系统维修费用进行分析,其中预防性维修阈值为固定值.文献[11 ]建立了多部件系统检测维修策略优化模型,其中部件故障间的相关性采用Gaussian Copula函数来描述.文献[12 ]采用Levy Copula函数描述系统二元退化间的相关关系,并在此基础上提出系统定期检测与机会维修策略,利用仿真方法获得了最优预防性维修阈值和机会维修概率阈值.在现有的基于Copula函数建立多元退化模型的文献中,很少研究在连续监测条件下的CBM决策优化问题,且没有考虑系统备件约束,而对一些价值较高的重要系统而言,很难保证时刻有可用备件,通常需要提前预定.在连续监测条件下,文献[13 ]研究了单一退化过程的预防性更换阈值优化问题.在此基础上,文献[14 ]考虑二元相关Levy退化过程,以系统长期运行条件下的期望维修费用率和可用度为准则,计算获得在不同系统结构下的最优预防性更换阈值.上述研究均是在维修延迟(维修资源的到达需要一定时间,此处的维修资源也可指备件)条件下,对预防性更换阈值进行决策,没有综合考虑备件订购阈值以及预防性维修阈值的优化问题. ...
Optimal maintenance planning for repairable multi-component systems subject to dependent competing risks
1
2015
... 对于有多个相关退化过程的系统而言,CBM决策的基础是对不同退化过程间的相关性进行建模描述.目前,常用的方法包括多维退化模型[4 ,5 ] 、退化率相关模型[6 ,7 ] 和基于Copula函数的退化相关模型[8 ,9 ] .文献[4 ,5 ]分别采用多元正态分布和二元Birnbaum-Saunders(B-S)分布来描述不同退化过程间的相关关系,而分布形式也限制了此类方法的应用.文献[6 ,7 ]将不同退化过程间的相关关系表示为一个部件退化水平对其余部件退化速率的影响.在工程实际中,这种影响关系及程度很难采用定量化方法描述和确定.相对而言,Copula函数在定量描述随机变量相关性方面具有很高的灵活性,并且提供了分析相关结构对系统可靠性影响的参数化方法.基于不同种类的Copula相关结构,文献[10 ]利用仿真方法对定期检测策略下的系统维修费用进行分析,其中预防性维修阈值为固定值.文献[11 ]建立了多部件系统检测维修策略优化模型,其中部件故障间的相关性采用Gaussian Copula函数来描述.文献[12 ]采用Levy Copula函数描述系统二元退化间的相关关系,并在此基础上提出系统定期检测与机会维修策略,利用仿真方法获得了最优预防性维修阈值和机会维修概率阈值.在现有的基于Copula函数建立多元退化模型的文献中,很少研究在连续监测条件下的CBM决策优化问题,且没有考虑系统备件约束,而对一些价值较高的重要系统而言,很难保证时刻有可用备件,通常需要提前预定.在连续监测条件下,文献[13 ]研究了单一退化过程的预防性更换阈值优化问题.在此基础上,文献[14 ]考虑二元相关Levy退化过程,以系统长期运行条件下的期望维修费用率和可用度为准则,计算获得在不同系统结构下的最优预防性更换阈值.上述研究均是在维修延迟(维修资源的到达需要一定时间,此处的维修资源也可指备件)条件下,对预防性更换阈值进行决策,没有综合考虑备件订购阈值以及预防性维修阈值的优化问题. ...
A condition-based maintenance policy for multi-component systems with Lévy copulas dependence
1
2016
... 对于有多个相关退化过程的系统而言,CBM决策的基础是对不同退化过程间的相关性进行建模描述.目前,常用的方法包括多维退化模型[4 ,5 ] 、退化率相关模型[6 ,7 ] 和基于Copula函数的退化相关模型[8 ,9 ] .文献[4 ,5 ]分别采用多元正态分布和二元Birnbaum-Saunders(B-S)分布来描述不同退化过程间的相关关系,而分布形式也限制了此类方法的应用.文献[6 ,7 ]将不同退化过程间的相关关系表示为一个部件退化水平对其余部件退化速率的影响.在工程实际中,这种影响关系及程度很难采用定量化方法描述和确定.相对而言,Copula函数在定量描述随机变量相关性方面具有很高的灵活性,并且提供了分析相关结构对系统可靠性影响的参数化方法.基于不同种类的Copula相关结构,文献[10 ]利用仿真方法对定期检测策略下的系统维修费用进行分析,其中预防性维修阈值为固定值.文献[11 ]建立了多部件系统检测维修策略优化模型,其中部件故障间的相关性采用Gaussian Copula函数来描述.文献[12 ]采用Levy Copula函数描述系统二元退化间的相关关系,并在此基础上提出系统定期检测与机会维修策略,利用仿真方法获得了最优预防性维修阈值和机会维修概率阈值.在现有的基于Copula函数建立多元退化模型的文献中,很少研究在连续监测条件下的CBM决策优化问题,且没有考虑系统备件约束,而对一些价值较高的重要系统而言,很难保证时刻有可用备件,通常需要提前预定.在连续监测条件下,文献[13 ]研究了单一退化过程的预防性更换阈值优化问题.在此基础上,文献[14 ]考虑二元相关Levy退化过程,以系统长期运行条件下的期望维修费用率和可用度为准则,计算获得在不同系统结构下的最优预防性更换阈值.上述研究均是在维修延迟(维修资源的到达需要一定时间,此处的维修资源也可指备件)条件下,对预防性更换阈值进行决策,没有综合考虑备件订购阈值以及预防性维修阈值的优化问题. ...
Maintenance policy for a continuously monitored deteriorating system
2
2003
... 对于有多个相关退化过程的系统而言,CBM决策的基础是对不同退化过程间的相关性进行建模描述.目前,常用的方法包括多维退化模型[4 ,5 ] 、退化率相关模型[6 ,7 ] 和基于Copula函数的退化相关模型[8 ,9 ] .文献[4 ,5 ]分别采用多元正态分布和二元Birnbaum-Saunders(B-S)分布来描述不同退化过程间的相关关系,而分布形式也限制了此类方法的应用.文献[6 ,7 ]将不同退化过程间的相关关系表示为一个部件退化水平对其余部件退化速率的影响.在工程实际中,这种影响关系及程度很难采用定量化方法描述和确定.相对而言,Copula函数在定量描述随机变量相关性方面具有很高的灵活性,并且提供了分析相关结构对系统可靠性影响的参数化方法.基于不同种类的Copula相关结构,文献[10 ]利用仿真方法对定期检测策略下的系统维修费用进行分析,其中预防性维修阈值为固定值.文献[11 ]建立了多部件系统检测维修策略优化模型,其中部件故障间的相关性采用Gaussian Copula函数来描述.文献[12 ]采用Levy Copula函数描述系统二元退化间的相关关系,并在此基础上提出系统定期检测与机会维修策略,利用仿真方法获得了最优预防性维修阈值和机会维修概率阈值.在现有的基于Copula函数建立多元退化模型的文献中,很少研究在连续监测条件下的CBM决策优化问题,且没有考虑系统备件约束,而对一些价值较高的重要系统而言,很难保证时刻有可用备件,通常需要提前预定.在连续监测条件下,文献[13 ]研究了单一退化过程的预防性更换阈值优化问题.在此基础上,文献[14 ]考虑二元相关Levy退化过程,以系统长期运行条件下的期望维修费用率和可用度为准则,计算获得在不同系统结构下的最优预防性更换阈值.上述研究均是在维修延迟(维修资源的到达需要一定时间,此处的维修资源也可指备件)条件下,对预防性更换阈值进行决策,没有综合考虑备件订购阈值以及预防性维修阈值的优化问题. ...
... E sup ( χ t M ) =max{ E ( X 1 ( t M ) ) Q L 1 , E ( X 2 ( t M ) ) Q L 2 ,…, E ( X n ( t M ) ) Q L n },系统各性能退化均为Gamma过程,在任意时刻退化量联合分布的边缘分布为Gamma分布,根据鞅停时定理[13 ] ,容易得到: ...
A preventive maintenance policy for a continuously monitored system with correlated wear indicators
4
2012
... 对于有多个相关退化过程的系统而言,CBM决策的基础是对不同退化过程间的相关性进行建模描述.目前,常用的方法包括多维退化模型[4 ,5 ] 、退化率相关模型[6 ,7 ] 和基于Copula函数的退化相关模型[8 ,9 ] .文献[4 ,5 ]分别采用多元正态分布和二元Birnbaum-Saunders(B-S)分布来描述不同退化过程间的相关关系,而分布形式也限制了此类方法的应用.文献[6 ,7 ]将不同退化过程间的相关关系表示为一个部件退化水平对其余部件退化速率的影响.在工程实际中,这种影响关系及程度很难采用定量化方法描述和确定.相对而言,Copula函数在定量描述随机变量相关性方面具有很高的灵活性,并且提供了分析相关结构对系统可靠性影响的参数化方法.基于不同种类的Copula相关结构,文献[10 ]利用仿真方法对定期检测策略下的系统维修费用进行分析,其中预防性维修阈值为固定值.文献[11 ]建立了多部件系统检测维修策略优化模型,其中部件故障间的相关性采用Gaussian Copula函数来描述.文献[12 ]采用Levy Copula函数描述系统二元退化间的相关关系,并在此基础上提出系统定期检测与机会维修策略,利用仿真方法获得了最优预防性维修阈值和机会维修概率阈值.在现有的基于Copula函数建立多元退化模型的文献中,很少研究在连续监测条件下的CBM决策优化问题,且没有考虑系统备件约束,而对一些价值较高的重要系统而言,很难保证时刻有可用备件,通常需要提前预定.在连续监测条件下,文献[13 ]研究了单一退化过程的预防性更换阈值优化问题.在此基础上,文献[14 ]考虑二元相关Levy退化过程,以系统长期运行条件下的期望维修费用率和可用度为准则,计算获得在不同系统结构下的最优预防性更换阈值.上述研究均是在维修延迟(维修资源的到达需要一定时间,此处的维修资源也可指备件)条件下,对预防性更换阈值进行决策,没有综合考虑备件订购阈值以及预防性维修阈值的优化问题. ...
... 由图5 可以看出,当迭代次数达到100次后,最优期望费用率值不再变化,且各优化过程均收敛于相同的最优值,对案例2也有类似的结果.当退化过程增多时,可适当增加算法中种群数量的设置,以提高算法收敛速度及精度.在此基础上,得到的系统最优维修及备件订购策略如表4 所示.此外,假如在这里采取文献[14 ]中的策略,即不单独考虑系统备件订购阈值,而是在预防性更换阈值( Q M 1 , Q M 2 ) 处订购备件,即 Q A i = Q M i . 在此策略下,所获得的系统最优维修策略见表4 . ...
... Optimal system maintenance and spare parts ordering policy
Tab.4 维修策略 案例1 案例2 $Q_{{A}_{1}}^{*}$ $Q_{{A}_{2}}^{*}$ $Q_{{M}_{1}}^{*}$ $Q_{{M}_{2}}^{*}$ $γ_{∞}^{*}$ $Q_{{A}_{1}}^{*}$ $Q_{{A}_{2}}^{*}$ $Q_{{M}_{1}}^{*}$ $Q_{{M}_{2}}^{*}$ $γ_{∞}^{*}$ 本文 5.92 5.92 8.03 8.03 8.922 6.07 5.35 7.96 6.50 8.649 文献[14 ] — — 6.06 6.06 9.586 — — 6.94 4.80 9.192
通过对比表4 中的数据可知,相对于文献[14 ]中的维修策略,综合优化备件订购阈值和预防性更换阈值能够有效降低系统期望维修费用率.在案例1和案例2下系统最优维修及备件订购策略对应的系统期望维修费用率 γ ∞ * 分别降低了6.9%和5.9%,这是由于同时优化调整 Q A i 和 Q M i 值可以在降低停机费用的同时,增加更新周期长度,使得系统寿命得到更为充分的利用. ...
... 通过对比表4 中的数据可知,相对于文献[14 ]中的维修策略,综合优化备件订购阈值和预防性更换阈值能够有效降低系统期望维修费用率.在案例1和案例2下系统最优维修及备件订购策略对应的系统期望维修费用率 γ ∞ * 分别降低了6.9%和5.9%,这是由于同时优化调整 Q A i 和 Q M i 值可以在降低停机费用的同时,增加更新周期长度,使得系统寿命得到更为充分的利用. ...
1
2004
... 式中:C 为Copula函数;c 为Copula密度函数;x =[x 1 x 2 … xn ],xi ≥0; F X i ( t - u ) ( xi )为第i 个退化过程在[u ,t ]内增量的分布函数;θ =[θ 1 θ 2 … θm ]为Copula函数中的参数向量,其取值直接影响着随机变量间相关关系强弱.对于Copula函数参数估计问题,可基于退化数据通过最大似然估计或边际推断(IFM)法分两步获得.在此基础上,可采用赤池信息准则(AIC)选择合适的Copula函数描述系统退化相关性.关于Copula函数的参数估计方法参考文献[15 ]. ...
基于随机参数Gamma过程的剩余寿命预测方法
1
2015
... 由于Gamma函数和不完全Gamma函数的影响, F X i ( t ) ( Q L i ) 的计算较为复杂,可采用B-S分布对其进行近似计算[16 ] . F X i ( t ) ( Q L i ) 可近似表示为 ...
基于随机参数Gamma过程的剩余寿命预测方法
1
2015
... 由于Gamma函数和不完全Gamma函数的影响, F X i ( t ) ( Q L i ) 的计算较为复杂,可采用B-S分布对其进行近似计算[16 ] . F X i ( t ) ( Q L i ) 可近似表示为 ...
A comparative study of Artificial Bee Colony algorithm
1
2009
... 为解决这一问题,采用ABC 算法搜索系统最优维修与备件订购策略.ABC 算法是一种基于蜂群搜索蜜源行为的群体智能优化算法,在智能优化算法中属于比较新的算法,该算法控制参数少、易于实现、全局收敛性能好,并且其优化性能优于其他一些传统智能优化算法[17 ] ,如遗传算法(GA )、粒子群优化(PSO )算法、差分进化(DE )算法等,特别是在非线性优化函数求解方面具有良好的性能.ABC 算法已广泛应用于解决各类优化问题,关于ABC 算法的具体机理和优化过程可参考文献[18 ],这里不再赘述.在案例分析中,ABC 算法采用MATLAB 实现. ...
A powerful and efficient algorithm for numerical function optimization: Artificial bee colony (ABC) algorithm
1
2007
... 为解决这一问题,采用ABC 算法搜索系统最优维修与备件订购策略.ABC 算法是一种基于蜂群搜索蜜源行为的群体智能优化算法,在智能优化算法中属于比较新的算法,该算法控制参数少、易于实现、全局收敛性能好,并且其优化性能优于其他一些传统智能优化算法[17 ] ,如遗传算法(GA )、粒子群优化(PSO )算法、差分进化(DE )算法等,特别是在非线性优化函数求解方面具有良好的性能.ABC 算法已广泛应用于解决各类优化问题,关于ABC 算法的具体机理和优化过程可参考文献[18 ],这里不再赘述.在案例分析中,ABC 算法采用MATLAB 实现. ...