Joint optimization of condition-based maintenance and inventory ordering based on status monitoring for multi-unit system

<!-- metadata commented in wiki content


<big>'''Joint Optimization of Condition-based Maintenance and Inventory Ordering Based on Status Monitoring for Multi-unit system'''</big>

HUAYANG DENG [1], QUAN SHI [1] AND YADONG WANG [1]

<span id='_Hlk63188127'></span><span id='_Hlk63188103'></span><span id='_Hlk63188144'></span><span id='_Hlk63188231'></span><sup>1 </sup>''Deparent of Equipment Command and Management, Army Engineering University of PLA, Shijiazhuang 050003, China''
-->
==Abstract==

In the modern industry, in order to reduce the inventory pressure, a variety of parts began to use unified kind spare parts for maintenance. However, highly integrated equipment is more difficult to use traditional RCM models, and researchers begin to steering based on state monitoring methods. Deepen a prediction of equipment failure. This paper mainly discussed the data-driven analysis method based on the Wiener process to predict the fault law of the same type. The joint model innovatively adopts the (''s''-1, ''s'') policy considering the industrial characteristic and multi-period resupply. In the end, we analyze (''s''-1, ''s'') policy in joint optimization by comparison to draw the optimal policy combination.

'''Keywords: '''Multi-unit System; Condition-based maintenance; Spare parts; Joint optimization; Wiener process; Resupply.

==1. Introduction==

In industrial, military, and other fields, with the emergence of high-precision, high-information, and highly integrated equipment, the traditional fault prediction based on RCM using classical fault law experience shows some limitations in the more complex systematic analysis [1]. With the development of sensor technology, many researchers began to establish equipment life degradation model based on data-driven analysis, so as to predict the law of equipment failure [2].

The main methods in estimating life degradation based on data drive are: Gamma process, Gaussian inverse process, and Wiener process. After the exploration of many researchers, the data-driven model can be adapted to more and more modeling environments and can be used flexibly by adjusting formula parameters and expanding formulas [3]. So far, it has been widely used in the field of reliability [4-6]. Compared with the characteristics of the three models, the Gamma process and Gaussian inverse process are only suitable for modeling monotone degradation paths. In contrast, the Wiener process is more suitable for non-monotonic degradation processes caused by repair, self-healing, or use intensity [7]. It is also an important prerequisite for choosing a proper inventory policy [8].

In the research of inventory policy, the application amount of inventory policies is very high [9]. As two kinds of the most basic and classical inventory policies, they are still widely used in the modern supply chain. From the optimization goal, the (<math> T </math>, <math> S </math>) policy focuses on optimizing the inspection cycle <math> T </math>, while the (<math> s </math>, <math> S </math>)’s is the number of inventory resupply (<math> s </math>). Condition-based maintenance (CBM) is the more complicated preventive maintenance policy. The (<math> s </math>, <math> S </math>) policy which pays more attention to the change of inventory is closer along with the characteristics of CBM which based on the state of parts [10]. However, due to the high importance of the equipment using the data-driven method, the inventory cost and the loss caused by failure are also high, so in industrial practice, the factories which have this kind of equipment mostly adopt the (<math> s -1</math>, <math> s </math>) policy. After a part is replaced, the resupply application will be sent to replenish the inventory level.

Secondly, most of the delay time of inventory is fixed in all resupply periods. This fixed model is feasible in a certain time range. But for some industrial processes with multiple periods and a wide range of time in each period, the delay time also fluctuates with the influence of the transportation environment, local policy, or human factors [11]. In view of this characteristic, this article models the delay time separately and chooses the Poisson distribution which is generally accepted in the current research field to discuss [12] That is why simulating the multi-period resupply process.

The rest of this paper is organized as follows. Section 2 is the fault analysis method based on Wiener process and the research on the joint model between maintenance and inventory ordering. Then, section 3 mainly introduces performance degradation, life prediction by Wiener process, and CBM policy. In section 4, according to the conditions and constraints, combined with the (<math> s -1</math>, <math> s </math>) policy, the maintenance-inventory joint policy under the same type of multi-component condition monitoring is modeled. Besides that, the model of multi-period resupply time is also analyzed. Section 5 adopts Monte Carlo simulation and improved particle swarm optimization (PSO) algorithm to optimize the joint model and the optimal decision combination by balancing the model parameters. In section 6, a numerical example is given. While the sensitivity analysis of the resupply parameter, the corresponding comparison between the two kinds of inventory policies shows the difference. Conclusions and future work are given in section 7.

==2. Literature review==

===2.1 Performance degradation process based on Wiener process===

The Wiener process is one of the commonly used methods for modeling product performance degradation process in reliability assessment based on monitoring data [7]. The degradation process is a random process, driven by BM and Brownian motion. It is usually used to establish a model in the nonlinear degradation process.

In the current research, this kind of model has gradually developed to general and nonlinear problems related to time or state [13]. Baraldi, Mangili, Zio ''et al.'' sort out the diversity of model research, including time diversity, unit diversity, measurement diversity, and structural diversity [14]. Another point of the study is to consider the impact of the working environment on the performance degradation process. Bian, Gebraeel, Bian, Gebraeel ''et al.'' add covariates in the Wiener process as the environmental factors in the formula [15].

The model based on the Wiener process has also been applied to the field of spare parts. Elwany, Gebraeel ''et al''. firstly combine the maintenance decision based on the Wiener process with the spare parts ordering to evaluate the optimal replacement and ordering time. They finally put forward the general framework of the joint policy as well [16]. After that, Zhang, Chen ''et al.'' tried to introduce the Bayesian model into the joint decision model, and update the model parameters with real-time performance data [17]. On the other hand, Zhang, Si ''et al.'' extend the decision-making process to propose a random order lead time for inventory [18].

===2.2 Inventory policy===

Maintenance policy leads to inventory policy [19]. Talking about all inventory policies, according to the optimization goal, it can be divided into two main aspects. One is mainly represented by the inventory policies whose research goal is check cycle (<math> T </math>), such as (<math> T </math>, <math> S </math>) and (<math> T </math>, <math> Q </math>) [20]. However, in the long-term research, experts found it has limitation in applying due to the uncertainty of actual demand [21]. Therefore, the rigid and fixed inventory policies are difficult to meet the needs of multiple maintenance policies. Researchers began to discuss another kind of inventory policy represented mainly by (<math> s </math>, <math> S </math>) and (<math> s </math>, <math> Q </math>) policy [10,22,23]. The corresponding situation of parts leads to the change of inventory. In this view, the current number of inventory (<math> s </math>) is the focus of optimization.

However, there is a reality that with the development of technology, especially sensors, the reliability of parts has been gradually changed. The coin has two sides. A large number of spare parts are piled up in the warehouse, so it will result in a substantial increase in the proportion of inventory cost to the total cost. Lots of spare parts in some smaller warehouses cannot be transferred. What's more, there could be a large number of faulty spare parts and losses resulting from it. Considering the developed Internet and logistics industry, the researchers put forward the (<math> s -1</math>, <math> s </math>) policy. This policy aims to resupplied spare parts. After the part is replaced, the resupply information will be sent out. Sherbrooke firstly developed the policy into METRIC model under the limited cost and minimum out-of-stock [24]. After that, Grave proposed a simple algorithm for this research, which can be used in calculation [25].

===2.3 Joint optimization===

In the joint optimization of maintenance and inventory policy, considering the diversity and complexity, the previous literature focuses on traditional failure laws, such as Exponential distribution, Poisson distribution, and so on [26]. In recent years, data driven fault analysis was gradually developed. The joint optimization research based on that has been begun to discuss, especially the Wiener process.

Wang [27] applied the Wiener process to the joint optimization of CBM and inventory. Using the method of remaining life predicted by monitoring data, he analyzed the optimal policy for single part system [28,29]. With the deepening of the research on the single-part model, the research on multiple parts has been paid more and more attention. Minou and Nguyen [30] combined the Markov process and genetic algorithm with multi-period joint optimization for multiple parts, which is based on the (<math> s </math>, <math> S </math>) policy. Reviewing the previous research, there are few studies in the joint optimization of CBM based on the Wiener process and the (<math> s -1</math>, <math> s </math>) policy. Apart from that, most of the joint optimization studies focus on the parameters of maintenance policy, but there are few studies on the influence of inventory policy parameters, especially the time model of multi-period inventory delay.

==3. Maintenance policy based on Wiener process==

===3.1 Notations===

[[#tab-1|Table 1]] shows a set the notations.

<div class="center" style="font-size: 75%;">'''Table 1'''. Notations list</div>

<div id='tab-1'></div>
{| class="wikitable" style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:auto;" 
|-style="text-align:center"
! Symbol !! Meaning !! Unit
|-
|  style="text-align: center;"|<math>k</math>
|  style="text-align: center;"|Index of customers  (<math>k = 1,2,3\cdots K</math>)
|  style="text-align: center;"|
|-
|  style="text-align: center;"|<math>\tau </math>
|  style="text-align: center;"|Index of check cycles  (<math>\tau = 1,2,3\cdots \Psi </math>)Ψ
|  style="text-align: center;"|
|-
|  style="text-align: center;"|<math>\tau_i</math>
|  style="text-align: center;"|Index of resupply periods (<math>\tau_i = 1,2,3\cdots \Phi</math>)
|  style="text-align: center;"|
|-
|  style="text-align: center;"|<math>L</math>
|  style="text-align: center;"|Failure threshold
|  style="text-align: center;"|
|-
|  style="text-align: center;"|<math>n_k</math>
|  style="text-align: center;"|Index of equipment owned by customer <math>k</math>''
|  style="text-align: center;"|
|-
|  style="text-align: center;"|<math>L_{nkp}</math>
|  style="text-align: center;"|Threshold of preventive maintenance <br> for the <math>n_k^{th}</math>part
|  style="text-align: center;"|
|-
|  style="text-align: center;"|<math>\mu</math>
|  style="text-align: center;"|Mean value
|  style="text-align: center;"|Use cycle
|-
|  style="text-align: center;"|<math>\sigma</math>
|  style="text-align: center;"|Magnetic moment
|  style="text-align: center;"|
|-
|  style="text-align: center;"|<math>E_{nk}^\tau</math>
|  style="text-align: center;"|Monitoring data of <math>n_k^{th}</math> part in period  <math>\tau</math>
|  style="text-align: center;"|
|-
|  style="text-align: center;"|<math>T_c</math>
|  style="text-align: center;"|Interval time between the two check cycles
|  style="text-align: center;"|Use cycle
|-
|  style="text-align: center;"|<math>T_{nk}^\tau</math>
|  style="text-align: center;"|Predicted residual life of <math>n_k^{th}</math> part in period <math>\tau</math> in theory
|  style="text-align: center;"|Use cycle
|-
|  style="text-align: center;"|<math>t^\tau</math>
|  style="text-align: center;"|Time of check cycle
|  style="text-align: center;"|Use cycle
|-
|  style="text-align: center;"|<math>C_c</math>
|  style="text-align: center;"|Unit check cost
|  style="text-align: center;"|''yuan''
|-
|  style="text-align: center;"|<math>C_{pf}</math>
|  style="text-align: center;"|Unit cost of preventive maintenance
|  style="text-align: center;"|''yuan''
|-
|  style="text-align: center;"|<math>C_{pr}</math>
|  style="text-align: center;"|Unit cost of collective maintenance
|  style="text-align: center;"|''yuan''
|-
|  style="text-align: center;"|<math>C_s</math>
|  style="text-align: center;"|Unit delay loss
|  style="text-align: center;"|''yuan''
|-
|  style="text-align: center;"|<math>C_k</math>
|  style="text-align: center;"|Unit inventory cost
|  style="text-align: center;"|''yuan''
|-
|  style="text-align: center;"|<math>C_b</math>
|  style="text-align: center;"|Unit ordering cost
|  style="text-align: center;"|''yuan''
|-
|  style="text-align: center;"|<math>S</math>
|  style="text-align: center;"|Maximum of inventory level
|  style="text-align: center;"|
|-
|  style="text-align: center;"|<math>s_k^\tau</math>
|  style="text-align: center;"|Current inventory level in customer <math>k</math> in period <math>\tau</math>
|  style="text-align: center;"|''yuan''
|-
|  style="text-align: center;"|<math>Q_k^{\tau i}</math>
|  style="text-align: center;"|Number of spare parts replenishments of customer <math>k</math> in resupply period <math>\tau_i</math>
|  style="text-align: center;"|
|-
|  style="text-align: center;"|<math>Q_k^{\tau i}</math>
|  style="text-align: center;"|Spare parts consumption of customer <math>k</math> in resupply period <math>\tau_i</math>
|  style="text-align: center;"|
|-
|  style="text-align: center;"|<math>t_k^{o\tau i}</math>
|  style="text-align: center;"|Start time of resupply in customer <math>k</math> in resupply period <math>\tau_i</math>
|  style="text-align: center;"|Use cycle
|-
|  style="text-align: center;"|<math>T_k^{o\tau i}</math>
|  style="text-align: center;"|Interval time of inventory delay in resupply period <math>\tau_i</math>
|  style="text-align: center;"|Use cycle
|-
|  style="text-align: center;"|<math>t_k^{r\tau i}</math>
|  style="text-align: center;"|Ending time of resupply in customer <math>k</math> in resupply period <math>\tau_i</math>
|  style="text-align: center;"|Use cycle
|-
|  style="text-align: center;"|<math>N_c</math>
|  style="text-align: center;"|Amount times of checking
|  style="text-align: center;"|
|-
|  style="text-align: center;"|<math>N_{pf}</math>
|  style="text-align: center;"|Amount times of preventive maintenance
|  style="text-align: center;"|
|-
|  style="text-align: center;"|<math>N_{pr}</math>
|  style="text-align: center;"|Amount times of collective maintenance
|  style="text-align: center;"|
|-
|  style="text-align: center;"|<math>N_o</math>
|  style="text-align: center;"|Amount of ordering
|  style="text-align: center;"|
|-
|  style="text-align: center;"|<math>N_s</math>(<math>t^\tau</math>)
|  style="text-align: center;"|Amount of delay parts from <math>t^{\tau^{-1}}</math> to <math>t^\tau</math>
|  style="text-align: center;"|
|-
|  style="text-align: center;"|<math>N_v</math> (<math>t^\tau</math>)
|  style="text-align: center;"|Amount of inventory at <math>t^\tau</math>
|  style="text-align: center;"|
|}

===3.2 Life distribution===

Suppose ''p''-dimensional random vector, as (1) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>E(t)=[E_1(t),E_2(t),E_3(t)\cdots E_p(t)]{}'</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (1)
|}


It is satisfying <math>p</math>-meta continuous stochastic processes {<math>E(t)</math>, <math>t\ge 0</math>} that concluded the following natures [31]:

(1) Assume initial quantity is <math>E(t) = (0,0,0\cdots 0)</math>. When <math>t = 0</math>, <math>E(t)</math> is continuous.

(2) <math>\Delta E = E(t + \Delta t) - E(t)</math> shows the changing amount between <math>t + \Delta t</math> and <math>t</math>. It obeys <math>p</math>-dimensional Normal distribution

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\Delta E\sim N(\Delta t\mu ,\Delta t\sigma )</math> 
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (2)
|}

where <math>\mu</math> is the mean value and <math>\sigma</math> is the standard deviation, as (3) and (4) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\mu {}'=(\begin{array}{cccc}
{\mu }_1\mbox{ } & {\mu }_2 & \mbox{ }\mbox{ }\mbox{ }\cdots  & {\mu }_\mbox{p}
\end{array}\mbox{ }\mbox{ }){}'\mbox{ }\mbox{ }\mbox{ }</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (3)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\mbox{ }\sigma =\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ }\mbox{ }\left(\mbox{ }\mbox{ }\mbox{ }\begin{array}{cccc}
{\sigma }_1^2 & {\sigma }_{12}^{} & \mbox{ }\mbox{ }\cdots  & {\sigma }_{1p}^{}\\
{\sigma }_{12}^{} & {\sigma }_2^2 & \cdots  & {\sigma }_{2p}^{}\\
\vdots  & \vdots  & \ddots  & {\sigma }_{3p}^{}\\
{\sigma }_{1p}^{} & {\sigma }_{2p}^{} & {\sigma }_{3p}^{} & {\sigma }_p^2
\end{array}\right)\mbox{ }\mbox{ }</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (4)
|}

Then regard

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>{\rho }_{xy}=\frac{{\sigma }_{xy}}{{\sigma }_x{\sigma }_y}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (5)
|}

as the correlation coefficient of <math>E_x(t)</math> and <math>E_y(t)</math>, (<math>x, y = 1,2\cdots p</math>, <math>x \neq y</math>).

(3) In any two non-adjacent time intervals [<math>t_1</math>, <math>t_2</math>] and [<math>t_3</math>, <math>t_4</math>] <math>(t_1</math>< <math>t_2</math>< <math>t_3</math>< <math>t_4)</math>, the respective increments <math>E(t_4) - E(t_3)</math> and <math>E(t_2) - E(t_1)</math> are independent of each other.

If satisfying, <math>E(t)</math> is regarded as a Wiener process. The Normal distribution satisfied by the increment contains the parameters of the Wiener process.

If <math>B(t)</math> is defined as a standard Brownian motion, that is, the univariate Wiener process can be expressed as shown in (6):

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>E(t)=\mu t+\sigma B(t)</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (6)
|}

where <math>\mu</math> is the drift parameter, <math>\sigma</math> is the diffusion parameter in this formula. Because of long-term use, the parts are bound to fail so that <math>\sigma</math> is higher <math>0 (\sigma > 0)</math>. From (6), <math>E(t)</math> is Brownian motion with linear drift [32].

According to the hypothesis and definition, each increment of the multivariate Wiener process obeys independent Normal distribution and belongs to non- monotonic Markov process. The mean value and variance of them are:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>E(E(t))=\mu t</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (7)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>Var(E(t))={\sigma }^2t</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (8)
|}

Because <math>E(t)</math> is a Wiener process, <math>E_x(t)</math> which is its arbitrary row vector is also a Wiener process. It is also subject to (6). If the performance degradation process in a product conforms to it and the degradation value reaches the failure threshold (<math>L</math>), the product is determined to be a failure. The life time (<math>T</math>) can be defined as the time when the first performance of that degrades to <math>L</math>:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>T=inf\lbrace t\vert E(t)=L,t\geq 0\rbrace </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (9)
|}

The reference [33] shows that the life distribution is subject to the inverse Gaussian distribution so as to obtain the probability density function (PDF) and life distribution function (LDF). As (10) and (11) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>f(t)=\frac{L}{\sqrt{2\pi {\sigma }^2t^3}}exp\left[- \frac{{\left(L-\mu t\right)}^2}{2{\sigma }^2t}\right]</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (10)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>F(t)=\Phi \left(\frac{\mu t-L}{\sigma \sqrt{t}}\right)+</math><math>exp\left(\frac{2\mu t}{{\sigma }^2}\right)\Phi \left(\frac{-L-\mu t}{\sigma \sqrt{t}}\right)</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (11)
|}

As a result, the expectation and variance of <math>T</math> are determined to be as (12) and (13):

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>E(T)=\frac{L}{\mu }</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (12)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>D(T)=\frac{L{\sigma }^2}{{\mu }^3}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (13)
|}

where <math>E(t)</math> is the expectation value and <math>D(t)</math> is the variance value.

Secondly, the parameters are determined. In references [29,30,31], <math>\mu</math> and <math>\sigma</math> are estimated by the method of maximum likelihood estimation. Suppose to get the data through monitoring the degradation process. Then, these can be calculated by preprocessing. According to the definition of LDF, what is the difference value between the current data and the former is subject to Normal distribution.

Likelihood function can be established by the former steps, as (14) follow:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>L(\mu ,{\sigma }^2)=\prod_{x=1}^{n_k}\prod_{y=1}^{m_k}\frac{1}{\sqrt{2{\sigma }^2\pi \Delta t_x^y}}exp\left[-\frac{{\left(\Delta E_x^y-\mu \Delta t_x^y\right)}^2}{2{\sigma }^2\pi \Delta t_x^y}\right]</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (14)
|}

The maximum likelihood estimation value of <math>\mu</math> and σ are as follow (15) and (16):

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\overset{\mbox{ˆ}}{\mu }=\frac{\displaystyle\sum_{x=1}^{n_k}E_x^y}{\displaystyle\sum_{x=1}^{n_k}t_x^y}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (15)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\hat{\sigma}=\left\{\frac1{\displaystyle\sum_{x=1}^{n_k}{m_k}}\left[\sum_{x=1}^{n_k}{\displaystyle\sum_{y=1}^{m_k}{\displaystyle\frac{\left(\Delta E_x^y\right)^2}{\Delta t_x^y}}}-\displaystyle\frac{\left(\displaystyle\sum_{x=1}^{n_k}{E_x^{m_k}}\right)^2}{\displaystyle\sum_{x=1}^{n_k}{t_x^{m_k}}}\right]\right\}^2</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (16)
|}

Following the former steps, the performance degradation function of parts based on the Wiener process can be established.

===3.3 Condition-based maintenance===

From the former, the system gets <math>E_{n_k}^{\tau}</math> at <math>t^{\tau}</math> by the continuous check. In modern industry, the economic cost caused by long-term reliance on corrective maintenance is getting higher and higher. It is generally accepted that the preventive maintenance policy can improve the availability of equipment. For the maintenance needs time, especially the delay loss in production is relatively high, the maintenance department generally replaces the part directly. There are two main ways in CBM, which are corrective and preventive maintenance.

By the collective maintenance, it would be happened in <math>t^{\tau}</math>, when the part is failed in the interval time from <math>t^{\tau^{-1}}</math> to <math>t^{\tau}</math>. For the preventive maintenance, the <math>n_k^{th}</math> part will be under that process when <math>E_{nk}^\tau</math> is up to <math>L_{nkp}</math>. Because the failure may lead to other concurrent failures, the associated parts should be thoroughly examined when the replacement is carried out. Resulting from it, <math>C_{pr}</math> will be higher than <math>C_{pf}</math> (<math>C_{pr} > C_{pf}</math>).

Secondly, in order to avoid delay loss as much as possible, the part will be executed the preventive maintenance when its degradation value is not up to <math>L</math> so that it can avoid collective maintenance in the turning process to next check cycle. That’s as (17) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>F(t^{\tau })=\Phi \left(\frac{\mu t^{\tau }-L_{n_k p}-E_{n_k}^{\tau }}{\sigma \sqrt{t^{\tau }}}\right)+</math><math>exp\left(\frac{2\mu t^{\tau }}{{\sigma }^2}\right)\Phi \left(\frac{-L_{n_kp}-E_{n_k}^{\tau }-\mu t^{\tau }}{\sigma \sqrt{t^{\tau }}}\right)>\alpha </math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (17)
|}


When <math>E_{n_k}^{\tau }</math> is up to the threshold of preventive maintenance or collective maintenance, the part will be replaced with a new one if there are spare parts in inventory.

However, it is likely to occur that the spare parts are in delay in reality. In this case, there is still a delay before replacing. The parts that cannot undergo preventive maintenance will continue to be used in order to decrease delay cost. But if the resupply is gotten, the maintenance will continue. The parts which are going to be collective maintenance will be replaced firstly. What’s more, the replacing time is so short than the lifetime that it is ignored when modeling.

Take the condition-based maintenance process for two types of parts as an example. The initial degradation value <math>E_{n_k}^1</math> is 0. [[#img-1|Figure 1]] shows the main process of condition-based maintenance.

<div id='img-1'></div>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[Image:Review_510195169132-image18.jpeg|center|600px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figure 1'''. Condition-based maintenance process of two parts 
|}


In [[#img-1|Figure 1]], <math>E_1^1 </math>  which is the degradation value of the first part at <math>t^1</math> is lower than <math>L_{1_p} (E_1^1 < L_{1_p}</math>) which is its threshold value of preventive maintenance. But the predicted residual life <math>T_1^{1'} </math> is less than <math>T_c </math> (<math>T_1^{1'} < T_c</math>). So, the part is under preventive maintenance at this time. Turn to <math>t^2</math> which is the second check cycle. The degradation value <math>E_1^2</math> which is the degradation value at the current cycle is also lower (<math>E_1^2  < L_{1_p}</math>). Different from roe former time, the predicted residual life <math>T_1^{2'} </math> is higher than <math>T_c </math> (<math>T_1^{2’} > T_c</math>). So, it needn’t be paired. However, turn to <math>t^3</math>, the situation is similar to that in <math>t^1</math>. The part is going to be under preventive maintenance as well.

Next to talking about the second part. At ''t''<sup>1</sup>, ''E''<sub>2</sub><sup>1</sup> which is the degeneration value of second part is lower than ''L''<sub>2</sub>''<sub>p</sub>'' (''E''<sub>2</sub><sup>1</sup> < ''L''<sub>2</sub>''<sub>p</sub>'') which is its threshold of preventive maintenance. According to former formulas, the predicted residual life ''T''<sub>2</sub><sup>1</sup>''’'' is less than ''T<sub>c</sub>'' (''T''<sub>2</sub><sup>1</sup>''’'' > ''T<sub>c</sub>''). So, it needn’t to brplacedweve. Hor, in the interval time between ''t''<sup>1</sup> and ''t''<sup>2</sup>, there is a failure in the part. It caused the breakdown loss. Apart from that, the part won’t be replaced because it isn’t reached the check cycle. It didn’t begin the collective maintenance until it is at ''t''<sup>3</sup>. Turning to the next check cycle, ''E''<sub>2</sub><sup>3</sup> which is the current degeneration value is lower than ''L''<sub>2</sub>''<sub>p</sub>'' (''E''<sub>2</sub><sup>3</sup> < ''L''<sub>2</sub>''<sub>p</sub>''). But ''T''<sub>2</sub><sup>3</sup>''’'' which is the current predicted residual life is less than ''T<sub>c</sub>'' (''T''<sub>2</sub><sup>3</sup>''’'' > ''T<sub>c</sub>''). The part is replaced at this time.

==4. Joint maintenance-inventory policy==

===4.1 Propose a joint policy===

(''s''-1, ''s'') policy is one of the new aspects in modern inventory policy. For some valuable parts, it is often adopted and requires monitoring of inventory. Assuming that in the beginning, the ''k''<sup>th</sup> inventory is up to the ''S''. The parts couldn’t degenerate in inventory. Once the replacement happens, the resupply applicant will be sent out.

It takes time for spare parts from application to arrival. During it, the parts on work are still performing. So, ''Q<sub>k</sub><sup>τi</sup>'' is not only determined by Δ''s<sub>k</sub><sup>τ</sup>'' which is the difference value in period ''τ'' between ''S'' and ''s<sub>k</sub><sup>τ</sup>'', but also ''Q<sub>k</sub><sup>τi</sup>’'' which is the consumptions in ''τ<sub>i</sub>''.

Still take the condition-based maintenance process of the two parts mentioned in Fig. 1 as an example. It supposes the ''k''<sup>th</sup> customer owns these two types of parts. The parts play an important role in each machine and consume the same kind of spare parts. So, the threshold of them is the same (''L''<sub>1</sub> ''='' ''L''<sub>2</sub> and ''L''<sub>1</sub>''<sub>p</sub>'' ''='' ''L''<sub>2</sub>''<sub>p</sub>''). We can make a decision in inventory ordering. That is as Fig. 2 following.
[[Image:Review_510195169132-image19.jpeg|center|600px]]<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Figure 2''' Analysis process in maintenance and inventory for the two kinds of parts</div>

At ''t''<sup>1</sup>, the first part has completed the preventive maintenance, but the degradation of the second part does not meet the replacement conditions. So, the inventory decreases 1. At the same time, the resupply application is sent. Turning to ''t''<sup>2</sup>, the first part has not replaced, and the second part is under collective maintenance. So, the inventory continues to decrease 1. However, the system can’t send the application because it has done in the former cycle. Time changes to ''t''<sup>3</sup>. The first part is under preventive maintenance. As for the second part, the condition of preventive maintenance has been reached. So, the inventory decreases 2. At the same time, the customer has gotten the first resupply so that the second resupply application can be sent.

==4.2. Resupply of Spare Parts==

<span id='_Hlk67870711'></span><span id='_Hlk67870737'></span>As Fig. 2 following, the first resupply application is sent to the system at ''t''<sup>1</sup>. When the application is sent, ''t''<sup>1</sup> is also ''t<sub>k</sub><sup>o</sup>''<sup>1</sup> which is the start time of the first resupply. Spending ''T<sub>k</sub><sup>o</sup>''<sup>1</sup> which is the delay time of the first inventory resupply, it goes to ''t<sub>k</sub><sup>r</sup>''<sup>1</sup> which is the ending time. At this time, the first resupply has finished. ''t<sub>k</sub><sup>r</sup>''<sup>1</sup> is equal to ''t<sub>k</sub><sup>o</sup>''<sup>1</sup> and ''T<sub>k</sub><sup>o</sup>''<sup>1</sup> (''t<sub>k</sub><sup>r</sup>''<sup>1</sup> = ''t<sub>k</sub><sup>o</sup>''<sup>1</sup> + ''T<sub>k</sub><sup>o</sup>''<sup>1</sup>). Turning to ''t''<sup>3</sup>, the second resupply application is sent. The process is the same as the former. ''t''<sup>3</sup> is also ''t<sub>k</sub><sup>o</sup>''<sup>2</sup> which is the start time of the second resupply. ''t<sub>k</sub><sup>r</sup>''<sup>1</sup> which is the ending time is determined by not only ''t<sub>k</sub><sup>o</sup>''<sup>2</sup>, but also ''T<sub>k</sub><sup>o</sup>''<sup>2</sup> which the delay time of the second inventory resupply (''t<sub>k</sub><sup>r</sup>''<sup>2</sup> = ''t<sub>k</sub><sup>o</sup>''<sup>2</sup> + ''T<sub>k</sub><sup>o</sup>''<sup>2</sup>). As a result, they can be obtained that the inventory delay time, start time and, ending time in every resupply period.

However, the spare parts which are adopted (''s''-1, ''s'') policy generally need long-term and multi-period resupply. Therefore, it is unreasonable for the previous fixed delay time model to continue to be used in this kind of spare parts. In this article, it is assumed that the delay time satisfies the Poisson distribution.

=5. Modeling in Joint Policy=

For the same kind of parts that adopt data monitoring to estimate the residual life. Too frequent inspection cannot provide a more accurate data basis for the life prediction. What’s more, it may affect the fitting speed of fault function based on the Wiener process because of the large amount. At the same time, this behavior will continue to increase the cost of testing, resulting in excessive cost.

On the contrary, it may increase the probability of part failure in the next inspection that ''T<sub>c</sub>'' is too large so that it could lose the significance of preventive maintenance policy. Secondly, if ''L<sub>nkp</sub>'' is too high, it will increase the times and the ordering cost. Apart from the two aspects, increasing maximum inventory could affect the inventory cost as well. What’s more, the decreasing inventory may cause the delay loss so that the availability of equipment will be affected.

In conclusion, the joint policy of maintenance and inventory ordering focuses on adjusting the decision variables which are ''T<sub>c</sub>'',  ''S,'' and ''L<sub>nkp</sub>'' to control the costs of preventive maintenance, collective maintenance, check, ordering even delay loss. Under the condition of certain equipment reliability, it may reduce the total cost as much as possible to obtain the optional joint combination of corresponding decision variables.

When calculating the joint model, there are mainly two methods. One is establishing an accurate set of formulas after a complete mathematical derivation. But it is difficult to use explicit expressions to clearly express the relationship between objective functions and decision variables. Then, the joint optimization for the same kind of multi-type parts in this article is a logical process triggered randomly at a discrete time by the change of decision variables. Another is the Marko simulation. It can solve this kind of problem well. By adjusting the decision parameters and analyzing a large number of simulations, the value distribution of objective function in this interval of time is obtained to get the optional joint policy.

The objective function adopted the expense cost rate (''EC'') to calculate the joint model. It is argued the target of ''EC'' is minimum in this model, as (18) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>minEC_{\infty }=minf_{EC}(T_c,S,L_{n_kp})</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (18)
|}


Because the monitoring object is the same kind of parts, their threshold values of preventive maintenance are all roughly the same, (18) can be adjusted to (19):

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>minEC_{\infty }=minf_{EC}(T_c,S,L_p)</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (19)
|}


'' 
{|
|-
| <math>s.t.</math>
| <math>t^{\tau }=\tau \times T_c,\tau \mbox{=}1,2,3...\psi </math>
|}
''

'' <math>T_c>0</math> ''

'' <math>S\in N^+</math> ''

'' <math>L_p>0</math> ''

According to the method in calculating ''EC<sub>∞</sub>'' from refence [34], it can be changed to (20):

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>EC_{\infty }=\underset{t_m\rightarrow \infty}{lim}\frac{C_c\cdot N_c+C_{pf}\cdot N_{pf}+C_{pr}\cdot N_{pr}+C_o\cdot N_o+C_s\cdot \sum_{\tau =1}^{\psi }N_s(t^{\tau })+C_v\cdot \sum_{\tau =1}^{\psi }N_v(t^{\tau })}{t_m\cdot N}</math>
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (20)
|}

However, it is a relatively large project to use Monte Carlo simulation to logically analyze the maintenance-ordering process. It spends so much time and achieved a little part for the whole process. In order to speed up the optimization rate and reduce the cost-effectiveness ratio of research, an intelligent optimization algorithm is used. It not only speeds up the calculation, but also provides ways to the real-time decision-making inventory for the same kind of spare parts in the future. The particle swarm optimization (PSO) algorithm is chosen among all intelligent optimization algorithms in this article

The PSO algorithm is one of the most widely used algorithms in multi-period model optimization [35]. It is paid attention to by researchers because of its easy implementation, high precision, fast convergence, and so on. It has demonstrated its superior algorithm characteristics in practical applications. What’s more,  the PSO algorithm is also a parallel algorithm. It can improve the efficiency of optimization by making full use of the parallel computing advantages of processors.

==5.1. PSO Algorithm==

In the field of algorithm research, there are many new algorithms emerging in endlessly. Many researchers combine the advantages of different algorithms to improve [36].

Whether it is the traditional algorithm or the new intelligent algorithm, it is mainly divided into two kinds: heuristics and meta-heuristics. It has been recognized by more and more researchers that the optimization policy can be obtained quickly by using the PSO algorithm. However, the convergence effect of the traditional PSO algorithm is not satisfactory because of the diverse and complex conditions in the joint model. Therefore, some researchers combine the algorithm with dynamic optimization [37], and some researchers introduce adaptive weight into the algorithm [38].

The improved PSO algorithm from reference [39] is adopted in this article to solve multi-period continuous optimization problems. It adjusts the convergence weight of particles in global and local optimization by combining the nonlinear dynamic inertia weight used to detect the global optimal point and the penalty factor of dynamic correction. It can adapt to the characteristics of multi-period iteration in this model to improve the convergence rate of calculation.

==5.2. Maintenance and Ordering==

1) ''Initialization''

Firstly, initialize ''E<sub>n</sub>''<sub>1</sub><sup>0</sup>, ''E<sub>n2</sub>''<sup>0</sup>, ''E<sub>n3</sub>''<sup>0</sup>... ''E<sub>nK</sub>''<sup>0</sup>. Those are the initial state of parts. And then, initialize the parameters by inputting, such as ''T<sub>c</sub>'', ''S'', ''L<sub>p</sub>'', ''μ'', ''σ, n<sub>K</sub>, ''Φ, ''K,'' ''L,'' ''σ, T''<sup>0</sup>''<sub>knk</sub>,'' ''C<sub>c</sub>,'' ''C<sub>pf</sub>,'' ''C<sub>pr</sub>,'' ''C<sub>o</sub>,'' ''C<sub>s</sub>,'' ''C<sub>v</sub>''.

Before turning to the general degradation process, the system will pay attention to whether there is resupply getting or not. As (21) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>OR_k^{\tau }=\lbrace \begin{array}{c}
\begin{array}{cc}
0 & \left(t_k^{rt_i}/T_c\leq \tau \right)\cup \left(OR_k^{\tau -1}=1\right)\\
OR_k^{\tau -1} & others
\end{array}\mbox{ }
\end{array}</math> 
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (21)
|}


Where ''OR<sub>k</sub><sup>τ</sup>'' is the resupply applicant remark. When there is resupply getting, it is updated to ‘''OR<sub>k</sub><sup>τ</sup>'' = 0’. If not, it will be the same as the former. In the beginning, there is no resupply. ‘''OR<sub>k</sub><sup>τ</sup>''’ will be updated to 1.

2) ''Maintenance for Delay Parts''

The system will replace the parts whose collective maintenance process is a delay. As (22) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>SPN_{n_k}^{\tau }=\lbrace \begin{array}{c}
\begin{array}{cc}
0 & \left(SPN_{n_k}^{\tau -1}\mbox{=}1\right)\cup \left(s_k^{\tau }>0\right)\\
SPN_{n_k}^{\tau -1} & others
\end{array}
\end{array}</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (22)
|}


Where ''SPN<sub>nk</sub><sup>τ</sup>'' is the delay mark in the current cycle, ''SPN<sub>nk</sub><sup>τ</sup>''<sup>-1</sup> is the mark in the last cycle. When the parts are replaced, ''SPN<sub>nk</sub><sup>τ</sup>'' will be changed to 0.

Because of the maintenance, it needs a spare part to replace. ''SPPF<sub>nk</sub><sup>τ</sup>'' and ''s<sub>k</sub><sup>τ</sup>'' are updated as (23) and (24) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>SPPF_{n_k}^{\tau }=1</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (23)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>s_k^{\tau }=s_k^{\tau }-1</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (24)
|}


3) ''Degradation Process''

From the former life prediction analysis, the process of ''n<sub>k</sub>''<sup>th</sup> degradation changing at a continuous point of time could be as (25) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>E_{n_k}^{\tau }=E_{n_k}^{\tau -1}\mbox{+}\Delta E_{n_k}^{\tau }</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (25)
|}


Where Δ''E<sub>nk</sub><sup>τ</sup>'' which obeys Normal distribution is a random increment at the interval time from ''τ''-1 to ''τ'', as (26) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\Delta E_{n_k}^{\tau }(T_c)=\mu \cdot T_c+\sigma \cdot B(T_c)</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (26)
|}


According to section 3, the state of ''n<sub>k</sub>'' in ''τ'' is as (27) and (28) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>SiF_{n_k}^{\tau }=\lbrace \begin{array}{cc}
2 & E_{n_k}^{\tau }\geq L\\
0 & E_{n_k}^{\tau }<L
\end{array}</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (27)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>SiF_{n_k}^{\tau }=\lbrace \begin{array}{cc}
1 & E_{n_k}^{\tau }\geq L_p\\
0 & E_{n_k}^{\tau }<L_p
\end{array}</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (28)
|}


:1) When the value of degradation reaches ''L'', the part is a failure. If there is inventory, it will be going to collective maintenance. It sets a remark which there is a part that need collective maintenance. As (27) following.

:2) When the degradation value is lower than ''L'', the part is normal. However, it may need preventive maintenance.

<span id='_Hlk67995631'></span><span id='_Hlk67995687'></span>
:3) When the degradation is lower than ''L<sub>p</sub>'', the part is normal. Except that, the system will predict the life degradation of part. If ''T<sub>nk</sub><sup>τ</sup>’'' is shorter than ''T<sub>c</sub>'' (''T<sub>nk</sub><sup>τ</sup>’'' < ''T<sub>c</sub>'') and there is inventory in ''k''<sup>th</sup> customer. it will be going to preventive maintenance. As (29) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>\mbox{ }SiF_{n_k}^{\tau }=\lbrace \begin{array}{cc}
1 & \left(\mbox{ }SiF_{n_k}^{\tau }=0\right)\cup \left(T_{n_k}^{\tau }{}'<T_c\right)\\
0 & \mbox{ }others
\end{array}</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (29)
|}


Where ''SiF<sub>nk</sub><sup>τ</sup>'' is the preparation mark of preventive maintenance.

<span id='_Hlk67994995'></span>
:4) When the degradation value is up to ''L<sub>p</sub>'', the part will be going to preventive maintenance. As (28) following.

In conclusion, (27) (28) and (29) can be sum up into to (30):

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>SiF_{n_k}^{\tau }=\lbrace \begin{array}{cc}
2 & E_{n_k}^{\tau }\geq L\\
1 & \left(\left(\mbox{ }E_{n_k}^{\tau }<L_p\mbox{ }\right)\cup \left(T_{n_k}^{\tau }{}'<T_c\right)\right)\cap \left(E_{n_k}^{\tau }\geq L_p\right)\\
0 & others
\end{array}</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (30)
|}


4) ''Inventory Changing''

Before replacing, the system will check the inventory. The amount of preventive and collective maintenance will be marked. For it, the replacement mark is set to mark the process. As (31) and (32) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>SPPF_{n_k}^{\tau }=\lbrace \begin{array}{cc}
1 & \left(SiF_{n_k}^{\tau }=1\right)\cup \left(s_k^{\tau }>0\right)\\
0 & others
\end{array}</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (31)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>SPRF_{n_k}^{\tau }=\lbrace \begin{array}{cc}
1 & \left(SiF_{n_k}^{\tau }=2\right)\cup \left(s_k^{\tau }>0\right)\\
0 & others
\end{array}</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (32)
|}


Where ''SPPF<sub>nk</sub><sup>τ</sup>'' is the preventive maintenance mark, ''SPRF<sub>nk</sub><sup>τ</sup>'' is the collective maintenance mark. When there is inventory, it will finish the corresponding maintenance. When there is little inventory, the replacement process will be delayed. So, ''SPN<sub>nk</sub><sup>τ</sup>'' is updated, as (33) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>SPN_{n_k}^{\tau }=\lbrace \begin{array}{c}
\begin{array}{cc}
1 & \left(SiF_{n_k}^{\tau }=2\right)\cup \left(s_k\mbox{=}0\right)\\
SPN_{n_k}^{\tau } & others
\end{array}
\end{array}</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (33)
|}

When there is little inventory, the process of parts that are going to collective maintenance will be delayed. It is marked ‘''SPN<sub>nk</sub><sup>τ</sup>'' = 1’. But the parts which are going to preventive maintenance will continue to be used because of no inventory. What's more, when inventory is insufficient, priority is given to the parts which are going to collective maintenance.

After maintenance, the number of inventories will be updated. As (34) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>s_k^{\tau }=\lbrace \begin{array}{c}
\begin{array}{cc}
s_k^{\tau }-1 & SPF_{n_k}^{\tau }=1\mbox{ }\\
s_k^{\tau }\mbox{ } & others
\end{array}\mbox{ }
\end{array}</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (34)
|}


Then the state of parts would be changed because of the maintenance. As (35) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>E_{n_k}^{\tau }{}'=\lbrace \begin{array}{cc}
E_{n_k}^{\tau } & others\\
0 & \mbox{ }\left(SPPF_{n_k}^{\tau }\mbox{=}1\right)\cap \left(SPRF_{n_k}^{\tau }\mbox{=}1\right)
\end{array}</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (35)
|}


Where ''E<sub>nk</sub><sup>τ</sup>’'' is the final value of degradation in this part in the current check cycle. When the part is replaced, the degradation will recalculate. If not, it will go on declining. There is one more word that the degradation process wouldn’t be affected by monitoring. Because of that, the start time of life for this part is determined whether to replace or not. The residual life is as (36) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>T_{kn_k}^{\tau }{}'=\lbrace \begin{array}{cc}
T_{kn_k}^0 & \left(SPPF_{n_k}^{\tau }\mbox{=}1\right)\cap \left(SPRF_{n_k}^{\tau }\mbox{=}1\right)\\
T_{kn_k}^{\tau }\mbox{ } & other
\end{array}</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (36)
|}


Where ''T<sub>knk</sub><sup>τ</sup>''<sup>’</sup> is the final value of residual life for this part, ''T''<sup>0</sup>''<sub>knk</sub>'' is the initial life. When the part is replaced, it will be recalculated. If not, the final residual life is equal to the theoretical value in current check period.

5) ''Ordering''

Once the inventory decreases, the resupply applicant would be sent at once. So, ''OR<sub>k</sub><sup>τ</sup>'' will be updated, as (37) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>~OR_k^{\tau }\mbox{=}\lbrace \begin{array}{cc}
1 & \left(s_k^{\tau }<S\right)\mbox{ }\cap \left(OR_k^{\tau }=0\right)\\
0 & others
\end{array}\mbox{ }\mbox{ }\mbox{ }\mbox{ }</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (37)
|}

The resupply period is started to be calculated. When the resupply applicant is sent, ''τ<sub>i</sub>'' will be updated as (38) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>~{\tau }_i\mbox{=}{\tau }_i\mbox{+}1</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (38)
|}


At the same time, the parameters of resupply would be marked, such as ''t<sub>k</sub><sup>oτi</sup>'', ''T<sub>k</sub><sup>oτi</sup>'', ''t<sub>k</sub><sup>rτi</sup>''. As (39) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>t_k^{o{\tau }_i}modT_c=0</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (39)
|}


''T<sub>k</sub><sup>oτi</sup>'' obeys Poisson distribution, as (40) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>T_k^{o{\tau }_i}=\frac{a^{{\tau }_i}}{{\tau }_i!}\cdot e^{-a}</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (40)
|}


Where ''τ<sub>i</sub>'' belongs to the positive integer (''τ<sub>i</sub>'' = 1, 2, 3…''Φ'').

From the former sections, ''t<sub>k</sub><sup>rτi</sup>'' is determined by delay time of inventory and the start time of resupply, as (41) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>t_k^{r{\tau }_i}=T_k^{o{\tau }_i}+t_k^{o{\tau }_i}</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (41)
|}

What's more, ''Q<sub>k</sub><sup>τi</sup>'' is made up of Δ''s<sub>k</sub>'' and ''Q<sub>k</sub><sup>τi</sup>’'', as (42) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>Q_k^{{\tau }_i}=\Delta s_k+Q_k^{{\tau }_i}{}'=S_k-</math><math>s_k^{\tau }+Q_k^{{\tau }_i}{}'</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (42)
|}


Δ''s<sub>k</sub>'' and ''Q<sub>k</sub><sup>τi</sup>’'' are determined by the state of parts. With the resupply getting, the inventory is updated. As (43) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>s_k^{\tau }=\lbrace \begin{array}{cc}
s_k^{\tau -1}+Q_k^{{\tau }_i} & OR_k^{\tau }=0\\
s_k^{\tau -1} & others
\end{array}</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (43)
|}

==5.3. Calculating EC<sub>∞</sub>==

When the decisions of maintenance and inventory policies are made in one check cycle, the values of ''N<sub>pf</sub>'', ''N<sub>pr</sub>'', ''N<sub>c</sub>'', ''N<sub>s</sub>'' and ''N<sub>v</sub>'' are updated. As the formulas following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>~N_{pf}=N_{pf}+\sum_{n_k=1}^{n_K}\mbox{ }SPPF_{n_k}^{\tau }</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (44)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>~N_{pr}=N_{pf}+\sum_{n_k=1}^{n_K}\mbox{ }SPRF_{n_k}^{\tau }</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (45)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>~N_c=N_c+n_K</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (46)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>~N_s=\sum_{n_k=1}^{n_K}SPN_{n_k}^{\tau }</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (47)
|}

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>~N_v=\sum_{n_k=1}^{n_K}s_k^{\tau -1}</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (48)
|}


It is necessary that the system should update the value of the check cycle. As (49) following:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>~\tau \mbox{=}\tau \mbox{+}1</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (49)
|}


When the system turns to the next check cycle, ''τ'' will be updated. The simulation process will back to the ‘initialization’ section and go along with the former process period by period. What’s more, it is the once simulation.

Before the simulation, the researching duration of the model l should be made clear. It would be divided into several check cycles to get discrete syncopation points. The process of maintenance and inventory ordering is simulated by the Monte Carlo method. With the iteration of the check cycle, the changes of parts and inventory in each cycle are deduced. Through the cumulative calculation of costs for all cycles, the expected cost rate of a simulation is obtained finally. This is once result of a simulation. The system will process multi-period simulation by adopting the adaptive improved PSO algorithm. The algorithm is as Fig. 3 following, which expands the optimization process of decision variables ''T<sub>c</sub>'', ''S'', ''L<sub>p</sub>'' in the function of ''EC<sub>∞</sub>''.

To sum up, the whole joint policy is following Fig. 4(a) (b) (c) (d) and (e) to get the optional decision combination (''T<sub>c</sub>''<sup>*</sup>, ''S''<sup>*</sup>, ''L<sub>p</sub>''<sup>*</sup>) under the minimum average expense rate. Where ''T<sub>c</sub>''<sup>*</sup> is the optional check cycle, ''S''<sup>*</sup> is the optional inventory, ''L<sub>p</sub>''<sup>*</sup> is the optional threshold of preventive maintenance. Fig. 4(a) is the part of simulation and maintenance for delay parts. Fig. 4(b) is the general degradation part. Fig. 4(c) is the part of inventory. Fig. 4(d) is the periodic intermediate process. Fig. 4(e) is the ordering part of resupply.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

[[Image:Review_510195169132-image57.jpeg|center|600px]]
'''Figure 3''' Iterative optimization process of the improve PSO algorithm</div>

=6. Numerical investigation=

==6.1. Case Description==

In automotive engines, because of the core position of the engine in the composition of the vehicle, especially racing cars, the fault monitoring of engine (internal combustion engine) and the subsequent maintenance procedure has always been the important contents that researchers pay close attention to. From the view of the maintenance department, the internal combustion engine has two main characteristics. First is its mechanical structure being complex, in order to ensure the tightness of the cylinder and the completeness of combustion. For example, the cylinder of a gasoline internal combustion engine is generally made by pouring with an integrated mold. Apart from it, there are lots of parts loading in the engine, such as bearings, pistons, flywheels.

The subtle hidden trouble caused by any one part may lead to the failure in the engine. Another one is that most of the main body of the automobile engine is generally embedded in the automobile interior so that the general self-test program cannot check whether there is any obvious fault on the outside of the engine. What's more, it’s too hard to check the interior of the engine. Therefore, the German Volkswagen Group, General Motors of the United States, Toyota of Japan, and other companies gradually began to use electronic monitoring and control technology to carry out real-time monitoring of engines and other key components. On the one hand, the central control system notifies the maintenance team belonging to the fleet. On the other hand, it may send the fault analysis report to the customer center belonging to the depot through the remote platform as well. When the vehicle enters the maintenance workshop, the maintenance personnel can download the monitoring data from the central control system for data analysis, so as to provide options for later maintenance and provide suggestions for the improvement of the car to the R & D department of the car company.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

[[Image:Review_510195169132-image58.jpeg|center|1910x1910px]]
'''Figure 4(a)''' The once process of simulation (part I)</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

[[Image:Review_510195169132-image59.jpeg|center|970x970px]]
</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Figure 4(b)''' The once process of simulation (part II)</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''
[[Image:Review_510195169132-image60.jpeg|center|1063x1063px]]
'''
'''Figure 4(c)''' The once process of simulation (part III)</div>

[[Image:Review_510195169132-image61.jpeg|center|228px]]
'''
'''

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Figure 4(d)''' The once process of simulation (part IV)</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

[[Image:Review_510195169132-image62.jpeg|center|1080x1080px]]
'''Figure 4(e)''' The once process of simulation (part V)</div>

According the data from a fleet which is in the interval of time from 2012 to 2017, the team is mainly responsible for the maintenance and modification of 10 vehicles in the fleet, including training, racing, and spare vehicles. All vehicles are monitored for the part of the engineer. According to reference [42], the data was preprocessed. Then the data is processed based on performance degradation modeling. The process does not belong to the main content of this article, so it will not be elaborated here. The result is as Table 2.

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''No.''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''g''1
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''g''2
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''|μ''<sub>1</sub>''|''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''|μ''<sub>2</sub>''|''
|-
|  style="text-align: center;"|1
|  style="text-align: center;"|0.4465
|  style="text-align: center;"|2.929
|  style="text-align: center;"|1.2532
|  style="text-align: center;"|0.722
|-
|  style="text-align: center;"|2
|  style="text-align: center;"|-0.3172
|  style="text-align: center;"|2.7418
|  style="text-align: center;"|1.0619
|  style="text-align: center;"|1.1577
|-
|  style="text-align: center;"|3
|  style="text-align: center;"|0.7363
|  style="text-align: center;"|3.0332
|  style="text-align: center;"|2.0687
|  style="text-align: center;"|1.7633
|-
|  style="text-align: center;"|4
|  style="text-align: center;"|0.2786
|  style="text-align: center;"|2.9975
|  style="text-align: center;"|0.7838
|  style="text-align: center;"|1.4075
|-
|  style="text-align: center;"|5
|  style="text-align: center;"|-2.445
|  style="text-align: center;"|2.6326
|  style="text-align: center;"|0.6725
|  style="text-align: center;"|2.0883
|-
|  style="text-align: center;"|6
|  style="text-align: center;"|0.0526
|  style="text-align: center;"|2.9214
|  style="text-align: center;"|0.1487
|  style="text-align: center;"|0.6422
|-
|  style="text-align: center;"|7
|  style="text-align: center;"|-0.082
|  style="text-align: center;"|2.7745
|  style="text-align: center;"|0.2304
|  style="text-align: center;"|0.8263
|-
|  style="text-align: center;"|8
|  style="text-align: center;"|-0.2252
|  style="text-align: center;"|2.7765
|  style="text-align: center;"|0.605
|  style="text-align: center;"|0.6049
|-
|  style="text-align: center;"|9
|  style="text-align: center;"|0.1047
|  style="text-align: center;"|2.6556
|  style="text-align: center;"|0.2942
|  style="text-align: center;"|2.0192
|-
|  style="border-bottom: 1pt solid black;text-align: center;"|10
|  style="border-bottom: 1pt solid black;text-align: center;"|-0.154
|  style="border-bottom: 1pt solid black;text-align: center;"|2.6817
|  style="border-bottom: 1pt solid black;text-align: center;"|0.4327
|  style="border-bottom: 1pt solid black;text-align: center;"|1.7516
|}


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Table 2''' Statistical results of data analysis</div>

Where column 2-3 are the third-order center distance and fourth-order center distance in order to check the value of ''k''<sup>th</sup> skewness and kurtosis, column 4-5 is the parameter value of statistical test.

Assume that ''α'' which is the significance level for determining whether the degradation process of that can be express by the Wiener process is equal to 0.05. If the value of ''|μ''<sub>1</sub>''|'' and ''|μ''<sub>2</sub>''|'' are both less than ''Z<sub>α</sub>''<sub>/4</sub>, the data can be used by this method. It is known that ''Z<sub>α</sub>''<sub>/4</sub> is equal to 2.25. So, the values of column 3 and 4 aren’t higher than 2.25. After checking, the degradation process can be expressed by Wiener process.

Maintenance and inventory costs have fluctuated over the past 6 years. Resulting from it, we considered the impact of economic purchasing power fluctuations on modeling in this article to process the cost data of 6 years. The calculated costs of every aspect are as Table 3 following.

<span id='_Hlk68630111'></span>According to the data, ''t'' which is the whole optimization time is equal to 20000 use cycle. ''C<sub>c</sub>'' which is the unit check cost is equal to 500 ''yuan''. When the degradation value is less than ''L'', it is assumed that it is a failure. After it was finished the collective maintenance, the system would cost 20000 ''yuan'' which is the unit cost of collective maintenance (''C<sub>pr</sub>''). It would turn back to the initial state until the replacement. However, there will be delay loss if no inventory. It may lead to the cost of personnel, equipment, and reputation because of the inability to carry out daily training, unnormal competitions and other planned schedule work. It is roughly 30000 ''yuan'' through calculating. In the meantime, the team will do some preventive maintenance while the value of degradation is less than ''L<sub>p</sub>''. It would cost 5000 ''yuan'' which is the unit cost of preventive maintenance.

The team adopted the (''s''-1, ''s'') policy to manager. The inventory cost in every check cycle includes the lease of the warehouse, the salary of managerment, and the arrangement cost of the inventory environment. It is roughly 8000 ''yuan'' which is the unit inventory cost. It would not send the resupply applicant until there is no inventory. The ordering cost of unit spare part (''C<sub>o</sub>'') is 8000 ''yuan''.

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''C<sub>c</sub>''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''C<sub>pf</sub>''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''C<sub>pr</sub>''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''C<sub>o</sub>''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''C<sub>s</sub>''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''C<sub>v</sub>''
|-
|  style="border-bottom: 1pt solid black;text-align: center;"|''fitness''

(Unit: ''yuan'')
|  style="border-bottom: 1pt solid black;text-align: center;"|500
|  style="border-bottom: 1pt solid black;text-align: center;"|5000
|  style="border-bottom: 1pt solid black;text-align: center;"|20000
|  style="border-bottom: 1pt solid black;text-align: center;"|8000
|  style="border-bottom: 1pt solid black;text-align: center;"|30000
|  style="border-bottom: 1pt solid black;text-align: center;"|8000
|}


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Table 3''' List of costs after processing</div>

==6.2. Modeling in Degradation==

According to section 2, the maximum parameter estimation method is used to calculate the above data. what set the initial value of state is equal to 75 (''E<sub>nk</sub>''<sup>1</sup> = 75). They are calculated that  <math>\overset{\frown}{\mu }</math> which is equal to -0.01475 and  <math>\overset{\frown}{\sigma }</math> which is equal to 0.39985. The degradation process of internal combustion engine can use (50) to follow:

{| class="formulaSCP" style="width: 100%; text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;" 
|-
| <math>~E(t)=75-0.01475t+0.39985B(t)</math> ''' '''
|}
| style="width: 5px;text-align: right;white-space: nowrap;" | (50)
|}


Set the value of ''L'' to 0. The whole process is simulated by Monte Carlo method.

It is obvious that the state degraded with the time going in Fig. 5. What's more, as time goes by, the offset increases gradually. Except that, it is possible that the predicted value of the next period is higher than the former. It’s the result in the deviation between the linear function and the linear function is getting larger and larger that the scope of Brownian motion is getting wider and wider with the time goes.

According to (12) and (13), the system would predict the residual life, as Fig. 6 following.

==6.3. Result Analysis==

From the equation of ''EC<sub>∞</sub>'', the decision parameters which need optimization conclude: ''T<sub>c</sub>'', ''S'', ''L<sub>p</sub>''. The improved PSO algorithm discussed in section 5 is used to optimize the joint policy.


In order to ensure the efficiency and accuracy of optimization, set the model parameters: ''N'' which is the number of species is equal to 100, ''c''<sub>1</sub> is equal to 1.6962 and ''c''<sub>2</sub> is equal to 1.8962 which are the learning factors, ''w''<sub>min</sub> is 0.6 and ''w''<sub>max</sub> is 0.8 which are inertia weight minimum and maximum, ''M'' which is the number of iterations is equal to 60. As for the resupply, the time of the first inventory delay is 300 use cycles. ''a'' which is the parameter of resupply model is equal to 110.3640 through calculating. Those will be adopted to the optimization process. After that, the system got the optimal decision combination which is [''T<sub>c</sub>''<sup>*</sup>, ''S''<sup>*</sup>, ''L<sub>p</sub>''<sup>*</sup>]. Where it is [1220,5,8]. The mean value of the expense cost rate is 133.835 ''yuan''. The convergence process is as Fig. 7 following.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">[[Image:Review_510195169132-image66.jpeg|center|1222x1222px]]
'''Figure 5''' Performance degradation curve obtained by Monte Carlo simulation</div>

From Fig. 7, the objective function curve tends to be smooth in the last stage of the convergence process, and the fluctuation range of the end decreases gradually. At the same time, the penalty function also tends to be flat, and the value of the function tends to 0. As Fig. 8 following, it’s a conclusion that the convergence effect is good so that the result is the optimal solution after optimization. The optional decision combination shows a fact that the fleet should adjust the check cycle to do one check after 1220 use cycles. On the one hand, the maintenance department should execute the policy of collective maintenance by monitoring. On the other hand, it should obey the policy of preventive maintenance when the value of degradation is less than 8. What’ 'sore, it’s a balance between maintenance efficiency and economic performance that the maximum inventory is set to 5. Once the replacement comes across, the resupply applicant will be sent to ensure the next maintenance.

In refence [40] and [41], the researchers analyzed the influence of parameters on optimization results, such as cost and fault threshold. There has been a comprehensive sensitivity analysis as far. However, the model is built for the multi-type parts of the same kind. The discussion of the new joint policy focuses on the resupply time model to study. And then, the parameters of resupply model are the next targets for discussion. As Table 4 following, the result of sensitivity analysis in ''a'' is gotten.

From Table 4, with the value of ''a'' is higher, the delay time of resupply will be longer. It causes the value of ''T<sub>c</sub>''<sup>*</sup>, ''S''<sup>*</sup>, and ''L<sub>p</sub>''<sup>*</sup> to get higher so that the expense cost rate is also higher. That's easy to be understood. Considering the availability of machines and delay cost, it will do some measures to solve it. One is up to inventory to avoid delay. Another one is up to the threshold of preventive maintenance to decrease the possibility of collective maintenance. However, those methods both result in the higher cost. The first causes the higher cost of inventory. The second cost of preventive maintenance is higher. On the other hand, when the value of ''a'' decreases, the delay time of inventory and the decision parameters don’t decrease as predicted. On the contrary, the threshold of preventive failure increases. In refence [42], the researcher has discussed the numerical analysis in the updating process of multi-period (''s'', ''S'') policy. It can be inferred that the probability of delay in spare parts is greatly reduced when the time of resupply is reducing. Because of it, the center of optimization changes to the maintenance policy from the inventory policy. The optional joint policy points out it that does well in reducing the total cost through raising the threshold of preventive maintenance appropriately. It reduced the cost of collective maintenance, while the cost of preventive maintenance is getting higher.

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

[[Image:Review_510195169132-image67.jpeg|center|1268x1268px]]
'''Figure 6''' Life prediction of part in performance degradation process</div>


[[Image:Review_510195169132-image68.jpeg|center]]

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Figure 7''' Objective function change process by improved PSO algorithm</div>

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">

[[Image:Review_510195169132-image69.jpeg|center|1268x1268px]]
'''Figure 8''' Penalty function change process by improved PSO algorithm</div>

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''No.''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''a''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''T<sub>c</sub>''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''S''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''L<sub>p</sub>''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''E''(''EC<sub>∞</sub>'')
|-
|  style="border-top: 1pt solid black;text-align: center;"|1
|  style="border-top: 1pt solid black;text-align: center;"|183.9400
|  style="border-top: 1pt solid black;text-align: center;"|1260
|  style="border-top: 1pt solid black;text-align: center;"|8
|  style="border-top: 1pt solid black;text-align: center;"|11
|  style="border-top: 1pt solid black;text-align: center;"|148.690
|-
|  style="text-align: center;"|2
|  style="text-align: center;"|147.1520
|  style="text-align: center;"|1240
|  style="text-align: center;"|6
|  style="text-align: center;"|10
|  style="text-align: center;"|144.185
|-
|  style="text-align: center;"|3*
|  style="text-align: center;"|110.3640
|  style="text-align: center;"|1220
|  style="text-align: center;"|5
|  style="text-align: center;"|8
|  style="text-align: center;"|133.835
|-
|  style="border-bottom: 1pt solid black;text-align: center;"|4
|  style="border-bottom: 1pt solid black;text-align: center;"|73.5759
|  style="border-bottom: 1pt solid black;text-align: center;"|1217
|  style="border-bottom: 1pt solid black;text-align: center;"|5
|  style="border-bottom: 1pt solid black;text-align: center;"|10
|  style="border-bottom: 1pt solid black;text-align: center;"|134.305
|}


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Table 4''' Sensitivity analysis in parameters of resupply model</div>

What's more, the model is adopted (''s''-1, ''s'') policy. The optimization effect has been gotten by compared with the (''s'', ''S'') policy. The result is as Table 5 following. Under the same condition, the value of ''S''<sup>*</sup> is smaller, while the values of ''T<sub>c</sub>''<sup>*</sup> and ''L<sub>p</sub>''<sup>*</sup> are higher in the (''s''-1, ''s'') policy with the comparison of the (''s'', ''S'') policy. Through analysis, the resupply applicant will be sent when the inventory begins to reduce. It is unnecessary that maintaining a relatively high inventory. It reduces the inventory cost, as Fig. 9(a) and (b) following. Where the blue square stands for the cost of (''s''-1, ''s'') policy, the red one is on behalf of the cost of (''s'', ''S'') policy. However, the joint policy adopted (''s'', ''S'') policy increases the check cycle and preventive maintenance threshold to avoid the cost of corrective maintenance, while the cost of collective maintenance is reduced.

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|Inventory policy
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''T<sub>c</sub>''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''L<sub>p</sub>''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''S''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''s''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|''E''(''EC<sub>∞</sub>'')
|-
|  style="border-top: 1pt solid black;text-align: center;"|(''s''-1, ''s'') policy
|  style="border-top: 1pt solid black;text-align: center;"|1220
|  style="border-top: 1pt solid black;text-align: center;"|8
|  style="border-top: 1pt solid black;text-align: center;"|5
|  style="border-top: 1pt solid black;text-align: center;"|/
|  style="border-top: 1pt solid black;text-align: center;"|135.484
|-
|  style="border-bottom: 1pt solid black;text-align: center;"|(''s'', ''S'') policy
|  style="border-bottom: 1pt solid black;text-align: center;"|1160
|  style="border-bottom: 1pt solid black;text-align: center;"|6
|  style="border-bottom: 1pt solid black;text-align: center;"|7
|  style="border-bottom: 1pt solid black;text-align: center;"|2
|  style="border-bottom: 1pt solid black;text-align: center;"|145.254
|}


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Table 5''' Comparison of the results of inventory strategy between (''s'', ''S'') and (''s''-1, ''s'') under the same conditions</div>


[[Image:Review_510195169132-image70.jpeg|center|1268x1268px]]

<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
'''Figure 9(a)''' The comparison of every aspect of costs between the (''s''-1, ''s'') policy and (''s'', ''S'') policy (The part of ''C'', ''C<sub>pr</sub>'', ''C<sub>o</sub>'', ''C<sub>s</sub>'', ''Cv'')</div>


<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
[[Image:Review_510195169132-image71.jpeg|center|1268x1268px]]
'''Figure 9(b)''' The comparison of every aspect of  costs between the (''s''-1, ''s'') policy and (''s'', ''S'') policy (The part of ''C<sub>c</sub>'', ''C<sub>pr</sub>'')</div>


The optional decision combination adopted (''s'', ''S'') policy needs the wider space for inventory either because of the larger number of maximum inventories. However, it is an important aspect of considering for some smaller factory or maintenance team, for example, the fleet discussed in this article. In conclusion, the joint policy adopted (''s''-1, ''s'') policy will ensure the minimum cost under the certain maintenance efficiency.

=7. Conclusion=

The center of this article focuses on the joint policy of maintenance and inventory ordering for the multi-style of the same kind. Performance degradation model is obtained by using widespread life prediction methods based on the Winer-based process.

There is a fact that the condition of the model is very complex. The maintenance budget is relatively low for some customers so that the inventory space is the same. The (''s''-1, ''s'') policy adopted in this model has some characteristics which can solve the problems, such as high-security efficiency, good economic benefits, and so on. Then, the resupply processes are multi-period and fluctuate changing. For this characteristic, the joint model has rebuilt the resupply process. We balanced to optimize the decision parameters which conclude check cycle, maximum inventory, and the threshold of preventive maintenance to get the optional decision combination. On the other hand, in order to reduce the cost-effectiveness ratio, the optimizing process is adopted the Monte Carlo method and the improved PSO algorithm.

Based on the joint model and optional decision combination, the effect of resupply model parameter on the optimization result is gotten. At the same time, the superiority of the optional policy is shown by the comparison with the result adopted by the (''s'', ''S'') policy.


Those are the whole research results in this article. While providing optimal solutions for industries, they can also provide a feasible reference for the multi-style part of the same kind in complex maintenance environments.



==Refences:==

<span id='ZOTERO_BREF_LAx9KFHLoLQz'></span><big>[1] G.T. Friedlob, F.J. Plewa. Understanding return on investment. Wiley. New York, May 1996.</big>

<big>[2] Z. Liu, L. Zhang. A review of failure modes, condition monitoring and fault diagnosis methods for large-scale wind turbine bearings. Measurement. 149: 107002, 2020.</big>

<big>[3] K. Le Son, M. Fouladirad, A. Barros, E. Levrat, B. Iung. Remaining useful life estimation based on stochastic deterioration models: A comparative study, Reliab. Eng. Syst. Saf. 112: 165–175, 2013.</big>

<big>[4] R.B. Chinnam. On-line reliability estimation for individual components using statistical degradation signal models. Qual. Reliab. Eng. Int. 18: 53–73, 2002.</big>

<big>[5] G.A. Whitmore, F. Schenkelberg. Modelling Accelerated Degradation Data Using Wiener Diffusion With A Time Scale Transformation. Lifetime Data Anal. 3: 27–45, 1997.</big>

<big>[6] M.E. Robinson. Bayesian methods for a growth-curve degradation model with repeated measures. Lifetime Data Anal. 6: 357–374, 2000. </big>

<big>[7] Guo Chiming, W. Wang, B. Guo, X. Si. A maintenance optimization model for mission-oriented systems based on Wiener degradation. Reliab. Eng. Syst. Saf. 111: 183–194, 2013.</big>

<big>[8] W. Wang. A stochastic model for joint spare parts inventory and planned maintenance optimization. Eur. J. Oper. Res. 216: 127–139, 2012.</big>

<big>[9] G.J. Feeney, C.C. Sherbrooke. The (''s'' − 1, ''s'' ) Inventory Policy Under Compound Poisson Demand. Manag. Sci. 12: 391–411, 1966. </big>

<big>[10] S. Bashyam, M.C. Fu. Optimization of (s, S) Inventory Systems with Random Lead Times and a Service Level Constraint. Manag. Sci. 44: S243–S256, 1998.</big>

<big>[11] T. Shibuya, T. Dohi, S. Osaki. Optimal continuous review policies for spare part provisioning with random lead times. Int. J. Prod. Econ. 55: 257–271, 1998.</big>

<big>[12] S.K. Srinivasan. Analysis of s, S inventory systems with general lead time and demand distribution and adjustable reorder size. Optimization. 19: 557–576, 1988.</big>

<big>[13] Z. Zhang, X. Si, C. Hu, Y. Lei. Degradation data analysis and remaining useful life estimation: A review on Wiener-process-based methods. Eur. J. Oper. Res. 271: 775–796, 2018.</big>

<big>[14] H. Wang, W. Song, E. Zio, A. Kudreyko, Y. Zhang. Remaining useful life prediction for Lithium-ion batteries using fractional Brownian motion and Fruit-fly Optimization Algorithm, Measurement. 161, 107904, 2020.</big>

<big>[15] L. Bian, N. Gebraeel. Stochastic modeling and real-time prognostics for multi-component systems with degradation rate interactions. IIE Trans. 46: 470–482, 2014.</big>

<big>[16] A.H. Elwany, N.Z. Gebraeel. Sensor-driven prognostic models for equipment replacement and spare parts inventory. IIE Trans. 40: 629–639, 2008.</big>

<big>[17] J.-L. Zhang, J. Chen. Bayesian solution to pricing and inventory control under unknown demand distribution. Oper. Res. Lett. 34: 517–524, 2006.</big>

<big>[18] X. Zhang, J. Zeng. Joint optimization of condition-based opportunistic maintenance and spare parts provisioning policy in multiunit systems. Eur. J. Oper. Res. 262: 479–498, 2017.</big> 

<big>[19] J. Lee, M. Ghaffari, S. Elmeligy. Self-maintenance and engineering immune systems: Towards smarter machines and manufacturing systems. Annu. Rev. Control. 35: 111–122, 2011.</big>

<big>[20] R.P. Covert, G.C. Philip. An EOQ Model for Items with Weibull Distribution Deterioration. E Trans. 5: 323–326, 1973. </big>

<big>[21] D. Petrovic, R. Roy, R. Petrovic. Supply chain modelling using fuzzy sets. Int. J. Prod. Econ. 59: 443–453, 1999.</big>

<big>[22] Q. Hu, J.E. Boylan, H. Chen, A. Labib. OR in spare parts management: A review. Eur. J. Oper. Res. 266: 395–414, 2018.</big>

<big>[23] U.S. Rao. Properties of the Periodic Review (''R'', ''T'') Inventory Control Policy for Stationary, Stochastic Demand. Manuf. Serv. Oper. Manag. 5: 37–53, 2003.</big>

<big>[24] C.C. Sherbrooke, Metric: A Multi-Echelon Technique for Recoverable Item Control, Oper. Res. 16 (1968) 122–141.</big> 

<big>[25] S.C. Graves. A Multi-Echelon Inventory Model for a Repairable Item with One-for-One Replenishment. Manag. Sci. 31: 1247–1256, 1985. </big>

<big>[26] T.S. Vaughan, Failure replacement and preventive maintenance spare parts ordering policy, Eur. J. Oper. Res. 161 (2005) 183–190.</big>

<big>[27] Zhaoqiang Wang, Wenbin Wang, Changhua Hu, Xuejuan Liu. An integrated decision model for critical component spare parts ordering and condition based replacement with prognostic information. Chem. Eng. Trans. 33: 1063–1068, 2013.</big>

<big>[28] Z.-Q. Wang, W. Wang, C.-H. Hu, X.-S. Si, W. Zhang. A Prognostic-Information-Based Order-Replacement Policy for a Non-Repairable Critical System in Service. IEEE Trans. Reliab. 64: 721–735, 2015.</big>

<big>[29] Z. Wang, C. Hu, W. Wang, X. Kong, W. Zhang. A prognostics-based spare part ordering and system replacement policy for a deteriorating system subjected to a random lead time, Int. J. Prod. Res. 53: 4511–4527, 2015.</big>

<big>[30] K.-A. Nguyen, P. Do, A. Grall. Condition-based maintenance for multi-component systems using importance measure and predictive information. Int. J. Syst. Sci. Oper. Logist. 1: 228–245, 2014.</big>

<big>[31] R.S. Chhikara, L. Folks. The inverse Gaussian distribution: theory, methodology, and applications. M. Dekker. New York, November 1989.</big>

<big>[32] I.I. Gihman, A.V. Skorokhod. The Theory of Stochastic Processes I. Springer Berlin Heidelberg. Berlin, Heidelberg, 2004.</big>

<big>[33] B. Jorgensen, V. Seshadri. The Inverse Gaussian Distribution: A Case Study in Exponential Families. J. Am. Stat. Assoc. 91: 432, 1996.</big>

<big>[34] D. Acharya, G. Nagabhushanam, S.S. Alam. Jointly Optimal Block-Replacement and Spare Provisioning Policy. IEEE Trans. Reliab. 35: 447–451, 1986.</big>

<big>[35] Wen-Fung Leong, G.G. Yen. PSO-Based Multiobjective Optimization With Dynamic Population Size and Adaptive Local Archives. IEEE Trans. Syst. Man Cybern. Part B Cybern. 38: 1270–1293, 2008.</big>

<big>[36] A.C. Nearchou. The effect of various operators on the genetic search for large scheduling problems. Int. J. Prod. Econ. 88: 191–203, 2004.</big>

<big>[37] O.R. Isah, A.D. Usman, A.M.S. Tekanyi. A Hybrid Model of PSO Algorithm and Artificial Neural Network for Automatic Follicle Classification. Int. J. Bioautomation. 21: 43-58, 2017. </big>

<big>[38] H. Shao, G. Zheng. Boundedness and convergence of online gradient method with penalty and momentum. Neurocomputing. 74: 765–770, 2011.</big>

<big>[39] Y. Wang, Q. Shi, Improved Dynamic PSO-Based Algorithm for Critical Spare Parts Supply Optimization Under (''T'', ''S'') Inventory Policy, IEEE Access. 7: 153694–153709, 2019.</big>

<big>[40] M.J. Zuo, Renyan Jiang, R.C.M. Yam. Approaches for reliability modeling of continuous-state devices. IEEE Trans. Reliab. 48: 9–18, 1999.</big>

<big>[41] Luochun Xiao. Joint optimization of condition-based maintenance and spare parts supply for the same type of multi-component system. Nanjing University of Aeronautics and Astronautics. Nanjing, 2016.</big>

<big>[42] Fei Zhao, Xuejuan Liu. Joint optimization of periodic inspection and spare parts strategy considering imperfect maintenance. Theory and practice of system Engineering. 37: 3201–3214, 2017.</big>