# Synchronization Reliability Evaluation Method for Mechanisms with Different Time Distribution

Authors: Huan Pang; Tianxiang Yu; Ning Wang
Journal: Mathematical Problems in Engineering
 (2018)

---

## Abstract

Synchronization reliability problems are common in the mechanical field; motion asynchrony will affect the performance or even lead to failure. As action time of mechanisms does not share the identical distribution, the existing synchronization reliability evaluation method has great limitations. Aimed at the problem, a novel synchronization reliability evaluation method is proposed. Starting from two mechanisms, the synchronization reliability ofN mechanisms can be obtained with recursive algorithm, where synchronization reliability of the mechanisms is expressed as a multiple integral, through dividing the integral domain into several independent domains, the uncertain integral limits are translated into certain integral limits, and then multiple integral can be solved. The numerical examples show that errors between the proposed methods and the MC simulation are very small, which proved that the methods are correct. Finally, synchronization reliability of folding wings is evaluated. The proposed methods make up for the limitations of the existing method and have good versatility.


---

## Body

## 1. Introduction

In the bookSynchronization of Mechanical Systems, Professor Alejandro Rodriguez Angeles of the University of Twente in Holland points out that “the development of today’s technology and the quality requirements of production processes have led to complex integrated production systems. In this case, the use of multiple combination systems is considered the most effective solution.” The composition of the combined system can be a number of identical or different units, which work together to perform tasks, and the “together” here refers to the synchronization or coordination of the work of each unit [1].In engineering fields, especially in aerospace field, it is common that several mechanisms that have the same structure work together to accomplish a function. For example, four folding wings deployed together after missile was launched to make the missile in control (Figure1(a)); several lock mechanisms work together to lock or unlock the landing gear doors (Figures 1(b) and 1(c)). And such examples widely exist.Figure 1
Illustration of the systems that have synchronization requirement ((a) missile with four folding wings; (b) landing gear door system with three locks; (c) cabin door system with eight locks).
(a) (b) (c)Motion asynchrony will affect the performance or even lead to failure [2–8]. For example, too large time difference of the four folding wings will lead to the missile being out of control or even lead to the missile crashing. Therefore, synchronicity should be clearly put forward in the design process.The theme of synchronization has attracted a lot of attention over the past few decades, which focuses on one of the basic issues in science, that is, how to implement the synchronization of units in multiple systems. In his doctoral thesis (1955), Blekhman of the Institute of Theoretical Engineering of the Academy of Sciences of Russian Academy of Sciences gives a clear definition of synchronization: the time consistency of two or more units, which can be expressed in certain specific relationships of the system. In addition, he expounded the synchronization problem in dynamics system in 1964 [9]. Later, Blekhman (1988) pointed out that synchronization is seen in such varied examples as weakly coupled pendulums, organ tubes, celestial bodies, electrical, electromagnetic, and quantum generators, vibrations in process equipment, turbine blades, populations of cells, swarms of fireflies, flocks of birds, schools of fish, and applauding or marching people [10]. Alejandro Rodriguez Angeles (2002) pointed out that there are several basic problems in system synchronization. First, the function of describing the desired synchronization targets must be formulated. For this reason, the type of system components and the interesting variables must be taken into account. For mechanism synchronization, the function can be defined as the norm of the difference between the variables of interest, such as the position and speed; secondly, the design of the coupling or the interconnection and feedback controllers must ensure synchronization behavior; finally, the condition that guarantees the synchronization target must be determined [1]. Further, Nagaev (2003) studied synchronization of inertial vibration exciters and dynamical objects of the general type, in which weak interaction of anisochronous and isochronous objects is considered [11].As to the synchronization performance analysis and reliability evaluation, early studies are focused on the fire-industrial products. The minimum movement synchronization precision is obtained by kinematic simulation [4]. The reliability of synchronous function time of detonator during long term storage is obtained by the Monte Carlo method [12]. From the above references, it is clear that acting time of fire-industrial products is primarily affected by the random factors such as material of the electrodes, geometry shapes, and surface conditions. Usually, acting time of them can be regarded as identically distributed (both distribution types and distribution parameters are identical). With the identical distribution assumption, Wang ZL proposed the synchronization reliability evaluation method for two products [5].Synchronization reliability evaluation method for three or more products has been studied by some scholars [13, 14], in which the action time of each product is also assumed to share the same normal distribution N(u,σ2). After sorting the variables according to the ascending order, x(1),x(2),…,x(n), the extreme difference is expressed as X=x(n)-x(1); according to the probability density function (PDF) of the extreme value distribute, the synchronization reliability of the n products can be expressed as R=Pxn-x1⩽ε. But this method is very difficult to solve, and we need to calculate the extreme value distribution of variables.Synchronization reliability problems for mechanical system are similar to these of fire-industrial products. Ideally, several mechanisms with the same structure should have the same performance, but the presence of many random factors makes the acting time different from each other. The random factors are the driving force, geometric parameters, physical parameters (e.g., coefficients of friction and damping), and the working load. Under normal circumstances, the driving force, geometrical parameters, and physical parameters of each mechanism can be regarded as sharing the same distribution, but the installation position and the attitude are usually different, which will lead to different working load distribution types or different parameters; as shown in Figure1(a), difference of the flight attitude makes aerodynamic load of the four wings different; as shown in Figures 1(b) and 1(c), installation position of the lock mechanism results in significant differences in their working load. Accordingly, the assumption that acting time distribution types or distribution parameters are identical no longer holds in sync reliability analysis of mechanisms, which means distribution type or parameters are often different, and their extreme value distribution cannot be obtained.Therefore, the existing synchronization reliability evaluation method cannot solve the synchronization reliability problems of mechanisms in engineering area. In order to solve the above engineering problems, this paper is going to research novel synchronization reliability evaluation method and the synchronization reliability of the folding wings will be evaluated using the proposed method.


## 2. Novel Methods for Synchronization Reliability Evaluation

For a system withn components, action time of each component is a random variable; it can be expressed by xi,i=1,…,n. As each component shares the same structure style and working principle, the distribution type of the action time of the components should be identical. However, the distribution parameters of each action time are often different because of the random factors such as manufacturing errors and the external load. Therefore, in the following text, probability density function of action time of ith component is expressed as fi(xi),i=1,…,n, where action time of all the components sharing the same distribution type and parameters is a particular case.Synchronization requires that the difference of the action time should be less thanε, namely, xi-xj<ε. Therefore, synchronization reliability can be expressed as R=P{xi-xj<ε}.Synchronization ofn-1 components is the premise of synchronization of n components. When the action time of the first component is determined, in order to satisfy the requirement of the synchronization, the action time of the second component should be among [x1-ε, x1+ε].In the same way, while action times of then-1 components are determined and they meet the requirement of the synchronization, sort the action time from small to large, as shown in Figure 2; the minimum action time is x(1) and the maximum action time is x(n-1), where, x(n-1)-x(1)<ε.Figure 2
Reliable domain analysis forxn.Then, it is obvious that when action time ofxn is among [xn-1-ε,x(1)+ε], the n components will satisfy the synchronization requirement. Therefore, synchronization reliability of n components can be expressed as(1)Rsyn=∫0∞f1x1∫x1-εx1+εf2x2∫max⁡x1,x2-εmin⁡x1,x2+εf3x3…∫max⁡x1,x2,…,xn-1-εmin⁡x1,x2,…,xn-1+εfnxndxn…dx3dx2dx1Equation (1) shows that synchronization reliability of n components can be expressed as aN-ple integral; it can be seen that the limitation of the most outer integral is 0 to ∞, because the action time should be positive number. And the limitation of the inner integral is determined by the order of the outside variables; the lower limitation is the maximum value minus ε and the upper limitation is the minimum value plus ε.However, (1) cannot be solved when the order of the outside variables is uncertain. Therefore, the integral space is divided into several independent domains to make the integral limitation determinately. Thus, the uncertain problem is transferred to certain problem. The dividing process is as follows.(  1) Integral Domain Dividing for Two Components.As to the synchronization reliability problem of two components, it is easy to see that when x1 is determined, integral upper limit and lower limit of x2 are [x1-ε,x1+ε], as shown in Figure 3.Figure 3
Reliable domain analysis forx2.Therefore, synchronization reliability of two components can be expressed as(2)Rsyn=∫0∞f1x1∫x1-εx1+εf2x2dx2dx1(  2) Integral Domain Dividing for Three Components.As to the synchronization reliability analysis for three components, since x1 and x2 are not determined, the integral domain of x3 cannot be determined and (1) cannot be solved. Therefore, the integral interval is divided into two sections: x1>x2 and x1<x2. When x1>x2, x3 should be among [x1-ε,x2+ε], while x1<x2,x3 should be among [x2-ε,x1+ε]. The sketch map is shown in Figure 4 and the independent domains are listed in Table 1.Table 1
Independent domains for three-component synchronization analysis.
Domain ofx2 Domain ofx3 x 1-ε<x2<x1 x 1-ε<x3<x2+ε x 1<x2<x1+ε x 2-ε<x3<x1+εFigure 4
Reliable domain analysis forx3.Therefore, synchronization reliability of the three components can be expressed as(3)Rsyn=∫0∞f1x1∫x1-εx1f2x2∫x1-εx2+εf3x3dx3dx2dx1+∫0∞f1x1∫x1x1+εf2x2∫x2-εx1+εf3x3dx3dx2dx1(  3) Integral Domain Dividing for Four Components. As to the synchronization reliability analysis for four components, since x1, x2, and x3 are not determined, the integral domain of x4 cannot be determined and (1) cannot be solved. Therefore, the integral interval is divided into six sections as shown in Figure 5 and the independent domains are listed in Table 2.Table 2
Independent domains for four-component synchronization analysis.
Independent Domains Independent Domains x 2 x 3 x 4 x 2 x 3 x 4 x 1-ε<x2<x1 x 1-ε<x3<x2 x 1-ε<x4<x3+ε x 1<x2<x1+ε x 2-ε<x3<x1 x 2-ε<x4<x3+ε x 2<x3<x1 x 1-ε<x4<x2+ε x 1<x3<x2 x 2-ε<x4<x1+ε x 1<x3<x2+ε x 3-ε<x4<x2+ε x 2<x3<x1+ε x 3-ε<x4<x1+εFigure 5
Reliable domain analysis forx4.It can be seen that the integral domain can be divided into six independent intervals when analyzing synchronization reliability for four components. So the synchronization reliability is the sum of the reliability of each independent domain as follows:(4)Rsyn=∫0∞f1x1∫x1-εx1f2x2∫x1-εx2f3x3∫x1-εx3+εf4x4dx4dx3dx2dx1+∫0∞f1x1∫x1-εx1f2x2∫x2x1f3x3∫x1-εx2+εf4x4dx4dx3dx2dx1+∫0∞f1x1∫x1-εx1f2x2∫x1x2+εf3x3∫x3-εx2+εf4x4dx4dx3dx2dx1+∫0∞f1x1∫x1x1+εf2x2∫x2-εx1f3x3∫x2-εx3+εf4x4dx4dx3dx2dx1+∫0∞f1x1∫x1x1+εf2x2∫x1x2f3x3∫x2-εx1+εf4x4dx4dx3dx2dx1+∫0∞f1x1∫x1x1+εf2x2∫x2x1+εf3x3∫x3-εx1+εf4x4dx4dx3dx2dx1(  4) Number Analysis of the Independent Integral Domains.According to the discussion of the integral domain, the number of independent domains for two components, three components, and four components is one, two, and six, respectively.As ton components, since integral limit of nth variable is determined by the order of then-1 components, the independent integral domain for n components isN=(n-1)!, and the relationship between number of the components and that of the independent integral domain is listed in Table 3.Table 3
Relationship between number of the components and that of the independent integral domains.
Total number of the components 2 3 4 5 … n Number of the independent domains 1 2 6 24 … (n-1)!Since the (n-1)! integral domains are independent, the synchronization reliability of n components is the sum of reliability in (n-1)! domains, and the reliability in each domain is an-ple integral, namely,(5)Rsyn=∑i=1n-1!Rn,iFor most of the distribution types, although the probability density functions are elementary functions, the integral cannot be expressed by elementary function. Therefore, numerical algorithms are often used to solveRn,i, such as numerical integral method and Monte Carlo simulation method [15–17].(  5) Flowchart of the Method.The synchronization reliability of N mechanisms can be obtained by recursive algorithm. That means beginning fromn=2 to obtain the synchronization reliability, and then adding a mechanism, untiln=N; the detailed steps are as follows.Step 1.
Define the number of the mechanismsN (N>2) and obtain the probability density function (PDF) fi(xi) of action time of each mechanism.Step 2.
Present the synchronization reliability of two mechanisms by (2).Step 3.
Keep the integral limits ofx1 to xn-1 unchanged, divide the integral domain of xn into n independent domains, and obtain the corresponding reliable domain of xn+1, as shown below:(6)Dn,1:x1-ε,x1andDn+1,1:x1-ε,xn+εDn,2:x1,x2Dn+1,2:xn-ε,x1+ε;……Dn,n:xn-1-ε,xn-1+εDn+1,n:xn-ε,x1+εStep 4.
According to the independent domains in Step3, replace ∫Dn-1,ifn(xn) in (2) by ∑i=1n∫Dn,ifn(xn)∫Dn+1,ifn+1(xn+1)dxn+1dxn; synchronization reliability of n mechanisms is expressed.Step 5.
Add a mechanism to the presentn, namely,n=n+1. Ifn<N, execute Steps 3 and 4, or solve the expression in Step 4, which is the synchronization reliability of N mechanisms Rsyn,N°


## 3. Validation of the Proposed Method

To validate the proposed methods, the following numerical examples are used. As shown in Table4, there are four components, whose action times are x1,x2,x3, and x4, respectively, and the distribution types and parameters are given in Table 4.Table 4
Distribution type and parameters of the action time of four components.
No. of the wing Deploy time(ms) Distribution type Mean value Standard deviation i=1 x 1 Normal 10 1 i=2 x 2 Normal 11 0.9 i=3 x 3 Normal 12 1.2 i=4 x 4 Normal 13 0.95According to the performance index, the difference in action time should be less thanε. In order to demonstrate the applicability of the proposed method, synchronization reliability of two, three, and four components is estimated using the proposed methods and the verification method, and the comparison results are listed in Table 5.Table 5
Results comparison of the reliability.
Number of components Time difference Synchronization reliabilityRsyn Proposed method MC method Error/% N=2(i=1,2) ε=2 0.7585 0.7587 -0.02636 ε=3 0.9300 0.930025 -0.00269 ε=4 0.9870 0.987063 -0.00638  N=3(i=1,2,3) ε=2 0.3829 0.38269 0.054875 ε=3 0.6831 0.683033 0.009809 ε=4 0.8831 0.883581 -0.05444  N=4(i=1,2,3,4) ε=2 0.1122603 0.112018 0.216305 ε=3 0.3710127 0.374236 -0.8613 ε=4 0.6897235 0.690952 -0.1778According to the examples, it can be seen that the error between the reliability obtained from MC method and that of the proposed method is very small, so the proposed methods are proved available.


## 4. Synchronization Reliability Evaluation for Folding Wings

### 4.1. Failure Mechanism Analysis

As shown in Figure1(a), a missile contains four folding wings, which are powered by four actuating cylinders, respectively. By screw pair, rectilinear motion is transferred into rotatory motion to achieve the wings movement.Considering four folding wings have the same structure, the movement out of sync is mainly caused by three aspects:(1) inherent factors: motion friction coefficient is different because of manufacture errors; (2) power factor: maximum pressure and maximum flow are different due to the dosage of the dynamite and the error of damping cavity; (3) external load factors: both of the magnitude and the direction of the aerodynamic load acting on each wing are quite different. Influence factors analysis for synchronization failure is shown in Figure 6.Figure 6
Influence factors analysis for synchronization failure.Performance index of the folding wing requires that the time difference should be less than 80 ms. Therefore, the synchronization reliability can be expressed as(7)R=Pti-tj<80msi,j=1,2,3,4


### 4.2. Dynamic Model and Agent Model

The dynamic model of the folding wing is built using LMS Virtual Lab, and the driving model is built using AMESim; then the cosimulation model is built combined with the dynamic model and the driving model.As the model contains a lot of contact collision units, the calculation speed is slow, and simulating large number of samples is time-consuming; therefore, agent model is needed to replace the simulation model.According to the failure influence factors analysis, six random factors affect the action time of the deploy wing. As a result, treat the six parameters as design variables, 50 samples are generated using the optimal Latin Hypercube Sampling (LHS) [18], and then RSM model [19], Kriging model [20], and RBF model [21] are used to establish the agent model of the four wings, respectively.In order to validate the three models, 13 new samples are generated and simulated, and the predicted values are compared with the simulated values. The agent model verification of the 1st wing is shown in Figure7, as those of the other three wings are similar to the 1st wing, which are not given in the paper.Figure 7
Comparison between predicted values and actual values of different agent models.
(a)
RSM Model (b)
Kriging Model (c)
RBF ModelFrom the comparison between predicted values and actual values of different agent models, it can be seen that RBF model is more accurate than the RSM model and Kriging model. Therefore, for subsequent analysis, RBF model is used instead of the simulation model.


### 4.3. Synchronization Reliability Evaluation for the Folding Wings

Flowchart of synchronization reliability evaluation for mechanism system is shown in Figure8, First, the dynamic simulation model of the folding wing mechanism is set up. Secondly, the agent model of the action time of the folding wing is constructed, and the model is accurate enough by iteration. Again, after obtaining the distribution type and distribution parameters of the action time of each wing, the proposed method is used to evaluate the synchronous reliability.Figure 8
Flowchart of synchronization reliability evaluation for mechanism system.According to the failure mechanism analysis, action time of the folding wings is influenced by six factors. The distribution types and parameters of each factor are shown in Table6.Table 6
Distribution types and parameters of the influence factors.
No. Variables Meaning Distribute type Mean value Standard deviation 1 P Max Pressure /MPa Normal 27 1 2 D Orificing radius /mm Normal 0.6 0.00667 3 f 1 Friction coefficient Normal 0.4 0.01 4 f 4 Friction coefficient Normal 0.15 0.01667 5 k F Force coefficient Normal 1 0.05 6 k M Moment coefficient Normal 1 0.1According to distribution types and parameters, 10000 random samples are generated, and then, with the four RBF models, action times of wings are obtained, through hypothesis testing [22]; all of the action times of four wings obey the normal distribution; the histogram is shown in Figure 9 and the distribution parameters are shown in Table 7.Table 7
Distribution types and parameters of the action time the four wings.
No. of the wing Action time(ms) Distribution type Mean value Standard deviation 1 t 1 normal 137 21 2 t 2 normal 118 19 3 t 3 normal 83 13 4 t 4 normal 86 16Figure 9
Histogram of the action time of the four wings.Applying the proposed method, the synchronization reliabilityRsyn is 0.9348.


## 4.1. Failure Mechanism Analysis

As shown in Figure1(a), a missile contains four folding wings, which are powered by four actuating cylinders, respectively. By screw pair, rectilinear motion is transferred into rotatory motion to achieve the wings movement.Considering four folding wings have the same structure, the movement out of sync is mainly caused by three aspects:(1) inherent factors: motion friction coefficient is different because of manufacture errors; (2) power factor: maximum pressure and maximum flow are different due to the dosage of the dynamite and the error of damping cavity; (3) external load factors: both of the magnitude and the direction of the aerodynamic load acting on each wing are quite different. Influence factors analysis for synchronization failure is shown in Figure 6.Figure 6
Influence factors analysis for synchronization failure.Performance index of the folding wing requires that the time difference should be less than 80 ms. Therefore, the synchronization reliability can be expressed as(7)R=Pti-tj<80msi,j=1,2,3,4


## 4.2. Dynamic Model and Agent Model

The dynamic model of the folding wing is built using LMS Virtual Lab, and the driving model is built using AMESim; then the cosimulation model is built combined with the dynamic model and the driving model.As the model contains a lot of contact collision units, the calculation speed is slow, and simulating large number of samples is time-consuming; therefore, agent model is needed to replace the simulation model.According to the failure influence factors analysis, six random factors affect the action time of the deploy wing. As a result, treat the six parameters as design variables, 50 samples are generated using the optimal Latin Hypercube Sampling (LHS) [18], and then RSM model [19], Kriging model [20], and RBF model [21] are used to establish the agent model of the four wings, respectively.In order to validate the three models, 13 new samples are generated and simulated, and the predicted values are compared with the simulated values. The agent model verification of the 1st wing is shown in Figure7, as those of the other three wings are similar to the 1st wing, which are not given in the paper.Figure 7
Comparison between predicted values and actual values of different agent models.
(a)
RSM Model (b)
Kriging Model (c)
RBF ModelFrom the comparison between predicted values and actual values of different agent models, it can be seen that RBF model is more accurate than the RSM model and Kriging model. Therefore, for subsequent analysis, RBF model is used instead of the simulation model.


## 4.3. Synchronization Reliability Evaluation for the Folding Wings

Flowchart of synchronization reliability evaluation for mechanism system is shown in Figure8, First, the dynamic simulation model of the folding wing mechanism is set up. Secondly, the agent model of the action time of the folding wing is constructed, and the model is accurate enough by iteration. Again, after obtaining the distribution type and distribution parameters of the action time of each wing, the proposed method is used to evaluate the synchronous reliability.Figure 8
Flowchart of synchronization reliability evaluation for mechanism system.According to the failure mechanism analysis, action time of the folding wings is influenced by six factors. The distribution types and parameters of each factor are shown in Table6.Table 6
Distribution types and parameters of the influence factors.
No. Variables Meaning Distribute type Mean value Standard deviation 1 P Max Pressure /MPa Normal 27 1 2 D Orificing radius /mm Normal 0.6 0.00667 3 f 1 Friction coefficient Normal 0.4 0.01 4 f 4 Friction coefficient Normal 0.15 0.01667 5 k F Force coefficient Normal 1 0.05 6 k M Moment coefficient Normal 1 0.1According to distribution types and parameters, 10000 random samples are generated, and then, with the four RBF models, action times of wings are obtained, through hypothesis testing [22]; all of the action times of four wings obey the normal distribution; the histogram is shown in Figure 9 and the distribution parameters are shown in Table 7.Table 7
Distribution types and parameters of the action time the four wings.
No. of the wing Action time(ms) Distribution type Mean value Standard deviation 1 t 1 normal 137 21 2 t 2 normal 118 19 3 t 3 normal 83 13 4 t 4 normal 86 16Figure 9
Histogram of the action time of the four wings.Applying the proposed method, the synchronization reliabilityRsyn is 0.9348.


## 5. Conclusions

In summary, in view of the synchronous reliability characteristics of the mechanism, we presented a synchronization reliability evaluation method for mechanisms whose action times are differently distributed. In the method, the synchronization reliability of mechanisms is obtained with recursive algorithm, and the proposed method is validated by a numerical example. In addition, synchronization reliability of the folding wings is analyzed by the proposed method, and the engineering practicability of this method is verified. However, the proposed method is not exhaustive and universal enough to cover the synchronization of electric generators, vibration exciters, and some other objects, in which the dynamic connections exist. In addition, load-sharing characteristic of multiple mechanisms leads to the correlation between action times, of which synchronization reliability will be considered in our future work.



---

*Source: 8756975-2018-09-25.xml*