# Scenario Reduction Method based on Output Performance for Condition-Based Maintenance Optimization.

1. IntroductionCondition-based maintenance (CBM) is a maintenance program that recommends maintenance decisions based on the information collected through condition monitoring [1]. CBM is one of the important maintenance methodologies for prognostics and health management (PHM) [2]. Since the degradation process for component\system is mainly nonlinear and non-stationary, forecasting of degradations is of uncertainty [3]. Therefore, it is essential to deal with CBM optimization incorporating stochastic prognostics [1]. Degradation forecasting is one of important prognostics for CBM optimization, which includes physical-based model [4], data-driven model [5, 6] and reliability model combing monitoring data [1,7].

For CBM incorporating with stochastic degradation, the current research work can be divided into two types in accordance of the planning horizon: infinite horizon and finite horizon [8]. CBM for a specified finite horizon can be more realistic than that for infinite horizon. In particular, the degradation possibility distribution can be dynamical and non-stationary [9]. Simulation-based CBM model can be applied to deal with this problem [10]. However, for the simulation-based CBM optimization, the computational complexity can increase in exponent with the number of monitoring conditions and stages. Therefore, it is essential to explore effective methods to reduce the computational complexity for CBM within the finite horizon. In other words, effective discretization methods for forecasting degradations can be considered a way to solve the bottle-neck in reducing computational complexity for CBM in finite horizon.

Uncertainty also exists in decision making problems for other application areas, as well as for CBM. Many computational methods for solving this kind of problem are based on approximating the original uncertainty distribution by discretization of the possibility model presenting uncertainty with finitely many scenarios [11]. It is natural to try to approximate the original scenarios set with much smaller subset that keeps key characteristics, and this process is called scenario reduction, as shown in Fig. 1. In general, fine discretizations can lead to numerically intractable problems, but coarse discretizations may lead to questionable model by neglecting important scenarios. Thus, scenario reduction has been an important topic in the past decades addressing uncertainty in decision making mainly for power system by stochastic programming [12-14].

The basic idea of scenario reduction is to select a few representative scenarios among a large number of scenarios (Fig. 1), and to assign probability to each selected scenario, so as to span possible output scenarios and capture wide range of outcomes in a way that is useful for decisionmaking [11]. According to the reduction mechanisms, current research of scenario reduction is mainly based on the input performance of uncertain parameters [13-15]. For scenario reduction based on input performance, the distance of scenarios is measured by the probability metric distance, and then the possibilities are redistributed among the selected scenarios by distance given in terms of mass transportation problems, to minimize the probability distance between initial and reduced scenarios. Heitsch and Romisch [15] proposed a kind of fast forward and backward reduction to reduce the computational complexity, applying the discretization of possibility distribution instead of the initial continuous probability distribution. A scenario optimal reduction technique, introduced by Dupacova et al [14], applied the Foret-Mourier distance and duality theory to compute the distance between two probability measures. Another variant of the scenario optimal reduction, introduced by De Oliveira et al. [13], combined with global and local (stage-wise) reduction, was to select a small set of sequences representing the stochastic process well enough.

These scenario reduction methods based on input performance focus mainly on the scenario parameters, but overlooking where the uncertainties appear in the problem modelling and their impacts on the optimal decision [16]. To deal with this problem, a novel scenario reduction method is proposed mainly based on the CBM model features. In other words, it is a kind of problem-dependent scenario reduction method, instead of isolating the scenario reduction from the problem.

The paper is organized as follows: Firstly, CBM optimization model and degradation model are introduced in Section 2 and 3, respectively. Secondly, in Section 4, we propose the scenario reduction based on the output-performance for forecasting degradations, aim to increase the computational efficiency of the CBM optimization. A case study based on degradation data from rotating machinery is presented in Section 5. Finally, conclusions and future works are outlined in Section 6.

2. Control-limit policy for CBM

For a single-unit system which degrades during the operation, the degradations can be detected by sensors at equally spaced discrete time epochs k, k=0, 1, 2...NT, with [DELTA]t being the constant time interval between two consecutive observations. The system is considered to fail whenever its degradation level reaches a given (i.e. not a decision variable) failure threshold [17]. The decision variable is preventive maintenance (PM) threshold [d.sup.i.sub.1] based on the monitoring degradation scenario. The expected maintenance cost rate including the costs entailed by corrective maintenance (CM), PM within the planning horizon is written as follows:

where [c.sub.F], [c.sub.P] are, respectively, CM and PM maintenance cost per time. We define [IF.sup.i.sub.t] and [IP.sub.t.sup.i] are, respectively, the indicator whether the system should be taken CM action, the indicator whether the system should be taken PM action, at time epoch t for degradation scenario i. The control-limit policy for CBM for a single-unit system is as follows [18]:

1) For each scenario i, if the degradations [z.sup.i.sub.t] exceeds the failure threshold, CM will be triggered. CM cost will be incurred and [IF.sup.i.sub.t] = 1. Otherwise, [IF.sup.i.sub.t] = 0.

2) PM actions will be triggered if the degradation is above PM threshold but it does not exceed the failure threshold. PM cost will be incurred as a result and [IF.sup.i.sub.t] = 1. Otherwise, [IF.sup.i.sub.t] = 0.

3) No maintenance actions are required if the degradation is below PM threshold. Therefore, no maintenance will be incurred, and both [IF.sup.i.sub.t] = 0 and [IP.sub.t.sup.i] = 0. For the CBM model, we assume that the system is as good as new after each maintenance action. As a result, the degradation of the system after each maintenance is equal to [W.sup.i.sub.1], if [IF.sup.i.sub.t], + [IP.sup.i.sub.t] > 1. Otherwise, the degradation will follow the forecasting model, i.e., [Z.sup.I.sub.t+1] = [W.sup.l.sub.t+1]. Furthermore, the degradation model can be referred to Section 3. For the control-limit policy of CBM within the finite horizon, the objective function is to minimize the expected maintenance cost rate for all degradation scenarios as follows:

[mathematical expression not reproducible] (2)

where [p.sub.i] is the possibility of each degradation scenario. From (1), (2) and equations in control-limit policy, optimal PM threshold can be obtained for each selected scenario.

3. Mathematical formulation for degradation forecasting model

Since exponential degradation model is common for rotating machinery and bearing application [9, 19], we use this type of degradation model which is expressed as follows:

[mathematical expression not reproducible], (3)

where [phi] is a known nonnegative constant, and [theta] is a lognormal random variable, [beta] is normal random variable, with mean [[micro].sub.1] and variance [[sigma].sup.2.sub.1] and [epsilon](t) is the error term which follows a Brownian with mean 0 and variance [[sigma].sup.2]. For simplification, the logarithm of degradation signal l(t) is formed as:

l(t) = ln(W(t)-[phi])=[theta]'+[beta]'t + [epsilon](t), (4)

where [theta]' and [beta]' are independent and follow normal distribution with mean [[micro].sub.0] and variance [[sigma].sup.2.sub.0], and with mean [[micro]'.sub.1] and variance [[sigma].sup.2.sub.1], respectively. Therefore, [theta]' and [beta]' can be expressed by [theta]' ~ N([[micro].sub.0], [[sigma].sub.0.sup.2]) and [beta]' ~ N([[micro]'.sub.1], [[sigma].sub.1.sup.2]), respectively. And [[micro]'.sub.1] = [[micro].sub.1] - [[sigma].sup.2] /2, [beta]' = [beta]- [sigma]2/2, [phi] is set to be 0.

PROPOSITION 1.The predictive distribution of the random variable L of stochastic degradation is normal with mean and variance [19]:

[micro](k, [I.sub.k]) = [l.sub.k] + [[micro]'.sub.[beta]] (k, [l.sub.k])[DELTA]t, (5)

[[sigma].sup.2](k) = [[sigma]'.sup.2.sub.[beta]](k)[DELTA][t.sup.2] + [[sigma].sup.2][DELTA]t. (6)

PROPOSITION 2.The posterior distribution of [beta]' , is normal with the following parameters [19]:

[mathematical expression not reproducible] (7)

[mathematical expression not reproducible] (8)

From Proposition 1 and 2, it can be seen that each time a signal is observed, it is used to update the predictive distribution of L. Therefore, it is a kind of non-stationary process for degradation, whose evolution is updated based on the monitoring condition for the system. For CBM, the optimal decision is made based on various forecasting degradation scenarios, since the degradation evolution process is stochastic within the planning horizon. Therefore, the initial scenarios can be generated through Proposition 1 and 2 by Monte Carlo simulation.

4. Scenario reduction method based on output performance

For CBM optimization combined with stochastic forecasting degradations, it is essential to make the scenario reduction to increase the calculation efficiency. Therefore, a kind of scenario reduction method based output performance is proposed, instead of isolating the scenario reduction from the problem. The procedures of this proposed method are shown in Fig. 2. It is based on the pattern recognition among the output performance (optimal PM threshold and the objective function value) and the forecasting degradation scenarios. Specific relationships can be derived from the analysis of the solutions for the CBM optimization with various degradation scenarios. The major procedure can be described by the following four steps.

Step 1. Monotonize the non-monotonic degradation processes.

Some degradation scenarios may be non-monotonic, and the non-monotonic degradation processes W(T) are monotonized to be the monotonic degradations Z(T) as follows:

Z(T) = max {w([tau]):[tau][less than or equal to]T}. (9)

Step 2. Identify feature points from the monotonic degradations.

1) For degradation scenario i, check whether it will exceed the failure threshold within the planning horizon. If it does, the first point which is [DELTA]Z higher than the failure threshold is set to be feature point, [Z.sub.1]. Otherwise, there is no feature point [Z.sub.1] for scenario i. [Z.sub.1] can be formulated as follows:

[T.sub.1] =min{[tau]:Z([tau]) + [DELTA]Z [greater than or equal to] [Z.sub.F]}, [Z.sub.1] = Z([T.sub.1]), (10)

where [DELTA]Z is the unit of discrete degradation evolution, and [T.sub.1] is the corresponding time of [Z.sub.1]. However, if there is no degradation over failure threshold, [T.sub.1] is set to be equal to NT.

2) For degradation scenario i, the first largest degradation point below the failure threshold [Z.sub.F] is set to be feature point [Z.sub.2], and alternatively, it can be formulated based on [T.sub.1] as follows:

[mathematical expression not reproducible] (11)

where [T.sub.2] is the corresponding time for feature point [Z.sub.2].

Meanwhile, N ([T.sub.2]) represents that the total number of PM

actions needed within the finite planning horizon if the first maintenance action is made at time [T.sub.2]. E.g. if [T.sub.2] = 0.4 NT, so it satisfies [T.sub.2] [member of] (NT/3, NT/2], then q = 2, and N ([T.sub.2]) = 2. That is, if the system is performed first maintenanceattime [T.sub.2] = 0.4 NT, then the second maintenance will be performed at time [T.sub.2] = 0.8 NT periodically and there will be two maintenance actions within the planning horizon. Moreover, it is assumed that the system is renewable after each maintenance action and it will follow the same degradation scenario as that before maintenance.

3) For degradation scenario i, search the first feature point [Z.sub.3] whose corresponding time [T.sub.3] satisfies that

[T.sub.3] = min{t: N (t) = N ([T.sub.2])]. Alternatively, [T.sub.3] can be formulated based on N ([T.sub.2]) as follows:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE], (12)

where Z ([T.sub.3]) is the degradation evolution process for scenario i during time period [T.sub.3]. Meanwhile, [T.sub.3] is round number.

4) For degradation scenario i, find the last feature point [Z.sub.4] which satisfies that [T.sub.4] = max{t: N(t) = N([T.sub.2])-1}. Alternatively, [Z.sub.4] can be formulated based on [T.sub.3] as follows:

[T.sub.4] = [T.sub.3] - 1,[Z.sub.4] = Z([T.sub.4]). (13)

Step 3. Identify characteristic patterns among solutions and scenario feature points.

After identifying these four feature points from the degradation scenarios, there are four kinds of modes can be indentified with different relationships among these feature points, as shown in Fig. 3.

For Mode 1, no degradations are over failure threshold, so there is no feature point [Z.sub.1]. The lower and upper bound of the optimal threshold are, respectively, [d.sub.1L] = [Z.sub.2] + [DELTA]Z and [d.sub.1U] = [Z.sub.F] - [DELTA]Z. Since no degradation points are over failure threshold, the number of PMs is zero within the planning horizon. For Mode 2, [Z.sub.1] exists and [Z.sub.2] = [Z.sub.3]. So, [d.sub.1L] = Z([T.sub.2] - 1) + [DELTA]Z and [d.sub.1U] = [Z.sub.2]. For Mode 3, [Z.sub.1] exists, [Z.sub.2] > [Z.sub.3] and [Z.sub.3] = [Z.sub.4]. So, [d.sub.1L] = [Z.sub.3] + [DELTA]Z and [d.sub.1U] = [Z.sub.2]. For Mode 4, [Z.sub.1] exists, [Z.sub.2] > [Z.sub.3] and [Z.sub.3] > [Z.sub.4]. So, [d.sub.1L] = [Z.sub.4] + [DELTA]Z and [d.sub.1U] = [Z.sub.2].

For Mode 2, 3 and 4, the number of PM is N([T.sub.2]) within the planning horizon.

Step 4. Objective function analysis.

From the relationship between modes and optimal solutions for different degradation scenarios, [N.sub.PM] can be determined by the specific modes of the degradation scenarios. The objective function is proportional to [N.sub.PM], since the objective is to minimize the total maintenance cost. As shown in Fig. 3, for Mode 2, 3, 4, the total number of PMs is determined by the time when the first PM occurs. That is, it is determined by position of the feature point [Z.sub.2] in the degradation scenarios. For each degradation scenario, the optimal threshold is an interval-value data, instead of single-value data. Therefore, the scenarios share the same objective function value can be reduced to one representative scenario for the proposed method.

After the scenario reduction by the proposed method, optimal PM thresholds for CBM can be solved based on the selected degradation scenarios.

5. Case study

5.1. Original data

For the application of proposed scenario reduction method based on output performance (short for PSR), we utilize the test data of vibration-based signals of the paper by Elwany, Gebraeel et al. [9]. According to the monitoring signals, exponential degradation model can be used to model the bearing degradation. The prior means, variances, and the variance parameter of the Brownian error terms are [[micro].sub.0] = -6.031, [[micro]'.sub.1] = 8.061 x [10.sup.-3], [[sigma].sup.2.sub.0] = 0.346, [[sigma].sup.2.sub.1] = 1.034 x [10.sup.-5], [[sigma].sup.2] = 0.0073 respectively. The failure threshold is defined by the root mean square (RMS) of all over acceleration, and it is set to be [Z.sub.F] = 0.025[V.sub.rms] (root mean square volts) [9]. Cost data in this application including [c.sub.F]=$120 and [c.sub.P]=$10. The planning horizon is NT = 300 min, and the time unit of degradations is set to be minute. For simplification, minute is the default time unit in section 5.2 and 5.3 unless otherwise stated.

5.2. Scenario generation for degradation forecasting

The models of randomness with their finite, discrete realizations are called scenario generators [20]. We generate [10.sup.-3] degradation scenarios according to the exponential degradation model as mentioned in Section 3. Firstly, the initial signal value is initialized ([l.sub.0] = log (0.001) = -6.9078). Secondly, the posterior distribution of [beta]' can be obtained based on the signal value at time t = 0. The forecasting signal value of t = 1 can be stochastically chosen. Then the posterior distribution of [beta]' of t = 2 can be obtained accordingly, and so forth. Then the forecasting degradation state, [Z.sub.t.sup.i], i = 1, 2,..., NT, can be simulated accordingly.

5.3. Results and discussions

The initial [10.sup.3] degradation scenarios are reduced to 10 representative scenarios by PSR (Fig. 4), according to the objective function values. For each kind of degradation scenarios, which share the same value of objective function value, are replaced by average degradation scenario as one representative. Meanwhile, the possibility of the new representative scenario is added by that of the same kind of degradation scenarios. According to the proposed scenario reduction method, the probability of the reduced scenarios can be added to the representative scenario.

In this section, one of scenario reduction method based on input performance, fast forward selection (short for FFS) [21], is chosen to make the comparison study to the proposed scenario reduction method, PSR. The initial degradation scenarios are reduced by FFS, and the relative possibility distance increases as the number of reduced scenarios decrease. The relative possibility distance is around 0.3 as the number of scenarios is reduced to 10. The corresponding possibility and optimal solutions of each selected scenario are listed in Table 1. Some of these 10 selected scenarios share the same number of PMs, instead of completely different. E.g. for scenario #426, #823, #910 and #1000, [N.sub.PM] = 1. Meanwhile, [N.sub.PM] for the rest of the selected scenarios are 0, 2, 3, 5, 6, and 8, respectively.

For the control-limit policy of CBM for single-component system, the optimal threshold is interval-value data, instead of single-value data. Therefore, the output performance of CBM can be analysed based on the objective function value (Q), optimal PM threshold ([d.sub.1]), the lower bound and upper bound for the optimal threshold ([d.sub.1L] and [d.sub.1U]).

The comparison results are shown in Fig. 5 and Table 2. Since objective function value is proportional to [N.sub.PM], the output performance of objective function value can be analysed by [N.sub.PM], as shown in Fig. 5, a. The probability distribution of [N.sub.PM] for selected scenarios (FFS, PSR) and initial scenarios (short for OPTI) are shown in Fig. 5, a. The probability distance is 0.276 between the probability distributions of [N.sub.PM] for selected scenarios by FFS and by OPTI. Meanwhile, the probability distance is 0.000 between the probability distributions of [N.sub.PM] for selected scenarios by PSR and by OPTI. Compared to FFS, the output performance, objective function value Q, of PSR is much closer to that of OPTI.

The probability distributions of [d.sub.1] for selected scenarios (FFS, PSR) and initial scenarios (OPTI) are shown in Fig. 5, b. The probability distance is 0.067 between the probability distributions of [d.sub.1] for selected scenarios by FFS and by OPTI. Meanwhile, the probability distance is 0.059 between the probability distributions of [d.sub.1] for selected scenarios by PSR and by OPTI. Compared to FFS, the output performance, [d.sub.1], of PSR is much closer to that of OPTI. The probability distributions of [d.sub.1U] for selected scenarios (FFS, PSR) and initial scenarios are shown in Fig. 5, d. The probability distance is 0.402 between the probability distributions of [d.sub.1U] for selected scenarios by FFS and by OPTI. Meanwhile, the probability distance is 0.396 between the probability distributions of [d.sub.1u] for selected scenarios by PRS and by OPTI. Compared to FFS, the output performance distance, [d.sub.1U], of PSR is much closer to that of OPTI.

However, the probability distributions of [d.sub.1L] for selected scenarios and initial scenarios are shown in Fig. 5,. The probability distance is 0.235 between the probability distributions of [d.sub.1L] for selected scenarios by FFS and by OPTI. Meanwhile, the probability distance is 0.336 between the probability distributions of [d.sub.1L] for selected scenarios by PRS and by OPTI. Compared to FFS, the output performance distance, [d.sub.1L], of PSR is a little bit greater to that of OPTI. It is the pay for the distance reduction of objective function value and [d.sub.1] by the proposed scenario reduction method, PSR.

From the comparison of the output performance for selected scenarios and initial scenarios, it shows that both the objection value (Q or [N.sub.PM]) and optimal threshold ([d.sub.1]) distribution distance between selected scenarios by PSR and initial scenarios are much closer than that between selected scenarios by FFS and initial scenarios. The reason is that the proposed scenario reduction method can consider the role of the uncertainty in the CBM model. Moreover, it can capture the characteristics, which have influence on the output performance. Therefore, the scenarios with same kind of features can be reduced to one representative scenario.

6. Conclusions

To improve the efficiency of scenario reduction for CBM problem involving uncertainties, such as stochastic degradation, a kind of scenario reduction method based on the output performance is proposed. This scenario reduction method can utilize the influence of the scenarios on their output performance, instead of input performance of scenarios. Therefore, the proposed method can increase the efficiency of scenario reduction so as to reduce computational complexity for CBM in finite horizon.

Further extension of this work will focus on investigating the feature modes for CBM policy of multi-component system that based on the proposed scenario reduction method. In addition, we assume that the state after maintenance is as good as mew. It is also possible to consider the case that the maintenance effect is partial repair or stochastic.

Acknowledgments

This work is supported by National Natural Science Foundation of China (Grant Nos. 71501148, 51275366, 51505349). The research work is also supported by National Science and Technology Major Project of China (Grant No. 2016ZX05038-002-LH002).

References

[1.] Jardine, A.K.S.; Lin, D.M.; Banjevic, D. 2006. A review on machinery diagnostics and prognostics implementing condition-based maintenance, Mechanical Systems and Signal Processing 20(7): 1483-1510. http://dx.doi.org/10.1016/jymssp.2005.09.012.

[2.] Lee, J.;Wu, F.; Zhao, W.; Ghaffari, M.; Liao, L.; Siegel, D. 2014. Prognostics and health management design for rotary machinery systems-Reviews, methodology and applications, Mechanical Systems and Signal Processing 42(1-2): 314-334. http://dx.doi.org/10.1016/jymssp.2013.06.004.

[3.] Gorjian, N.; Ma, L.; Mittinty, M.; Yarlagadda, P.; Sun, Y. 2009. A review on degradation models in reliability analysis. In Proceedings of the 4th World Congress on Engineering Asset Management. Athens. http://dx.doi.org/10.1007/978-0-85729-320-6_42.

[4.] Li, C.J.; Lee, H. 2005. Gear fatigue crack prognosis using embedded model, gear dynamic model and fracture mechanics, Mechanical systems and signal processing 19(4): 836-846. http://dx.doi.org/10.1016/j.ymssp.2004.06.007.

[5.] Sun, J.; Zuo, H.; Wang, W.; Pecht, M.G. 2014. Prognostics uncertainty reduction by fusing on-line monitoring data based on a state-space-based degradation model, Mechanical Systems and Signal Processing 45(2): 396-407. http://dx.doi.org/10.1016/jymssp.2013.08.022.

[6.] Benmounah, A. 2014. Optimization of maintenance intervals preventive and repair of gas turbines. Case of Algerian gas pipelines, Mechanika 20(3): 305-309. http://dx.doi.org/10.5755/j01.mech.20.3.7394.

[7.] Dragomir, O. E.; Gouriveau, R; Dragomir, F.; Minca, E.; Zerhouni, N. 2009. Review of prognostic problem in condition-based maintenance, European Control Conference, ECC'09., Budapest: Hungary. https://hal.archives-ouvertes.fr/hal-00418761.

[8.] Nakagawa, T.; Mizutani, S. 2009. A summary of maintenance policies for a finite interval, Reliability Engineering & System Safety 94(1): 89-96. http://dx.doi.org/10.1016/j.ress.2007.04.004.

[9.] Elwany, A.H.; Gebraeel, N.Z.; Maillart, L.M. 2011. Structured Replacement Policies for Components with Complex Degradation Processes and Dedicated Sensors, Operations Research 59(3): 684-695. http://dx.doi.org/10.1287/opre. 1110.0912.

[10.]Alrabghi, A.; Tiwari, A. 2015. State of the art in simulation-based optimisation for maintenance systems, Computers & Industrial Engineering 82: 167-182. http://dx.doi.org/10.1016/j.cie.2014.12.022.

[11.] Li, Z.; Floudas, C.A. 2014. Optimal scenario reduction framework based on distance of uncertainty distribution and output performance: I. Single reduction via mixed integer linear optimization Computers & Chemical Engineering 70(0): 50-66. http://dx.doi.org/10.1016/j.compchemeng.2014.03.019.

[12.] Feng, Y.; Ryan, S.M. 2013. Scenario construction and reduction applied to stochastic power generation expansion planning, Computers & Operations Research 40(1): 9-23. http://dx.doi.org/10.1016/j.cor.2012.05.005.

[13.] De Oliveira, W.L.; Sagastizabal, C; Penna, D.D.J.; Maceira, M.E.P.; Damazio, J.M. 2010. Optimal scenario tree reduction for stochastic streamflows in power generation planning problems, Optimisation Methods & Software 25(6): 917-936. http://dx.doi.org/10.1080/10556780903420135.

[14.] Dupacova, J.; Growe-Kuska, N.; Romisch, W. 2003. Scenario reduction in stochastic programming: an approach using probability metrics, Mathematical programming 95(3): 493-511. http://dx.doi.org/10.1007/s10107-002-0331-0.

[15.]Heitsch, H.; Romisch, W. 2009. Scenario tree modeling for multistage stochastic programs, Mathematical Programming 118(2): 371-406. http://dx.doi.org/10.1007/s10107-007-0197-2.

[16.] Sumaili, J.; Keko, H.; Miranda, V.; Zhou, Z.; Botterud, A.; Wang, J. 2011. Finding representative wind power scenarios and their probabilities for stochastic models. In Intelligent System Application to Power Systems (ISAP), 2011 16th International Conference on. IEEE. http://dx.doi.org/10.1109/ISAP.2011.6082195.

[17.] Zhao, X.; Fouladirad, M.; Berenguer, C; Bordes, L. 2010. Condition-based inspection/replacement policies for non-monotone deteriorating systems with environmental covariates, Reliability Engineering & System Safety 95(8): 921-934. http://dx.doi.org/10.1016/j.ress.2010.04.005.

[18.] Tian, Z.; Liao, H. 2011. Condition based maintenance optimization for multi-component systems using proportional hazards model, Reliability Engineering & System Safety 96(5): 581-589. http://dx.doi.org/10.1016/j.ress.2010.12.023.

[19.] Gebraeel, N.Z.; Lawley, M.A.; Li, R; Ryan, J.K 2005. Residual-life distributions from component degradation signals: A Bayesian approach HE Transactions 37(6): 543-557. http://dx.doi.org/10.1080/07408170590929018.

[20.] Di Domenica, N.; Mitra, G.; Valente, P.; Birbilis, G 2007. Stochastic programming and scenario generation within a simulation framework: an information systems perspective, Decision Support Systems 42(4): 2197-2218. http://dx.doi.org/10.1016/j.dss.2006.06.013.

[21.] Growe-Kuska, N.;Heitsch, H.; Romisch, W. 2003. Scenario reduction and scenario tree construction for power management problems. In 2003 IEEE Bologna Power Tech Conference Proceedings. http://dx.doi.org/10.1109/PTC.2003.1304379.

Summary

Condition-based maintenance (CBM) has been widely applied to maintenance policy. Due to the stochastic forecasting degradation, scenario reduction method has been developed to improve the efficiency of CBM. However, most existing scenario reduction methods focus mainly on the input performance of the forecasting degradation without considering the significant output performance characteristic based on the model.

In order to warrant the CBM optimization precision while reducing stochastic degradation scenarios efficiently, a new scenario reduction method is formulated that the scenarios with same objective function values can be reduced to one representative scenario. As a result, the reduced scenarios by the proposed method can maintain the probability distributions of objective values, while keeping optimal thresholds close to that of initial scenarios. Finally, the method is applied to select the representative degradation scenarios for CBM optimization model by utilizing vibration-based degradation signals from a rotating machinery application. Compared to the traditional scenario reduction method, the proposed method further improves accuracy and reduction efficiency of CBM optimization.

Keywords: condition-based maintenance, degradation, scenario reduction, output performance, probability, preventive maintenance, optimization.

Received July 01, 2016

Accepted October 13, 2017

Xinbo QIAN (*), Qiuhua TANG (**)

(*) Key Laboratory of Metallurgical Equipment and Control Technology of Ministry of Education, Hubei Key Laboratory of Mechanical Transmission and Manufacturing Engineering, Wuhan University of Science and Technology, Wuhan, China, 430081, E-mail: xinboqian@wust.edu.cn

(**) Department of Industrial Engineering, Wuhan University of Science and Technology, Wuhan, China, 430081, E-mail: tangqiuhua@wust. edu. cn

http://dx.doi.org/10.5755/j01.mech.23.5.15435

Table 1 Ten selected scenarios by FFS Scenario #96 #426 #823 #828 #848 Pi 0.001 0.094 0.211 0.001 0.079 NPM 5 1 1 8 2 Scenario #866 #910 #923 #990 #1000 Pi 0.017 0.249 0.033 0.286 0.029 NPM 3 1 6 0 1 Table 2 Output performance distance comparison Comparison for [N.sub.PM] [d.sub.1] [d.sub.1U] [d.sub.1L] OPTI -FFS 0.276 0.067 0.402 0.235 OPTI-PSR 0.000 0.059 0.396 0.336

Printer friendly Cite/link Email Feedback | |

Author: | Qian, Xinbo; Tang, Qiuhua |
---|---|

Publication: | Mechanika |

Date: | Sep 1, 2017 |

Words: | 4909 |

Previous Article: | Study on Clearance Optimization of All-metal Screw Pumps: Experiment and Simulation. |

Next Article: | A Numerical Study on Performance of Dental Air Turbine Handpieces. |