Abstract

The redundancy optimization problem is formulated for an uncertain parallel-series system with warm standby elements. The lifetimes and costs of elements are considered uncertain variables, and the weights and volumes of elements are random variables. The uncertain measure optimization model (UMOM), the uncertain optimistic value optimization model (UOVOM), and the uncertain cost optimization model (UCOM) are developed through reliability maximization, lifetime maximization, and cost minimization, respectively. An efficient simulation optimization algorithm is provided to calculate the objective values and optimal solutions of the UMOM, UOVOM, and UCOM. A numerical example is presented to illustrate the rationality of the models and the feasibility of the optimization algorithm.

1. Introduction

The primary goal of reliability design is to improve the reliability of a system. To maintain the reliability to a higher level, the redundancy allocation is an effective method in the system design phase. While improving system reliability by a redundancy method, the cost, weight, and volume also increase. Thus, it is an important topic for system decision-makers to determine the optimal number of redundant elements under certain system constraints.

In the traditional redundancy optimization problem, various kinds of optimization models have been proposed under the assumption that the lifetimes of the elements are random variables. Due to imprecision of data for element lifetimes in certain situations, fuzzy redundancy optimization models [14] are then developed based on fuzzy set theory [5, 6]. Furthermore, Zhao and Liu [7] proposed three redundancy optimization models under the assumption that the lifetimes of the elements are presented as fuzzy variables. Wang and Watada [8] developed two fuzzy random redundancy allocation models for a parallel-series system when the lifetimes of the elements are treated as fuzzy random variables. Recently, some researchers have addressed reliability optimization designs of some systems by considering interval-valued component reliability in an uncertain environment. Roy et al. [9] applied the symmetrical form of interval numbers by interval-valued parametric functional form to evaluate the optimum system reliability and system cost of the redundancy allocation problem. Zhang and Chen [10] investigated an interval multiobjective optimization problem for reliability redundancy allocation of a series-parallel system. Moreover, some researchers concentrated on some hybrid uncertainty optimization problems for system reliability [11, 12].

The probability, interval, and fuzzy theories have been widely used to handle the high level of uncertainty in various real-world applications. With the development of the research on the uncertainty phenomena, the mathematical model based on the probability, interval, and fuzzy theories is not enough to solve all problems, especially when we have no available samples but belief degree from the experts. Belief degree function is a type of distribution function for indeterminate quantity. Since it usually deviates far from the frequency, using probability theory may lead to counterintuitive results. In this case, we should use uncertainty theory. The uncertainty theory provides a useful tool to study reliability modeling and optimization problem of systems with human uncertainty phenomena. The basic uncertainty theory was founded by Liu [13] in 2007. It was refined by Liu [14] in 2010 based on normality axiom, duality axiom, subadditivity axiom, and product axiom. Nowadays, uncertainty theory has become a branch of mathematics for modeling human uncertainty. Many theories and applications have been done based on uncertainty theory, for example, uncertain statistics [15, 16], uncertain programming [17, 18], uncertain logic [1921], uncertain inference [2224], uncertain process [2527], uncertain differential equations [28, 29], and uncertain graph [30].

The system reliability via uncertainty measure was first studied by Liu [31]. Afterwards, Liu [32] investigated reliability of redundant systems with cold and warm redundant elements based on uncertainty theory. Liu et al. [33] studied the reliability and MTTF of unrepairable systems with uncertain lifetimes. Wen and Kang [34] analyzed an uncertain random system based on chance theory which is a generalization of both probability theory and uncertainty theory. Gao et al. [35] studied the reliability of -out-of- systems with uncertain random lifetimes. Gao et al. [36] studied the reliability of the -out-of- system with uncertain weights. Zeng et al. [37] defined belief reliability as an uncertain measure due to the explicit representation of epistemic uncertainty and investigated the belief reliability for coherent systems based on minimal cut sets.

Making use of an uncertain variable as a tool to characterize the lifetimes and costs of elements, we will discuss the redundancy optimization problems for a parallel-series system with warm standby elements in this paper. In this work, three uncertain optimization models are developed, and an efficient simulation optimization algorithm is given to solve these models. In Section 2, some basic concepts and theorems concerning uncertainty theory are presented. The problem formulation of an uncertain parallel-series system is considered in Section 3. Section 4 shows the three uncertain optimization models and gives a solution approach to these models. A numerical example is provided in Section 5, and Section 6 presents a general conclusion.

2. Preliminaries

Definition 1 (see [13, 38]). Let be a -algebra on a nonempty set . A set function is called an uncertain measure if it satisfies the following axioms:

Axiom 1 (normality axiom). for the universal set .

Axiom 2 (duality axiom). for any event .

Axiom 3 (subadditivity axiom). For any countable sequence of events , we have the triple is called an uncertainty space.

Axiom 4 (product axiom). Let be the uncertainty space for . Then, the product uncertain measure is an uncertain measure satisfying where are arbitrarily chosen events from for , respectively.

Definition 2 (see [13]). An uncertain variable is a measure function from an uncertainty space to the set of real numbers; that is, for any Borel set of real numbers, the set is an event.

Definition 3 (see [13]). The uncertainty distribution of an uncertain variable is defined as for any real number .

Definition 4 (see [13]). An uncertainty distribution is said to be regular if it is a continuous and strictly increasing function with respect to at which , and In addition, the inverse function is called the inverse uncertainty distribution of .

Definition 5 (see [17]). A variable is said to be an -optimistic value if is an uncertain variable and where .

Definition 6 (see [13]). An uncertain variable is said to be linear if it has a linear uncertainty distribution. which is denoted by . Apparently, the linear uncertain variable is regular and has an inverse uncertainty distribution .

Definition 7 (see [13]). An uncertain variable is said to be lognormal if has a normal uncertainty distribution. denoted by , where and are real numbers with . The uncertain variable is regular, and its inverse uncertainty distribution is

Theorem 1 (see [13]). Assume that are independent uncertain variables with regular uncertainty distributions , respectively. If the function is strictly increasing with respect to , then has an inverse uncertainty distribution .

Theorem 2 (see [13]). Let be an uncertain variable with a regular uncertainty distribution . Then, .

Theorem 3 (see [13]). Let and be independent uncertain variables with finite expected values. Then, for any real numbers and , we have

Theorem 4 (see [13]). Assume that are independent uncertain variables with uncertainty distributions , respectively. Then, , , and have uncertainty distributions , , and , respectively.

3. Problem Formulation of an Uncertain Parallel-Series System

Consider a warm standby redundant parallel-series system composed of subsystems , and subsystem consists of components connected in series, as shown in Figure 1. The component in subsystem contains one original element and warm standby redundant elements, and .

Throughout the paper, we assume the following for the warm standby redundant parallel-series system under consideration. (1)The system, components, and elements are only in one of two states (up state or down state) at any time.(2)All the elements are independent.(3)The system starts to work at time 0.(4)For each component, there is one element available.(5)The system or element is nonrepairable.(6)The standby redundant elements may deteriorate in the standby period.(7)The lifetime of each element is an uncertain variable, and the standby deterioration rates of all the elements in the standby period are the same constant.(8)Conversion switches of the standby system are absolutely reliable, and the conversion is instantaneous.(9)The costs of all the elements are independent uncertain variables, and the weights and volumes of all the elements are independent random variables.

3.1. Component Lifetime

According to the assumptions, the lifetime of the component in the subsystem can be expressed as where where represents the number of the elements in the component for the subsystem and represents the working lifetime of the th element in the component for the subsystem , and .

For simplicity, it is assumed that the deterioration rates of all the warm standby elements are the same constant, denoted by (very small number). According to [32], the working lifetime of the th element in the component for the subsystem is expressed as where represents the uncertain lifetime of the th element in the component for the subsystem from time 0 and (i.e., the first element starts to work at time 0). Thus, the working lifetime of all elements in the component for the subsystem can be expressed as follows:

By (11) and (14), the lifetime of component in the subsystem can be expressed as

3.2. System Lifetime

The lifetime of subsystem can be obtained by

The lifetime of the system can be determined as

Using (15), (16), and (17), can be written as

4. The System Optimization Model

The purpose of redundancy optimization is to find the optimal solution (the numbers of elements in each component for each subsystem) for improving the system performances under certain constraints. In this section, three different standby redundancy optimization models are considered for different management purposes. In these models, the costs of all the elements are presented as uncertain variables and the weights and volumes of all the elements are presented as random variables.

4.1. The Uncertain Measure Optimization Model

The uncertain measure is used to define the system reliability, and the optimization objective is to maximize the uncertainty measure that the system lifetime is greater than or equal to the given time under the expected cost, expected weight, and expected volume constraints. In this way, the general form of the uncertain measure optimization model (UMOM) is as follows: where , , and represent the uncertain cost, random weight, and random volume of each element in the component for the subsystem , respectively.

Assume that the uncertain lifetimes have the same uncertainty distribution . According to Theorem 4 and [32], the uncertain lifetime distribution functions of the subsystems , denoted by , respectively, can be obtained as where is the uncertain lifetime distribution function of the component in subsystem . The lifetime of the system has an uncertainty distribution function.

According to the system reliability definition and duality of uncertain measure, we have

According to Theorem 3 and the linear operation of the expected value of the random variable, we have

Therefore, the UMOM (19) is equivalent to

4.2. The Uncertain Optimistic Value Optimization Model

If we use the system lifetime to measure the system performance, then the optimization objective is to maximize the -optimistic value of the system lifetime when is equal to the given value under the same constraints as model (19). In this way, the form of the uncertain optimistic value optimization model (UOVOM) can be expressed as where represents the time related to decision vector . By Theorem 1, the inverse uncertainty distribution of the system lifetime is where is the inverse uncertainty distribution of the lifetime of the component in the subsystem . By using the duality of uncertain measure, can be replaced by which is equivalent to . According to the monotonicity of , can reach its maximum when . Therefore, the UOVOM (25) is equivalent to

4.3. The Uncertain Cost Optimization Model

In the following, we give an uncertain cost optimization model with system reliability, expected weight, and expected volume constraints. Since the cost of the system is an uncertain variable, it cannot be directly minimized. Based on the expected value of the system cost, we may minimize its expected value. Then, we have the uncertain cost optimization model (UCOM) as follows:

The UCOM (28) is equivalent to

4.4. Solution Approach for the Optimization Model

It is easy to see that the models (24), (27), and (29) which are equivalent to the UMOM, UOVOM, and UCOM, respectively, are crisp nonlinear integer programming models. An efficient simulation optimization algorithm is provided to calculate the objective values and optimal solutions of the UMOM, UOVOM, and UCOM. The main steps of the proposed algorithm are as follows:

Step 1. Set the range of the integer decision variable to be , , .

Step 2. Randomly generate integer vectors from the interval , and construct a matrix with rows and columns. where the row vector is the possible solution vector.

Step 3. Substitute each row of the matrix into and for the UMOM, UOVOM, and UCOM, respectively. We can obtain the feasible solutions and can also construct feasible solution matrixes for the three different models.

Step 4. Calculate , and based on the row vectors of the feasible solution matrixes in Step 3 for the three different models, respectively. The feasible solutions that make , and the minimum are the best solutions for the three different models, respectively.

Step 5. Change in Step 2, repeat Step 2 to Step 4, observe the changes of the objective function values obtained under different simulation times, and determine the most appropriate number of simulation .

Step 6. Return the best solutions obtained under the value in Step 5 as the final optimal solutions for the three different models.

5. Numerical Example

In this section, we give an illustration of the optimization models and their solutions for a warm standby redundant parallel-series system shown in Figure 2. The detailed data used in this example are those given in Table 1.

In the UMOM, we take , , , , and . The lower and upper bounds of the number of the redundant elements are and , respectively. Then, the model can be expressed as

By using an enumeration algorithm in MATLAB, the exact solution of the model (34) is (4,7,8,10,3) and the corresponding objective function value is 0.9173. That is, the components 11, 12, 21, 22, and 23 have 3, 6, 7, 9, and 2 redundant elements, respectively. They are solved in about 100 s. It is very difficult to use the enumeration algorithm to solve the optimization problem when the decision variables are very large. However, the proposed simulation optimization algorithm in Section 4.4 is very efficient. The effectiveness of the algorithm is illustrated below. Firstly, the simulation times are increased by a step size of 100. Table 2 shows the change of reliability, optimal solution, running time, and error with different simulation times, where the error represents the absolute value of the difference between the real value and the optimal value obtained by our algorithm. We can see from Table 2 that the fluctuation of the error is very obvious and the running time increases as the simulation times increase. Secondly, the simulation times are increased by a step size of 1000, and the data obtained are shown in Table 3. It can be seen from Table 3 that the system reliability is stable at 0.9173 and the error fluctuation is zero. In addition, we can see that completing the computation takes very little time by using the simulation optimization algorithm. Figure 3 illustrates the objective function value (system reliability) versus simulation times. The optimal solutions for different deterioration rates are given in Table 4.

In the UOVOM, we take , , , , and . The lower and upper bounds of the number of the redundant elements are and , respectively. Then, the model can be expressed as

The relationship between the simulation times and the -optimistic value of system lifetime is illustrated in Figure 4 when the step size of simulation times is 5000.

It can be seen that the -optimistic value of system lifetime has a common upper bound, whose value is the optimal value, and the corresponding solution is the optimal solution. The optimal solution is (6,2,1,19,1), and the maximum lifetime of the system is 228.7692. That is, the components 11, 12, 21, 22, and 23 have 5, 1, 0, 18, and 0 redundant elements, respectively. The optimal solutions and the maximum -optimistic value of system lifetime at different degradation rates are obtained in Table 5.

In the UCOM, we take , , , , , , and , then the model can be expressed as

The relationship between the simulation times and system expected cost is illustrated in Figure 5 when the step size of simulation times is 1000. From Figure 5, we can see that the system expected cost under different simulation times has common lower bound, whose value is the minimum expected cost of the system. The optimal solution is (3,1,1,8,1), and the corresponding minimum expected cost is 144.1000. That is, the components 11 and 22 have 2 and 7 redundant elements, respectively. The components 12, 21, and 23 have no redundant elements. The optimal solutions and the minimum expected costs at different degradation rates are given in Table 6.

6. Conclusion

In this paper, an uncertain parallel-series system with warm standby elements was investigated. Under the assumption that the element lifetime and cost are uncertain variables and the weight and volume of an element are random variables, the optimization model of warm standby redundancy for the uncertain parallel-series system was proposed. We formulated three different optimization models—UMOM, UOVOM, and UCOM—based on reliability maximization, lifetime maximization, and cost minimization, respectively. An efficient simulation optimization algorithm was designed to calculate the objective values and optimal solutions of the UMOM, UOVOM, and UCOM. The numerical example showed the rationality of the proposed models and the effectiveness of the simulation optimization algorithm. In the future, the redundancy optimization model with priority will be constructed according to the preference of the decision-maker. Also, the uncertain parallel-series system considered in this paper is a binary system; how to model the uncertain multistate parallel-series system is another future research direction.

Data Availability

In this paper, the authors present a numerical example (the data of the elements are assumed to be uncertain or random variables in the numerical example).

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This research work was supported in part by the Hebei Province Natural Science Foundation of China (Grant no. A2018203088), the National Natural Science Foundation of China (Grant no. 11601469), and the Science Research Project of Hebei Province Education Department (nos. ZD2017079 and QN2016016).