Research Article | Open Access

Volume 2015 |Article ID 937565 | https://doi.org/10.1155/2015/937565

Tong Li, Qingxuan Jia, Gang Chen, Hanxu Sun, "Motion Reliability Modeling and Evaluation for Manipulator Path Planning Task", Mathematical Problems in Engineering, vol. 2015, Article ID 937565, 13 pages, 2015. https://doi.org/10.1155/2015/937565

# Motion Reliability Modeling and Evaluation for Manipulator Path Planning Task

Accepted02 Mar 2015
Published17 Mar 2015

#### Abstract

Motion reliability as a criterion can reflect the accuracy of manipulator in completing operations. Since path planning task takes a significant role in operations of manipulator, the motion reliability evaluation of path planning task is discussed in the paper. First, a modeling method for motion reliability is proposed by taking factors related to position accuracy of manipulator into account. In the model, multidimensional integral for PDF is carried out to calculate motion reliability. Considering the complex of multidimensional integral, the approach of equivalent extreme value is introduced, with which multidimensional integral is converted into one dimensional integral for convenient calculation. Then a method based on the maximum entropy principle is proposed for model calculation. With the method, the PDF can be obtained efficiently at the state of maximum entropy. As a result, the evaluation of motion reliability can be achieved by one dimensional integral for PDF. Simulations on a particular path planning task are carried out, with which the feasibility and effectiveness of the proposed methods are verified. In addition, the modeling method which takes the factors related to position accuracy into account can represent the contributions of these factors to motion reliability. And the model calculation method can achieve motion reliability evaluation with high precision and efficiency.

#### 1. Introduction

As a kind of complicated multichain structure, manipulator can achieve various operations and has good environmental adaptability, which makes it widely used in industrial manufacturing, medical and aerospace fields, and so forth. In order to guarantee the positioning accuracy requirement of various tasks and complicated environment, the motion safety and reliability of manipulator attract extensive attention . As a result, motion reliability  which concentrates on the motion accuracy is used to give a quantitative description for the motion performance of manipulator. Generally, the motion reliability is defined as the probability that the position guarantees the accuracy requirement under the affection of various factors . Since manipulator achieves various tasks via positioning operation, the operation status of manipulator can be reflected in numerical by evaluation of motion reliability.

In order to evaluate motion reliability, the model of motion reliability should be established firstly, during which factors related to position accuracy should be considered. Firstly, the relationship between factors and position accuracy should be derived. Zhuang et al.  and Chen et al.  established the relationship between position error and parameters deviation of manipulator, respectively. Wu [13, 14] analyzed uncertain factors influencing the mechanical positioning accuracy with interval method. The methods mentioned above can be referenced in establishing motion reliability model. Pandey and Zhang  discussed the motion reliability of manipulator considering the clearance in joint independently based on probability theory. The clearance is transformed into errors of joint angles to establish the relationship with motion reliability. However, single factor considered is not enough for the modeling of motion reliability. Rao and Bhatti  evaluated the kinematics reliability and dynamics reliability of a two-link structure. Unfortunately, factors related to motion reliability were not taken into account. Moreover, the reliability evaluation for manipulator is more complicated than two-link structure, which needs further discussion.

Based on the analysis above, factors related to motion reliability are not given enough attention. And a universal model which reflects the relationship between factors and motion reliability is not achieved. The model of motion reliability should take the contribution of factors to position accuracy into account. Actually, factors related to position accuracy are various, such as clearance , friction , and wear . Many of them are unobservable and uncontrollable and even have coupling relationship between each other, which makes it difficult to establish a model including all factors. In order to reflect the influence on position accuracy caused by clearance, friction, wear, and so forth into model, kinematics parameters and joint angles can be taken as intermediate variables to deliver the influence, which has been implemented in . In this way, the motion reliability model can be established based on the deviations of kinematics parameters and joint angles.

After the establishment of motion reliability model, the model should be calculated to achieve evaluation. For this purpose, Kim et al.  used the first order reliability method (FORM) to calculate reliability. Generally, FORM can only solve reliability problem related to a few performance functions. The performance functions which are established with position accuracy related to various factors are too many for FORM to be applied. Kumar et al.  and Sharma et al.  used genetic algorithms and fuzzy methodology to analyze the reliability, respectively. Dashuang et al.  used the traditional Monte Carlo method to calculate the dynamic precision reliability of six degrees of freedom (DOF) mechanism. And the authors  have analyzed the motion reliability of a 8-DOF modular robot with Monte Carlo method based on a simple model. However, the recursive relationship between factors and motion reliability is not systematically derived in the model. And too many samples are needed in motion reliability calculation with Monte Carlo method. Actually, the fuzzy methodology and Monte Carlo method are quite complex, which leads to high computation cost during reliability computation. Novi Inverardi and Tagliani  calculated the probability density function with the principle of maximum entropy, which greatly decreased the samples in computation compared with Monte Carlo method. The method with maximum entropy shows a new way to achieve motion reliability calculation with high efficiency. Considering the practical application, a simple and feasible method for calculating motion reliability model with high precision should be proposed.

Since many tasks depend on the positioning operation of manipulator, path planning  takes an important role in achieving various operations. The paper aims at establishing the motion reliability models of path planning task and achieving evaluation. The motion reliability model of path planning task represents the probability that the position of manipulator guarantees the accuracy threshold during the task. The modeling method proposed in the paper takes the contributions of various factors to position accuracy into account. And, in order to reflect the influence of various factors, kinematics parameters and joint angles are taken as intermediate variables to deliver the influence. As a basis, the motion reliability model can be established based on the devotions of kinematics parameters and joint angles. Considering the complex of multidimensional integrals in the model of continuous trajectory tracking, the approach of equivalent extreme value is introduced to achieve model simplification. As a result, multidimensional integrals for PDF are converted into one dimensional integral. Since huge samples and computation time are needed in traditional method for PDF calculation, a method based on the maximum entropy principle is proposed to achieve calculation with high precision and efficiency. So far, the evaluation of motion reliability can be achieved, and an intuitionistic expression for the performance of manipulator can be obtained with the evaluation result. Meanwhile, calibration and control can be devoted to improving the motion reliability of manipulator for further discussion.

In conclusion, the paper is organized as follow. In Section 2, the mathematical relationship between motion reliability and position accuracy is firstly derived. The influence of the factors related to the position accuracy is delivered to the deviations of kinematics parameters and joint angles, based on which the motion reliability models of path planning task are established. In Section 3, the approach of equivalent extreme value is introduced to simplify the motion reliability model, and method based on the maximum entropy principle is proposed to achieve model calculation. In Section 4, simulations are carried out to verify the effectiveness and correctness of the proposed modeling and evaluation method. And strategies for improving the motion reliability are discussed further based on the evaluation result. The last part is summary.

#### 2. Motion Reliability Models of Path Planning Task

Path planning tasks are always divided into point-to-point path planning and continuous trajectory tracking, and their motion reliability models have different focuses. For the model of continuous trajectory tracking, position accuracy of the entire trajectory should be considered, while, for the model of point-to-point path planning, position accuracy of target point is paid more attention. Aiming at the two kinds of path planning, motion reliability models are established based on the analysis of factors related to position accuracy.

##### 2.1. The Mathematical Expression for Motion Reliability of Path Planning

The motion reliability for a single trajectory point is discussed. Set position accuracy threshold (namely, the maximum acceptable value of the deviation between the actual and desired position) as at each direction of the point. Then an enveloping space is formed as a sphere. The desired position is set as center and is set as radius. As a result, the Cartesian space at the end-effector can be divided into three subspaces shown in Figure 1. The actual position locating in the sphere is considered to be motion reliable, while it is considered to be motion unreliable when actual position locates outside the sphere. The actual positions locating on the shell are critical positions, which can either be reliable or unreliable at the next movement. Based on the analysis above, the mathematical expression when path planning task is motion reliable can be achieved.

(a) Point-to-Point Path Planning. and represent the actual and desired values of the target point, respectively. represents factors related to position accuracy. represents the position accuracy threshold; then the expression when point-to-point path planning task is motion reliable can be expressed as

Suppose as the PDF related to position accuracy of the target point, and make ; then (1) can be turned into

(b) Continuous Trajectory Tracking. When continuous trajectory tracking is motion reliable, the entire trajectory should guarantee the position accuracy threshold. Define and as the actual and desired position at point . Assume position accuracy threshold for every trajectory point is the same, which is expressed as . The number of trajectory points is . Then the mathematical expression when continuous trajectory tracking is motion reliable can be expressed as

Define the PDF related to the position accuracy of the entire trajectory points as , and make ; then (3) can be turned into

##### 2.2. Analysis of Factors Related to Position Accuracy

In (2) and (4), represents the factors related to position accuracy. In practice, there are too many factors influencing position accuracy, such as clearance, friction, wear, deformation, assembly tolerances, and flexibility. However, many of them are unobservable and uncontrollable; the accurate model between these factors and position accuracy is hard to obtain. Furthermore, not all the factors have direct influence on the position accuracy in Cartesian space. For example, clearance and friction have influence on the accuracy of joint angles, while assembly tolerances and deformation mainly have influence on geometry parameters. Therefore, the influence caused by these factors can be represented by the deviation of kinematics parameters and joint angles, which is considered as intermediate variables to deliver the influence to the position accuracy at the end-effector. Besides, the deviations of kinematics parameters and joint angles are easy to be measured and adjusted. In this way, deviations of DH parameters are used to form the factor . And motion reliability can be modeled based on the relationship between deviations of DH parameters and position accuracy.

For manipulator, the transformation matrix between th and th coordinate system is described based on parameters , , :

and represent the actual and nominal transformation matrix between th and th coordinate system, respectively. The deviation between them is expressed as

Due to (plug (5) into (6)), it can be obtained that

Combining with (7), the following equation is obtained:wherein and are the position deviation and orientation deviation of th coordinate system relative to th coordinate system, respectively. Meanwhile, . Consider

For the entire manipulator, the nominal transformation matrix from the coordinate system of the end-effector to the base coordinate system can be expressed as

and represent the actual and nominal transformation matrix from the coordinate of the end-effector to the base coordinate system, respectively. Then,wherein is the differential of the transformation matrix from the coordinate of the end-effector to the base coordinate system. First-order approximation is used for (11). Then, according to ,

Define . Combining with , it is concluded as

and represent the position and orientation deviation; they can be expressed as

Combine (8) with (14); the following relationship is concluded:wherein and , , , , which can be expressed in detail as

Then the position deviation of the end-effector can be obtained from (15):

So far, relationship between deviations of kinematics parameters and position error is derived. Meanwhile, the contribution of each parameter to the position accuracy is also obtained as matrix . In fact, the deviations of kinematics parameters are the intermediate variables to reflect the influence of factors related to position accuracy. Then motion reliability can be obtained via the integral for the PDF of position accuracy. The relationships are shown in Figure 2.

Figure 2 gives a schematic view of the influence relationship between factors and motion reliability. When different factors are considered, their effect can be delivered via kinematics parameters. Thus, each kinematics parameters can be expressed as a function, such aswherein is the th factor which can deliver the influence on the position accuracy to parameter . Other kinematics parameters can be expressed in a similar form. If (18) can be linearized, the contribution weight of each factor can be obtained as which is listed in Figure 2. However, many problems exist in the derivation and linearization of functions like (18), which will not be discussed further in the paper.

##### 2.3. Motion Reliability Models of Path Planning Task

In summary, combining (17) with (2) and (4), the motion reliability models for point-to-point path planning and continuous trajectory tracking can be obtained, respectively. The motion reliability model for point-point path planning is expressed as

The motion reliability model for continuous trajectory tracking is expressed as

So far, models of motion reliability for path planning task are obtained. According to (19) and (20), the models mean to obtain the probability that the position deviation guarantees the position accuracy threshold. Meanwhile, deviations of kinematics parameters are taken into account, which successfully introduces the influence of factors related to position accuracy into the models. Then the models should be calculated to achieve motion reliability evaluation, which will be discussed in the next section.

#### 3. Simplification and Calculation for Motion Reliability Model

In order to calculate motion reliability models of path planning task, samples of position deviations should be firstly obtained by introducing errors into factors. As a basis, PDF related to the samples can be computed and the motion reliability can be calculated by integral for PDF. For point-to-point path planning, motion reliability can be calculated according to (19) by one dimensional integral for PDF. However, it is more complicated to calculate the motion reliability of continuous trajectory tracking because there are so many trajectory points. The PDF in the task is concerned with position deviations of entire trajectory points; multidimensional integrals for PDF should be carried out in motion reliability evaluation, which brings about incredible computation cost simultaneously. In order to solve model (20) simply, approach of equivalent extreme value is introduced to simplify the model.

##### 3.1. Simplification for Motion Reliability Model

Define theorem as follows.

Theorem 1. Suppose are random variables; define the maximum value of as ; then it goes that .

Proof. Firstly, prove .
Define as the PDF related to ; then,Because the maximum value of , can be expressed asit can be derived thatConsider is proved.
The maximum value of is expressed as . Define , , and then .
Combing with , it goes that

According to the theorem, define as the maximum position deviation of entire trajectory points; then motion reliability model of (20) can be turned intowherein represents the PDF related to the maximum position deviation of the manipulator. Thus, the multidimensional integrals for the PDF related to position deviation of the entire trajectory points can be turned into one dimensional integral for PDF related to the maximum position deviation. This method is called the approach of equivalent extreme value, with which the motion reliability model of continuous trajectory tracking can be simplified and easily solved.

##### 3.2. Calculation for PDF Based on Maximum Entropy Principle

After the model is simplified, the core of motion reliability evaluation is the calculation of PDF. Generally, Monte Carlo method can be used to solve any problems about probability distribution. However, massive samples are needed and computation cost is high. In order to achieve fast calculation of PDF with high precision, the information entropy is introduced and a method based on the principle of maximum entropy is proposed. The principle of maximum entropy  was firstly proposed by Jaynes, which aimed at choosing the probability distribution at the maximum state of information entropy. The distribution obtained at maximum entropy state is considered to be the most realistic. Meanwhile, the probability distribution based on the maximum entropy principle has high computation precision, which has been verified in  by comparison with Monte Carlo method. Therefore, the method based on the maximum entropy principle has huge advantages in fast and precise computation.

When evaluating the motion reliability of continuous trajectory tracking, the PDF can be calculated as follows. Suppose as the maximum position deviation of the entire trajectory, and represents the PDF of . Then the information entropy of can be expressed as

Equation (26) shows the relationship between the entropy and PDF. The PDF of continuous trajectory tracking task should meet the constraints shown in (27a) and (27b), wherein constraint (27a) means that the integral for PDF equals 1. Constraint (27b) is a group of constraints, which denote the moment of . The number of constraints (27b) depends on the moment . During constraint (27b), the moment can be either integer or fraction:

As the maximum entropy principle describes, the PDF is the most realistic when the information entropy is maximum. In order to find the maximum value of entropy, Lagrange multiplier method  is introduced. Suppose as Lagrange multipliers; the Lagrange function related to information entropy is established as

Take partial derivative with respect to , and make . Then,

The expression of at the maximum entropy state can be derived:

Meanwhile, can be obtained by plugging (27a) into (30):

In order to calculate the unknown constants in , Kullback-Leibler (K-L) divergence  is introduced, which represents the deviation between and the actual value. The expression is shown as follows:

Substitute (26), (27a), (27b), (30), and (31) into (31):

When is the closest to the actual value, the K-L divergence is minimal. Suppose and . Since is unrelated to and , is minimal, meaning the following formula is minimal:

According to (27b), the moment of equals . can be approximately expressed as based on moment estimation, where is the number of samples. Therefore, the calculation of the PDF is turned into

Set the initial value of ; then and can be iteratively obtained via search algorithm until is minimum. Substituting and into (30) and (31), the value of and the expression of can be achieved. Finally, the motion reliability of continuous trajectory tracking can be calculated by one dimensional integral for .

##### 3.3. Process for Calculating Motion Reliability Models of Path Planning Task

Based on the analysis above, the process for calculating the motion reliability of path planning task can be concluded, which is also shown in Figure 3.(a)Motion reliability model is firstly established for a particular path planning task. For model of continuous trajectory tracking, the approach of equivalent extreme value is introduced. As a result, multidimensional integrals for PDF related to the entire trajectory points are turned into one dimensional integral for PDF related to the maximum position deviation.(b)Introduce errors into factors related to position accuracy of the model, and simulate the path planning task for multiple times. For point-to-point path planning, position deviation at the target point of each simulation is obtained. For continuous trajectory tracking, position deviations at each trajectory point are obtained.(c)For point-to-point path planning, the position deviation at the target point of each simulation is selected as samples. For continuous trajectory tracking, the maximum position deviation of each simulation is selected as samples.(d)Lagrange function about information entropy of position deviation samples is established, and the expression of PDF related to position deviation samples is derived at the state of maximum entropy.(e)The unknown coefficients of the PDF are calculated by using K-L divergence. Then the motion reliability evaluation is achieved by one dimensional integral for PDF.

#### 4. Simulation on Motion Reliability Evaluation of Path Planning Task

A 7-DOF manipulator is used for motion reliability evaluation of path planning task. The coordinate systems are defined as shown in Figure 4, while the nominal DH parameters are listed in Table 1. A particular continuous trajectory task is defined for simulation as follows.

 Link (°) (m) (m) (°) 1 0 0.6 0 90 2 90 0.5 0 −90 3 0 0.5 5 0 4 0 0.5 5 0 5 0 0.5 0 90 6 −90 0.5 0 −90 7 0 0.6 0 0

The initial configuration of manipulator is set as in joint space, and the position of the end-effector at this configuration is denoted by point . The target pose is set as , which is denoted by point . Continuous trajectory tracking task is carried out from to . Task cycle is 20 s, and control cycle is 0.05 s. As a result, 400 trajectory points are planned during the task.

##### 4.1. Motion Reliability Evaluation of a Particular Path Planning Task

During the motion reliability evaluation of task, errors are introduced into kinematics parameters. Errors following normal distribution are introduced into , , and errors following normal distribution are introduced into , . The task is simulated for 500 times and position deviations of the entire trajectory points are collected for each simulation . The maximum position deviation is selected from each simulation to obtain 500 samples .

With the samples of maximum position deviation, PDF of the continuous trajectory tracking task is calculated. The initial value of is set as , and the order of moment is set as . Search algorithm is carried out according to (35). , are obtained as and when is minimum. As a basis, the curve of PDF related to the position accuracy threshold is obtained as in Figure 5. Via one dimensional integral for the PDF, the motion reliability can be calculated at different accuracy threshold, which is shown in Figure 6.

Meanwhile, the PDF and probability curves calculated by Monte Carlo method are also provided in Figures 5 and 6. With Monte Carlo method, the planning is simulated for 10000 times to obtain results. By comparison, the PDF and probability curves calculated with maximum entropy principle have the same accuracy with the results obtained with Monte Carlo method. At the same accuracy threshold, the deviation of calculated probability with the two methods is only about 1%. For example, when accuracy threshold is  m, the motion reliability is 0.9167 calculated by maximum entropy principle, while it equals 0.9064 with Monte Carlo method. It is verified that the evaluation result with maximum entropy principle has high precision.

Besides, huge computation cost and samples requirement limit Monte Carlo method in computing the PDF of trajectory tracking task. Nearly 3000 seconds are used to handle the 10000 samples to achieve the PDF. However, the method based on the maximum entropy principle needs fewer samples (just 500 samples), which greatly reduces the computation cost. In conclusion, the motion reliability evaluation method based on the maximum entropy principle is feasible with high efficiency and computation precision.

##### 4.2. Motion Reliability Evaluation with Different Trajectory Points

When evaluating motion reliability of continuous tracking task, the number of trajectory points will have influence on the motion reliability. If the number is small, the characteristics of the trajectory cannot be reflected completely, which leads to a credible evaluation result. In addition, with the number increase, the computation cost will become higher. Therefore, proper number of trajectory points should be selected by taking both computation efficiency and precision into account.

The trajectory mentioned above is taken as an example, and the motion reliability of the task is evaluated with different trajectory points. The position accuracy threshold is set as  m. Motion reliability varies with the number of trajectory points, which is shown in Figure 7. It is obvious that when the number of trajectory points is small, the motion reliability changes greatly with the number of trajectory points, because small numbers of trajectory points cannot reflect the characteristics of the trajectory completely. After the number is larger than 200, the reliability becomes stable. After the number is larger than 400, the reliability is nearly constant. The increased trajectory points only bring about higher computation cost in motion reliability evaluation. In summary, 400–500 is the proper number for the trajectory , which can completely reflect the characteristics of the trajectory and show good operational efficiency.

Besides, trajectory points are determined by task cycle and control cycle. The motion reliability evaluation will benefit from the proper number of trajectory points, which also indicates that the arrangement for task cycle and control cycle of the task is reasonable. In this way, reasonable arrangement for trajectory points is the basis to achieve fast and precise motion reliability evaluation.

##### 4.3. Contributions of Factors Related to Motion Reliability

Factors influencing the position accuracy of manipulator can be included into the model of motion reliability. Via introducing errors into them, motion reliability of a particular path planning task can be evaluated, which is achieved in Section 4.1. In this section, the contributions of factors related to the motion reliability are discussed. , and , are differed into two groups for different units. Errors following the same normal distribution are introduced into each group, respectively. is introduced into , , and the decrease in motion reliability caused by , is shown in Figure 8(a). is introduced into , ; the decrease in motion reliability caused by , is shown in Figure 8(b).

It is obvious that decreases caused in motion reliability are different within different position accuracy threshold for one factor. Compared with the matrix in (17) which shows the contribution of factor at a moment in the task, Figure 8 can comprehensively reflect the influence of the factor on motion reliability during the entire task cycle. Further, the decrease in reliability caused by is larger than within any accuracy threshold, while the decreases caused by and are basically the same. The decrease values in motion reliability caused by different factors at position accuracy threshold are also given out in the figure. Thus, the contributions can be expressed as , with the same errors. Meanwhile, the factors have coupling effect on motion reliability. The decreases in reliability caused by factors are listed at the accuracy threshold 0.004 m in Table 2.

 Angular factors Stand. deviation of error (°) Decrease of reliability Coupling effect 0.01 0.08733 0.3809 0.01 0.1197 Length factors Stand. deviation of error (m) Decrease of reliability Coupling effect 0.001 0.1174 0.3608 0.001 0.1121

In Table 2, when factors with the same unit act together, the decrease of reliability is greater than the sum of the decrement when they act separately. For example, decrement caused by and is larger than the sum of decrement caused by them separately. Besides, when the standard deviations of introduced errors are different, the PDF and decrease in motion reliability are shown as in Figure 9 by analyzing factor . It can be concluded that the decrease in motion reliability depends on the position accuracy threshold and standard deviations of errors, and the variation of which can be expressed by quadric surface as the figure shows. At the same position accuracy threshold, with the increase of standard deviation in errors, the decrease in motion reliability becomes larger. Variations of other factors can be expressed in similar ways, which are not given out here.

Thus the contributions of kinematics parameters to motion reliability are analyzed, and the coupling effect of parameters is discussed, too. Meanwhile, the variation according to the position accuracy threshold and standard deviation of errors is derived. As a result, the importance of each kinematics parameter in improving motion reliability can be obtained. When improving the motion reliability of manipulator, factors with high contribution should be given high priority to handle. According to the different contributions of parameters, the strategy to achieve motion reliability improvement is discussed in the next section.

##### 4.4. Discussion of Trajectory Optimization Based on Maximizing Motion Reliability

Since , , , and are easy to be observed and handled during task, the motion reliability can be improved by reducing the errors in these factors based on the contributions to motion reliability. For example, since contributes more than , , and , based on decreasing errors with kinematics calibration, the accuracy of should be improved as much as possible. Thus, aiming at maximizing motion reliability, a trajectory optimization method is discussed via analyzing the contribution of .

Errors following normal distribution with different standard deviation are introduced into . Motion reliability curves are achieved with the changes of position accuracy thresholds, which are shown in Figure 10. With the increase of errors, the motion reliability decreases obviously at the same accuracy threshold. At accuracy threshold 0.005 m, when increases from to , the motion reliability decreases . Thus, the motion reliability can be improved signally by reducing errors.

Because of the time-varying characteristics of gear clearance, transmission error, and friction, the motion errors at different joint angles are different, which can be established as . Generally, every joint is limited into a certain range, which is determined by the characteristics of joint motor and configuration requirement of manipulator. For redundant manipulator, infinite feasible trajectories exist between and . When position accuracy is guaranteed well in the entire trajectory, the motion reliability is high. Considering function , a proper trajectory should be chosen to arrange the joint angle at the range with smaller introduction of errors. Thus, strategy aiming at maximizing motion reliability can be described as follows.

Task constraints should be fulfilled firstly; namely, the position accuracy should be guaranteed in the entire trajectory, which can be expressed as . Then threshold constraints of parameters should be fulfilled during task, such as joint velocity and acceleration, which have the expression of . The objective function is established as the motion reliability of the trajectory tracking task. The output is motion reliability; the input is the joint angle range related to the selected trajectory. Then the expression of the trajectory optimization can be concluded:

In trajectory , the corresponding range for joint is arranged as , during which the introduced joint angel errors are minimum. The objective function represents the motion reliability of the task, during which represents the PDF related to factors .

#### 5. Conclusion

Motion reliability as a criterion can reflect the motion performance of manipulator synthetically. With the modeling and evaluation of motion reliability, the status of systems can be obtained and monitored in numerical during task operation. This paper establishes a universal model for motion reliability of a particular path planning task. The model which takes the factors related to position accuracy (such as clearance, wear, and friction) into account can reflect the contribution of these factors to motion reliability. Meanwhile, motion reliability can be evaluated by calculating the model. For model calculation, a model simplification method based on the approach of the equivalent extreme value is introduced. With this method, the multidimensional integral for PDF can be converted into one dimensional integral for PDF. Further, a method based on the maximum entropy principle is proposed to obtain the PDF. This method needs less samples and computation cost than the traditional Monte Carlo method, which makes it successful to achieve motion reliability evaluation with high precision and efficiency.

In the paper, the authors mean to establish a motion reliability model which can reflect the influence of factors related to position accuracy, such as clearance, wear, and friction,. Since it is difficult to derive the relationship between these factors and position accuracy, kinematics parameters are considered to be the intermediate variables to deliver the influence of the factors. Thus, the relationship between deviations of kinematics reliability and position accuracy is derived, based on which the model is established. However, the relationships between factors and kinematics parameters are not discussed in the paper, which should be devoted to further research.

From the model and motion reliability evaluation results, the contribution of each kinematics parameter to motion reliability can be concluded. According to the contribution, the importance of each kinematics parameter in improving motion reliability can be decided. When improving the motion reliability of manipulator, factors with high contribution should be given priority to handle. For example, in Section 4.4, because the joint angle contributes more than other parameters, the errors should be decreased firstly in operations. According to the motion reliability model and evaluation result, trajectory optimization or control strategy aiming at improving motion reliability of manipulator is a new optimization problem, which is meaningful for further research.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

This work is supported by the National Key Basic Research Program of China (2013CB733000), the National Natural Science Foundation of China (61175080), and Specialized Research Fund for the Doctoral Program of Higher Education (20120005120004).

1. B. S. Dhillon, A. R. M. Fashandi, and K. L. Liu, “Robot systems reliability and safety: a review,” Journal of Quality in Maintenance Engineering, vol. 8, no. 3, pp. 170–212, 2002. View at: Publisher Site | Google Scholar
2. J. Xia, Z. N. Jiang, H. Liu, H. Cai, and G. Wu, “A novel hybrid safety-control strategy for a manipulator,” International Journal of Advanced Robotic Systems, vol. 11, no. 1, pp. 58–67, 2014. View at: Publisher Site | Google Scholar
3. S. Rubrecht, V. Padois, P. Bidaud, M. de Broissia, and M. da Silva Simoes, “Motion safety and constraints compatibility for multibody robots,” Autonomous Robots, vol. 32, no. 3, pp. 333–349, 2012. View at: Publisher Site | Google Scholar
4. J. Xia, Z. Jiang, H. Liu, H. Cai, and G. Wu, “A manipulator's safety control strategy based on fast continuous collision detection,” in Proceedings of the IEEE International Conference on Robotics and Biomimetics (ROBIO '13), pp. 2380–2385, IEEE, December 2013. View at: Publisher Site | Google Scholar
5. C. Zhang and G. Bai, “Reliability analysis of a two-link flexible robot manipulator,” in Proceedings of the 2nd International Conference on Information Science and Engineering (ICISE '10), pp. 2269–2272, Hangzhou, China, December 2010. View at: Publisher Site | Google Scholar
6. X. Huang and Y. Zhang, “Probabilistic approach to system reliability of mechanism with correlated failure models,” Mathematical Problems in Engineering, vol. 2012, Article ID 465853, 11 pages, 2012. View at: Publisher Site | Google Scholar
7. C. Carreras and I. D. Walker, “Interval methods for improved robot reliability estimation,” in Proceedings of the Annual Reliability and Maintainability Symposium, pp. 22–27, January 2000. View at: Google Scholar
8. J. Li, J.-B. Chen, and W.-L. Fan, “The equivalent extreme-value event and evaluation of the structural system reliability,” Structural Safety, vol. 29, no. 2, pp. 112–131, 2007. View at: Publisher Site | Google Scholar
9. E. Taufer, S. Bose, and A. Tagliani, “Optimal predictive densities and fractional moments,” Applied Stochastic Models in Business and Industry, vol. 25, no. 1, pp. 57–71, 2009. View at: Publisher Site | Google Scholar | MathSciNet
10. C. T. J. Xun, C. Zhang, and Y. Jiang, Electrical System Reliability Engineering, Science Press, Beijing, China, 2010.
11. H. Zhuang, L. K. Wang, and Z. S. Roth, “Error-model-based robot calibration using a modified CPC model,” Robotics and Computer Integrated Manufacturing, vol. 10, no. 4, pp. 287–299, 1993. View at: Publisher Site | Google Scholar
12. G. Chen, Q. Jia, T. Li, and H. Sun, “Calibration method and experiments of robot kinematics parameters based on error model,” Jiqiren/Robot, vol. 34, no. 6, pp. 680–688, 2012. View at: Publisher Site | Google Scholar
13. W. Wu and S. S. Rao, “Uncertainty analysis and allocation of joint tolerances in robot manipulators based on interval analysis,” Reliability Engineering & System Safety, vol. 92, no. 1, pp. 54–64, 2007. View at: Publisher Site | Google Scholar
14. W. Wu and S. S. Rao, “Interval approach for the modeling of tolerances and clearances in mechanism analysis,” Journal of Mechanical Design, Transactions of the ASME, vol. 126, no. 4, pp. 581–592, 2004. View at: Publisher Site | Google Scholar
15. M. D. Pandey and X. Zhang, “System reliability analysis of the robotic manipulator with random joint clearances,” Mechanism and Machine Theory, vol. 58, pp. 137–152, 2012. View at: Publisher Site | Google Scholar
16. S. S. Rao and P. K. Bhatti, “Probabilistic approach to manipulator kinematics and dynamics,” Reliability Engineering and System Safety, vol. 72, no. 1, pp. 47–58, 2001. View at: Publisher Site | Google Scholar
17. H. P. Jawale and H. T. Thorat, “Positional error estimation in serial link manipulator under joint clearances and backlash,” Journal of Mechanisms and Robotics, vol. 5, no. 2, Article ID 021003, 2013. View at: Publisher Site | Google Scholar
18. H. Abdellatif and B. Heimann, “On compensation of passive joint friction in robotic manipulators: Modeling, detection and identification,” in Proceedings of the Computer Aided Control System Design, IEEE International Conference on Control Applications, IEEE International Symposium on Intelligent Control, pp. 2510–2515, Munich, Germany, October 2006. View at: Publisher Site | Google Scholar
19. Z. F. Bai, Y. Zhao, and J. Chen, “Dynamics analysis of planar mechanical system considering revolute clearance joint wear,” Tribology International, vol. 64, pp. 85–95, 2013. View at: Publisher Site | Google Scholar
20. J. Kim, W.-J. Song, and B.-S. Kang, “Stochastic approach to kinematic reliability of open-loop mechanism with dimensional tolerance,” Applied Mathematical Modelling, vol. 34, no. 5, pp. 1225–1237, 2010. View at: Publisher Site | Google Scholar | MathSciNet
21. N. Kumar, J.-H. Borm, and A. Kumar, “Reliability analysis of waste clean-up manipulator using genetic algorithms and fuzzy methodology,” Computers & Operations Research, vol. 39, no. 2, pp. 310–319, 2012. View at: Publisher Site | Google Scholar
22. S. P. Sharma, D. Kumar, and A. Kumar, “Reliability analysis of complex multi-robotic system using GA and fuzzy methodology,” Applied Soft Computing Journal, vol. 12, no. 1, pp. 405–415, 2012. View at: Publisher Site | Google Scholar
23. L. Dashuang, W. Kesheng, Z. Longlong et al., “Dynamic precision reliability analysis for six degrees of freedom micro-displacement mechanism of reticle stage in lithography machine based on Monte-Carlo,” in Proceedings of the International Conference on Quality, Reliability, Risk, Maintenance, and Safety Engineering (QR2MSE '13), pp. 250–252, 2013. View at: Google Scholar
24. T. Li, Q. X. Jia, G. Chen et al., “Motion reliability assessment of robot for trajectory tracking task,” Systems Engineering and Electronics, vol. 36, no. 12, pp. 2556–2561, 2014. View at: Google Scholar
25. P. L. Novi Inverardi and A. Tagliani, “Maximum entropy density estimation from fractional moments,” Communications in Statistics. Theory and Methods, vol. 32, no. 2, pp. 327–345, 2003. View at: Publisher Site | Google Scholar | MathSciNet
26. G. Chen, L. Zhang, Q. X. Jia, M. Chu, and H. Sun, “Repetitive motion planning of free-floating space manipulators,” International Journal of Advanced Robotic Systems, vol. 10, pp. 253–263, 2013. View at: Publisher Site | Google Scholar
27. E. Jaynes, “Review of ‘principles of statistical mechanics—the information theory approach’ (Katz, A.; 1967),” IEEE Transactions on Information Theory, vol. 14, no. 4, pp. 611–612, 1968. View at: Publisher Site | Google Scholar
28. X. Li, The Method Study about Probability Distribution Based on the Principle of Maximum Entropy, North China Electric Power University, Beijing, China, 2008.

#### More related articles

We are committed to sharing findings related to COVID-19 as quickly as possible. We will be providing unlimited waivers of publication charges for accepted research articles as well as case reports and case series related to COVID-19. Review articles are excluded from this waiver policy. Sign up here as a reviewer to help fast-track new submissions.