Reliability Modeling of Electromechanical System with Meta-Action Chain Methodology

To establish a more flexible and accurate reliability model, the reliability modeling and solving algorithm based on the metaaction chain thought are used in this thesis. Instead of estimating the reliability of the whole system only in the standard operating mode, this dissertation adopts the structure chain and the operating action chain for the system reliability modeling. The failure information and structure information for each component are integrated into the model to overcome the given factors applied in the traditional modeling. In the industrial application, there may be different operating modes for a multicomponent system. The meta-action chain methodology can estimate the system reliability under different operating modes by modeling the components with varieties of failure sensitivities. This approach has been identified by computing some electromechanical system cases. The results indicate that the process could improve the system reliability estimation. It is an effective tool to solve the reliability estimation problem in the system under various operating modes.


Introduction
During recent decades, the electromechanical systems have been widely used in the industrial field.Therefore, the importance of the reliability of electromechanical systems has increased accordingly.On the other hand, the growing complexity of systems has made the job of reliability modeling more challenging.Many qualitative and quantitative techniques have been proposed to analyze the reliability.Traditional methods are effective under the hypothesis that system reliability can be described by macroscopic stochastic models while pure probability models only describe the stochasticity of failures in system level without considering failure mechanism.This paper is proposed to identify the connection between macroscopic system failure probability and microscopic failure cause.System operation information and structure information are integrated into the model to describe the state changes for a given system.In fact, macroscopic system operation is the input of the microscopic failure processes of each component in the system.On the other hand, the microscopic failure processes make up the macroscopic stochastic model.Herein, meta-action structure chain and operating series chain are proposed in this paper to combine these two aspects.
System reliability modeling has been widely used in researches in recent years [1][2][3][4][5][6][7].In the previous research, the failure modes and effects analysis (FMEA) has been often used to determine the criticality and frequency of failure.Yang et al. [8] presented cubic transformed functions of FMEA.Kim et al. [9] summarized the limitations of the FMEA methods.Besides, fault tree analysis (FTA) is always utilized to calculate the probability of failures [10][11][12][13][14][15][16].Majdara and Wakabayashi [17] put forward a component-based modeling of systems for the automated fault tree generation.Lindhe et al. [18] optimized the FTA with the approximate calculations.
In addition, the fuzzy set theory and Bayesian approach have been widely used especially in the cases that the systems have small sample [19][20][21].Marquez et al. [19] utilized Bayesian to propose improved reliability.Bae et al. [20] utilized the Bayesian approach to model the two-phase degradation.Sun et al. [21] adopted the fuzzy fault tree analysis to make fault diagnosis for the servo system.These static reliability models have shown advantages in reliability analysis for the electromechanical systems in the industrial fields.However, these methods cannot model the multistate nature of the systems.If system operation does not satisfy the rated conditions, these stochastic models would not be effective.In recent years, the research on dynamic reliability has made a considerable contribution to the system reliability modeling in terms of overcoming the mentioned shortcoming.Wang et al. [22] defined the system status as reliability block diagrams to model different operating regimes of a parallel system.Chiacchio et al. [23] proposed a stochastic automaton model and developed a Simulink environment to integrate the features of a dynamic system.This model considers the multistate operational and failure nature of a system.Compared with static methods, dynamic reliability modeling is much more suitable for industrial application, having a vast research potential.
In this paper, system operation is expressed as state transitions versus time dimension and the state of the system is expressed with all the status of meta-actions.State transitions of the specific meta-action are utilized to combine the information of failure modes.Traditional reliability modeling estimates the system reliability with system operating time.In fact, the operating time of each component in the system has an implicit function to the total operating time of the whole system and some failures are even not sensitive to the operating time (e.g., sensitive to rotation numbers and total travel).System structure information, operating information, and failure information should be integrated into one reliability model to make reliability estimation more practical.
Thus, the impact of the paper can be summarized as follows: (1) The meta-action chain method is presented, which can be used to integrate the structure information and operating information to estimate the reliability for the electromechanical system.The definition and calculation of the formula are given to support the modeling.
(2) The reliability solving algorithm is designed to solve the Markov transfer matrices and the reliability curve.The discretization approach is used to solve the discontinuous time interval.
(3) The proposed methods in this thesis establish the relationship between the microscopic failure modes and the macroscopic reliability estimation.More information from the accelerated loading experiments should be taken into consideration in the modeling, making the reliability curve more reasonable and precise.
The rest of the thesis is organized as follows: Section 2 introduces some basic theories used in the modeling; the meta-action chain methodology (MCM) is presented in Section 3; Section 4 introduces the reliability estimation for an electromechanical system; In Section 5, an example is presented to illustrate the proposed approach and some comparisons with the existing approaches are discussed; Conclusion is made in the final section.

The Meta-Action Theory and Markov Model
The Pedigree-Function-Movement-Action (PFMA) structural decomposition shown in Figure 1 is proposed to study the operation of the electromechanical system [24][25][26][27].Herein, the action is defined as the geometric position variation of a component in the system.The definition of metaaction is a structure-independent action to realize the specific function.Meta-action cannot be further decomposed.The motions of the electromechanical system may be complicated while the meta-actions must be simple.Components in the electromechanical only have two types of meta-action.One is rotation, such as rotations of motor, screw, gear, and shaft.The other is translation, like the translation of nut and piston.This decomposition can establish the motion model for the electromechanical systems.
The meta-action theory has been utilized in the study of precision control and assembly reliability analysis.This method focuses on the minimum fault analytical unit to precisely control the reliability of the electromechanical system as shown in Figure 2 and Figure 3. Furthermore, the meta-action decomposition has established the relationship between the faults and component actions.Compared with traditional modeling, not only the structural information of the system but also the functions of the system are considered in metaaction modeling.Such approach shows great advantage in the reliability analysis and assembly precision analysis for complex electromechanical systems.

The Meta-Action Chain Methodology
As meta-action focuses on the basic moving unit of the system, this method is used to analyze the system precisely.However, the reliability and precision are controlled by optimizing the sensitive parameters of the single essential meta-action in the former research.The failure mechanism and work mode effect are not included in the traditional modeling.On the one hand, the meta-action chain methodology (MCM) is put forward in this chapter to solve these problems.On the other hand, the traditional Markov chain methods in the reliability research only focus on the failure probability in the  given states of the system.The relationship between the states selection and system operation has not been established.In the MCM model proposed by this study, the system operating states are represented by the meta-actions matrix and the mapping matrix.The meta-actions in each state have their own computing chains to ensure the consistency with the failure mechanism.

Definition of Meta-Action Chain.
The definition of a meta-action chain is a series of ordered meta-actions that is driven by the same power input to describe the states of specific function in given condition and operating mode.Two kinds of information are considered in the meta-action chain.One is structure information and the other is operating state information.Meta-action chain includes two parts, the structure chain and operating series chain.The operating series chain can be mapped into the Markov chain for different working modes.
In structure chain, there are  power inputs (usually motors) denoted as  and meta-actions are denoted as   .The structure chain is the universal set of meta-actions and it contains all the elements used in operating chain.The outputs of each line is the actuator denoted as .The number of outputs  for each input power   is l  .For example, the number of outputs for the first input power is  1 .As the transmission chain may have two or more outputs, the metaaction chain would have branches, shown as   1 ,   1 , and   1 in Figure 4.This structure decomposition is able to express the operation for electronic systems clearly.
To realize specific operation cycle, electromechanical systems work orderly action by action as shown in Figure 5. From function 1 to function , system states change from  11 to   .Here,  11 means the first state in function 1 and  12 means the second state in function 1.Then,   means the   th state in function .  refers to the number of states in function .In state  11 , there are  11 meta-actions operating together.During the operation cycle, all the meta-actions compose a meta-action operating series chain to describe the states that the system changes into. 11,1 is another expression of meta-action and refers to the first meta-action in the first state of function 1.For example,  11,1 =    11 means that the mentioned meta-action is    11 in the structure.

The Reliability Calculation of MCM.
To determine the reliability of each state for the system, the operating metaaction chain is expressed by the Boolean variable matrices.
Herein,  is the number of the structure meta-actions in the th power input.Subsequently, the state of the system can be expressed as a one-dimensional matrix.
1, the meta-action is operating 0, the meta-action is not operating. .The electromechanical system changes from one state to another in the given operating series.If the state changes into the wrong state, the system is considered to have a fault.The reliability of the system is the probability that all the given state series of the system can be reached successfully.Figure 6 provides an example of four states to show the state transition of the system.The system state changes from  11 to  22 .If the system changes into the wrong state, a failure will occur.   means the probability that system state changes from state  to state  in  steps.In different running modes and different failure mechanisms, the reliability calculations of meta-actions are different.
Take one column as an example in Figure 7; it represents the states' change process of the specific meta-action.
In different operating modes, the state changes of the meta-actions are different.The operating mode information can be used in the modeling just as the failure mode information shown in Figure 8.
The common failure modes include degradation, switch fault, and alternating stress.The failure mechanism can be reflected on the meta-actions as clearly shown in the figure .If the meta-action is designed to always run during the operation, the column is full of 1 as the normal operation column in the figure.When the degradation failures happen, the continuous state changes from 1 to 0. The second type of failure is the switch fault.The reliability of such kind of metaaction depends on the probability of successful action.The failure occurs when the state is not able to change.The third kind of failure is the alternating stress failure.The state change frequency of such type of meta-action is high, which may cause additional degradation to the function units.Thus, the reliability calculations for the three kinds of meta-actions are different.
In the degradation process, the reliability of meta-action depends on the operating time.The operating time is expressed as the discrete variable.In the given time interval , the meta-action state remains unchanged. is the sampling interval and is also the resolution of the calculation.The operating time is calculated by the number of time intervals   .The switch fault is sensitive to the number of the switch times.The frequency of action is   and the failure rate  is considered as a constant.The alternating stress may cause additional degradation to the function unit.The correction coefficient   is used to modify the formula.
PLP (Power Law Process) is utilized to model the failure density function in the degradation process.
(3)  is the size parameter,  is the shape parameter,  de () is the degradation failure probability, and  is the time between failures.
The failure rate of switch failure  is considered as constant during the operation as the statistic result shows that it is low and changes little.
LND (Logarithmic Normal Distribution) is used to model the probability density function in the alternating stress process.
The correction coefficient   depends on the standard cycle times  0 , the index of stress state , and the cycle times Actuator O 1l I ,      .Usually,  0 = 10 7 and  = 9.  is logarithm to the standard cycle time  0 and  is variation.
On the other hand, the reliability of the whole system at time  is expressed as the product of reliability of all the meta-actions.Herein, the reliability model of the system is considered as a series model.Different failure modes happen independently.
Reliability transfer matrix  of each meta-action at time  is calculated in the Markov chain form.

𝑝 𝐴 𝑗𝑖
V refers to the probability that meta-action   translates from state  to state V.
A flow chart is shown in Figure 9 to illustrate the calculation algorithm.

Reliability Estimation for the Electromechanical System with MCM
The subsystem of a CNC machine tool is used to illustrate the method proposed by this thesis.The structure is presented in Figure 10.The subsystem is used to realize the indexing rotation function of the rotary table.And the meta-action structure chain is established in Figure 11.
A typical work circle of the system concludes 2 functions, including the exchanging worktable and indexing table.The operating series chains of the functions are shown in Figure 12.
Figure 13 is a timing sequence diagram of the system.This figure reflects the time variation of the rotation angle of the worktable in a given working condition. is the rotation angle of the worktable. is the load mass of the system. is the hydraulic pressure. is the rotational speed.
The operation time of this cycle is 1 hour (3600 seconds) and the worktable needs to be exchanged before each cycle in 5 minutes (300 seconds).Thus, the sampling interval  in 1 second is selected and a 3900 × 12 operating states matrix .  = 3900 is created as shown in Figure 14.
The failure modes of each meta-action can be determined by the operating states matrix and the previous failure data.Table 1 presents the failure modes of each meta-action.Failure on meta-action

Numerical Example
The original failure data of the given system are obtained from experiments, customer survey, and customer service data.The parameters of reliability models are determined by the fitting probability density function.They are listed in Table 1.
Herein, time  is calculated by the hour.

The Reliability Estimation for the Standard Operating
Mode.At the beginning of the curve, the reliability of the system decreases smoothly as the main failure mode is degradation.After 3000 hours, the number of alternating stress cycles reaches the critical value.Thus, the reliability of the system decreases sharply for the reason that the main failure mode changes into the alternating stress.This curve is smooth and the turning point is clear.The solving accuracy is acceptable.
The traditional reliability estimation models the system reliability curve with Weibull distribution as shown in Figure 15.The test result under the significance level of 5% is ksstat = 0.0740 < cv = 0.1112, indicating that the hypothesis is not rejected.The shape parameter  is 1.8035 and the scale parameter  is 8125.8331.
Compared with the Weibull distribution, the MCM curve contains more failure modes information.In fact, in the given working condition, alternating stress failures occur more frequently than degradation failures after 3000 hours.On the other hand, the Weibull distribution curve is still smooth after the occurrence of the critical value.The traditional modeling methods cannot describe such kind of trend.If the MCM method is set as the benchmark, a cumulative relative error of Weibull distribution is shown as follows: This value shows the difference between the two methods directly.Generally, the two methods can reflect the overall trend of the failure rate of the given system.According to some specific details, the maximum error of these two methods is observed in  max = 0.0581.And the proposed MCM method result shows more realistic behavior than the traditional Weibull distribution.

Reliability Estimation for Different Operating Modes.
The system reliability in different operating modes is estimated to show the failure effects from each operating parameter.In other operating modes, the timing sequence diagrams of the system are different from the standard operating mode.Hence, the operating states matrices are also different.Three typical operating modes are utilized to illustrate the MCM model as Figure 16 presents.
In operating mode 2, the degradation time of the system decreases into half of the standard and the number of alternating stress cycle numbers decreases into quarter.In operating mode 3, the number of the alternating stresses decreases into half while in operating mode 4, the alternating stress cycle number increases to 2.5 times.

Reliability Estimation for Other Electromechanical Systems.
Figure 18 is the structure of a lifting part of tray automatic exchanging device.
As the function of this system is not complex (only exchanging tray), its failure rate is sensitive to the exchange frequency.According to the former failure data, only three  meta-action failures were observed during the service life as shown in Table 5.
Both the MCM and Weibull distribution are utilized to estimate the system under the y 8 min −1 , 10 min −1 , and 12 min −1 exchange frequency.The results are presented in Figure 19.
Cumulative relative errors and maximum errors are presented in Table 6.
As the failure mechanisms of this example are simpler than the former one, there exists the difference between the two methods, as shown in Table 6.With the increment of the exchange frequency, the sensitivity of system failure changes from time dimension to the combination of time and exchange numbers.At a low exchange frequency, there is little difference between the MCM and Weibull distribution for the reason that system failure is sensitive to the time dimension.Both of the two methods can reflect such kind of trend.However, as the exchange frequency increases, the traditional Weibull distribution needs to be refreshed frequently to adjust the changeable failure mechanism.For instance, in the frequency of 10 min −1 , at 2000 h, the shape parameter  is 1.6133 and the scale parameter  is 8625.2314,while at the 3000 h shape parameter,  is 1.8221 and the scale parameter  is 7636.5287.On the other hand, MCM is much more stable as the independent input data of each meta-action are more than the data of the system operation.Besides, system failure data and meta-action data are combined as the input of MCM instead of the only system failure.For these reasons, MCM can estimate the system dynamic reliability precisely.

Conclusions
This paper presents the meta-action chain method to identify the connection between microscopic failure mechanism and microscopic system reliability calculation.The main scientific novelties are the modeling technique of the meta-action structure chain and operating series chain, reliability calculation modeling with Markov transferring matrices, and solving algorithm for industrial application.The meta-action chain modeling technique defines system operation according to the states of each meta-action instead of the traditional reliability block diagram.The transitions of these states cause failures directly and the failure mechanism could be researched independently as failure sensitivity is much clearer.On the other hand, this method remains consistent with the traditional reliability modeling methods.Existing researches could be further developed and expanded by using microscopic failure data to obtain more precise and practical results.
Reliability calculation modeling uses the discretization approach to solve the discontinuous time interval.Solving

Figure 3 :
Figure 3: Meta-action of gear shaft rotation assembly unit.

Figure 7 :
Figure 7: States change of specific meta-action.

Figure 15 :
Figure 15: Reliability estimation for the NC rotary table with MCM and Weibull distribution.
The universal set of the metaactions    is

Table 1 :
Failure modes of each meta-action.

Table 2 ,
Table 3, and Table 4 present distribution function parameters of each meta-action.

Table 2 :
Degradation failure distribution function parameters of each meta-action.

Table 3 :
Switch failure rate of each meta-action.

Table 4 :
Alternating stress failure distribution function parameters of each meta-action.

Table 5 :
Failure modes of the lifting part.

Table 6 :
Cumulative relative errors and maximum errors.