«上一篇
文章快速检索     高级检索
下一篇»
  哈尔滨工程大学学报  2017, Vol. 38 Issue (6): 950-955  DOI: 10.11990/jheu.201603079
0

引用本文  

周炳海, 刘玉旺. 考虑设备劣化及随机需求的最优生产周期模型[J]. 哈尔滨工程大学学报, 2017, 38(6), 950-955. DOI: 10.11990/jheu.201603079.
ZHOU Binghai, LIU Yuwang. An optimal production period model for a buffered system with machine deteriorating upstream and stochastic demand downstream[J]. Journal of Harbin Engineering University, 2017, 38(6), 950-955. DOI: 10.11990/jheu.201603079.

基金项目

国家自然科学基金项目(71471135)

通信作者

周炳海, E-mail:bhzhou@tongji.edu.cn

作者简介

周炳海(1965-), 男, 教授, 博士生导师

文章历史

收稿日期:2016-03-23
网络出版日期:2017-05-04
考虑设备劣化及随机需求的最优生产周期模型
周炳海, 刘玉旺    
同济大学 机械与能源工程学院, 上海 201804
摘要:针对生产过程中由于设备劣化引起产品缺陷率增大以及顾客需求波动问题,构建了综合考虑设备劣化和需求随机的最优生产周期模型。利用伽马过程对上游生产单元的劣化过程进行建模,并假设当其劣化量达到中间某一水平时,生产活动从“受控”转变为“失控”状态;以此为基础,在假设下游需求随机的条件下,以最小化单位产品的总成本为优化目标建立了数学模型;借助数值实例演示了本模型的有效性,并对相关参数作了敏感性分析。分析结果表明:本文所提模型对于描述设备劣化程度对产品质量的影响,降低企业运行成本(库存、返工、维护等)有一定的指导意义。
关键词设备劣化    伽马过程    随机需求    缓冲    优化    
An optimal production period model for a buffered system with machine deteriorating upstream and stochastic demand downstream
ZHOU Binghai, LIU Yuwang    
School of Mechanical and Energy Engineering, Tongji University, Shanghai 201804, China
Abstract: During production, machine deterioration and stochastic demand can cause an increased probability of producing defective items and a shortage of WIPs, respectively. Hence, an optimal production period model considering machine deterioration upstream and stochastic demand downstream was formulated. First, a gamma process was introduced to model the degradation of the upstream machine; after the degradation state exceeded a predetermined threshold, the production state transformed from "in-control" to "out-of-control." Based on the abovementioned descriptions, a mathematical model with an objective function for minimizing the total expected cost per unit was derived, assuming the downstream demand rate was stochastic. Finally, a numerical example was provided to illustrate the application of the proposed model, and sensitivity analysis of the model with respect to some key parameters was performed. Results show that the proposed model has certain instructive significance in quality control and for cost reduction of imperfect production systems.
Key words: machine    deterioration    Gamma    process    stochastic    demand    buffer    optimization    

在实际制造系统中,随着生产设备的运行,其劣化程度逐渐加剧,使得加工产品的不合格率也随之增大[1-3]。此外,顾客的需求模式逐渐从以往的“大批量、单一化”向“多品种、小批量”转变,这一转变对如何保障制造过程的平稳性提出了新的挑战。合理地决策设备的运行周期以及缓冲区的容量是解决上述两个问题的有效途径。

许多文献将生产系统简化为一个两设备一缓冲系统,以便将设备的劣化过程、维护策略以及缓冲区的容量控制同时纳入一个数学模型中综合考虑。Chelbi等假设对上游设备进行定期预防性维护,建立了缓冲区库存S和预防性维护间隔T的联合优化模型[4]。在此基础上,Rezg等又考虑一个生产周期内缓冲开始建立的时刻A,将模型扩展为(T, S, A)的三维联合优化模型[5]。成国庆等利用几何过程来刻画上游生产设备的退化特征,研究了修复非新条件下的带缓冲区的生产系统维修更换策略[6]。Karamatsoukis等假设上游设备有m个工作状态,1个失效状态,为其最优的预防性维护建立了离散时间的Markov决策模型,即所谓的“控制限”策略[7-8]。之后Karamatsoukis等又将该模型扩展至上下游设备均随时间逐渐劣化并均可能发生失效的情形[9]。Sana[1]考虑到生产过程的不完美性(即存在不合格品),在假设对系统进行定期预防性维护的条件下确定了使单位产品总成本最小的最优库存水平[3]。Radhoui同样考虑一个不完美生产系统,为其质量控制和预防性维护作了联合优化[10]。文献[11-12]则是考虑一类存放时随时间逐渐劣化的产品,对经典的EPQ模型进行了扩充,两者的不同在于分别假设维护时间为定值和随机量。

上述文献虽然都考虑到了上游设备的劣化、失效及对应的维护策略(TBM[13]或CBM[14]),但它们大都假设生产出的产品均是完美的(文献[3, 9]除外),并未将设备的劣化状态与产品质量两者联系起来,并且在对缓冲区容量进行优化时都假设下游的需求速率恒定。本文以生产-缓冲-装配系统为研究对象,假设随着上游生产单元劣化量的增大,产品的不合格率也随之增大,并且下游装配单元的需求速率是一随机量,以单位产品总成本最小为目标函数,建立并求解了数学模型。

1 问题描述及模型构建 1.1 问题域描述

考虑一个由生产单元(P)、装配单元(A)和中间缓冲区(B)组成的系统,生产单元生产出的产品配件存放于缓冲区,装配线以随机的需求速率X从缓冲区取用配件进行装配(如图 1)。

图 1 生产-缓冲-装配系统示意图 Fig.1 Structure of the production-buffer-assembly system

为进一步描述研究问题域,作如下假设:

1) 在t=0时刻,生产单元以速率p开始生产,装配单元以随机的需求速率X(p0Xp)取用配件,pX,则缓冲区库存以速率pX进行累积;经过生产周期T1之后,生产单元停止生产,装配单元仍以速率X消耗配件,直至时刻T2缓冲耗尽。

2) 随着生产活动的进行,生产单元会发生劣化,假设其劣化过程为连续过程。当tτ时,生产活动处于受控状态,产品不合格率为a;当τttf时,生产活动处于失控状态,产品不合格率为b,且ba(如图 2)。

图 2 生产单元的劣化过程及生产过程的状态转变示意图 Fig.2 Schematic illustration of the degradation process and the transition of production state

3) 对生产单元进行连续监测,在生产过程中其不会发生失效(即生产过程中不会发生X(t)>DF)。

4) 在正常生产周期T1结束后,对不合格品统一进行返工。

5) 存放在缓冲区的产品不发生劣化或报废。

1.2 生产过程的不完美性建模

考虑到劣化过程的不可逆性,而伽马过程有非负的独立增量特性,因此我们选用定态伽马过程(形状参数α(t)=αt,尺度参数β)来刻画设备的劣化量随时间逐渐增大这一特性。设X(t)表示t时刻设备的劣化量,则其概率密度函数可表示为[15]

$Ga\left( {x\left| {\alpha t,\beta } \right.} \right) = \frac{1}{{{\beta ^{\alpha t}}\mathit{\Gamma }\left( {\alpha t} \right)}}{x^{\alpha t - 1}}{{\rm{e}}^{ - x/\beta }}$ (1)

其中,$\mathit{\Gamma }\left( \alpha \right) = \int_0^\infty {{z^{a - 1}}{{\rm{e}}^{ - z}}{\rm{d}}z\left( {\alpha > 0} \right)} $表示α的伽马函数。如图 2,设备在时刻τ劣化量达到状态转变阈值DP,则τ的累积分布函数和概率密度函数分别为

$\begin{array}{*{20}{c}} {{F_\tau }\left( t \right) = \Pr \left[ {\tau \le t} \right] = \Pr \left[ {X\left( t \right) \ge {D_p}} \right] = }\\ {\int_{{D_p}}^\infty {G{a_{X\left( t \right)}}\left( x \right){\rm{d}}x} = \frac{{\mathit{\Gamma }\left( {\alpha t,{D_p}/\beta } \right)}}{{\mathit{\Gamma }\left( {\alpha t} \right)}}} \end{array}$ (2)
$\begin{array}{*{20}{c}} {{f_\tau }\left( t \right) = \frac{\partial }{{\partial t}}\left[ {\frac{{\mathit{\Gamma }\left( {\alpha t,{D_P}/\beta } \right)}}{{\mathit{\Gamma }\left( {\alpha t} \right)}}} \right] = \frac{\partial }{{\partial \tilde \alpha }}\left[ {\frac{{\mathit{\Gamma }\left( {\bar \alpha ,{D_P}/\beta } \right)}}{{\mathit{\Gamma }\left( {\bar \alpha } \right)}}} \right]\alpha '\left( t \right) = }\\ {\frac{\alpha }{{\mathit{\Gamma }\left( {\alpha t} \right)}}\int_{{D_P}/\beta }^\infty {\left\{ {\log z - \psi \left( {\alpha t} \right)} \right\}{z^{\alpha t - 1}}{{\rm{e}}^{ - z}}{\rm{d}}z} } \end{array}$ (3)

其中,α=αtψ(α)是伽马函数对数的导数:

$\psi \left( \alpha \right) = \frac{{\mathit{\Gamma '}\left( \alpha \right)}}{{\mathit{\Gamma }\left( \alpha \right)}} = \frac{{\partial \log \mathit{\Gamma }\left( \alpha \right)}}{{\partial \alpha }},\alpha > 0$ (4)
1.3 随机需求条件下的成本分析 1.3.1 库存持有成本

根据生产周期长度与生产单元状态转变时刻的相对大小关系,可分为以下两种情况。

情形1τT1,生产过程全部处于受控状态,产品不合格率为a,生产结束后对设备进行预防性维护,并假设修复如新。任一时刻合格品的库存水平可通过下式表达:

$\begin{array}{l} \frac{{{\rm{d}}{I_{g11}}\left( t \right)}}{{{\rm{d}}t}} = \left( {1 - a} \right)p - X,\\ 0 < t \le {T_1}\;且\;{I_{g11}}\left( 0 \right) = 0 \end{array}$
$\begin{array}{l} \frac{{{\rm{d}}{I_{g12}}\left( t \right)}}{{{\rm{d}}t}} = - X,{T_1} \le t \le {T_2}\;且\\ {I_{g12}}\left( {{T_1}} \right) = \left[ {\left( {1 - a} \right)p - X} \right]{T_1},{I_{g12}}\left( {{T_2}} \right) = 0 \end{array}$

利用上述边界条件,可以求得

${I_{g11}}\left( t \right) = \left[ {\left( {1 - a} \right)p - X} \right]t,0 \le t \le {T_1}$ (5)
${I_{g12}}\left( t \right) = - Xt + \left( {1 - a} \right)p{T_1},{T_1} \le t \le {T_2}$ (6)

Ig12(T2)=0,可以得到

${T_2} = \frac{{\left( {1 - a} \right)p{T_1}}}{X}$ (7)

则合格品的库存持有成本为

$C_{H1}^G = c_h^g\left[ {\int_0^{{T_1}} {{I_{g11}}\left( t \right)dt} + \int_{{T_1}}^{{T_2}} {{I_{g12}}\left( t \right){\rm{d}}t} } \right]$ (8)

化简得

$C_{H1}^G = c_h^g\left[ { - \frac{{pT_1^2}}{2} + \frac{{apT_1^2}}{2} + \frac{{{p^2}T_1^2}}{{2X}} - \frac{{a{p^2}T_1^2}}{X} + \frac{{{a^2}{p^2}T_1^2}}{{2X}}} \right]$

不合格品的库存水平可通过下式表达:

${I_{r11}}\left( t \right) = apt,0 \le t \le {T_1}$ (9)

则不合格品的库存持有成本为

$C_{H1}^G = c_h^r\int_0^{{T_1}} {{I_{r11}}\left( t \right){\rm{d}}t} = c_h^r\frac{{apT_1^2}}{2}$ (10)

情形2τT1tf,生产过程先是处于受控状态,产品不合格率为a;在时刻τ之后转入失控状态,不合格率为b(ba)。在生产活动结束后,通过对设备进行预防性维护使其恢复至初始状态。任一时刻合格品的库存水平可由下式表达:

$\frac{{{\rm{d}}{I_{g21}}\left( t \right)}}{{{\rm{d}}t}} = \left( {1 - a} \right)p - X,0 < t \le \tau \;且\;{I_{g21}}\left( 0 \right) = 0$
$\begin{array}{l} \frac{{{\rm{d}}{I_{g22}}\left( t \right)}}{{{\rm{d}}t}} = \left( {1 - b} \right)p - X,\tau \le t \le {T_1}\;且\\ {I_{g22}}\left( \tau \right){\rm{ = }}\left[ {\left( {1 - a} \right)p - X} \right]\tau \end{array}$
$\begin{array}{l} \frac{{{\rm{d}}{I_{g23}}\left( t \right)}}{{{\rm{d}}t}} = - X,{T_1} \le t \le {T_2}\;且\\ {I_{g23}}\left( {{T_1}} \right) = {I_{g22}}\left( {{T_1}} \right),{I_{g23}}\left( {{T_2}} \right) = 0 \end{array}$

利用上述边界条件,可以求得

${I_{g21}}\left( t \right) = \left[ {\left( {1 - a} \right)p - X} \right]t,\;\;\;0 \le t \le \tau $ (11)
$\begin{array}{*{20}{c}} {{I_{g22}}\left( t \right) = \left[ {\left( {1 - b} \right)p - X} \right]t + \left( {b - a} \right)p\tau ,}\\ {\tau \le t \le {T_1}} \end{array}$ (12)
$\begin{array}{*{20}{c}} {{I_{g23}}\left( t \right) = - Xt + \left( {1 - b} \right)p{T_1} + \left( {b - a} \right)p\tau ,}\\ {{T_1} \le t \le {T_2}} \end{array}$ (13)

Ig23(T2)=0,可以得到

${T_2} = \frac{{\left( {1 - b} \right)p{T_1} + \left( {b - a} \right)p\tau }}{X}$ (14)

则合格品的库存持有成本为

$C_{H2}^G = c_h^g\left[ {\int_0^\tau {{I_{g21}}\left( t \right){\rm{d}}t} + \int_\tau ^{{T_1}} {{I_{g22}}\left( t \right){\rm{d}}t} + \int_{{T_1}}^{{T_2}} {{I_{g23}}\left( t \right){\rm{d}}t} } \right]$ (15)

化简得

$\begin{array}{l} C_{H2}^G = c_h^g\left( { - \frac{{pT_1^2}}{2} + \frac{{bpT_1^2}}{2} + \frac{{{p^2}T_1^2}}{{2X}} - \frac{{b{p^2}T_1^2}}{X} + } \right.\\ \;\;\;\;\;\;\;\;\;\;\frac{{{b^2}{p^2}T_1^2}}{{2X}} - \frac{{a{p^2}{T_1}\tau }}{X} + \frac{{b{p^2}{T_1}\tau }}{X} + \frac{{ab{p^2}{T_1}\tau }}{X} - \\ \;\;\;\;\;\;\;\;\;\;\frac{{{b^2}{p^2}{T_1}\tau }}{{2X}} + \frac{{ap{\tau ^2}}}{2} - \frac{{bp{\tau ^2}}}{2} + \frac{{{a^2}{p^2}{\tau ^2}}}{{2X}} - \\ \;\;\;\;\;\;\;\;\;\;\left. {\frac{{ab{p^2}{\tau ^2}}}{X} + \frac{{{b^2}{p^2}{\tau ^2}}}{{2X}}} \right) \end{array}$

不合格品的库存水平可通过下式表达:

$\frac{{{\rm{d}}{I_{r21}}\left( t \right)}}{{{\rm{d}}t}} = ap,0 < t \le \tau \;且\;{I_{r21}}\left( 0 \right) = 0$
$\frac{{{\rm{d}}{I_{r22}}\left( t \right)}}{{{\rm{d}}t}} = bp,\tau \le t \le {T_1}\;且\;{I_{r22}}\left( \tau \right) = ap\tau $

利用边界条件,可求得

${I_{r21}}\left( t \right) = apt,0 < t \le \tau $ (16)
${I_{r22}}\left( t \right){\rm{ = }}bpt + \left( {a - b} \right)p\tau ,\tau \le t \le {T_1}$ (17)

则不合格品的库存持有成本为

$C_{H2}^R = c_h^r\left[ {\int_0^\tau {{I_{r21}}\left( t \right){\rm{d}}t} + \int_\tau ^{{T_1}} {{I_{r22}}\left( t \right){\rm{d}}t} } \right]$ (18)

化简得

$C_{H2}^R = c_h^r\left( { - \frac{{ap{\tau ^2}}}{2} + \frac{{bp{\tau ^2}}}{2} + ap{T_1}\tau - bp{T_1}\tau + \frac{{bpT_1^2}}{2}} \right)$

综合以上两种情形,合格品和缺陷品的总期望库存持有成本分别为

$C_H^G = \int_{{p_0}}^p {\left[ {\int_{{T_1}}^\infty {C_{H1}^Gf\left( \tau \right){\rm{d}}\tau } + \int_0^{{T_1}} {C_{H2}^Gf\left( \tau \right){\rm{d}}\tau } } \right]g\left( x \right){\rm{d}}x} $ (19)
$C_H^R = \int_{{p_0}}^p {\left[ {\int_{{T_1}}^\infty {C_{H1}^Rf\left( \tau \right){\rm{d}}\tau } + \int_0^{{T_1}} {C_{H2}^Rf\left( \tau \right){\rm{d}}\tau } } \right]g\left( x \right){\rm{d}}x} $ (20)
1.3.2 其他成本

两种情形分别对应的不合格品的总数目:

$N = \left\{ \begin{array}{l} ap{T_1},\tau \ge {T_1}\\ ap\tau + bp\left( {{T_1} - \tau } \right),\tau < {T_1} \end{array} \right.$ (21)

其期望值

$\begin{array}{l} E\left( N \right) = ap{T_1}\int_{{T_1}}^\infty {f\left( \tau \right){\rm{d}}\tau } + ap\int_0^{{T_1}} {\tau f\left( \tau \right){\rm{d}}\tau } + \\ \;\;\;\;\;\;\;\;\;\;\;\;bp\int_0^{{T_1}} {\left( {{T_1} - \tau } \right)f\left( \tau \right){\rm{d}}\tau } = \\ \;\;\;\;\;\;\;\;\;\;\;\;ap{T_1} + ap\int_0^{{T_1}} {\left( {\tau - {T_1}} \right)f\left( \tau \right){\rm{d}}\tau } + \\ \;\;\;\;\;\;\;\;\;\;\;\;bp\int_0^{{T_1}} {\left( {{T_1} - \tau } \right)f\left( \tau \right){\rm{d}}\tau } = \\ \;\;\;\;\;\;\;\;\;\;\;\;ap{T_1} + \left( {b - a} \right)p\int_0^{{T_1}} {\left( {{T_1} - \tau } \right)f\left( \tau \right){\rm{d}}\tau } \end{array}$ (22)

于是返工成本为

${C_{RE}} = {c_{re}}\left[ {\alpha p{T_1} + \left( {\beta - \alpha } \right)p\int_0^{{T_1}} {\left( {{T_1} - \tau } \right)f\left( \tau \right){\rm{d}}\tau } } \right]$ (23)

一个生产周期内的生产成本:

${C_P} = {c_p}p{T_1}$ (24)

随着生产单元劣化程度的增大,维护所需付出的人力物力成本也随之增大,此处假设对生产单元的预防维护成本与其劣化量成正比(CMX(T1)),即下列等式成立:

$\begin{array}{l} {C_M} = {\gamma _0}E\left[ {X\left( {{T_1}} \right)} \right] = {\gamma _0}\int_0^\infty {xGa\left( {\left| {\alpha {T_1},\beta } \right.} \right){\rm{d}}x} = \\ \;\;\;\;\;\;\;\;{\gamma _0}\int_0^\infty {x\frac{1}{{{\beta ^{\alpha T}}\Gamma \left( {\alpha {T_1}} \right)}}{x^{\alpha {T_1} - 1}}{{\rm{e}}^{ - x/\beta }}{\rm{d}}x} \end{array}$ (25)

式中:γ0为一常数,γ>0。

ω(ξ)表示缺货时间,则

$\omega \left( \xi \right) = \left\{ \begin{array}{l} 0,\;\;\;\;\xi \le {T_d}\\ \xi - {T_d},\;\;\;\xi > {T_d} \end{array} \right.$

其中缓冲耗尽所需时间Td可通过下式表达:

${T_d} = E\left[ {\frac{{\left( {p - X} \right){T_1}}}{X}} \right] = \int_{{p_0}}^p {\left( {\frac{p}{x} - 1} \right){T_1}{g_x}\left( x \right){\rm{d}}x} $ (26)

考虑到维护时间会随着生产单元劣化量的增大而增大,假设ξ服从参数为$\left( {\lambda \propto \frac{1}{{E\left[ {X\left( {{T_1}} \right)} \right]}}} \right)$的指数分布,则ξ的概率密度函数为

$\varphi \left( \xi \right) = \lambda {{\rm{e}}^{ - \lambda \xi }}$ (27)
$\lambda = \frac{\gamma }{{E\left[ {X\left( {{T_1}} \right)} \right]}} = \frac{{{\gamma _1}}}{{\int_0^\infty {x\frac{1}{{{\beta ^{\alpha T}}\mathit{\Gamma }\left( {\alpha {T_1}} \right)}}{x^{\alpha {T_1}}}{{\rm{e}}^{ - x/\beta }}{\rm{d}}x} }}$ (28)

式中:γ1为一常数,γ1>0。

于是,缺货成本的期望值为

${C_s} = {c_s}E\left( {\xi - {T_d}} \right) = {c_s}\int_{{T_d}}^\infty {\left( {\xi - {T_d}} \right)\varphi \left( \xi \right){\rm{d}}\xi } $ (29)

设生产准备成本为K,由式(19)~(20),(23)~(25) 及(29),可以得到一个生产周期内单位产品的总成本为

${\rm{CR = }}\frac{{K + {C_p} + C_H^G + C_H^R + {C_{RE}} + {C_M} + {C_s}}}{{P{T_1}}}$ (30)
2 数值实例

本节通过一个数值实例来对所提模型进行演示,模型中各参数的设置见表 1。其中,伽马过程的参数αβ采用最大似然法进行估计[16],其他参数参考实际生产中的相关数据确定。

表 1 模型中的各项参数 Tab.1 Parameters of the machines

首先,假设下游装配单元的需求速率分别服从区间[1, 5]的均匀分布和参数为0.5的指数分布,分别计算这两种情形下产品的总成本率随生产周期T1的变化而发生的变化情况,如图 3

图 3 产品的总成本率CR随生产周期T1的变化 Fig.3 Production period versus totatl expected costs per unit

图 3可知,若下游需求服从区间[1, 5]的均匀分布,则最优生产周期T1*=26,此时单位产品的总成本最小,为57.127;若下游需求服从参数为0.5的指数分布,可求得最优生产周期T1*=20,对应的最小单位产品总成本为58.05。

下面对模型中的相关参数作敏感性分析,旨在分析这些参数的变化对上述各项成本及最优生产周期所造成的影响(以X~U[1, 5]为例)。

Dp的变化对T1*及各项成本的影响见表 2

表 2 Dp的变化对T1*及各项成本的影响 Tab.2 Threshold of state transition versus optimal production period and various costs

表 2可以得到以下结论:

1) 最优生产周期T1*Dp的变化并不敏感。当Dp在20上下以2个单位递减或递增时,T1*始终保持为26;当Dp增大到40时,T1*才增大为28。这表明在Dp的增幅较大时,T1*会随着Dp的增大而增大;

2) 当T1*=26时,维护成本CM和缺货成本CS的值均保持不变,这是因为假设维护时间与生产单元的劣化量成正比,根据伽马过程的性质,在时刻T1,生产单元的期望劣化量E[X(T1)]=αT1,于是可以得到:维护时间与生产周期T1也成正比。同理可证:在T1不变时,缺货成本CS也是一定值。

3) 当T1的值固定时,随着Dp的增大,合格品的总库存持有成本CHG逐渐增大,非合格品的总库存持有成本CRH和返工成本CRE逐渐减小,总成本率逐渐减小。原因在于:Dp的大小决定了生产单元由“受控”状态转变到“失控”状态的快慢(后者的不合格品率更高)。Dp值越小,则生产单元就越快地转变为“失控”状态,在同等的生产周期T1下生产出的不合格品数目越多,导致不合格品的库存成本和返工成本越大。

4) 对比表 2中的最后两行,可以看到:Dp=40对应的不合格品库存成本大于Dp=30对应的不合格品库存成本,但前者的返工成本却小于后者。这表明前者的不合格品数实际上小于后者,而前者的不合格品库存成本较大的原因在于其对应的生产周期T1=28>26(库存成本不仅取决于库存数量,还取决于持有时间)。

K的变化对T1*和CR的影响见表 3,由表 3可知:最优生产周期T1*对生产准备K的变化较为敏感。随着K的增大,最优生产周期T1*和单位产品总成本CR均随之增大。

表 3 K的变化对T1*和CR的影响 Tab.3 Set-up costs versus optimal production period and expected costs per unit

缺货时间随T1的变化情况如图 4

图 4 缺货时间随T1的变化情况 Fig.4 production period versus shortage time

图 4可以看到:随着生产周期T1的增大,缺货时间逐渐增大,两者大致呈线性关系。这是因为尽管T1的增大使积累的合格品库存变大,在下游需求模式未发生改变的情况下可维持装配单元继续工作的时间越长,但是由于生产周期变长引起的维护时间的增大量更多,导致两者的差值(即缺货时间)越来越大。

3 结论

1) 以一个两单元一缓冲系统为研究对象,假设上游生产单元的劣化过程服从定态伽马过程,当其劣化量达到一定水平时,生产过程由“受控”转变为“失控”状态(后者的缺陷品率较大),这一假设更好地符合实际生产情形,为生产企业及时开展设备维护、产品质量控制提供了一种新的思路;

2) 模型假设下游装配单元的的需求速率不恒定,为一随机变量,也很好契合了实际中客户需求难以预测这一特点。

3) 在以上两条假设的基础上,以生产周期长度为决策变量,构建了最小化单位产品总成本的最优生产周期模型,并通过数值实例验证了所提模型的有效性。分析结果表明:本文优化模型对描述产品不合格率随生产设备劣化量的增大而增大、以及需求随机条件下的库存、维护等成本分析是切实有效的,对于企业合理安排生产计划、降低运行成本具有一定指导意义。

关于本文今后的研究工作有以下几个方面:假设在每一生产周期结束后对生产单元进行的预防维护只能“修复非新”,决策变量可以增加为对生产周期和生产速率的联合优化等。

参考文献
[1] 周炳海, 刘子龙. 考虑质量损失的退化系统维护建模[J]. 浙江大学学报:工学版, 2016, 50(12): 2270-2276.
ZHOU Binghai, LIU Zilong. Maintenance modelling for deteriorating system considering quality loss[J]. Journal of Zhejiang University:engineering science, 2016, 50(12): 2270-2276. (0)
[2] ZHOU B H, LIU Z L. Optimizing preventive maintenance:a deteriorating system with buffers[J]. Industrial management & data systems, 2016, 116(8): 1719-1740. (0)
[3] SANA S S. Preventive maintenance and optimal bufferinventory for products sold with warranty in an impe-rfect production system[J]. International Journal of production research, 2012, 50(23): 6763-6774. DOI:10.1080/00207543.2011.623838 (0)
[4] CHELBI A, AIT-KADI D. Joint optimal buffer inventory and preventive maintenance strategy for a randomly failing production unit[J]. Journal of decision systems, 2003, 12(1): 21-30. DOI:10.3166/jds.12.21-30 (0)
[5] REZG N, DELLAGI S, CHELBI A. Joint optimal inventory control and preventive maintenance policy[J]. International Journal of Production Research, 2008, 46(19): 5349-5365. DOI:10.1080/00207540802273942 (0)
[6] 成国庆, 周炳海, 李玲. 考虑缓冲区库存的退化系统最优维修更换策略[J]. 计算机集成制造系统, 2015, 21(6): 1593-1600.
CHENG Guoqing, ZHOU Binghai, LI Ling. Optimal maintenance strategy for deteriorating systems with intermediate buffers[J]. Computer integrated manufacturing systems, 2015, 21(6): 1593-1600. (0)
[7] KARAMATSOUKIS C C, KYRIAKIDIS E G. Optimal maintenance of a production inventory system with-idle periods[J]. European journal of operational research, 2009, 196: 744-751. DOI:10.1016/j.ejor.2008.04.012 (0)
[8] PAVITSOS A, KYRIAKIDIS E G. Markov dicision models for the optimal maintenance of a production unitwith an upstream buffer[J]. Computers & Operations Research, 2009, 36: 1993-2006. (0)
[9] KARAMATSOUKIS C C, KYRIAKIDIS E G. Optimal maintenance of two stochastically deteriorating machines with an intermediate buffer[J]. European journal of operational research, 2010, 207: 297-308. DOI:10.1016/j.ejor.2010.04.022 (0)
[10] RADHOUI M, REZG N, CHELBI A. Joint quality control and preventive maintenance strategy for imperfect production processes[J]. Journal of intelligent manufacturing, 2010, 21(2): 205-212. DOI:10.1007/s10845-008-0198-x (0)
[11] LIN G C, GONG D C. On a production-inventory system of deteriorating items subject to random machinebreakdowns with a fixed repair time[J]. Mathematical and computer modelling, 2006, 43: 920-932. DOI:10.1016/j.mcm.2005.12.013 (0)
[12] WIDYADANA G A, WEE H M. Optimal deterioratingitems production inventory models with random machine breakdown and stochastic repair time[J]. Applied mathematical modelling, 2011, 35: 3495-3508. DOI:10.1016/j.apm.2011.01.006 (0)
[13] CHEN W J. Minimizing total flow time in the singlemachine scheduling problem with periodic maintenance[J]. Journal of the operational research society, 2006, 57(4): 410-415. DOI:10.1057/palgrave.jors.2601998 (0)
[14] ZEQUEIRA R I, BERENGUER C. Periodic imperfect preventive maintenance with two categories of competing failure modes[J]. Reliability engineering & system safety, 2006, 91(4): 460-468. (0)
[15] 陶红玉, 周炳海. 基于随机退化的串行生产系统机会维护模型[J]. 上海交通大学学报, 2013, 47(12): 1911-1917.
TAO Hongyu, ZHOU Binghai. Opportunistic maintenance model for series production systems based on stochastic degradations[J]. Journal of Shanghai Jiao Tong University, 2013, 47(12): 1911-1917. (0)
[16] 翟子青. 基于随机过程的蒸汽发生器传热管腐蚀失效寿命分析[D]. 上海: 上海交通大学, 2011: 34-38.
ZHAI Zi Qing. Stochastic modeling for the degradation of steam generator tubes due to pitting corrosion[D]. Shanghai:Shanghai Jiao Tong University, 2011:34-38. http://cdmd.cnki.com.cn/Article/CDMD-10248-1011268171.htm (0)