Reconfiguration point decision method based on dynamic complexity for reconfigurable manufacturing system (RMS)

To address the problem of how to identify the best time to implement reconfiguration for the reconfigurable manufacturing system (RMS), a dynamic complexity-based RMS reconfiguration point decision method is proposed. This method first identifies factors that affect RMS dynamic complexity (including both positive and negative complexity) at the machine tool and manufacturing cell levels. Next, based on information entropy theory, a quantitative model for RMS dynamic complexity is created, which is solved via state probability analysis for processing capability and the processing function. This model is combined with cusp catastrophe theory to establish an RMS reconfiguration decision model. Both positive and negative complexity are control variables for cusp catastrophe. Cusp catastrophe’s state condition is used to identify RMS state catastrophe at the final stage of production. This catastrophe point is the RMS reconfiguration point. Finally, the case study result shows that this method can effectively identify the RMS state catastrophe moment so that system reconfiguration is implemented promptly to improve RMS’s responsiveness to the market.


Introduction
In 2010, the German government proposed the concept of Industry 4.0 in its report "High-Tech Strategy (2010)". That report marked the prologue to the fourth Industrial Revolution, which is dominated by intelligent manufacturing. Lee Jie, Professor and Distinguished Chair at The University of Cincinnati (US) and director of the National Science Foundation's Intelligent Maintenance System Industry/University Cooperative Center, believes that there are three pillars in the development of Industry 4.0. One of them is that "manufacturing systems, including the machine tool itself, are able to make adjustments automatically and promptly according to differences in the processed product and changes in processing condition, and possess the so-called 'self-reflection' capability" (Lee et al. 2015). The RMS reconfiguration capability satisfies the requirement of intelligent manufacturing and is a critical component in future intelligent manufacturing environments. Mehrabi et al. (2000aMehrabi et al. ( , 2002 provide an outlook for RMS research focus and future development and suggest that RMS is a desirable next step in the evolution of production systems. The reconfiguration moment decision is essential for RMS to be the backbone of intelligent manufacturing; it is also a critical and fundamental problem waiting for a solution on the road to intelligent manufacturing. RMS is a manufacturing system with fast response capability. Its architecture, hardware and software can be reassembled and adjusted to rapidly change the system function to respond to impacts from market fluctuation, technology innovation and policy change (Koren et al. 1999;Koren 2013;Mehrabi et al. 2000b). Because of its huge market potential, researchers have studied various aspects of RMS, including reconfigurability (Wang 2000), system layout planning (Goyal et al. 2012;Wu et al. 2007), system performance analysis (Cai 2004), part-family formation (Wang et al. 2016a;Goyal et al. 2013;Gupta et al. 2013;Hasan et al. 2014) and reconfiguration schemes evaluation (Wang et al. 2016b). Zhao et al. (2000aZhao et al. ( , b, 2001a provide a stochastic RMS model that includes a framework, optimal configuration, optimal selection policy and performance measurement. Bi et al. (2008) provides an overview of RMS. Renzi et al. (2014) elaborate the advantages of an RMS in terms of cost and efficiency and conduct an in-depth analysis of key technologies in RMS. Although RMS has excellent architecture that can respond to rapid changes in market demand (Wang 2000), the issue of when to implement such a change is a complex and difficult decision problem. The decision about when to implement RMS reconfiguration should consider various production factors such as cost, time and order; a conclusion based on the analysis of any single factor is not convincing. Therefore, this paper has analyzed the RMS reconfiguration point from the system's perspective. Researchers have conducted comprehensive analyzing manufacturing systems via system state. Rao (2006) provide a quantitative description of manufacturing system state complexity from the perspective of information theory; they analyze the manufacturing system's static complexity and dynamic complexity and create an information entropy-based static and dynamic complexity measurement model. To address deficiencies in the manufacturing system architecture complexity modeling method, Duan (2012) analyze the effect of the intermediate buffer zone state on the manufacturing system state; they also employ the universal generating function and information entropy theory to create and evaluate a model for non-serial manufacturing system complexity. Smart et al. (2013) study the manufacturing system's dynamic complexity based on improved information entropy theory; their main focus is an analysis of equipment and queue state in the production process and an elaboration of the relationship between system complexity and system operation. Zhang (2011) analyze the manufacturing system dissipation structure; create a manufacturing system entropy variation model based on information entropy theory; analyze the relationship between the manufacturing system control variable and the state variable based on catastrophe theory; and perform a quantitative analysis of system dissipation structure formation. The aforementioned manufacturing state analysis is only based on a part state or the machine tool state, whereas the manufacturing system state is determined by interaction between parts and machine tools. The part type demands a specific processing function from the manufacturing system; the part quantity demands a corresponding processing capability from the manufacturing system. As provider of processing function and processing capability, the machine tool provides one or more processing functions and processing capabilities of specific quantities. To summarize, the analysis of manufacturing system complexity in terms of processing function and processing capability fully describes the interaction between parts and machine tools, reflects the nature of the manufacturing system state, and thus provides evidence for the RMS reconfiguration decision. Therefore, this paper proposes a dynamic complexitybased RMS reconfiguration point decision method. Based on an analysis of the RMS reconfiguration mechanism, the effect of system processing function and processing capability on system complexity is investigated and an information entropy-based RMS complexity quantitative model is created. This model is combined with a cusp catastrophe to analyze system state variation and identify the RMS state catastrophe moment.

RMS reconfiguration mechanism and problem analysis
After RMS construction is completed, a running-in period is required; that period is called the ramp-up time (Rösiö and Säfsten 2013). During this period, system stability is relatively low, and problems such as machine failure and substandard products are likely to emerge. After a period of tuning and running-in, problems in the new system are gradually fixed; the system operates in a highly efficient and stable state and is capable of high-quality, high-yield and low-cost production. This period is called the stable production period. As RMS continues its operation, because of random internal and external factors such as machine tool failure and new production order queue-jumping, the system fails to provide sufficient processing functions and capability. The leads to delayed order delivery, more difficult scheduling, declined productivity and chaos in production, which prompts system reconfiguration. At this moment, the system runs into the end stage of production. RMS stops production at the right moment and enters the reconfiguration period. Machine tools and components are added, rearranged, removed and adjusted to reconfigure the RMS processing function and processing capability. When reconfiguration is completed, the RMS processing function and processing capability are updated; system capability is able to meet new market demand; and RMS starts a new round of production cycles. In this cyclic process, the quantitative accumulation of machine tool blocks and failures leads to a qualitative change in the system state, or state catastrophe. After a catastrophe, system-level reconfiguration is required and scheduling alone is unable to solve system problems. The key to this process is to identify the catastrophe moment, which is also this paper's focus. Every RMS reconfiguration updates system capability and thus always maintains a relatively high market response capability and productivity. Figure 1 shows the RMS reconfiguration mechanism and implementation process. That diagram shows that after a temporary ramp-up time, RMS can provide processing function and processing capability for current production requirements, and the system is experiencing stable production. At the end stage of production, because of the combined effect of system internal and external factors, three problems arise in the RMS processing function and processing capability system: (1) the processing-function requirement is met, but the processing capability is inadequate; (2) processing capability is adequate, but the processing function is inadequate; and (3) both processing function and processing capability are inadequate. At this moment, system operation efficiency declines and system production capacity gradually declines, leading to delayed order delivery. At a certain moment, processing function and processing capability deficiency leads to manufacturing system state catastrophe; the system enters a reconfiguration period to rebuild processing function and processing capability and then moves to the next production cycle. The RMS state is the result of interaction between processing function and processing capability driven by orders. Factors such as order fluctuation and machine tool failure lead to variation in processing function and processing capability, resulting in system state fluctuation. RMS complexity measures fluctuation in a system's processing function and processing capability. System complexity variation reflects variations in system state stability. In a normal situation, as the system processing function gradually increases, the processing function and processing capability fluctuation are relatively stable, system complexity gradually increases, and system state stability gradually increases. At the end stage of production, the system processing function gradually becomes saturated, the processing function and processing capability fluctuation increase; and system complexity continues to increase. However, system state stability demonstrates a trend of accelerated decline. Figure 2 shows the Complexity Stability Fig. 2 Relationship between RMS complexity and stability relationship between system complexity and system state stability (Qin 2004).
When the RMS processing function or processing capability fails to meet order demand, if reconfiguration is implemented immediately, the normal result is a relatively high cost attributable to production down time. When there are a large number of delayed orders, system reconfiguration is more difficult and costly. The decision about the reconfiguration moment becomes a critical step in RMS implementation. In other words, at the end stage of production, system complexity variation is analyzed to identify the catastrophe point of the system stable state. When system stable state catastrophe occurs, continued production normally results in a dramatic increase in cost. Moreover, if production continues until the system is in an extremely unstable state, system reconfiguration requires significantly higher cost and a much longer time. Therefore, when a system stable state catastrophe occurs, it is the best time to implement system reconfiguration. This reconfiguration moment is the RMS reconfiguration point. In this paper, based on information entropy, system complexity variation is analyzed quantitatively and is combined with cusp catastrophe to complete the decision about the RMS reconfiguration point.

RMS complexity analysis
Manufacturing system complexity variation should be considered in the RMS reconfiguration decision process. There are two categories of manufacturing system complexities: static complexity and dynamic complexity. Static complexity is primarily concerned with the production state of the manufacturing system according to its scheduling plan. However, in actual production, because of uncertain factors such as market fluctuation and machine failure, production normally deviates from the scheduling plan. Dynamic complexity describes a manufacturing system's actual operational state, which represents system state variation in actual production. Therefore, in this paper, system complexity is simplified as dynamic complexity and RMS state analysis is based on dynamic complexity.
System complexity is a measurement of system uncertainty. Information entropy (Czinner and Mena 2016;Duan 2012;Rao 2006;Tsujimura and Gen 1999;Shannon 1948) represents the amount of information contained in a system state. Information entropy accurately reflects the system complexity variation. When system complexity is very low, system uncertainty is low, which means that the system is in an ordered state. In this case, less information is required to describe the system state, and system information entropy is relatively small. However, when system complexity is very high, the system state is difficult to predict, which means that the system is in a chaotic state. In this case, more information is required to describe the system, which features relatively large information entropy. Therefore, a quantitative analysis of RMS complexity is based on information entropy theory, whose mathematical description is shown in formula (1) (Shannon 1948).
where X represents a system; E(X ) represents information contained in system X or the information entropy of system X ; p i (i = 1, 2, . . . , n) represents the probability of system X in the i-th symbol; and lb represents logarithmic to the base 2.
Based on the RMS definition, RMS's rapid-response capability is determined by each manufacturing cell's ability to adjust processing capability and processing function; i.e., RMS reconfiguration is driven by adjusting the processing capability and processing function. Therefore, the RMS processing capability and processing function are chosen as indexes to describe system complexity. At end stage of each production period, RMS processing capability and processing function state probability are analyzed to create the RMS complexity information-entropy model. In the analysis of the system processing capability state, machine tool productivity and the buffer zone state are chosen as analysis indexes; in the analysis of the system processing function state, the machine tool processing state is chosen as the analysis index.

Information entropy-based RMS dynamic complexity model
RMS is composed of several manufacturing cells. Each manufacturing cell contains numerous machine tools. Each type of machine tool is classified as the same type of processing function. The processing route establishes the relationship between processing functions. Part quantity is a test of a machine tool's processing capability; i.e., that affects a machine tool's productivity. When RMS is in operation, the machine tool is in specific state; machine tool quantity and machine tool state in the cellt directly determine cell state. Similarly, overall states of cells in the system determine the system state. Therefore, the machine tool state in the cell should be identified.
When the machine tool is in normal operation and contributes to system stability, it is in an operative/non-blocked state. When the machine tool can complete its processing task but its efficiency is declining and compromises system stability, that tool is in an operative/blocked state. When the machine tool fails to complete its production task because of environmental disturbance (such as a new order or equipment failure) and causes production to be put on hold, it is in an inoperative state. Combined fluctuation in machine tool quantity, machine tool state and cell state determines system complexity, whereas machine tool state type determines complexity type. Based on the effect of the machine tool state on system stability, system complexity is classified as positive complexity and negative complexity. This classification is combined with information entropy theory to create the system dynamic complexity quantitative model shown in formulae (2), (3) and (4): where, E X represents RMS dynamic complexity; E p represents RMS positive complexity; and E n represents RMS negative complexity. RMS dynamic complexity is weighted sum of positive complexity and negative complexity. Positive complexity is information entropy for the operative/nonblocked state; negative complexity is the sum of information entropy for the operative/blocked and inoperative states. N represents the quantity of cells in RMS; G i represents types of parts in the i-th cell; and S j represents the processing route for the j-th type of part in the i-th cell. For instance, when a manufacturing cell contains 2 types of parts, processing routes for parts 1 and 2 are {a b c} and {a c d}, respectively; next, there is where a, b, c and drepresent processing functions. Each type of processing function can contain multiple machine tools. |S j | represents processing route length; p 1 ijk represents the probability of the operative/non-blocked state when the j-th type of part in the i-th cell calls for the k-th processing function in processing route S j ; p 2 ijk represents the probability of the operative/blocked state when the j-th type of part in the i-th cell calls for the k-th processing function in processing route S j ; p 3 ijk represents the probability of the inoperative state when the j-th type of part in the i-th cell calls for the k-th processing function in processing route S j .
The machine tool occupies three states: operative/nonblocked, operative/blocked and inoperative. The operative and inoperative states are related to the processing function. The non-blocked and blocked states are related to processing capability. Because processing function and processing capability are independent from each other, operative/inoperative states and unblock/blocked states are also independent from each other. Machine tool state probability is treated as the combined probability of two types of states: processing function and processing capability. In other words, operative/non-blocked states probability is the product of the probability of the machine tool in its operative state and the probability of the machine tool in its non-blocked state; operative/blocked state probability is the product of the probability of the machine tool in its operative state and the probability of the machine tool in its blocked state. To summarize, the solution for the probability of the machine tool's three states is shown in formula (5): where p K ijk represents the probability of the operative state when the j-th type of part in the i-th cell calls for the k-th processing function in processing route S j; p C ijk represents the probability of the non-blocked state when the j-th type of part in the i-th cell calls for the k-th processing function in processing route S j; p D ijk represents the probability of the blocked state when the j-th type of part in the i-th cell calls for the k-th processing function in processing route S j; and p B ijk represents the probability of the inoperative state when the j-th type of part in the i-th cell calls for the k-th processing function in processing route S j.

Analysis of processing function probability
Each machine tool may provide multiple processing functions; i.e., the relationship between machine tool and processing function is 1:1 or 1:n. The same machine tool may provide different processing functions during the different processing stages of a part. Therefore, machine tools with similar processing functions in a cell are classified as the same type of processing function. The relationship between machine tools and processing functions is shown in formula (6).
where N i represents the total number of machine tools in the i-th manufacturing cell in RMS; n represents the types of processing functions in the manufacturing cell; and D i j represents the quantity of machine tools with the j-th type of processing function in the i-th manufacturing cell. The machine tool is in one of the four operation states: idle state, processing state, adjustment state (maintenance state) and failure state. This section only focuses on whether or not a cell satisfies the processing-function requirementi.e., whether the corresponding machine tool provides this function. Because a machine tool in the adjustment state or failure state fails to satisfy the operation condition, the two are classified as a single type: the failure state. Therefore, the machine tool state in this paper includes the processing state, the failure state and the idle state. During operation, when the machine tool is in the failure state, that means all of this machine tool's processing functions are in the inoperative state. In other words, when all of the machine tools providing a specific type of processing function are in the failure state, that processing function is inoperative. When the machine tool is in the operation state, it only provides a 0 1 h

Input transition rate
Output transition rate single processing function and the other processing functions of the machine tool are all in the inoperative state. Therefore, for parts, even when machine tools with the corresponding processing function in the processing route are in the operative state, if the provided processing function is not what is required by this part, then the processing function required by this part is inoperative. To summarize, the probability of whether or not the processing function is operative is calculated via formulae (7) and (8): where p G ijkm represents failure probability of the m-th machine tool with the k-th type of processing function for the j-th type of part in the i-th cell; pJ ijkm represents the probability of the m-th machine tool with the k type of processing function for the j-th type of part in the i-th cell that is in the operative state but fails to provide the processing function required by this type of part; and M k represents the quantity of machine tools with the k-th processing function. Both P G ijk and PJ ijk are represented as state time frequency. Machine tool state duration in a specific period is collected, and the frequency of each state is calculated as ratio of state duration to overall time.
To summarize, the probability analysis for the processing function is completed, which includes the probability of the processing function in the operative state and the probability of the processing function in the inoperative state.

Probability analysis of processing capability
According to the analysis in "Analysis of processing function probability" section, when the processing function is in the operative state, the non-blocked and blocked states of the relevant machine tool (currently, the machine tool for part processing is determined; i.e., the relationship between the processing function and the part is limited to 1:1) should be further analyzed to obtain the processing capability state. When the buffer zone before the machine tool has stock and the buffer zone after the machine tool has vacancy, the machine tool is in the non-blocked state; when the buffer zone before the machine tool has no stock or the buffer zone after the machine tool has no vacancy, the machine tool is in the blocked state. The buffer zone state represents the machine tool state; therefore, states of buffer zones before and after the machine tool are analyzed to calculate the probability of processing capability, as shown in formulae (9) and (10): where p Q ijk represents the probability that the buffer zone before the machine tool has no stock; and p H ijk represents the probability that the buffer zone after the machine tool has no vacancy.
Buffer zone probability is analyzed further. Assume that the buffer zone capacity is h. Then, from empty to full stock, there are h+1 states. When the machine tool before the buffer zone completes the processing of a part, the buffer zone state moves backward one cell; i.e., the previous machine tool's productivity is the buffer zone's input transition rate. When the machine tool after the buffer zone starts to process a part, the buffer zone state moves forward one cell; i.e., the subsequent machine tool's productivity is the buffer zone's output transition rate. The buffer zone state transition process is shown in Fig. 3.
Based on the above analysis, the state transition equation is constructed as shown in formula (11): where p l ijk represents the probability of the buffer zone in the l-th state; when l = 0, p 0 ijk = p Q ijk ; i.e., the probability of a buffer zone with no stock; when l = h, p h ijk = p H ijk ; i.e., the probability of a buffer zone with no vacancy; w R ijk represents the productivity of the machine tool before the buffer zone; and w C ijk represents the productivity of the machine tool after the buffer zone; i.e., the productivity of the k-th processing function for the j-th part in the i-th cell.
The buffer zone state transition equation set (11) is solved to obtain the probability of the buffer zone state, as shown in formula (12): When the buffer zone before the machine tool has no stock, there is l = 0 in formula (12); the probability of a buffer zone before the machine tool with no stock is shown in formula (13).
When the buffer zone after the machine tool has no vacancy, there is l = h in formula (12); the probability of a buffer zone after a machine tool with no vacancy is shown in formula (14).
To summarize, the probability analysis of processing capability is complete, and both non-blocked and blocked states are calculated and analyzed.

Cusp catastrophe model-based RMS reconfiguration point identification
Various complexity components control system operation and determine system stability. When positive complexity of a maintenance system's stability is predominant, the system has relatively high stability and productivity. When the negative complexity component that leads to system instability is predominant, system stability significantly declines, which could result in system state catastrophe and trigger system reconfiguration. Analysis shows that during operation, RMS will experience state catastrophe because of internal and external factors such as new orders and machine tool failure. Catastrophe theory, which was proposed by French mathematician THOM in 1972, is a universal method used to investigate transition, discontinuity and catastrophe. In catastrophe theory, the system potential function is the study object, which comprises the state variable and external control parameter; the critical point of the system balance state is calculated from the potential function (Barunik and Vosvrda 2009;Chow et al. 2012;Dou and Ghose 2006;Hu and Xia 2015;Sethi and King 1998;Thom et al. 1975;Saunders 1980). Therefore, in this paper, RMS dynamic complexity is analyzed and combined with catastrophe theory to calculate and analyze the reconfiguration point. In research by THOM, when there are no more than four control variables, no more than two types of states, and a maximum of seven types of basic catastrophe models. Four common types of catastrophe models are the folded catastrophe, the cusp catastrophe, the fork-tailed catastrophe and the butterfly catastrophe. The cusp catastrophe state variation is shown in Fig. 4. The diagram shows a folded catastrophe manifold in which the surface top lobe and bottom lobe are in stable balanced states, where the intermediate lobe is in an unstable balanced state. This diagram shares the same characteristic as the RMS implementation process in Fig. 1, in which each production period contains a stable operation period and an unstable ramp-up time. In an ideal situation, the RMS state varies stably along AB, and the system is in a dynamic balanced state. In actual production, because there are uncontrollable factors such as order change and machine tool failure, the system state varies along CD. After a period of stable operation, RMS state catastrophe occurs; i.e., catastrophe occurs at point E and the system moves from the production period to the reconfiguration period. After the system reconfiguration is completed, the new production period starts; after running-in during the ramp-up time, the new stable operation period begins. This process is called the RMS cycle. System dynamic complexity dynamically reflects this cycle; positive complexity and negative complexity are two control variables that balance each other and maintain the system state. When positive complexity dominates, the system is in stable operation period; when negative complexity dominates, system state catastrophe occurs; when positive complexity and negative complexity are balanced, the system is either in a ramp-up time or at the end stage of production. If the system is in a ramp-up time, positive complexity gradually dominates; if it is at the end stage of production, negative complexity gradually dominates.
Based on the above analysis, as production progresses, the gradual wearing of the machine tool results in system state change. The functional relation formula for the machine tool state and system state is expanded in the Taylor series up to the 5th power, as shown in formula (15). Because the highest power of cusp catastrophe's state function is 4th, the Taylor series up to 5th power is enough to describe state catastrophe of RMS.
where Srepresents the RMS state; M represents the machine tool state; and a i (i = 0, 1, . . ., 5) are a set of data related to RMS state complexity. Derivation is performed on both sides of formula (15) to obtain the system state's potential function, as shown in formula (16): where S represents the potential function of RMS state.
Because the cusp catastrophe model has no cubic coefficient item, formula (16) is transformed. Let M = n − m and m = a 4 /(5a 5 ). The cubic coefficient item in the potential function is eliminated to obtain the standard format of the cusp catastrophe potential function, as shown in formula (17).
where V represents the standard format of cusp catastrophe function; b i are the matrix transformation result of a i (i = 0, 1, . . . , 5). Let n = 4 1 b 4 θ , b 4 > 0 . This is substituted into formula (17) to obtain the new system state potential function V , as shown in formula (18).
where θ is the RMS state variable and u and v are two control variables in RMS. Based on the previous analysis, the value of the control variable is determined by positive complexity and negative complexity. The value of the control variable is calculated from the ratio of the complexity component with a different effect to overall system complexity. v represents the positive complexity control variable; u represents the negative complexity control variable. Because negative complexity has a negative effect on system stability, the value for the negative complexity control variable is negative, i.e., u < 0. b 0 is a constant, which does not change the potential function's catastrophe characteristic. Generally, the origin of the potential function is changed to eliminate constant item b 0 . Therefore, the potential function is expressed as formula (19) (Saunders 1980).
where V represents RMS's standard format potential function; u and v are the control variable of RMS based on information entropy theory. Derivation is performed on both ends of potential function V in formula (19). The derived function is set to 0 to obtain a balanced surface equation, as shown in formula (20).
Formulae (20) and (21) are combined to eliminate θ and obtain the catastrophe characteristic equation, as shown in formula (22).
= 8u 3 + 27v 2 where represents state catastrophe decision variable of RMS. When >0, there is 27v 2 > −8u 3 ; i.e., positive complexity is predominant, RMS is in the dynamic stable state and reconfiguration is not needed. When = 0, there is 27v 2 = −8u 3 . Positive complexity and negative complexity are in a dynamic balanced state; RMS is at a critical point. Slight interference may lead to RMS state catastrophe and the need to reconfigure. When <0, there is 27v 2 < −8u 3 . Negative complexity is predominant. RMS stable state catastrophe occurs, and immediate reconfiguration is needed.

Test verification
An RMS plant in a research institute is used as an example to verify the method proposed in this paper. When there is order change, machine tool failure or processing route change, operations such as adding, rearranging and removing machine tools are required. In this paper, these operations are treated as the processing function and processing capability state change. When the current RMS processing function and processing capability are unable to meet demand, problems such as order delay and system instability arise; RMS reconfiguration is needed to reconfigure the system processing function and processing capability to respond to market fluctuation. The issue of when to implement system reconfiguration has a significant impact on reconfiguration cost, reconfiguration time and production benefit. Therefore, the system state should be analyzed. RMS dynamic complexity variation is analyzed to determine whether system catastrophe occurs and then to decide whether to implement RMS reconfiguration. Therefore, the case study is based on the state data of the machine tool and parts in the RMS plant. Processing function and processing capability variation in an order production process are analyzed. Based on information entropy theory, state probability is converted to system complexity data, which are combined with a cusp catastrophe to analyze the RMS reconfiguration decision. When RMS has a new order or implements a new processing route, this algorithm is reapplied in the system state analysis.
This paper focuses on system-level reconfiguration. In production, different cells are in different states; states of multiple cells comprise the RMS state. Moreover, cross-cell operation in system reconfiguration requires the analysis of multiple cells. Therefore, RMS in this case contains three processing cells. Cell 1 contains four machine tools: M 11 , M 12 , M 3 and M 4 . Cell 2 contains five machine tools: M 3 , M 5 , M 6 , M 7 and M 8 . Cell 3 contains three machine tools: M 1/2 , M 5 and M 8 . Among them, M 11 and M 12 mean that these two machine tools provide processing function 1; machine tool M 1/2 means that this machine tool provides processing functions 1 and 2 and can switch processing function according to requirements. Other machine tools only provide a single processing function, and each machine tool can only process a single part. The processing rule is first come, first serve. If different processing functions of the same machine tool are required, a part coming late automatically enters the waiting queue. There is one buffer zone before and one after each machine tool; the buffer zone capacity is 5. The RMS processing function and processing capability are shown in Fig. 6. The order contains six types of parts. The corresponding processing route is planned according to plant processing function and processing capability. The production task is listed in Table 1. In the table, the part number represents the part category; the processing route represents the required processing function codes to complete the processing of the corresponding part; the quantity represents the processing capability required for each part. For instance, part 03 follows processing route 6, 7, 9. This means that processing this part requires machine tools M 6 , M 7 and M 9 . Based on the processing capability of existing machine tools in the plant, the overall time to process each type of part and the processing time for each step are listed in Table 2. Each line represents the steps and processing time for each part; each column represents the part type processed by each machine tool and the corresponding processing time. Processing task  is assigned according to order information. Parts 01 and 02 are processed in cell 1; parts 03 and 04 are processed in cell 2; and parts 05 and 06 are processed in cell 3. Three processing cells operate simultaneously and the processing in each cell follows the first come, first serve (FCFS) rule. The data in Table 2 are the fixed procedure time defined in the process design. However, failure in production is random, and the production process in this paper is a stochastic process. Based on the processing time distribution in Table 2, the machine tool and part state in the system are collected every 10 working days and probabilities of the processing function and processing capability are calculated. Because there is a large amount of data, machine tool M 1/2 is used as an example to demonstrate the probability calculation, which includes the processing function's operative probability p K ijk and inoperative probability p B ijk , the failure probability p G ijkm , the non-required function processing probability pJ ijkm , the processing capability's non-blocked state probability p C ijk and blocked state probability p D ijk , the probability of a previous buffer zone p Q ijk , the probability of a subsequent buffer zone p H ijk , the machine tool's operative/non-blocked state probability p 1 ijk , the operative/blocked state probability p 2 ijk , and the inoperative probability p 3 ijk , as listed in Table 3. The machine tool state is monitored and the duration of each state is collected every 10 working days. Machine tool state probabilities p G ijkm and pJ ijkm are calculated from the ratio of state duration to overall time. Based on formulae (7) and (8), p K ijk and p B ijk are calculated. The throughput of the machine tool in operation is monitored to obtain the average productivity of the corresponding machine tool. Based on formulae (13) and (14), p Q ijk and p H ijk , are calculated. Next, based on formulae (9) and (10), p C ijk and p D ijk are calculated. Based on data statistics and analysis, the operative/nonblocked state probability p 1 ijk , the operative/blocked state probability p 2 ijk , and the inoperative probability p 3 ijk for each machine tool in each cell are calculated as listed in Tables 4,  5 and 6. RMS complexity and RMS cusp catastrophe reconfiguration decision data are calculated from the above data via the algorithm proposed in this paper, which includes RMS complexity, positive complexity, negative complexity, decision index and system control variables, as listed in Table 7. RMS state variation trend diagram in Fig. 7 is based on data in Table 7. The system state trend diagram in Fig. 7 shows the trend of system state detection index : after a relatively stable production period, the system experiences a state catastrophe, which triggers system reconfiguration. In the first 40 working days, when is above 0 and stays around 0.8, fluctuation is relatively small. This shows that the system is in a stable operation period and reconfiguration is not needed. During the 40th ∼ 50th working days, decreases dramat- ically from 0.77 to −0.427. When is <0, the critical point is crossed and the cusp catastrophe condition is met. This shows that system state catastrophe occurs and immediate RMS reconfiguration is required. RMS monitoring is continued. Figure 5 shows that on the 60th working day, the system state continues the downward trend. This further confirms the necessity of reconfiguration.

Conclusions
In this paper, a dynamic complexity-based RMS reconfiguration point decision method is proposed. First, the effect of processing function and processing capability on RMS is analyzed and the relationship between RMS complexity and RMS state is elaborated. This shows that the moment of an RMS state catastrophe is the best time to implement RMS reconfiguration. Next, based on information entropy theory, a system complexity quantitative model is created to analyze the effect of the system complexity component on the system state. That is then combined with the cusp catastrophe theory to analyze system state variation under the effect of various complexity components. The cusp catastrophe determinant condition helps identify the state catastrophe point and decide the RMS reconfiguration point. The test shows that information entropy-based system complexity analysis provides a quantitative description of system complexity. System complexity variation is analyzed in terms of processing function and processing capability to reveal the nature of the system state. Application of the cusp catastrophe helps identify the system state catastrophe moment. This provides evidence for the decision maker to decide the reconfiguration point, promptly implement system reconfiguration and maintain RMS vitality. Subsequent work focuses on an indepth investigation of the RMS state catastrophe mechanism, which includes an analysis of the production factors that lead to system state catastrophe.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.