Fuzzy sets membership functions integrated with logistic map as the chaos generator were used to create reliability bifurcations diagrams of the system with redundancy of the components. This paper shows that increasing in the number of redundant components results in a postponement of the moment of the first bifurcation which is considered as most contributing to the loss of the reliability. The increasing of redundancy also provides the shrinkage of the oscillation orbit of the level of the system’s membership to reliable state. The paper includes the problem statement of redundancy optimization under conditions of chaotic behavior of influencing parameters and genetic algorithm of this problem solving. The paper shows the possibility of chaos-tolerant systems design with the required level of reliability.
1. Introduction
The classical reliability theory [1, 2] is based on the probabilistic approach. Essential limitations of this approach are connected with “the problem of the source data” which depend on many factors and which may not correspond to the real conditions of the system’s functioning. Besides, the statistical data used in the probabilistic reliability models fix only the facts of real failures and do not contain the information about the causes of these failures. Whereas the causes of these failures are connected with the elements’ variables (temperature, humidity, tension, etc.), which become more (or less) than a certain critical level. So, we can affirm that the probabilistic theory [1, 2] models the reliability in the space of events effects (i.e., failures) and suits badly for the reliability modelling in the space of events causes (i.e., variables).
The alternative for the probabilistic modeling of the reliability is the approach based on the fuzzy logic [3] and related possibility theory [4]. In this case the classical “failure probability” is replaced by “failure possibility” which is modeled by the membership function of the system (or the element) variables to the reliable state (Figure 1).
Relationship of the probability theory and fuzzy logic in reliability estimation.
The explicit dependence of the membership function on the variables (failure causes) makes convenient the integration of the fuzzy model of reliability with the technique of time series [5], which allows observing the change of the reliability level in the real time.
The chaos theory is a new approach to the analysis of nonlinear time series [6]. It uses the conceptual apparatus of the theory of nonlinear oscillations [7] and purposes to study the phase portrait of the dynamical system with its intrinsic states of stability (attractors) and bifurcations, that is, “jumps” between stable states. Unlike the classical oscillation theory [7] where the phase portrait is formed on the base of the system description by means of differential equations, the chaos theory [6] offers the methods of the phase-portrait extraction from the experimental data, that is, directly from the time series.
The integration of fuzzy reliability model with the phase portraits of variables creates preconditions for the construction of phase portrait reflecting the system reliability dynamics. The pattern of bifurcations which may be interpreted as failure instants is of particular interest.
The works on the fuzzy reliability theory began in 90s of the last century. The first specialized collection of papers in this field is the edited volume [8]. The first monographs containing the approaches to the construction of fuzzy reliability theory are [9–11]. In the works [12, 13] there is a fuzzy algorithmic approach to reliability modeling based on the algebra of regular algorithms [14] and expert assessments of the performance correctness of operators and conditions by means of membership functions [3]. The problem of reliability optimization of the control resources based on fuzzy algorithmic approach is solved in [15].
An approach to the online reliability evaluation based on the integration of fuzzy logic and forecasting methods of time series (exponential smoothing and Kalman filtering) is proposed in [16].
The idea of chaos theory application in the reliability modeling appears in [17]. Two real data bases about software failures are processed by the methods of chaos theory in the paper [18]. It was shown that the deterministic model of failures is more adequate to the experimental data than the traditional stochastic models, for example, the modified Poisson’s law, and so forth. The results of the work [18] can be considered as a new approach (alternative to the statistical) to the data processing about the failures on the level of elements. We do not know any publications about the chaos theory applications to the reliability modeling of the system taking into account its structure.
In this connection there are following the questions.
How does the system structure influence the phase portrait of reliability?
Is it possible to solve the redundancy optimization problem with deterministic (chaotic) order of the occurrences of failures?
As far as the failures are connected with the oscillations of variables, then the answers for these questions should be searched on the base of integration of the fuzzy logic and chaos theory.
In this paper we use the fuzzy algorithmic approach to parametric reliability modeling proposed in [12, 13], and the simple generator of chaotic oscillations of variables in the form of logistic function [19]. The further description is organized in the following way.
Section 2 describes the principles of reliability dynamics modelling by means of the composition of membership function and a generator of chaos.
Section 3 examines the fuzzy reliability model of an element with a multiple redundancy. We consider the results of computer experiment on the analysis of reliability level bifurcations depending on the multiplicity of redundancy.
Section 4 considers the redundancy optimization problem under the conditions of chaotic oscillations of the parameters of the elements of systems.
2. Basic Principles
The fuzzy chaotic approach to the reliability dynamics modeling is based on the following principles.
2.1. Fuzzy Correctness
According to this principle introduced in [12, 13], there is not a crisp boundary between “correct” (1) and “incorrect” (0) results of the functioning of a system and its elements. For the formal evaluation of the correctness level it is used the multidimensional (by the number of variables) membership function μ1(x1,x2,…) which depends on the measured parameters (input variables). The correctness of each variable is determined by the membership function μ1(xi) of the variable xi to the correct value.
The function μ1(xi) can be interpreted as the correctness distribution of the variable xi: extreme cases μ1(xi)=1(0) correspond to the maximal (minimal) level of the correctness of the variable xi. Pay attention that the correctness distribution μ1(xi) satisfies the axioms of fuzzy sets theory [3], in contrast to the probabilistic distributions used in the classical reliability theory [1, 2].
The typical correctness distributions (membership functions) are represented in Figure 2. They correspond to three possible cases of fuzzy boundaries between “correct” (1) and “incorrect” (0):
correct (1)—incorrect (0),
incorrect (0)—correct (1)—incorrect (0),
incorrect (0) correct (1).
The typical correctness distributions.
2.2. Integration of Membership Functions and Time Series
It is assumed that for the variable x it is known the time series of its values (x1,x2,x3,…) in discrete moments of time (t1,t2,t3,…). Putting these values in the membership function μ1(x), we receive the dynamics of the correctness level of the variable x in the form of the function μ1[x(t)], Figure 3.
Integration of membership functions and time series.
2.3. Chaos Generator
Chaos means the oscillations which seem random but in truth they are generated by the deterministic nonlinear model. In [6] about 40 models of the chaos generators are described. Each model contains variables whose values must be fitted on the base of the experimental data. Algorithm of the chaos generation is explained by means of iterative Lamerey diagram, widely used in the classical theory of nonlinear oscillations [7].
It is assumed that there is a known function xm+1=f(xm), connecting two neighboring elements of time series: xm and xm+1. Iterative diagram consists of this function and the bisector xm=xm+1 (Figure 4). Choosing the initial point x0 by means of vertical and horizontal lines we obtain the points on the axis xi as follows:
(1)x1=f(x0),x2=f(x1),x3=f(x2),….
Iterative Lamerey diagram.
The most popular generator of chaos is given by the logic map [19] as follows:
(2)xm+1=f(xm)=axm(1-xm),0≤xm≤1,
where 0≤a≤4 is the control parameter determining the nature of chaotic orbits.
Using the iterative equation (2) we can generate the consequence x1,x2,… for the given parameter a and initial point x0.
For example, if x0 = 0.25 and a=0.5, then we find that
With the corresponding values of the parameter a it is possible to get different types of attractors (Figure 5) by means of iteration algorithm (Figure 4):
Different types of logistic map attractors: (a) stable focus, (b) stable orbit, (c) double orbit, and (d) chaotic orbit.
Increasing gradually the parameter a, it is possible to observe the moments of bifurcations, that is, transitions from one type of the attractor to another. Figure 6 shows that in the moment a=3 there is a jump from one stable state to two other stable states. In the moment a=3.449… the number of stable states is doubled, and so forth. The moment a=4 corresponds to the complete chaos [19]. The described chaos generator will be used further for the integration with the fuzzy reliability model.
Bifurcation diagram of logistic map.
3. Fuzzy Chaotic Reliability Model
We consider the simple system with redundancy which is modeled by fuzzy algorithmic approach proposed in [12, 13] and logistic function (2).
3.1. Element with Redundancy
The element with redundancy is presented in Figure 7 in the form of the parallel circuit where the primary element (A) has y of redundant elements (y=0,1,2,…). All the elements are supposed to be homogeneous.
Element with redundancy.
The quality of the element A functioning depends on the variable x which varies during the time: x=x(t). To evaluate the reliability of the element A it is used: μA1(x) is a membership function which determines correctness distribution of variable x during the functioning of A.
The parallel circuit (Figure 7) assumes that the failure of the system requires the failure of all (y+1) elements similar to A. That is why the correctness of the element with redundancy functioning is evaluated by the following formula:
(3)μB1(y)=1-[1-μA1(x)]y+1.
3.2. Reliability Bifurcation
The model (3) allows observing the dynamics of the system reliability level, that is, of the correctness function μB1(y) during the chaotic oscillations of the variable x according to the logistic map (2).
The purpose of the computer experiment consisted of the research of bifurcations of the correctness level μB1(y) with different correctness distributions of the element A and different redundancy rates (y=1,3,5,7).
The experiment is carried out with two correctness distributions shown in Figure 8: triangle (a) and threshold (b). During the chaos generation the parameter a of the logistic map (2) was changed in the range from 2.5 to 4. For each distribution (Figure 8) we obtained 4 bifurcation diagrams, each of them corresponds to different redundancy rates (y). The results are represented in Figures 9 and 10, where the horizontal axis is the chaos parameter (a), and vertical axis is the reliability level μB1(y).
Correctness distributions of variable x in the experiment.
Reliability bifurcations for triangle correctness distribution.
Reliability bifurcations for threshold correctness distribution.
Figures 9 and 10 show that in spite of the chaos growth (parameter a) by the increasing of redundancy rate (y), it is possible
to postpone the moment of the first bifurcation which is associated with the reliability loss and
to decrease the diameter of an orbit around which there are oscillations of the level of system’s membership to the stable state.
That is why we can consider a redundancy optimization problem under chaotic oscillations of the parameters of elements.
4. Redundancy Optimization under Chaosof Parameters
We consider a sequential system where each element has some level of redundancy. This system is described by the series-parallel structure (Figure 11), where Ai is a component which depends on the variable xi, yi is the redundancy rate of the component Ai, and Y=(y1,y2,…,yn) is the vector of the redundancy rates yi, i=1,2,…,n.
Sequential system with redundancy of elements.
It is supposed to be known that
μAi1(xi) is correctness distribution of variable xi during the element Ai functioning,
[xi_,xi¯] is the range of possible values of the variable xi and
ci is mean cost of one redundant component like Ai.
For the system in the Figure 11 taking into account (3) we have
(4)μB1(Y)=∏i=1n{1-[1-μAi1(xi)]yi+1}(5)CB(Y)=∑i=1nciyi,
where μB1(Y) is the correctness of system’s functioning with redundancy rates given by vector Y, and CB(Y) is the total expenses for the redundancy.
It is supposed that the variable xi during the functioning of the component Ai makes the chaotic oscillations according to the law (2) as follows:
(6)xi(m+1)=axi(m)[1-xi(m)],
where xi(m) is the value of the variable xi on the mth step of the logistic map.
One of the practically important problems of the redundancy optimization under chaotic oscillation of parameters can be formulated as follows:
to find the vector Y=(y1,y2,…,yn) which provides
(7)CB(Y)→min,
subject to
(8)μB1(Y)∈[μB*,1],a=a*,
where μB* is the minimal acceptable level of the correctness of system’s functioning with redundancy of the elements.
a* is the parameter of the logistic map (6) determining the level of chaos.
The problem was solved for n=5, ci=1, i=1,2,…,5 and the correctness distributions from the Table 1. It was supposed that the variables xi for all the components Ai have the same level of chaos (parameter a*).
Correctness distribution of elements.
In order to solve the problem of optimization we used the genetic algorithm from the MATLAB presented in the Appendix.
Some of the optimal solutions are represented in Tables 2 and 3 for two minimum acceptable levels of the system functioning correctness (0.8 and 0.9) and four levels of chaos (a*).
Optimal redundancy rates (μB*=0.8).
Chaos parameter (a*)
Y
Criteria
y1
y2
y3
y4
y5
CB(Y)
μB(Y)
2.9
2
1
0
1
2
6
0.8167
3.2
5
4
4
5
3
21
0.8065
3.5
5
5
7
9
9
35
0.8043
3.6
6
5
8
13
13
45
0.8027
Optimal redundancy rates (μB*=0.9).
Chaos parameter (a*)
Y
Criteria
y1
y2
y3
y4
y5
CB(Y)
μB(Y)
2.9
3
1
0
2
2
8
0.9017
3.2
4
4
3
6
7
24
0.9015
3.5
8
6
7
12
11
44
0.9006
3.6
8
8
13
13
16
58
0.9002
Tables 2 and 3 show that by increasing redundancy rate under conditions of chaotic oscillations of parameters of elements it is possible to achieve the stabilization of system functioning correctness level in the required interval [μB*,1]. In the worst case system design the chaos parameter a* can be chosen as the biggest.
It is important to mention that the classical redundancy optimization problem [20] was solved without any statistic data about the reliability of elements. Instead of such a data we used expert information about fuzzy set membership functions describing the correctness distributions of parameters which influence the reliability.
5. Conclusion
The integration of fuzzy logic and chaos theory is a convenient methodology of the system’s reliability dynamics observation in the space of parameters connected with a failure. Using parameters correctness distribution in the form of fuzzy sets membership function, and a chaos generator in the form of logistic map we obtained the bifurcation diagrams of the system reliability functioning with redundancy of the elements.
It is shown that increasing the number of redundant elements postpones the moment of the first bifurcation connected with a loss of reliability, and decreases the orbit size around which there is an oscillation of the level of system’s membership to the reliable state.
It was proposed the problem statement of redundancy optimization under chaotic oscillations of parameters, and the example of its solution by genetic algorithm is provided.
The complex systems reliability dynamics modelling by integration of fuzzy logic and chaos theory can be considered as a perspective direction of the future research.
Appendix
Genetic algorithms are well-known methodology for the reliability optimization problem solving [21]. It is based on a definition and coding of initial population of the solutions and application to them the operations of crossover, mutation, and selection according to the given fitness function. The convenient tool of application of this methodology is MATLAB which provide the automation of coding and basic operation. To apply the MATLAB tools we need to define the fitness function.
P={Yj∣j=1,2,…,L} is a population, where Yj=(y1j,y2j,…,ynj) is a jth chromosome solution, and L is a size of population.
Operations of crossover and mutation were implemented in MATLAB. Mutation coefficient was chosen experimentally and was equal to 0.1. Fitness function FF(Yj) was chosen with the respect to the total expenses for the redundancy (5) and penalty function PF(Yj) as follows:
(A.1)FF(Yj)=∑i=1nci·yij+PF(Yj).
The use of penalty function allows eliminating solutions, which do not satisfy condition (8), by assigning them the value (1000), which is much greater than first member of summation in (A.1). This number 1000 was chosen experimentally. So we have
(A.2)PF(Yj)={0,ifcondition(8)issatisfied1000,otherwise.
As an example, Figure 12 shows the results redundancy optimization by genetic algorithm in MATLAB when μB*=0.8 and a*=2.9.
Solution found with genetic algorithm.
The upper part of the Figure 12 shows the dynamics of fitness function changing with the growth of the generations number. The lower part shows the current optimal solution Y=(2,1,0,1,2). It took 300 generations to find the optimal solution with population size of L=250.
GnedenkoB. V.BelyaevY. K.SolovyevA. D.1969New York. NY, USAAcademic PressBarlowR. E.ProshanF.1965London, UKWileyZadehL.The concept of a linguistic variable and its application to approximate reasoning—I–III197583199251no. 4, pp. 301–357, 1975—no. 9, pp 43–80, 1976ZadehL. A.Fuzzy sets as a basis for a theory of possibility1978113282-s2.0-49349133217BoxG. E. P.JenkinsG. M.ReinselG. C.19943rdEnglewood Cliffs, NJ, USAPrentice HallSprottJ. C.2003Oxford University PressButeninN. V.NejmarkY. I.FufaevN. A.1987Moscow, RussiaNaukaOnisawaT.KasprzykJ.Reliability and safety analysis under fuzziness19954Phisika, A Springer CompanyCaiK. Y.1996Boston, Mass, USAKluwer Academic PublishersRotshteinA.ShtovbaS.1997Vinnitsa, UkraineContinent-PRIMUtkinL.ShubinskyI.2000St. Petersburg, RussiaLubavich Publishing HouseRotshteinA.Fuzzy reliability analyses of human's algorithms activity2007212318RotshteinA.Fuzzy-algorithmic analysis of complex systems reliability2011476919931GlushkovV. M.Automata theory and formal microprogram transformations196515182-s2.0-3425050754210.1007/BF01071417KatelnikovD. I.RotshteinA. P.Fuzzy algorithmic simulation of reliability: control and correction resource optimization20104969679712-s2.0-7995174508010.1134/S1064230710060146KolarikW. J.WoldstadJ. C.LuS.LuH.Human performance reliability: on-line assessment using fuzzy logic20043654574672-s2.0-184273900510.1080/07408170490426161ZouF. Z.LiC. X.Chaotic model for software reliability20012432812912-s2.0-0035284597DickS.BethelC. L.KandelA.Software-reliability modeling: the case for deterministic behavior20073711061192-s2.0-3384574534310.1109/TSMCA.2006.886364MayR. M.Simple mathematical models with very complicated dynamics197626155604594672-s2.0-0017185443KozlovB. A.UshakovI. A.1975Moscow, RussiaSoviet RadioGenM.ChengR.1997John Wiley & Sons