Multi-Objective Optimum Design and Maintenance of Safety Systems: An In-Depth Comparison Study Including Encoding and Scheduling Aspects with NSGA-II
Abstract
:1. Introduction
- In this work, seven encoding alternatives are thoroughly explored: Real encoding (with Simulated Binary Crossover), Binary encoding (with 1 point Crossover), Binary encoding (with 2 point Crossover), Binary encoding (with Uniform Crossover), Gray encoding (with 1 point Crossover), Gray encoding (with 2 point Crossover) and Gray encoding (with Uniform Crossover). Their performances are compared using the Hypervolume indicator and statistical significance tests.
- Additionally, three accuracy levels on time are explored for the binary encoding; hours, days and weeks, in order to analyse the effect of chromosome length in the evolutionary search and final non-dominated set of solutions. Their performances are compared using the Hypervolume indicator and statistical significance tests.
- The methodology is applied in an industrial test case, obtaining an improved non-dominated front of optimum solutions that could be considered as both a benchmark case and reference solution.
2. Literature Review
2.1. Redundancy Allocation of Reliability Systems Design Optimisation
2.2. Preventive Maintenance Strategy Optimisation
2.3. Redundancy Allocation and Preventive Strategy Optimisation
3. Methodology and Description of the Model
3.1. Extracting Availability and Cost from Functionability Profiles
- n is the total number of operation times,
- is the i-th operation time in hours (Time To Failure or Time To Start a Preventive Maintenance activity),
- m is the total number of recovery times,
- is the j-th recovery time in hours (Time To Repair or Time To Perform a Preventive Maintenance activity).
- C is the system operation cost quantified in economic units,
- q is the total number of corrective maintenance activities,
- is the cost due to the i-th corrective maintenance activity,
- p is the total number of preventive maintenance activities,
- is the cost due to the j-th preventive maintenance activity.
3.2. Building Functionability Profiles to Evaluate the Objective Functions
- System mission time must be defined and then, the process continues for all devices.
- The device Functionability Profile (PF) must be initialised.
- The Time To Start a Preventive Maintenance activity (TP) proposed by the Multi-objective Evolutionary Algorithm is extracted from the individual (candidate solution) that is being evaluated and a Time To Perform a Preventive Maintenance activity (TRP) is randomly generated, within the limits previously fixed.
- With reference to the failure probability density function related to the device, a Time To Failure (TF) is randomly generated, within the limits previously fixed.
- If TP < TF, the preventive maintenance activity is performed before the failure. In this case, as many operating times units as TP units followed by as many recovery times units as TRP units are added to the device Functionability Profile. Each time unit represented in this way (both operating times and recovery times) is equivalent to one hour, day or week of real time, depending on the accuracy level chosen.
- If TP > TF, the failure occurs before carrying out the preventive maintenance activity. In this case, taking into consideration the repair probability density function related to the device, the Time To Repair after the failure (TR) is randomly generated, within the limits previously fixed. Then, as many operating times units as TF units followed by as many recovery times units as TR units are added to the device Functionability Profile. Each time unit represented in this way (both operating times and recovery times) is equivalent to one hour, day or week of real time, depending on the accuracy level chosen.
- Steps 4 to 6 have to be repeated until the end of the device mission time.
- Steps 2 to 7 have to be repeated until the construction of the Functionability Profiles of all devices.
- After building the Functionability Profiles of the devices, the system Functionability Profile is built by referring to the series (AND) or parallel (OR) distribution of the system devices.
3.3. Multi-Objective Optimisation Approach
4. The Case Study
- The number of redundant devices is limited as shown in Figure 1,
- two states are considered for each device: operation or failed state,
- the devices are independent of each other,
- a repair starts just after the failure of the device,
- a repair returns the device to the as-good-as-new state.
- It is necessary to establish the optimum period to perform a preventive maintenance activity for the system devices, and
- It is necessary to decide whether to include redundant devices such as P2 and/or V4 by evaluating design alternatives. Including redundant devices will improve the system Availability but it will also increase the system operation Cost.
5. Description of the Experiments Carried Out
5.1. Comparing Encodings
- Real encoding: This is formed by strings of real numbers (with 0 as a minimum value and 1 as a maximum value) following the shape []. The presence of redundant devices, P2 and V4, is defined by and , respectively, and the optimum Time To Start a Preventive Maintenance activity in relation to each system device is represented by to . The values of the decision variables must be scaled and rounded, i.e., and are rounded to the nearest integer (0 implies that the respective device is not selected whereas 1 implies the opposite). to are scaled among the respective and (depending on the type of device) and rounded to the nearest integer using Equation (5), where TP is the true value of the Time To Start a Preventive Maintenance activity, measured in hours (e.g., the decision variable represents the Time To Start a Preventive Maintenance for the valve V1 whose and has a value of 8760 h and 35,040 h, respectively. If the value of the decision variable is 0.532, the value of the Time To Start a Preventive Maintenance activity will be 8760 + 0.532 · (35,040 − 8760) ≈ 22,741) h.
- Binary encoding: This is formed by strings of binary numbers that vary between 0 and 1, where the total number of bits is 103 and they are:
- : This denotes the presence of the pump P2 in the system design (0 implies that the respective device is not selected whereas 1 implies the opposite).
- : This denotes the presence of the valve V4 in the system design (0 implies that the respective device is not selected whereas 1 implies the opposite).
- to : These denote the Time To Start a Preventive Maintenance activity to the valve V1. A binary scale that allows representation of the numbers from to is needed. has a value of 8760 h and has a value of 35,040 h so steps needed, where the step 0 represents a time of 8760 h and the step 26,279 represents a time of 35,040 h. A binary scale with at least 26,280 steps involves using 15 bits (as 26,280 steps lies between and . Since 26,280 steps are needed and 32,768 are possible on the scale, an equivalent relation must be used. Each step in the scale of 32,768 steps represents steps in the scale of 26,768 steps. Therefore, it is possible to calculate the true Time To Start a Preventive Maintenance activity (in hours) using the scale change shown by Equation (6), where B represents the decimal value of the binary string to (e.g., if the values of the decision variables in binary encoding are 1 0 1 1 0 1 1 0 0 0 1 1 1 0 1, the decimal value in the scale of 32,768 steps will be 23,325. If 26,768 steps are scaled, the number achieved is steps. Therefore, the true Time To Start a Preventive Maintenance activity amounts to 18,707 + 8760 = 24,467 h).
- to : These denote the Time To Start a Preventive Maintenance activity to the pump P2. A binary scale that allows representation of the numbers from to in needed. has a value of 2920 h and has a value of 8760 h so steps needed, where the step 0 represents the time of 2920 h and the step 5839 represents the time of 8760 h. A binary scale with at least 5840 steps involves using 13 bits (as 5840 steps lies between and . Since 5840 steps are needed and 8142 are possible on the scale, an equivalent relation must be used. Each step in the scale of 8142 represents steps on a scale of 5840 steps. Therefore, it is possible to calculate the true Time To Start a Preventive Maintenance activity (hours) using the scale change shown by Equation (7), where B represents the decimal value of the binary string to (e.g., if the values of the decision variables in binary encoding are 1 0 1 1 0 1 1 0 0 0 1 1 1, the value on a scale of 8192 steps will be 5831. If 5840 steps are scaled, the number achieved is steps. Therefore, the true Time To Start a Preventive Maintenance activity amounts to 4157 + 2920 = 7077 h).
- to : These denote the Time To Start a Preventive Maintenance activity to the pump P3. The behaviour of its encoding is similar to the behaviour explained for the pump P2.
- to : These denote the Time To Start a Preventive Maintenance activity to the valve V4. The behaviour of its encoding is similar to the behaviour explained for the valve V1.
- to : These denote the Time To Start a Preventive Maintenance activity to the valve V5. The behaviour of its encoding is similar to the behaviour explained for the valve V1.
- to : These denote the Time To Start a Preventive Maintenance activity to the valve V6. The behaviour of its encoding is similar to the behaviour explained for the valve V1.
- to : These denote the Time To Start a Preventive Maintenance activity to the valve V7. The behaviour of its encoding is similar to the behaviour explained for the valve V1.
- Gray encoding: When binary encoding is used, close numbers could bring big scheme modifications (e.g., the code for 15 is 0 1 1 1 while the code for 16 is 1 0 0 0, which represents changes in four (all) bits). Conversely, very similar numbers can represent numbers that are very apart (e.g., the code for 0 is 0 0 0 0 while the code for 16 is 1 0 0 0). In addition to standard binary encoding, Gray encoding is used. Gray code is a binary code where the difference among neighbouring numbers always differs by one bit [64,65,66].
5.2. Comparing Accuracies
- Binary encoding—Days: This is formed by strings of binary numbers that vary between 0 and 1, where the total number of bits is 73 and they are:
- : This denotes the presence of the pump P2 in the system design (0 implies that the device is not selected whereas 1 implies the opposite).
- : This denotes the presence of the valve V4 in the system design (0 implies that the device is not selected whereas 1 implies the opposite).
- to : These denote the Time To Start a Preventive Maintenance activity to the valve V1. A binary scale that allows representation of the numbers from to expressed in days as a time unit is needed. has a value of 8760 h (equivalent to 365 days) and has a value of 35,040 h (equivalent to 1460 days) so steps needed, where the step 0 represents the time of 365 days and the step 1094 represents the time of 1460 days. A binary scale with at least 1095 steps involves using 11 bits (due to the fact that 1095 steps lies between and ). Since 1095 steps are needed and 2048 are possible on the scale, an equivalent relation must be used. Each step on the scale of 2048 steps represents steps on a scale of 1095 steps. Therefore, it is possible to achieve the true Time To Start a Preventive Maintenance activity (days) using the scale change shown by Equation (8), where B represents the decimal value of the binary string to (e.g., if the values of the decision variables in binary encoding are 1 0 1 1 0 1 1 0 0 0 1, the decimal value on the scale of 2048 steps will be 1457. Scaling to a scale of 1095 steps, the number achieved is steps. Therefore, the true Time To Start a Preventive Maintenance activity amounts to 779 + 365 = 1144 days).
- to : These denote the Time To Start a Preventive Maintenance activity to the pump P2. A binary scale that allows representation of the numbers from to expressed in days as a time unit is needed. has a value of 2920 h (equivalent to 122 days) and has a value of 8760 (equivalent to 365 days) so steps needed, where the step 0 represents the time of 122 days and the step 242 represents the time of 365 days. A binary scale with at least 243 steps involves using 8 bits (as 243 steps lies between and . Since 243 steps are needed and 256 are possible on the scale, an equivalent relationship must be used. Each step in the scale of 256 represents steps in the scale of 243 steps. Therefore, it is possible to achieve the true Time To Start a Preventive Maintenance activity (days) using the scale change shown by Equation (9), where B represents the decimal value of the binary string to (e.g., if the values of the decision variables in binary encoding are 1 0 1 1 0 1 1 0, the value on the scale of 256 steps will be 182). Scaling to a scale of 243 steps, the number achieved is steps. Therefore, the true Time To Start a Preventive Maintenance activity amounts to 173 + 122 = 295 days).
- to : These denote the Time To Start a Preventive Maintenance activity to the pump P3. The behaviour of its encoding is similar to the behaviour explained for the pump P2.
- to : These denote the Time To Start a Preventive Maintenance activity to the valve V4. The behaviour of its encoding is similar to the behaviour explained for the valve V1.
- to : These denote the Time To Start a Preventive Maintenance activity to the valve V5. The behaviour of its encoding is similar to the behaviour explained for the valve V1.
- to : These denote the Time To Start a Preventive Maintenance activity to the valve V6. The behaviour of its encoding is similar to the behaviour explained for the valve V1.
- to : These denote the Time To Start a Preventive Maintenance activity to the valve V7. The behaviour of its encoding is similar to the behaviour explained for the valve V1.
- Binary encoding—Weeks: It is formed by strings of binary numbers that vary between 0 and 1, where the total number of bits is 54 and they are:
- : This denotes the presence of the pump P2 in the system design (0 implies that the device is not selected whereas 1 implies the opposite).
- : This denotes the presence of the valve V4 in the system design (0 implies that the device is not selected whereas 1 implies the opposite).
- to : These denote the Time To Start a Preventive Maintenance activity to the valve V1. A binary scale that allows representation of the numbers from to expressed in weeks as a time unit is needed. has a value of 8760 h (equivalent to 52 weeks) and has a value of 35,040 h (equivalent to 209 weeks) so steps needed, where step 0 represents a time of 52 weeks and step 156 represents a time of 209 weeks. A binary scale with at least 157 steps involves using 8 bits (as 157 steps lies between and ). Since 157 steps are needed and 256 are possible on the scale, an equivalent relationship must be used. Each step on a 256-steps scale represents steps on the 157-steps scale. Therefore, it is possible to achieve the true Time To Start a Preventive Maintenance activity (weeks) using the scale change shown by Equation (10), where B represents the decimal value of the binary string to (e.g., if the values of the decision variables in binary encoding are 1 0 1 1 0 1 1 0, the decimal value in the scale of 256 steps will be 182. Working with a scale of 157 steps, the number achieved is steps. Therefore, the true Time To Start a Preventive Maintenance activity amounts to 112 + 52 = 164 weeks).
- to : These denote the Time To Start a Preventive Maintenance activity to the pump P2. A binary scale that allows representation of the numbers from to expressed in weeks as a time unit is needed. has a value of 2920 h (equivalent to 17 weeks) and has a value of 8760 (equivalent to 52 weeks) so steps needed, where step 0 represents the time of 17 weeks and step 34 represents the time of 52 weeks. A binary scale with at least 35 steps involves using 6 bits (as 35 steps lies between and . Since 35 steps are needed and 64 are possible on the scale, an equivalent relationship must be used. Each step on the scale of 64-steps scale represents steps in the 35-steps scale. Therefore, it is possible to achieve the true Time To Start a Preventive Maintenance activity (weeks) using the scale change shown by Equation (11), where B represents the decimal value of the binary string to (e.g., if the values of the decision variables in binary encoding are 1 0 1 1 0 1, the value in the scale of 64 steps will be 45. Scaling on the 35-steps scale, the number achieved is steps. Therefore, the true Time To Start a Preventive Maintenance activity amounts to 45 + 17 = 62 weeks).
- to : These denote the Time To Start a Preventive Maintenance activity to the pump P3. The behaviour of its encoding is similar to the behaviour explained for the pump P2.
- to : These denote the Time To Start a Preventive Maintenance activity to the valve V4. The behaviour of its encoding is similar to the behaviour explained for the valve V1.
- to : These denote the Time To Start a Preventive Maintenance activity to the valve V5. The behaviour of its encoding is similar to the behaviour explained for the valve V1.
- to : These denote the Time To Start a Preventive Maintenance activity to the valve V6. The behaviour of its encoding is similar to the behaviour explained for the valve V1.
- to : These denote the Time To Start a Preventive Maintenance activity to the valve V7. The behaviour of its encoding is similar to the behaviour explained for the valve V1.
5.3. NSGA-II Configuration
- Crossover: Type of crossover during the evolutionary process. The Simulated Binary Crossover (SBX) is used for real encoding while one point (1PX), two point (2PX) and uniform crossover (UX) are used for binary and Gray encodings.
- Population size (N): The population sizes used are 50, 100 and 150 individuals.
- Mutation Probability (PrM): This is the expectation of the number of genes mutating. The central value is equivalent to 1/decision variables (for the case study, the number of decision variables is 9 using real encoding, 103 using standard and Gray encoding with the hour as a time unit, 73 using standard binary encoding with the day as a time unit and 54 using standard binary encoding with the week as a time unit). Two more probabilities, one above and the other below the central value (1.5/decision variables and 0.5/decision variables, respectively) are set.
- Mutation Distribution (disM): This is the mutation distribution index when the Simulated Binary Crossover is used. It is set to the common value of 20.
- Crossover Probability (PrC): The probability of applying a crossover operator is set to 1 in all cases.
- Crossover Distribution (disC): The crossover distribution index when the Simulated Binary Crossover is used. It is set to the common value of 20.
- The scale factor used to compute the Cost was 1700 economic units.
- The scale factor used to compute the system Unavailability was 0.003.
6. Results
6.1. Encoding Experiment
6.1.1. Real Encoding
6.1.2. Standard Binary Encoding
6.1.3. Gray Encoding
6.1.4. Comparing Encoding Configurations
6.2. Accuracy Experiment
6.2.1. Standard Binary Encoding (Days)
6.2.2. Standard Binary Encoding (Weeks)
6.2.3. Comparing Accuracy-Level Configurations
6.3. Accumulated Non-Dominated Set of Designs
7. Conclusions
Author Contributions
Funding
Institutional Review Board Statement
Informed Consent Statement
Data Availability Statement
Acknowledgments
Conflicts of Interest
Appendix A. Detailed Information of the Case Study Parameters
- Life Cycle. System mission time, expressed in hours.
- Corrective Maintenance Cost. The cost involved in developing a repair activity to recover the system following a failure, expressed in economic units per hour.
- Preventive Maintenance Cost. The cost involved in developing a Preventive Maintenance activity, expressed in relation to the Corrective Maintenance Cost.
- Pump . Minimum operation Time To Failure for a pump without Preventive Maintenance, expressed in hours.
- Pump . Maximum operation Time To Failure for a pump without Preventive Maintenance, expressed in hours.
- Pump . Failure rate for a pump, which follows an exponential failure distribution, expressed in hours raised to the power of minus six.
- Pump . Minimum Time To Repair or duration of a Corrective Maintenance activity for a pump, expressed in hours.
- Pump . Maximum Time To Repair or duration of a Corrective Maintenance activity for a pump, expressed in hours.
- Pump . Mean for the normal distribution followed for the Time To Repair assumed for a pump, expressed in hours.
- Pump . Standard deviation for the normal distribution followed for the Time To Repair assumed for a pump, expressed in hours.
- Pump . Minimum operation Time To Start a scheduled Preventive Maintenance activity for a pump, expressed in hours.
- Pump . Maximum operation Time To Start a scheduled Preventive Maintenance activity for a pump, expressed in hours.
- Pump . Minimum Time To Perform a Preventive Maintenance activity for a pump, expressed in hours.
- Pump . Maximum Time To Perform a Preventive Maintenance activity for a pump, expressed in hours.
- Valve . Minimum operation Time To Failure for a valve without Preventive Maintenance, expressed in hours.
- Valve . Maximum operation Time To Failure for a valve without Preventive Maintenance, expressed in hours.
- Valve . Failure rate for a valve, which follows an exponential failure distribution, expressed in hours raised to the power of minus six.
- Valve . Minimum Time To Repair or duration of a Corrective Maintenance activity for a valve, expressed in hours.
- Valve . Maximum Time To Repair or duration of a Corrective Maintenance activity for a valve, expressed in hours.
- Valve . Mean for the normal distribution followed for the Time To Repair assumed for a valve, expressed in hours.
- Valve . Standard deviation for the normal distribution followed for the Time To Repair assumed for a valve, expressed in hours.
- Valve . Minimum operation Time To Start a scheduled Preventive Maintenance activity for a valve, expressed in hours.
- Valve . Maximum operation Time To Start a scheduled Preventive Maintenance activity for a valve, expressed in hours.
- Valve . Minimum Time To Perform a Preventive Maintenance activity for a valve, expressed in hours.
- Valve . Maximum Time To Perform a Preventive Maintenance activity for a valve, expressed in hours.
References
- Misra, K.B. Reliability Engineering: A Perspective. In Handbook of Performability Engineering; Misra, K.B., Ed.; Springer Limited: London, UK, 2008; Chapter 19; pp. 253–289. ISBN 978-1-84800-130-5. [Google Scholar]
- Andrews, J.D.; Moss, T.R. Reliability and Risk Assessment, 2nd ed.; ASME Press: New York, NY, USA, 2002; p. 540. ISBN 0-7918-0183-7. [Google Scholar]
- Simon, D. Evolutionary Optimization Algorithms. Biologically Inspired and Population-Based Approaches to Computer Intelligence; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2013; p. 742. ISBN 978-0-470-93741-9. [Google Scholar]
- Kuo, W.; Prasad, R.; Tillman, F.A.; Mwang, C.L. Optimal Reliability Design: Fundamentals and Applications, 1st ed.; Cambridge University Press: Cambridge, UK, 2001; p. 389. ISBN 978-0-521-78127-5. [Google Scholar]
- Gao, Y.; Feng, Y.; Zhang, Z.; Tan, J. An optimal dynamic interval preventive maintenance scheduling for series systems. Reliab. Eng. Syst. Saf. 2015, 142, 19–30. [Google Scholar] [CrossRef]
- Cacereño, A.; Greiner, D.; Galván, B. Solving Multi-objective Optimal Design and Maintenance for Systems Based on Calendar Times Using NSGA-II. In Advances in Evolutionary and Deterministic Methods for Design, Optimization and Control in Engineering and Sciences; Gaspar-Cunha, A., Periaux, J., Giannakoglou, K.C., Gauger, N.R., Quagliarella, D., Greiner, D., Eds.; Springer Nature: Cham, Switzerland, 2021; pp. 245–259. [Google Scholar]
- Knezevic, J. Mantenibilidad, 1st ed.; Isdefe: Madrid, Spain, 1996; p. 209. ISBN 84-89338-08-6. [Google Scholar]
- Deb, K.; Pratap, A.; Agarwal, S.; Meyarivan, T. A Fast and Elitist Multiobjective Genetic Algorithm: NSGA-II. IEEE Trans. Evol. Comput. 2002, 6, 182–197. [Google Scholar] [CrossRef] [Green Version]
- Busacca, P.G.; Marseguerra, M.; Zio, E. Multiobjective optimization by genetic algorithms: Application to safety systems. Reliab. Eng. Syst. Saf. 2001, 72, 59–74. [Google Scholar] [CrossRef]
- Marseguerra, M.; Zio, E.; Podofillini, L.; Coit, D.W. Optimal Design of Reliable Network Systems in Presence of Uncertainty. IEEE Trans. Reliab. 2005, 54, 243–253. [Google Scholar] [CrossRef]
- Tian, Z.; Zuo, M.J. Redundancy allocation for multi-state systems using physical programming and genetic algorithms. Reliab. Eng. Syst. Saf. 2006, 91, 1049–1056. [Google Scholar] [CrossRef]
- Huang, H.Z.; Qu, J.; Zuo, M.J. Genetic-algorithm-based optimal apportionment of reliability and redundancy under multiple objectives. IIE Trans. 2009, 41, 287–298. [Google Scholar] [CrossRef]
- Zoulfaghari, H.; Hamadani, A.Z.; Ardakan, M.A. Bi-objective redundancy allocation problem for a system with mixed repairable and non-repairable components. ISA Trans. 2014, 53, 17–24. [Google Scholar] [CrossRef]
- Taboada, H.A.; Espiritu, J.F.; Coit, D.W. MOMS-GA: A Multi-Objective Multi-State Genetic Algorithm for System Reliability Optimization Design Problems. IEEE Trans. Reliab. 2008, 57, 182–191. [Google Scholar] [CrossRef]
- Taboada, H.A.; Baheranwala, F.; Coit, W.D.; Wattanapongsakorn, N. Practical solutions for multi-objective optimization: An application to system reliability design problems. Reliab. Eng. Syst. Saf. 2007, 92, 314–322. [Google Scholar] [CrossRef]
- Greiner, D.; Galván, B.; Winter, G. Safety Systems Optimum Design by Multicriteria Evolutionary Algorithms. In Evolutionary Multi-Criterion Optimization. Lecture Notes in Computer Science, Proceedings of the Second International Conference, EMO 2003, Faro, Portugal, 8–11 April 2003; Fonseca, C.M., Fleming, P.J., Zitzler, E., Deb, K., Thiele, L., Eds.; Springer: Berlin, Germany, 2003; pp. 722–736. [Google Scholar]
- Salazar, D.; Rocco, C.M.; Galván, B.J. Optimization of constrained multiple-objective reliability problems using evolutionary algorithms. Reliab. Eng. Syst. Saf. 2006, 91, 1057–1070. [Google Scholar] [CrossRef]
- Limbourg, P.; Kochs, H.D. Multi-objective optimization of generalized reliability design problems using feature models—A concept for early design stages. Reliab. Eng. Syst. Saf. 2008, 93, 815–828. [Google Scholar] [CrossRef]
- Kumar, R.; Izui, K.; Yoshimura, M.; Nishiwaki, S. Multi-objective hierarchical genetic algorithms for multilevel redundancy allocation optimization. Reliab. Eng. Syst. Saf. 2009, 94, 891–904. [Google Scholar] [CrossRef]
- Chambari, A.; Rahmati, S.H.A.; Najafi, A.A.; Karimi, A. A bi-objective model to optimize reliability and cost of system with a choice of redundancy strategies. Comput. Ind. Eng. 2012, 63, 109–119. [Google Scholar] [CrossRef]
- Safari, J. Multi-objective reliability optimization of series-parallel systems with a choice of redundancy strategies. Reliab. Eng. Syst. Saf. 2012, 108, 10–20. [Google Scholar] [CrossRef]
- Ardakan, M.A.; Hamadani, A.Z.; Alinaghian, M. Optimizing bi-objective redundancy allocation problem with a mixed redundancy strategy. ISA Trans. 2015, 55, 116–128. [Google Scholar] [CrossRef]
- Ghorabaee, M.K.; Amiri, M.; Azimi, P. Genetic algorithm for solving bi-objective redundancy allocation problem with k-out-of-n subsystems. Appl. Math. Model. 2015, 39, 6396–6409. [Google Scholar] [CrossRef]
- Amiri, M.; Khajeh, M. Developing a bi-objective optimization model for solving the availability allocation problem in repairable series–parallel systems by NSGA-II. J. Ind. Eng. Int. 2016, 12, 61–69. [Google Scholar] [CrossRef] [Green Version]
- Jahromi, A.E.; Feizabadi, M. Optimization of multi-objective redundancy allocation problem with non-homogeneous components. Comput. Ind. Eng. 2017, 108, 111–123. [Google Scholar] [CrossRef]
- Kayedpour, F.; Amiri, M.; Rafizadeh, M.; Nia, A.S. Multi-objective redundancy allocation problem for a system with repairable components considering instantaneous availability and strategy selection. Reliab. Eng. Syst. Saf. 2017, 160. [Google Scholar] [CrossRef]
- Sharifi, M.; Moghaddam, T.A.; Shahriari, M. Multi-objective Redundancy Allocation Problem with weighted-k-out-of-n subsystems. Heliyon 2019, 5, e02346. [Google Scholar] [CrossRef] [Green Version]
- Chambari, A.; Azimi, P.; Najafi, A.A. A bi-objective simulation-based optimization algorithm for redundancy allocation problem in series-parallel systems. Expert Syst. Appl. 2021, 173, 114745. [Google Scholar] [CrossRef]
- Zhao, J.H.; Liu, Z.; Dao, M.T. Reliability optimization using multiobjective ant colony system approaches. Reliab. Eng. Syst. Saf. 2007, 92, 109–120. [Google Scholar] [CrossRef]
- Chiang, C.H.; Chen, L.H. Availability allocation and multi-objective optimization for parallel–series systems. Eur. J. Oper. Res. 2007, 180, 1231–1244. [Google Scholar] [CrossRef]
- Elegbede, C.; Adjallah, K. Availability allocation to repairable systems with genetic algorithms: A multi-objective formulation. Reliab. Eng. Syst. Saf. 2003, 82, 319–330. [Google Scholar] [CrossRef]
- Khalili-Damghani, K.; Abtahi, A.R.; Tavana, M. A new multi-objective particle swarm optimization method for solving reliability redundancy allocation problems. Reliab. Eng. Syst. Saf. 2013, 111, 58–75. [Google Scholar] [CrossRef]
- Jiansheng, G.; Zutong, W.; Mingfa, Z.; Ying, W. Uncertain multiobjective redundancy allocation problem of repairable systems based on artificial bee colony algorithm. Chin. J. Aeronaut. 2014, 27, 1477–1487. [Google Scholar]
- Samanta, A.; Basu, K. An attraction based particle swarm optimization for solving multi-objective availability allocation problem under uncertain environment. J. Intell. Fuzzy Syst. 2018, 35, 1169–1178. [Google Scholar] [CrossRef]
- Muñoz, A.; Martorell, S.; Serradell, V. Genetic algorithms in optimizing surveillance and maintenance of components. Reliab. Eng. Syst. Saf. 1997, 57, 107–120. [Google Scholar] [CrossRef]
- Marseguerra, M.; Zio, E.; Podofillini, L. Condition-based maintenance optimization by means of genetic algorithms and Monte Carlo simulation. Reliab. Eng. Syst. Saf. 2002, 77, 151–166. [Google Scholar] [CrossRef]
- Gao, J.; Gen, M.; Sun, L. Scheduling jobs and maintenances in flexible job shop with a hybrid genetic algorithm. J. Intell. Manuf. 2006, 17, 493–507. [Google Scholar] [CrossRef]
- Sánchez, A.; Carlos, S.; Martorell, S.; Villanueva, J.F. Addressing imperfect maintenance modelling uncertainty in unavailability and cost based optimization. Reliab. Eng. Syst. Saf. 2009, 94, 22–32. [Google Scholar] [CrossRef]
- Wang, Y.; Pham, H. A multi-objective optimization of imperfect preventive maintenance policy for dependent competing risk systems with hidden failure. IEEE Trans. Reliab. 2011, 60, 770–781. [Google Scholar] [CrossRef]
- Ben Ali, M.; Sassi, M.; Gossab, M.; Harrath, Y. Simultaneous scheduling of production and maintenance tasks in the job shop. J. Intell. Manuf. 2011, 49. [Google Scholar] [CrossRef]
- An, Y.; Chen, X.; Zhang, J.; Li, Y. A hybrid multi-objective evolutionary algorithm to integrate optimization of the production scheduling and imperfect cutting tool maintenance considering total energy consumption. J. Clean Prod. 2020, 268, 121540. [Google Scholar] [CrossRef]
- Bressi, S.; Santos, J.; Losa, M. Optimization of maintenance strategies for railway track-bed considering probabilistic degradation models and different reliability levels. Reliab. Eng. Syst. Saf. 2021, 207, 107359. [Google Scholar] [CrossRef]
- Martorell, S.; Villanueva, J.F.; Carlos, S.; Nebot, Y.; Sánchez, A.; Pitarch, J.L.; Serradell, V. RAMS+C informed decision-making with application to multi-objective optimization of technical specifications and maintenance using genetic algorithms. Reliab. Eng. Syst. Saf. 2005, 87, 65–75. [Google Scholar] [CrossRef]
- Oyarbide-Zubillaga, A.; Goti, A.; Sanchez, A. Preventive maintenance optimisation of multi-equipment manufacturing systems by combining discrete event simulation and multi-objective evolutionary algorithms. Prod. Plan. Control. 2008, 19, 342–355. [Google Scholar] [CrossRef]
- Berrichi, A.; Amodeo, L.; Yalaoui, F.; Châtalet, E.; Mezghiche, M. Bi-objective optimization algorithms for joint production and maintenance scheduling: Application to the parallel machine problem. J. Intell. Manuf. 2009, 20, 389–400. [Google Scholar] [CrossRef]
- Moradi, E.; Fatemi Ghomi, S.M.T.; Zandieh, M. Bi-objective optimization research on integrated fixed time interval preventive maintenance and production for scheduling flexible job-shop problem. Expert Syst. Appl. 2011, 38, 7169–7178. [Google Scholar] [CrossRef]
- Hnaien, F.; Yalaoui, F. A bi-criteria flow-shop scheduling with preventive maintenance. IFAC Proc. Vol. 2013, 46, 1387–1392. [Google Scholar] [CrossRef]
- Wang, S.; Liu, M. Multi-objective optimization of parallel machine scheduling integrated with multi-resources preventive maintenance planning. J. Manuf. Syst. 2015, 37, 182–192. [Google Scholar] [CrossRef]
- Piasson, D.; Bíscaro, A.A.P.; Leão, F.B.; Sanches Mantovani, J.R. A new approach for reliability-centered maintenance programs in electric power distribution systems based on a multiobjective genetic algorithm. Electr. Power Syst. Res. 2016, 137, 41–50. [Google Scholar] [CrossRef] [Green Version]
- Sheikhalishahi, M.; Eskandari, N.; Mashayekhi, A.; Azadeh, A. Multi-objective open shop scheduling by considering human error and preventive maintenance. Appl. Math. Model. 2019, 67, 573–587. [Google Scholar] [CrossRef]
- Boufellouh, R.; Belkaid, F. Bi-objective optimization algorithms for joint production and maintenance scheduling under a global resource constraint: Application to the permutation flow shop problem. Comput. Oper. Res. 2020, 122, 104943. [Google Scholar] [CrossRef]
- Zhang, C.; Yang, T. Optimal maintenance planning and resource allocation for wind farms based on non-dominated sorting genetic algorithm-II. Renew. Energy 2021, 164, 1540–1549. [Google Scholar] [CrossRef]
- Berrichi, A.; Yalaoui, F.; Amodeo, L.; Mezghiche, M. Bi-Objective Ant Colony Optimization approach to optimize production and maintenance scheduling. Comput. Oper. Res. 2010, 37, 1584–1596. [Google Scholar] [CrossRef]
- Suresh, K.; Kumarappan, N. Hybrid improved binary particle swarm optimization approach for generation maintenance scheduling problem. Swarm Evol. Comput. 2013, 9, 69–89. [Google Scholar] [CrossRef]
- Li, J.Q.; Pan, Q.K.; Tasgetiren, M.F. A discrete artificial bee colony algorithm for the multi-objective flexible job-shop scheduling problem with maintenance activities. Swarm Evol. Comput. 2014, 38, 1111–1132. [Google Scholar] [CrossRef]
- Galván, B.; Winter, G.; Greiner, D.; Salazar, D.; Méndez, M. New Evolutionary Methodologies for Integrated Safety System Design and Maintenance Optimization. In Computational Intelligence in Reliability Engineering: Evolutionary Techniques in Reliability Analysis and Optimization; Levitin, G., Ed.; Springer: Berlin/Heidelberg, Germany, 2007; pp. 151–190. [Google Scholar]
- Okasha, N.M.; Frangopol, D.M. Lifetime-oriented multi-objective optimization of structural maintenance considering system reliability, redundancy and life-cycle cost using GA. Struct. Saf. 2009, 31, 460–474. [Google Scholar] [CrossRef]
- Adjoul, O.; Benfriha, K.; El Zant, C.; Aoussat, A. Algorithmic Strategy for Simultaneous Optimization of Design and Maintenance of Multi-Component Industrial Systems. Reliab. Eng. Syst. Saf. 2021, 208. [Google Scholar] [CrossRef]
- Lins, I.D.; Droguett, E.A.L. Multiobjective optimization of availability and cost in repairable systems design via genetic algorithms and discrete event simulation. Pesquisa Oper. 2009, 29, 43–66. [Google Scholar] [CrossRef] [Green Version]
- Lins, I.D.; López, E. Redundancy allocation problems considering systems with imperfect repairs using multi-objective genetic algorithms and discrete event simulation. Simul. Model. Pract. Theory 2011, 19, 362–381. [Google Scholar] [CrossRef]
- SINTEF. OREDA—Offshore Reliability Data Handbook, 5th ed.; Det Norske Veritas: Hovik, Norway, 2009; p. 796. ISBN 978-82-14-04830-8. [Google Scholar]
- Emmerich, M.; Deutz, A. A tutorial on multiobjective optimization: Fundamentals and evolutionary methods. Nat. Comput. 2018, 17, 585–609. [Google Scholar] [CrossRef] [Green Version]
- Coello, C.A. Multi-objective Evolutionary Algorithms in Real-World Applications: Some Recent Results and Current Challenges. In Advances in Evolutionary and Deterministic Methods for Design, Optimization and Control in Engineering and Sciences; Computational Methods in Applied Sciences; Greiner, D., Galván, B., Périaux, J., Gauger, N., Giannakoglou, K., Winter, G., Eds.; Springer International Publishing AG: Cham, Switzerland, 2015; pp. 3–18. [Google Scholar]
- Whitley, D.; Rana, S.; Heckendorn, R. Representation Issues in Neighborhood Search and Evolutionary Algorithms. In Genetic Algorithms and Evolution Strategies in Engineering and Computer Science; Quagliarella, D., Périaux, J., Poloni, C., Winter, G., Eds.; John Wiley & Sons: Chichester, UK, 1997; pp. 39–57. [Google Scholar]
- Savage, C. A Survey of Combinatorial Gray Codes. SIAM Rev. 1997, 39, 605–629. [Google Scholar] [CrossRef] [Green Version]
- Greiner, D.; Winter, G.; Emperador, J.M.; Galván, B. Gray Coding in Evolutionary Multicriteria Optimization: Application in Frame Structural Optimum Design. In Evolutionary Multi-Criterion Optimization. EMO 2005. Lecture Notes in Computer Science; Coello, C.A., Hernández Aguirre, A., Zitzler, E., Eds.; Springer: Berlin/Heidelberg, Germany, 2005; Volume 3410. [Google Scholar] [CrossRef]
- Tian, Y.; Cheng, R.; Zhang, X.; Jin, Y. PlatEMO: A MATLAB Platform for Evolutionary Multi-Objective Optimization [Educational Forum]. IEEE Comput. Intell. Mag. 2017, 12, 73–87. [Google Scholar] [CrossRef] [Green Version]
- Zitzler, E.; Thiele, L.; Laumanns, M.; Fonseca, C.M.; da Fonseca, V.G. Performance Assessment of Multiobjective Optimizers: An Analysis and Review. IEEE Trans. Evol. Comput. 2003, 7, 117–132. [Google Scholar] [CrossRef] [Green Version]
- Fonseca, C.M.; Paquete, L.; López-Ibáñez, M. An Improved Dimension-Sweep Algorithm for the Hypervolume Indicator. In Proceedings of the IEEE International Conference on Evolutionary Computation, Vancouver, BC, Canada, 16–21 July 2006; pp. 1157–1163. [Google Scholar] [CrossRef] [Green Version]
- Benavoli, A.; Corani, G.; Mangili, F. Should We Really Use Post-Hoc Tests Based on Mean-Ranks? J. Mach. Learn. Res. 2016, 17, 1–10. [Google Scholar]
Parameter | Value | Source |
---|---|---|
Life Cycle | 700,800 h | - |
Corrective Maintenance Cost | 0.5 units | Machinery & Reliability Institute |
Preventive Maintenance Cost | 0.125 units | Machinery & Reliability Institute |
Pump | 1 h | Machinery & Reliability Institute |
Pump | 70,080 h | Machinery & Reliability Institute |
Pump | 159.57 × 10 h | OREDA 2009 |
Pump | 1 h | Machinery & Reliability Institute |
Pump | 24.33 h | + 4 |
Pump | 11 h | OREDA 2009 |
Pump | 3.33 h | ( – TRmin)/3 |
Pump | 2920 h | Machinery & Reliability Institute |
Pump | 8760 h | Machinery & Reliability Institute |
Pump | 4 h | Machinery & Reliability Institute |
Pump | 8 h | Machinery & Reliability Institute |
Valve | 1 h | Machinery & Reliability Institute |
Valve | 70,080 h | Machinery & Reliability Institute |
Valve | 44.61 × 10 h | OREDA 2009 |
Valve | 1 h | Machinery & Reliability Institute |
Valve | 20.83 h | + 4 |
Valve | 9.5 h | OREDA 2009 |
Valve | 2.83 h | ( – TRmin)/3 |
Valve | 8760 h | Machinery & Reliability Institute |
Valve | 35,040 h | Machinery & Reliability Institute |
Valve | 1 h | Machinery & Reliability Institute |
Valve | 3 h | Machinery & Reliability Institute |
Method | Encoding | Crossover | Time Unit | Population | PrM | disM | PrC | disC |
---|---|---|---|---|---|---|---|---|
NSGA-II | Real | SBX | Hour | 50-100-150 | 0.5-1-1.5 | 20 | 1 | 20 |
Standard Binary | 1 point (1PX) | Hour | 50-100-150 | 0.5-1-1.5 | - | 1 | - | |
Standard Binary | 2 Point (2PX) | Hour | 50-100-150 | 0.5-1-1.5 | - | 1 | - | |
Standard Binary | Uniform (UX) | Hour | 50-100-150 | 0.5-1-1.5 | - | 1 | - | |
Gray | 1 Point (1PX) | Hour | 50-100-150 | 0.5-1-1.5 | - | 1 | - | |
Gray | 2 Point (2PX) | Hour | 50-100-150 | 0.5-1-1.5 | - | 1 | - | |
Gray | Uniform (UX) | Hour | 50-100-150 | 0.5-1-1.5 | - | 1 | - | |
Standard Binary | 1 Point (1PX) | Day | 50-100-150 | 0.5-1-1.5 | - | 1 | - | |
Standard Binary | 2 Point (2PX) | Day | 50-100-150 | 0.5-1-1.5 | - | 1 | - | |
Standard Binary | Uniform (UX) | Day | 50-100-150 | 0.5-1-1.5 | - | 1 | - | |
Standard Binary | 1 Point (1PX) | Week | 50-100-150 | 0.5-1-1.5 | - | 1 | - | |
Standard Binary | 2 Point (2PX) | Week | 50-100-150 | 0.5-1-1.5 | - | 1 | - | |
Standard Binary | Uniform (UX) | Week | 50-100-150 | 0.5-1-1.5 | - | 1 | - |
Encoding | Time Unit | Average Time | Sequential Time |
---|---|---|---|
Real SBX | Hour | 2 days, 18 h and 12 min. | 8 months, 4 days, 23 h and 22 min. |
Binary 1PX | Hour | 2 days, 18 h and 57 min. | 8 months, 7 days, 18 h and 48 min. |
Binary 2PX | Hour | 3 days, 2 h and 39 min. | 9 months, 6 days, 5 h and 52 min. |
Binary UX | Hour | 3 days and 39 min. | 8 months, 29 days, 3 h and 8 min. |
Gray 1PX | Hour | 2 days, 19 h and 35 min. | 8 months, 10 days, 3 h and 23 min. |
Gray 2PX | Hour | 2 days, 19 h and 56 min. | 8 months, 11 days, 11 h and 29 min. |
Gray UX | Hour | 2 days, 19 h and 40 min. | 8 months, 10 days, 10 h and 34 min. |
Total | 2 days, 21 h and 6 min. | 4 years, 11 months, 19 days, 8 h and 36 min. |
Identifier | Configuration | Average | Median | Max. | Min. | St. D. | Friedman’s Test |
---|---|---|---|---|---|---|---|
Av. Rank | |||||||
ID1 | N = 50; PrM = 0.5 | 2.2821 | 2.2857 | 2.3079 | 2.2531 | 0.0179 | 5.900 |
ID2 | N = 100; PrM = 0.5 | 2.2886 | 2.2831 | 2.3258 | 2.2682 | 0.0170 | 4.800 |
ID3 | N = 150; PrM = 0.5 | 2.2938 | 2.2883 | 2.3187 | 2.2730 | 0.0144 | 4.300 |
ID4 | N = 50; PrM = 1 | 2.2927 | 2.2932 | 2.3113 | 2.2710 | 0.0117 | 3.999 |
ID5 | N = 100; PrM = 1 | 2.2863 | 2.2916 | 2.3064 | 2.2584 | 0.0148 | 5.200 |
ID6 | N = 150; PrM = 1 | 2.2890 | 2.2885 | 2.3157 | 2.2619 | 0.0170 | 5.100 |
ID7 | N = 50; PrM = 1.5 | 2.2857 | 2.2788 | 2.3228 | 2.2701 | 0.0171 | 6.100 |
ID8 | N = 100; PrM = 1.5 | 2.2821 | 2.2826 | 2.3051 | 2.2561 | 0.0162 | 5.600 |
ID9 | N = 150; PrM = 1.5 | 2.2968 | 2.2962 | 2.3307 | 2.2641 | 0.0211 | 4.000 |
p-value | 0.5788 |
Encoding | Identifier | Configuration | Average | Median | Max. | Min. | St. D. | Friedman’s Test |
---|---|---|---|---|---|---|---|---|
Av. Rank | ||||||||
B1PX | ID1 | N = 50; PrM = 0.5 | 2.2936 | 2.2919 | 2.3297 | 2.2738 | 0.0169 | 4.600 |
ID2 | N = 100; PrM = 0.5 | 2.2919 | 2.2953 | 2.3057 | 2.2703 | 0.0107 | 4.500 | |
ID3 | N = 150; PrM = 0.5 | 2.2863 | 2.2915 | 2.2984 | 2.2554 | 0.0132 | 5.400 | |
ID4 | N = 50; PrM = 1 | 2.2929 | 2.2918 | 2.3151 | 2.2566 | 0.0171 | 4.400 | |
ID5 | N = 100; PrM = 1 | 2.2865 | 2.2822 | 2.3200 | 2.2674 | 0.0151 | 5.799 | |
ID6 | N = 150; PrM = 1 | 2.2932 | 2.2926 | 2.3216 | 2.2761 | 0.0150 | 4.700 | |
ID7 | N = 50; PrM = 1.5 | 2.2803 | 2.2775 | 2.2994 | 2.2642 | 0.0136 | 7.000 | |
ID8 | N = 100; PrM = 1.5 | 2.2970 | 2.3008 | 2.3126 | 2.2708 | 0.0134 | 3.100 | |
ID9 | N = 150; PrM = 1.5 | 2.2884 | 2.2871 | 2.3394 | 2.2616 | 0.0209 | 5.500 | |
p-value | 0.1228 | |||||||
B2PX | ID1 | N = 50; PrM = 0.5 | 2.2808 | 2.2800 | 2.3173 | 2.2401 | 0.0237 | 6.500 |
ID2 | N = 100; PrM = 0.5 | 2.3013 | 2.3051 | 2.3260 | 2.2714 | 0.0172 | 3.300 | |
ID3 | N = 150; PrM = 0.5 | 2.2916 | 2.2875 | 2.3627 | 2.2630 | 0.0298 | 5.400 | |
ID4 | N = 50; PrM = 1 | 2.2941 | 2.2926 | 2.3156 | 2.2724 | 0.0152 | 4.600 | |
ID5 | N = 100; PrM = 1 | 2.2908 | 2.2929 | 2.3056 | 2.2755 | 0.0089 | 4.800 | |
ID6 | N = 150; PrM = 1 | 2.2921 | 2.2905 | 2.3215 | 2.2644 | 0.0160 | 4.999 | |
ID7 | N = 50; PrM = 1.5 | 2.2924 | 2.2873 | 2.3483 | 2.2694 | 0.0226 | 5.100 | |
ID8 | N = 100; PrM = 1.5 | 2.2931 | 2.2930 | 2.3326 | 2.2606 | 0.0205 | 4.900 | |
ID9 | N = 150; PrM = 1.5 | 2.2883 | 2.2849 | 2.3107 | 2.2738 | 0.0145 | 5.399 | |
p-value | 0.4762 | |||||||
BUX | ID1 | N = 50; PrM = 0.5 | 2.2841 | 2.2828 | 2.3023 | 2.2611 | 0.0122 | 5.500 |
ID2 | N = 100; PrM = 0.5 | 2.2954 | 2.3000 | 2.3141 | 2.2647 | 0.0163 | 3.300 | |
ID3 | N = 150; PrM = 0.5 | 2.2883 | 2.2880 | 2.3265 | 2.2576 | 0.0199 | 5.200 | |
ID4 | N = 50; PrM = 1 | 2.2959 | 2.2939 | 2.3497 | 2.2702 | 0.0223 | 3.400 | |
ID5 | N = 100; PrM = 1 | 2.2848 | 2.2866 | 2.3046 | 2.2660 | 0.0140 | 5.700 | |
ID6 | N = 150; PrM = 1 | 2.2830 | 2.2850 | 2.2955 | 2.2687 | 0.0085 | 5.799 | |
ID7 | N = 50; PrM = 1.5 | 2.2893 | 2.2866 | 2.3112 | 2.2622 | 0.0163 | 5.100 | |
ID8 | N = 100; PrM = 1.5 | 2.2808 | 2.2800 | 2.2994 | 2.2579 | 0.0138 | 6.100 | |
ID9 | N = 150; PrM = 1.5 | 2.2857 | 2.2855 | 2.3036 | 2.2657 | 0.0121 | 4.899 | |
p-value | 0.2132 |
Encoding | Identifier | Configuration | Average | Median | Max. | Min. | St. D. | Friedman’s Test |
---|---|---|---|---|---|---|---|---|
Av. Rank | ||||||||
G1PX | ID1 | N = 50; PrM = 0.5 | 2.2833 | 2.2838 | 2.2929 | 2.2710 | 0.0064 | 5.500 |
ID2 | N = 100; PrM = 0.5 | 2.2990 | 2.3010 | 2.3556 | 2.2640 | 0.0252 | 3.600 | |
ID3 | N = 150; PrM = 0.5 | 2.2815 | 2.2850 | 2.2951 | 2.2626 | 0.0099 | 6.200 | |
ID4 | N = 50; PrM = 1 | 2.2865 | 2.2834 | 2.3017 | 2.2762 | 0.0083 | 5.200 | |
ID5 | N = 100; PrM = 1 | 2.2989 | 2.2986 | 2.3347 | 2.2652 | 0.0186 | 3.100 | |
ID6 | N = 150; PrM = 1 | 2.2812 | 2.2791 | 2.3043 | 2.2611 | 0.0117 | 6.000 | |
ID7 | N = 50; PrM = 1.5 | 2.2882 | 2.2916 | 2.3171 | 2.2526 | 0.0187 | 4.200 | |
ID8 | N = 100; PrM = 1.5 | 2.2786 | 2.2804 | 2.2992 | 2.2608 | 0.0137 | 6.100 | |
ID9 | N = 150; PrM = 1.5 | 2.2874 | 2.2820 | 2.3180 | 2.2706 | 0.0155 | 5.100 | |
p-value | 0.0943 | |||||||
G2PX | ID1 | N = 50; PrM = 0.5 | 2.2947 | 2.2971 | 2.3192 | 2.2757 | 0.0128 | 3.700 |
ID2 | N = 100; PrM = 0.5 | 2.2802 | 2.2814 | 2.2953 | 2.2592 | 0.0121 | 6.500 | |
ID3 | N = 150; PrM = 0.5 | 2.2856 | 2.2895 | 2.2978 | 2.2519 | 0.0136 | 4.299 | |
ID4 | N = 50; PrM = 1 | 2.2912 | 2.2868 | 2.3186 | 2.2659 | 0.0192 | 4.900 | |
ID5 | N = 100; PrM = 1 | 2.2832 | 2.2835 | 2.2951 | 2.2690 | 0.0070 | 5.600 | |
ID6 | N = 150; PrM = 1 | 2.2899 | 2.2913 | 2.3256 | 2.2640 | 0.0193 | 4.100 | |
ID7 | N = 50; PrM = 1.5 | 2.2866 | 2.2880 | 2.3132 | 2.2617 | 0.0169 | 4.800 | |
ID8 | N = 100; PrM = 1.5 | 2.2809 | 2.2813 | 2.3140 | 2.2585 | 0.0148 | 6.500 | |
ID9 | N = 150; PrM = 1.5 | 2.2920 | 2.2912 | 2.3364 | 2.2739 | 0.0169 | 4.600 | |
p-value | 0.2164 | |||||||
GUX | ID1 | N = 50; PrM = 0.5 | 2.2862 | 2.2892 | 2.3088 | 2.2540 | 0.0180 | 4.300 |
ID2 | N = 100; PrM = 0.5 | 2.2828 | 2.2813 | 2.2967 | 2.2696 | 0.0079 | 5.100 | |
ID3 | N = 150; PrM = 0.5 | 2.2885 | 2.2922 | 2.3077 | 2.2666 | 0.0136 | 4.100 | |
ID4 | N = 50; PrM = 1 | 2.2907 | 2.2939 | 2.3134 | 2.2683 | 0.0156 | 4.000 | |
ID5 | N = 100; PrM = 1 | 2.2879 | 2.2858 | 2.3165 | 2.2611 | 0.0187 | 4.699 | |
ID6 | N = 150; PrM = 1 | 2.2835 | 2.2816 | 2.3048 | 2.2701 | 0.0120 | 5.400 | |
ID7 | N = 50; PrM = 1.5 | 2.2862 | 2.2852 | 2.3056 | 2.2608 | 0.0143 | 5.600 | |
ID8 | N = 100; PrM = 1.5 | 2.2852 | 2.2851 | 2.3068 | 2.2524 | 0.0153 | 5.100 | |
ID9 | N = 150; PrM = 1.5 | 2.2755 | 2.2700 | 2.3064 | 2.2638 | 0.0142 | 6.700 | |
p-value | 0.4572 |
Identifier | Description | Configuration | Average | Median | Max. | Min. | St. D. | Friedman’s Test |
---|---|---|---|---|---|---|---|---|
Av. Rank | ||||||||
ID1 | Real | N = 50; PrM = 1 | 2.2927 | 2.2932 | 2.3113 | 2.2710 | 0.0117 | 4.600 |
ID2 | Binary 1P | N = 100; PrM = 1.5 | 2.2970 | 2.3008 | 2.3126 | 2.2708 | 0.0134 | 4.000 |
ID3 | Binary 2P | N = 100; PrM = 0.5 | 2.3013 | 2.3051 | 2.3260 | 2.2714 | 0.0172 | 3.000 |
ID4 | Binary U | N = 100; PrM = 0.5 | 2.2954 | 2.3000 | 2.3141 | 2.2647 | 0.0163 | 4.199 |
ID5 | Gray 1P | N = 100; PrM = 1 | 2.2989 | 2.2986 | 2.3347 | 2.2652 | 0.0186 | 3.500 |
ID6 | Gray 2P | N = 50; PrM = 0.5 | 2.2947 | 2.2971 | 2.3192 | 2.2757 | 0.0128 | 4.000 |
ID7 | Gray U | N = 50; PrM = 1 | 2.2907 | 2.2939 | 2.3134 | 2.2683 | 0.0156 | 4.699 |
p-value | 0.5979 |
Encoding | Hypervolume Accumulated Value |
---|---|
Real | 2.3943 |
Binary 1 Point Crossover | 2.4142 |
Binary 2 Point Crossover | 2.4298 |
Binary Uniform Crossover | 2.3984 |
Gray 1 Point Crossover | 2.4011 |
Gray 2 Point Crossover | 2.3982 |
Gray Uniform Crossover | 2.3829 |
Global | 2.4553 |
Encoding | Time Unit | Average Time | Sequential Time |
---|---|---|---|
Binary 1PX | Hour | 2 days, 18 h and 57 min | 8 months, 7 days, 18 h and 48 min. |
Binary 2PX | Hour | 3 days, 2 h and 39 min | 9 months, 6 days, 5 h and 52 min. |
Binary UX | Hour | 3 days and 39 min | 8 months, 29 days, 3 h and 8 min. |
Binary 1PX | Day | 2 days, 22 h and 3 min | 8 months, 19 days, 9 h and 14 min. |
Binary 2PX | Day | 2 days, 19 h and 13 min | 8 months, 8 days, 18 h and 18 min. |
Binary UX | Day | 2 days, 19 h and 13 min | 8 months, 8 days, 18 h and 55 min. |
Binary 1PX | Week | 2 days, 18 h and 3 min | 8 months, 4 days, 9 h and 39 min. |
Binary 2PX | Week | 2 days, 19 h and 21 min | 8 months, 9 days, 6 h and 55 min. |
Binary UX | Week | 2 days, 19 h and 20 min | 8 months, 9 days, 4 h and 40 min. |
Total | 2 days, 20 h and 53 min | 6 years, 4 months, 12 days, 22 h and 41 min. |
Encoding | Identifier | Configuration | Average | Median | Max. | Min. | St. D. | Friedman’s Test |
---|---|---|---|---|---|---|---|---|
Av. Rank | ||||||||
B1PX-D | ID1 | N = 50; PrM = 0.5 | 2.2834 | 2.2781 | 2.3430 | 2.2655 | 0.0220 | 5.800 |
ID2 | N = 100; PrM = 0.5 | 2.2942 | 2.2911 | 2.3198 | 2.2788 | 0.0133 | 4.100 | |
ID3 | N = 150; PrM = 0.5 | 2.2785 | 2.2834 | 2.2934 | 2.2563 | 0.0126 | 6.600 | |
ID4 | N = 50; PrM = 1 | 2.2819 | 2.2795 | 2.3127 | 2.2622 | 0.0152 | 5.800 | |
ID5 | N = 100; PrM = 1 | 2.2896 | 2.2870 | 2.3188 | 2.2621 | 0.0184 | 4.500 | |
ID6 | N = 150; PrM = 1 | 2.2964 | 2.2941 | 2.3193 | 2.2804 | 0.0103 | 3.200 | |
ID7 | N = 50; PrM = 1.5 | 2.2922 | 2.2932 | 2.3239 | 2.2707 | 0.0177 | 4.400 | |
ID8 | N = 100; PrM = 1.5 | 2.2776 | 2.2793 | 2.3057 | 2.2530 | 0.0171 | 6.800 | |
ID9 | N = 150; PrM = 1.5 | 2.2939 | 2.2947 | 2.3107 | 2.2758 | 0.0135 | 3.800 | |
p-value | 0.0246 | |||||||
B2PX-D | ID1 | N = 50; PrM = 0.5 | 2.2838 | 2.2795 | 2.3043 | 2.2681 | 0.0123 | 5.400 |
ID2 | N = 100; PrM = 0.5 | 2.2783 | 2.2785 | 2.2947 | 2.2609 | 0.0106 | 6.600 | |
ID3 | N = 150; PrM = 0.5 | 2.2867 | 2.2880 | 2.3043 | 2.2757 | 0.0084 | 4.300 | |
ID4 | N = 50; PrM = 1 | 2.2930 | 2.2865 | 2.3173 | 2.2746 | 0.0176 | 4.200 | |
ID5 | N = 100; PrM = 1 | 2.2874 | 2.2831 | 2.3188 | 2.2652 | 0.0192 | 5.399 | |
ID6 | N = 150; PrM = 1 | 2.2916 | 2.2937 | 2.3150 | 2.2712 | 0.0152 | 4.100 | |
ID7 | N = 50; PrM = 1.5 | 2.2877 | 2.2849 | 2.3249 | 2.2693 | 0.0161 | 4.600 | |
ID8 | N = 100; PrM = 1.5 | 2.2841 | 2.2817 | 2.3136 | 2.2658 | 0.0133 | 5.300 | |
ID9 | N = 150; PrM = 1.5 | 2.2887 | 2.2826 | 2.3372 | 2.2651 | 0.0216 | 5.100 | |
p-value | 0.5612 | |||||||
BUX-D | ID1 | N = 50; PrM = 0.5 | 2.2918 | 2.2948 | 2.3152 | 2.2643 | 0.0168 | 4.200 |
ID2 | N = 100; PrM = 0.5 | 2.2897 | 2.2865 | 2.3356 | 2.2602 | 0.0198 | 4.699 | |
ID3 | N = 150; PrM = 0.5 | 2.2844 | 2.2866 | 2.3104 | 2.2522 | 0.0179 | 5.600 | |
ID4 | N = 50; PrM = 1 | 2.2798 | 2.2732 | 2.3078 | 2.2601 | 0.0160 | 6.400 | |
ID5 | N = 100; PrM = 1 | 2.2897 | 2.2893 | 2.3127 | 2.2729 | 0.0151 | 4.800 | |
ID6 | N = 150; PrM = 1 | 2.2923 | 2.2907 | 2.3333 | 2.2666 | 0.0219 | 4.700 | |
ID7 | N = 50; PrM = 1.5 | 2.2894 | 2.2880 | 2.3201 | 2.2691 | 0.0138 | 5.100 | |
ID8 | N = 100; PrM = 1.5 | 2.2923 | 2.2901 | 2.3279 | 2.2637 | 0.0183 | 5.000 | |
ID9 | N = 150; PrM = 1.5 | 2.2908 | 2.2863 | 2.3461 | 2.2606 | 0.0249 | 4.500 | |
p-value | 0.8007 |
Comparison | p-Value | Conclusion |
---|---|---|
ID3 versus ID6 | 0.0059 | Significant difference found |
ID4 versus ID6 | 0.0371 | Significant difference found |
ID6 versus ID8 | 0.0371 | Significant difference found |
ID1 versus ID6 | 0.0840 | The null hypothesis cannot be rejected |
ID5 versus ID6 | 0.3223 | The null hypothesis cannot be rejected |
ID2 versus ID6 | 0.3750 | The null hypothesis cannot be rejected |
ID6 versus ID7 | 0.4922 | The null hypothesis cannot be rejected |
ID6 versus ID9 | 0.6250 | The null hypothesis cannot be rejected |
Encoding | Identifier | Configuration | Average | Median | Max. | Min. | St. D. | Friedman’s Test |
---|---|---|---|---|---|---|---|---|
Av. Rank | ||||||||
B1PX-W | ID1 | N = 50; PrM = 0.5 | 2.2957 | 2.2935 | 2.3430 | 2.2561 | 0.0250 | 4.300 |
ID2 | N = 100; PrM = 0.5 | 2.2872 | 2.2873 | 2.3127 | 2.2665 | 0.0142 | 5.100 | |
ID3 | N = 150; PrM = 0.5 | 2.2897 | 2.2808 | 2.3252 | 2.2623 | 0.0212 | 5.100 | |
ID4 | N = 50; PrM = 1 | 2.2812 | 2.2805 | 2.2977 | 2.2520 | 0.0157 | 5.799 | |
ID5 | N = 100; PrM = 1 | 2.2907 | 2.2893 | 2.3007 | 2.2812 | 0.0060 | 3.999 | |
ID6 | N = 150; PrM = 1 | 2.2825 | 2.2825 | 2.3004 | 2.2679 | 0.0101 | 6.100 | |
ID7 | N = 50; PrM = 1.5 | 2.2890 | 2.2899 | 2.3110 | 2.2597 | 0.0164 | 4.900 | |
ID8 | N = 100; PrM = 1.5 | 2.2874 | 2.2848 | 2.3066 | 2.2690 | 0.0123 | 5.000 | |
ID9 | N = 150; PrM = 1.5 | 2.2888 | 2.2892 | 2.3176 | 2.2634 | 0.0155 | 4.700 | |
p-value | 0.7979 | |||||||
B2PX-W | ID1 | N = 50; PrM = 0.5 | 2.2917 | 2.2911 | 2.3101 | 2.2764 | 0.0121 | 4.000 |
ID2 | N = 100; PrM = 0.5 | 2.2858 | 2.2858 | 2.3013 | 2.2671 | 0.0122 | 5.300 | |
ID3 | N = 150; PrM = 0.5 | 2.2892 | 2.2869 | 2.3237 | 2.2567 | 0.0250 | 4.900 | |
ID4 | N = 50; PrM = 1 | 2.2797 | 2.2743 | 2.3014 | 2.2703 | 0.0108 | 6.299 | |
ID5 | N = 100; PrM = 1 | 2.2840 | 2.2838 | 2.2965 | 2.2665 | 0.0085 | 5.200 | |
ID6 | N = 150; PrM = 1 | 2.2793 | 2.2780 | 2.3066 | 2.2502 | 0.0147 | 6.200 | |
ID7 | N = 50; PrM = 1.5 | 2.2921 | 2.2817 | 2.3465 | 2.2709 | 0.0246 | 4.899 | |
ID8 | N = 100; PrM = 1.5 | 2.2895 | 2.2880 | 2.3198 | 2.2546 | 0.0168 | 4.001 | |
ID9 | N = 150; PrM = 1.5 | 2.2918 | 2.2854 | 2.3130 | 2.2710 | 0.0145 | 4.200 | |
p-value | 0.4439 | |||||||
BUX-W | ID1 | N = 50; PrM = 0.5 | 2.2894 | 2.2847 | 2.3336 | 2.2659 | 0.0197 | 5.500 |
ID2 | N = 100; PrM = 0.5 | 2.2925 | 2.2922 | 2.3144 | 2.2756 | 0.0113 | 4.399 | |
ID3 | N = 150; PrM = 0.5 | 2.2911 | 2.2932 | 2.3216 | 2.2664 | 0.0156 | 4.399 | |
ID4 | N = 50; PrM = 1 | 2.2858 | 2.2836 | 2.3057 | 2.2660 | 0.0128 | 5.499 | |
ID5 | N = 100; PrM = 1 | 2.2826 | 2.2834 | 2.3144 | 2.2631 | 0.0140 | 6.700 | |
ID6 | N = 150; PrM = 1 | 2.2854 | 2.2857 | 2.3122 | 2.2645 | 0.0164 | 5.499 | |
ID7 | N = 50; PrM = 1.5 | 2.2893 | 2.2887 | 2.3094 | 2.2681 | 0.0125 | 4.899 | |
ID8 | N = 100; PrM = 1.5 | 2.3009 | 2.2942 | 2.3316 | 2.2870 | 0.0145 | 3.400 | |
ID9 | N = 150; PrM = 1.5 | 2.2921 | 2.2862 | 2.3163 | 2.2774 | 0.0143 | 4.700 | |
p-value | 0.3128 |
Id. | Description | Configuration | Average | Median | Max. | Min. | St. D. | Friedman’s Test |
---|---|---|---|---|---|---|---|---|
Av. Rank | ||||||||
ID1 | Binary 1P (hour) | N = 100; PrM = 1.5 | 2.2970 | 2.3008 | 2.3126 | 2.2708 | 0.0134 | 4.700 |
ID2 | Binary 2P (hour) | N = 100; PrM = 0.5 | 2.3013 | 2.3051 | 2.3260 | 2.2714 | 0.0172 | 3.699 |
ID3 | Binary U (hour) | N = 100; PrM = 0.5 | 2.2954 | 2.3000 | 2.3141 | 2.2647 | 0.0163 | 4.500 |
ID4 | Binary 1P (day) | N = 150; PrM = 1.0 | 2.2964 | 2.2941 | 2.3193 | 2.2804 | 0.0103 | 5.000 |
ID5 | Binary 2P (day) | N = 150; PrM = 1.0 | 2.2916 | 2.2937 | 2.3150 | 2.2712 | 0.0152 | 5.800 |
ID6 | Binary U (day) | N = 50; PrM = 0.5 | 2.2918 | 2.2948 | 2.3152 | 2.2643 | 0.0168 | 5.300 |
ID7 | Binary 1P (week) | N = 100; PrM = 1.0 | 2.2907 | 2.2893 | 2.3007 | 2.2812 | 0.0060 | 6.200 |
ID8 | Binary 2P (week) | N = 50; PrM = 0.5 | 2.2917 | 2.2911 | 2.3101 | 2.2764 | 0.0121 | 5.900 |
ID9 | Binary U (week) | N = 100; PrM = 1.5 | 2.3009 | 2.2942 | 2.3316 | 2.2870 | 0.0145 | 3.900 |
p-value | 0.4053 |
Encoding | Time Unit | Hypervolume Accumulated Value |
---|---|---|
Binary 1 Point Crossover | Hour | 2.4142 |
Binary 2 Point Crossover | Hour | 2.4298 |
Binary Uniform Crossover | Hour | 2.3984 |
Binary 1 Point Crossover | Day | 2.4007 |
Binary 2 Point Crossover | Day | 2.3942 |
Binary Uniform Crossover | Day | 2.3984 |
Binary 1 Point Crossover | Week | 2.4047 |
Binary 2 Point Crossover | Week | 2.4285 |
Binary Uniform Crossover | Week | 2.4037 |
Global | 2.4646 |
Id | Q | Cost [eu] | Time Unit | V1 [h] | P2 [h] | P3 [h] | V4 [h] | V5 [h] | V6 [h] | V7 [h] |
---|---|---|---|---|---|---|---|---|---|---|
1 | 0.002720 | 823.38 | Hours | 25,408 | 0 | 8633 | 0 | 34,179 | 34,903 | 31,386 |
2 | 0.002713 | 835.75 | Hours | 29,225 | 0 | 8633 | 0 | 27,070 | 33,454 | 33,690 |
3 | 0.002591 | 986.88 | Days | 1285 | 0 | 353 | 902 | 1441 | 967 | 1386 |
4 | 0.002576 | 988.00 | Weeks | 127 | 0 | 47 | 165 | 206 | 172 | 203 |
5 | 0.002295 | 992.75 | Days | 1435 | 0 | 350 | 830 | 1088 | 1459 | 1454 |
6 | 0.001495 | 1360.50 | Hours | 34,746 | 8239 | 8408 | 0 | 32,676 | 31,769 | 31,484 |
7 | 0.001334 | 1363.75 | Days | 1394 | 360 | 315 | 0 | 1125 | 1301 | 1026 |
8 | 0.001276 | 1424.25 | Weeks | 204 | 50 | 37 | 0 | 174 | 173 | 188 |
9 | 0.001189 | 1431.75 | Hours | 31,040 | 8617 | 8103 | 0 | 34,787 | 31,445 | 29,929 |
10 | 0.001174 | 1449.00 | Days | 178 | 51 | 50 | 0 | 182 | 171 | 195 |
11 | 0.001112 | 1496.12 | Weeks | 112 | 48 | 48 | 0 | 176 | 196 | 144 |
12 | 0.001056 | 1512.50 | Days | 880 | 358 | 354 | 1186 | 1020 | 1346 | 1098 |
13 | 0.001039 | 1524.00 | Weeks | 162 | 45 | 42 | 195 | 147 | 176 | 189 |
14 | 0.001002 | 1524.62 | Hours | 28,042 | 7228 | 6549 | 30,777 | 20,951 | 31,779 | 34,982 |
15 | 0.000939 | 1528.88 | Hours | 29,939 | 7028 | 7904 | 21,690 | 25,711 | 34,814 | 34,791 |
16 | 0.000913 | 1530.38 | Weeks | 194 | 47 | 49 | 103 | 137 | 145 | 179 |
17 | 0.000835 | 1567.75 | Days | 1028 | 363 | 340 | 1280 | 1430 | 986 | 1326 |
18 | 0.000813 | 1611.38 | Hours | 22,516 | 5955 | 6298 | 17,237 | 29,568 | 27,075 | 22,757 |
19 | 0.000776 | 1667.25 | Hours | 32,355 | 7185 | 8384 | 27,489 | 26,949 | 32,700 | 32,770 |
20 | 0.000771 | 1741.00 | Hours | 27,443 | 7556 | 8520 | 30,003 | 30,056 | 30,579 | 21,639 |
21 | 0.000742 | 1764.75 | Hours | 22,776 | 8696 | 8506 | 12,256 | 34,763 | 33,006 | 29,706 |
22 | 0.000725 | 1770.12 | Hours | 30,813 | 7371 | 8453 | 29,958 | 16,345 | 30,776 | 25,358 |
Experiment | Hypervolume Accumulated Value |
---|---|
Encoding | 2.4553 |
Accuracy | 2.4646 |
Total | 2.4651 |
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
Share and Cite
Cacereño, A.; Greiner, D.; Galván, B.J. Multi-Objective Optimum Design and Maintenance of Safety Systems: An In-Depth Comparison Study Including Encoding and Scheduling Aspects with NSGA-II. Mathematics 2021, 9, 1751. https://doi.org/10.3390/math9151751
Cacereño A, Greiner D, Galván BJ. Multi-Objective Optimum Design and Maintenance of Safety Systems: An In-Depth Comparison Study Including Encoding and Scheduling Aspects with NSGA-II. Mathematics. 2021; 9(15):1751. https://doi.org/10.3390/math9151751
Chicago/Turabian StyleCacereño, Andrés, David Greiner, and Blas J. Galván. 2021. "Multi-Objective Optimum Design and Maintenance of Safety Systems: An In-Depth Comparison Study Including Encoding and Scheduling Aspects with NSGA-II" Mathematics 9, no. 15: 1751. https://doi.org/10.3390/math9151751
APA StyleCacereño, A., Greiner, D., & Galván, B. J. (2021). Multi-Objective Optimum Design and Maintenance of Safety Systems: An In-Depth Comparison Study Including Encoding and Scheduling Aspects with NSGA-II. Mathematics, 9(15), 1751. https://doi.org/10.3390/math9151751