The Korean Society of Marine Engineering
[ Original Paper ]
Journal of the Korean Society of Marine Engineering - Vol. 38, No. 8, pp.1004-1009
ISSN: 2234-7925 (Print) 2234-8352 (Online)
Print publication date Oct 2014
Received 15 Sep 2014 Revised 27 Oct 2014 Accepted 29 Oct 2014
DOI: https://doi.org/10.5916/jkosme.2014.38.8.1004

An efficient metaheuristic for multi-level reliability optimization problem in electronic systems of the ship

JangKil-Woong1 ; KimJae-Hwan
1Department of Data Information, Korea Maritime and Ocean University

Correspondence to: Department of Data Information, Korea Maritime and Ocean University, 727 Taejong-ro, Yeongdo-gu, Busan 606-791, Korea, E-mail: jhkim@kmou.ac.kr, Tel: 051-410-4374

Copyright © The Korean Society of Marine Engineering

The redundancy allocation problem has usually considered only the component redundancy at the lowest-level for the enhancement of system reliability. A system can be functionally decomposed into system, module, and component levels. Modular redundancy can be more effective than component redundancy at the lowest-level because in modular systems, duplicating a module composed of several components can be easier, and requires less time and skill. We consider a multi-level redundancy allocation problem in which all cases of redundancy for system, module, and component levels are considered. A tabu search of memory-based mechanisms that balances intensification with diversification via the short-term and long-term memory is proposed for its solution. To the best of our knowledge, this is the first attempt to use a tabu search for this problem. Our tabu search algorithm is compared with the previous genetic algorithm for the problem on the new composed test problems as well as the benchmark problems from the literature. Computational results show that the proposed method outstandingly outperforms the genetic algorithm for almost all test problems.

Keywords:

Multi-level redundancy allocation problem, Tabu search, Genetic algorithm

Acronyms

RAP   Redundancy allocation problem

MRAP   Multi-level RAP

GA     Genetic algorithm

TS      Tabu search

ACO     Ant colony optimization

SA      Simulated annealing

VNS      Variable neighborhood search

PSO       Particle swarm optimization

HM     Hybrid metaheuristic

TSMRAP Our tabu search for MRAP


Nomenclature

n    The number of basic items

m    The number of constraints

xi    The number of redundancy allocated to the ith basic item

R (x)      The system reliability

(x)    The penalty function of the system reliability

gi (x)   The jth constraint function

ri     The reliability of the ith basic item

ci     The cost of the ith basic item

bj     The amount of allowable resource i

li      The lower bound of xi

ui     The upper bound of xi


1. Introduction

Reliability is considered to be one of most important design measures in various systems. Suppose that we want to maximize the system reliability by allocating redundancy units under constraints on cost, volume, weight and other variables in the system design phase. In addition, we have various information to components, parts, some modules and sub-systems which can be used to design our system. The basic problem is that it is difficult to satisfy the required system reliability with basic modules and components. Thus, we should consider redundant units to improve system reliability as a parallel structure. Redundancy allocation has been used mainly to enhance system reliability. The RAP ([1][2]) which involves selecting redundancy levels at each subsystem in order to maximize system reliability under several resource constraints, is a well-known combinatorial optimization problem. The problem arises frequently in system designs such as semi-conductor integrated circuits, nanotechnology, and most electronic systems of the ship. The RAP has various system structures such as series, series-parallel, complex network, k-out-of-n, and so on [3]. In this paper, we restrict ourselves to series-parallel system.

Solutions for the series-parallel RAP have been suggested by many authors. Fyffe et al. [4] originally set up the problem and suggested a solution algorithm utilizing a dynamic programming approach. Nakagawa and Miyazaki [5] developed 33 variations of Fyffe’s problem, where the weight constraint varied its value from 159 to 191. Coit & Liu [6] proposed zero-one integer programming for small size of problems. They constrained the solution space so that only the identical component type can be allowed for each subsystem. On the contrary, Coit & Smith [7] extended Fyffe’s problem in such a way that the parallel system could be more flexible. They allowed a mixing of component types within a subsystem and employed a GA to obtain optimal solutions. Better solutions have been presented by TS (Kulturel-Konak et al. [8][9]), ACO (Liang & Smith [10]), VNS (Liang & Chern [11]), and HM (Nahas et al. [12], Ouzineb et al. [13]). The above all metaheuristics proposed for only best solutions for 14-subsystem problems without referring to their global optimal solutions. However, Kim & Kim [14] suggested the global optimal solutions for these problems by the transformation of binary integer programming, and also showed that all solutions suggested by TS [9] for the 14-subsystem problems were the global optimum. For lager problems with up to 56-subsystem, Bae et al. [15] proprosed SA as a solution method and compared metaheuristics with global optimal solutions.

In the meanwhile, Yun & Kim [16] proposed a new kind of RAP, that is, MRAP. The traditional RAP is to consider only the component redundancy at the lowest-level. MRAP is to consider all cases of the redundancy for system, module, and component levels. A system can be functionally decomposed into system, module, and component levels. Modular redundancy can be more effective (see Boland & EL-Neweihi [17]) than component redundancy at the lowest-level because in modular systems, duplicating a module composed of several components can be easier, and requires less time and skill, than duplicating each component. Therefore the lower the level of the redundant item and the more spare parts added, the higher the cost of redundancy. Kumar et al. [18] studied a similar problem in which more complex structures are also included and proposed a hierarchical GA. Yun et al. [19] considered the extended problem (multiple MRAP) and proposed a GA. Kim & Jang [20] proposed a modified tabu search for the RAP of complex systems.

In this paper, we consider a MRAP in which all cases of redundancy for system, module, and component levels are considered. A TSMRAP of memory based mechanisms that balances intensification with diversification via the short-term and long-term memory is proposed for its solution. To the best of our knowledge, this is the first attempt to use a TS for MRAP. Our TSMRAP is compared with the previous GA [16] for MRAP on the new composed test problems as well as the benchmark problems from the literature. Computational results show that the TSMRAP outstandingly outperforms the GA for almost all test problems.

This paper is organized as follows. In Section 2, we explain the formulation of the presented model. The developed metaheuristic method, TSMRAP, is illustrated in Section 3. An example of MRAP’s superior solution to the traditional RAP is studied and computational results from several numerical experiments are discussed in Section 4. We conclude in Section 5 with some suggestions for further research.


2. Problem formulation

In MRAP, the number of available redundancy structure increases exponentially as the size of problem becomes large. In this paper, two assumptions are set up to exclude impossible redundancy structures (Yun et al. [19]).

Assumptions:

1) The combination of the basic items for redundancy should satisfy the function at the system level. If a basic item is used, all its sibling items should be used or its function should be satisfied by corresponding child items.

2) The basic items for redundancy should be used in parallel at one combination.

The problem objective is to maximize system reliability, R (x), given constraints on the system, only for system cost in this paper. The system is configured as a series-parallel system. The MRAP optimization model can be generally formulated as the following nonlinear integer programming problem:

This problem was proven to be a NP-hard problem (Chern [21]). We refer interested readers to Yun & Kim [16] for the detailed formulation of MRAP.


3. TS algorithm

In this paper, we propose a TS called the TSMRAP which is based on the TS [8] algorithm for the series-parallel RAP.

3.1 The general steps of TSMRAP

The steps of TSMRAP can be briefly expressed as follows:

Step 0: Generate random initial feasible solutions.

Step 1: Explore all possibilities of the neighborhood through the defined move methods.

Step 2: If the best solution of the step 1 is in the tabu list and is not better than the current best, then repeat Step 1 to find the next best. Otherwise, select the solution as the next best move.

Step 3: If the stopping criterion is satisfied, then stop. Otherwise, go to Step 1.

3.2 Random initial solutions

The initial solutions of TSMRAP are randomly generated, and the scheme of randomly constructing the initial solutions is nearly identical to the general scheme, except that it allows either the module or the component for redundancy, due to the characteristics of our problem. In our experiments, TSMRAP was applied 5 times with different initial solutions for each problem.

3.3 Tabu list

The tabu list is one of the mechanism to prevent cycling and guide the search toward unexplored regions of the solution space. Kultruel-Konak et al. [8] showed that the dynamic size of a tabu list plays an important role in finding the better solutions for RAP. For the long-term memory, the size of a tabu list is reset every 20 iterations to the value of between [n, 3n] uniformly distributed. Once the list is full, the oldest element of the tabu list is removed as a new one is added.

3.4 Penalty function

The TS generally adopts the penalty function to allow to explore the search towardthe promising infeasible regions. The penalty function of our TSMRAP is employed as the same that of TS [8] . The penalty function for only one cost constraint is as follows:

Where Rall is the unpenalized (feasible or infeasible) system reliability of the best solution found so far, Rfeas is the system reliability of the best feasible solution found so far, and ∆c represents the magnitude of the constraint violation. The initial value of NFTc is set to 1% of the cost constraints limit C and k is set to 1, though computational results are not sensitive to this value.

3.5 The structure of generating the neighborhood Solutions

In TSMRAP, neighborhood solutions are generated by two moves, that is, the first and second move. The first move is similar to the TS [8] for the series-parallel system, that is, to change the allocated redundant number of the module or component by adding or subtracting one. The second move is somewhat more complicated than the first move due to the characteristics of the problem. The scheme of the second move is to replace the current solution with the redundant value of the module or component. That is, if the current solution within a subsystem is assigned to the module, then it is replaced with a series of the component. Reversely, if the current solution within a subsystem is assigned to the component, then it is replaced by the module.

For example, let us consider the following system structure with 2 modules and 5 components in Table 1. The cost function is as follows.

An example with two modules

The cost limit is set to 95. Suppose that the current solution during the iterations of TSMRAP to obtain the optimal solution for our problem to be given in Table 2, that is, the encoding status of (2, 0, 0, 0, 0, 2, 1, 0, 0). From the current solution, all the neighborhood solutions generated by the two moves are shown in Table 2. The number of the neighborhood solutions by the first move (adding or subtracting one) are 5 cases shown in Table 2. In the first move, notice that when the redundant value of component (122) equals to 1, it is not allowed to subtract one from the component (122).

All the neighborhood solutions by two moves

The scheme of the second move is to replace the current solution with the redundant value of the module or component. The number of the neighborhood solutions generated by the second move is 4 cases. As shown in Table 2, when the value of module (11) equals to 2, it is replaced with two cases of components (111, 112, 113), that is, (1, 1, 1) and (2, 2, 2). Namely, it is allowed to assign a series of the redundant value of components up to the value of the redundant number of module (11), that is, 2. Reversely, when the value of component (121, 122) equals to (2, 1), it is replaced with two cases of module (12), that is, 1 and 2. We assign the redundant value of module up to the maximum value of the redundant number among components (121) and (122), that is, 2.

In Table 2, the system reliability for each neighborhood solution is evaluated by the penalty function of Eq. (2). For example, (2, 0, 0, 0, 0, 2, 2) has the penalty value of 0.796362 instead of 0.895469. Among 9 candidates of neighborhood solutions, (2, 2, 0, 0, 0, 0, 0), which has the maximum value of the system reliability (0.852996), is selected for the next best move of TSMRAP.


4. An example and computational results

In this section, we conduct three experiments. First, we mention the possibility of having MRAP’s superior solutions to the traditional RAP. Second, we compare the proposed TSMRAP and the previous GA for three cases of levels. Finally, we evaluate the performance of metaheuristics for additional 120 test problems for moderate size of system structure.

4.1 An example of MRAP’s superior solutions to the traditional RAP

We compare the traditional RAP of considering only the component redundancy for the lowest-level with MRAP of allowing the modular and the component redundancy. We consider the 3-level system. Table 3 gives the values of input variables for reliability, price, and additive cost parameter. Test problems are generated by varying the cost limit from 150 to 340. A comparison between the traditional RAP and MRAP is shown in Table 4. As the results in Table 4 indicate, MRAP provides the superior solutions to the traditional RAP for 19 cases of 20 test problems. For only one case (cost limit of 190), the traditional RAP and MRAP have the same result of 0.8878.

The input data for 3-level system

A comparison of the traditional RAP and MRAP

Specifically, for the cost limit of 170, both of the system structuresare shown in Figure 1. In this case, the optimal solution of the traditional RAP and MRAP are given by 0.8455 and 0.8511, respectively. This result indicates that MRAP considering the modular redundancy as well as the component redundancy at the lowest-level provides the possibility of having the better system reliability than the traditional RAP of considering only the component redundancy.

Figure 1:

The system structure

4.2 Computational results

To additionally evaluate the performance of the previous GA[16] and the proposed TSMRAP, the new test problems for moderate size of system structure are designed. They are composed of 4 sets of 10 test problems for each level, that is, totally 120 test problems. For the case of 3-level, we replicated the data in Table 3 by h-fold, where h ranges from 4, 6, 8, to 10, that is, the number of subsystems of each problem became n=3 × h (h=4, 6, 8, and 10). Similarly, for the cases of 4-level, the number of subsystems of each problem became n=2 × h (h =6, 9, 12, and 15). For each set, 10 test problems were generated according to increase the cost limit of the constraint by 100. We increase the cost limit(C) of the constraints appropriately for each problem according to the size of problems in order to assign the reasonable value (not to be very low) of system reliability for each problem. In our numerical experiments, the stopping criterion of TSMRAP was defined as 200 iterations without finding an improvement in the best feasible solution. TSMRAP was applied 5 times with different starting initial solutions for each test problem. The computational results for evaluating the performance between GA and TSMRAP are summarized in Table 5. Two algorithms are coded in C/C++, and experiments are performed on a Pentium IV 3.0 GHz PC. Performances of GA and TSMRAP are assessed in terms of average relative error(A), maximum relative error(M), optimality rate(O) and average execution time(sec.) of 10 problems(T) defined as follows.

O = the number of cases in which each method yields the best solution among 10 problems.

Rj = the system reliability obtained by each method for test problem j.

Rj*= the best system reliability obtained by both of GA and TSMRAP.

Computational results of GA and TSMRAP

As the results in Table 5 indicate, TSMRAP outstandingly outperformed the GA. TSMRAP obtained a higher system reliability in 115 of 120 test cases. TSMRAP failed to obtain the best solution for only two cases in the case of n=30 for 4-level. Even for these cases, the average relative error (A) of GA is very poor, scoring the higher value of 0.0356 than 0.0016 of TSMRAP. GA obtained totally the best solution in only 5 cases, in which two cases are the higher system reliability than TSMRAP and 3 cases are the same system reliability of TSMRAP.

In terms of computing time, GA and TSMRAP are nearly equal for the case n=12 and 18 of 3-level and 4-level. For the cases (n=24 and 30) of 3-level and 4-level, TSMRAP is generally faster than GA. Specifically, when n=30 of 3-level and 4-level, the computing time of GA is almost 3 times of TSMRAP’s.


5. Conclusions

The RAP has usually considered only the component redundancy at the lowest-level for the enhancement of system reliability. A system can be functionally decomposed into system, module, and component levels. Modular redundancy can be more effective than component redundancy at the lowest-levelbecause in modular systems, duplicating a module composed of several components can be easier, and requires less time and skill. This paper dealt with a MRAP in which all cases of redundancy for system, module, and component levels are considered. A TSMRAP of memory-based mechanisms that balances intensification with diversification via the short-term and long-term memory was proposed for its solution. To the best of our knowledge, this is the first attempt to use a TS for MRAP. Our TSMRAP was compared with the previous GA for MRAP on the new composed test problems as well as the benchmark problems from the literature. Computational results showed that the TSMRAP outstandingly outperformed the GA for almost all test problems.

For the further research, we will consider various RAPs, develop efficient metaheuristics and compare their performance. Even if we consider the static RAP, the RAP in dynamic reliability cases may be a promising area and in that case, operations problem can also be considered together with RAP.

Notes

This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0), which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

References

  • W. Kuo, V. R. Prasad, F. A. Tillman, and C. L. Hwang, “Optimal reliability design: Fundamentals and Applications”, Cambridge University Press, (2001).
  • W. Kuo, and V. R. Prasad, “An annotated overview of system-reliability optimization”, IEEE Transactions on Reliability, 49(2), p176-187, (2000).
  • W. Kuo, C. L. Hwang, and F. A. Tillman, “A note on heuristic methods in optimal system reliability”, IEEE Transactions on Reliability, 27(5), p320-324, (1978).
  • D. E. Fyffe, W. W. Hines, and N. K. Lee, “System reliability allocation and a computational algorithm”, IEEE Transactions on Reliability, 17(2), p64-69, (1968). [https://doi.org/10.1109/TR.1968.5217517]
  • Y. Nakagawa, and S. Miyazaki, “Surrogate constraints algorithm for reliability optimization problems with two constraints”, IEEE Transactions on Reliability, 30(2), p175-180, (1981). [https://doi.org/10.1109/TR.1981.5221024]
  • D. E. Coit, and J. Liu, “System reliability optimization with k-out-of-n subsystems”, International Journal of Reliability, Quality and Safety Engineering, 7(2), p129-142, (2000). [https://doi.org/10.1142/S0218539300000110]
  • D. E. Coit, and A. E. Smith, “Reliability optimization of series-parallel systems using a genetic algorithm”, IEEE Transactions on Reliability, 45(2), p254-260, (1996). [https://doi.org/10.1109/24.510811]
  • S. Kulturel-Konak, A. E. Norman, and D. W. Coit, “Efficiently solving the redundancy allocation problem using tabu search”, IIE Transactions, 35(6), p515-526, (2003). [https://doi.org/10.1080/07408170304422]
  • S. Kulturel-Konak, B. A. Norman, D. W. Coit, and A. E. Smith, “Exploiting tabu search memory in constrained problems”, INFORMS Journal on Computing, 16(3), p241-254, (2004). [https://doi.org/10.1287/ijoc.1030.0040]
  • Y. C. Liang, and A. E. Smith, “An ant colony optimization algorithm for the reliability allocation problem problem (RAP)”, IEEE Transactions on Reliability, 53(3), p417-423, (2004).
  • Y. C. Liang, and Y. C. Chen, “Redundancy allocation of series-parallel systems using a variable neighborhood search algorithm”, Reliability Engineering and System Safety, 92(3), p323-331, (2007). [https://doi.org/10.1016/j.ress.2006.04.013]
  • N. Nahas, M. Nourelfath, and D. Ait-Kadi, “Coupling ant colony and the degraded ceiling algorithm for the redundancy allocation problem of series-parallel systems”, Reliability Engineering and System Safety, 92(2), p211-222, (2007). [https://doi.org/10.1016/j.ress.2005.12.002]
  • M. Ouzineb, M. Nourelfath, and M. Gendreau, “An efficient heuristic for reliability design optimization problems”, Computers and Operations Research, 37(2), p223-235, (2010). [https://doi.org/10.1016/j.cor.2009.04.011]
  • J. H. Kim, and J. S. Kim, “Globally solving the redundancy allocation problem for the case of series-parallel systems”, Proceedings of the 2nd Asian International Workshop of Advanced Reliability Modeling Ⅱ, p150-157, (2006). [https://doi.org/10.1142/9789812773760_0018]
  • C. O. Bae, H. G. Kim, J. H. Kim, J. Y. Son, and W. Y. Yun, “Solving the redundancy allocation problem with multiple component choices using metaheuristics”, International Journal of Industrial Engineering, 14(4), p315-323, (2007).
  • W. Y. Yun, and J. W. Kim, “Multi-level redundancy optimization in series systems”, Computers and Industrial Engineering, 46(2), p337-346, (2004). [https://doi.org/10.1016/j.cie.2003.12.013]
  • P. Boland, and E. EL-Neweihi, “Component redundancy vs. system redundancy in the hazard rate ordering”, IEEE Transactions on Reliability, 44(4), p614-619, (1995). [https://doi.org/10.1109/24.475980]
  • R. Kumar, K. Izui, Y. Masataka, and S. Nisiwaki, “Multilevel redundancy allocation optimization using hierarchical genetic algorithm”, IEEE Transactions on Reliability, 57(4), p650-661, (2008). [https://doi.org/10.1109/TR.2008.2006079]
  • W. Y. Yun, Y. M. Song, and H. G. Kim, “Multiple multi-level redundancy allocation in series systems”, Reliability Engineering and System Safety, 92(3), p308-313, (2007). [https://doi.org/10.1016/j.ress.2006.04.006]
  • J. H. Kim, and K. W. Jang, “A modified tabu search for redundancy allocation of complex systems of ships”, Journal of the Korean Society of Marine Engineering, 38(2), p225-232, (2014). [https://doi.org/10.5916/jkosme.2014.38.2.225]
  • M. S. Chern, “On the computational complexity of reliability redundancy allocation in a series system”, Operations Research Letters, 11(5), p309-315, (1992). [https://doi.org/10.1016/0167-6377(92)90008-Q]

Figure 1:

Figure 1:
The system structure

Table 1:

An example with two modules

Unit Parent unit Reliability Price λ
1(system) - 0.40029 72 2
11 1 0.72675 26 2
12 1 0.76500 19 3
111 11 0.90000 5 3
112 11 0.95000 6 4
113 11 0.85000 5 4
121 12 0.90000 6 4
122 12 0.85000 7 4

Table 2:

All the neighborhood solutions by two moves

x=(11, 12, 111, 112, 113, 121, 122) c(x) Rp (x)
current sol. ( 2, 0, 0, 0, 0, 2, 1) 95 0.778669
1st Move Adding ( 3, 0, 0, 0, 0, 2, 1) 125 0.647355
( 2, 0, 0, 0, 0, 3, 1) 149 0.438873
( 2, 0, 0, 0, 0, 2, 2) 114 0.796362
Subtracting ( 1, 0, 0, 0, 0, 2, 1) 67 0.611560
( 2, 0, 0, 0, 0, 1, 1) 77 0.707881
2nd Move ( 0, 0, 1, 1, 1, 2, 1) 66 0.611560
( 0, 0, 2, 2, 2, 2, 1) 112 0.727356
( 2, 1, 0, 0, 0, 0, 0) 78 0.707881
( 2, 2, 0, 0, 0, 0, 0) 103 0.852996

Table 3:

The input data for 3-level system

Unit Parent unit Reliability Price λ
1(system) - 0.40029 72 2
11 1 0.72675 26 2
12 1 0.76500 19 3
13 1 0.72000 21 2
111 11 0.90000 5 3
112 11 0.95000 6 4
113 11 0.85000 5 4
121 12 0.90000 6 4
122 12 0.85000 7 4
131 13 0.90000 8 3
132 13 0.80000 7 4

Table 4:

A comparison of the traditional RAP and MRAP

Cost limit Traditional RAP MRAP
Total Cost Rel. Total Cost Rel.
150 150 0.7687 149 *0.8057
160 150 0.7687 158 *0.8309
170 168 0.8455 169 *0.8511
180 168 0.8455 175 *0.8668
190 186 *0.8878 186 *0.8878
200 186 0.8878 199 *0.9010
210 209 0.8959 202 *0.9136
220 209 0.8959 215 *0.9272
230 212 0.8959 228 *0.9319
240 239 0.9052 228 *0.9319
250 241 0.9174 241 *0.9457
260 241 0.9174 241 *0.9469
270 264 0.9258 270 *0.9609
280 264 0.9258 270 *0.9609
290 290 0.9342 270 *0.9609
300 296 0.9354 270 *0.9609
310 296 0.9354 304 *0.9755
320 317 0.9439 304 *0.9755
330 317 0.9439 304 *0.9755
340 317 0.9439 304 *0.9755

Table 5:

Computational results of GA and TSMRAP

3-level 4-level
GA TSMRAP GA TSMRAP
n=12 A
M
O
T
0.0084
0.0208
1/10
2.7
0.0
0.0
10/10
4.0
0.0143
0.0241
1/10
28.9
0.0
0.0
10/10
27.3
n =18 A
M
O
T
0.0155
0.0446
0/10
10.7
0.0
0.0
10/10
10.6
0.0205
0.0349
0/10
54.5
0.0
0.0
10/10
59.9
n=24 A
M
O
T
0.0118
0.0354
1/10
44.7
0.0
0.0
10/10
21.9
0.0166
0.0384
0/10
163.6
0.0
0.0
10/10
104.5
n=30 A
M
O
T
0.0118
0.0302
0/10
90.1
0.0
0.0
10/10
34.9
0.0356
0.1231
2/10(*)
461.5
0.0016
0.0133
8/10(*)
169.4