Abstract
The kinematic reliability of robots, defined as the probability that the end-effector falls inside the specified safe boundary, is of great significance in predicting the accuracy achieved in reality. This work selects the 7 degrees-of-freedom (7-DOF) redundant robot as an example to conduct reliability analysis by utilizing the envelope method against time-related issues in this work. Since variables in industrial robots are very small relative to their means, the motion error functions are commonly linearized by the first-order Taylor’s formula to simplify calculation, and the failure models in all directions and attitude angles are then established through the probability method over the entire input interval. As a result, the actual accuracy of the robot in each pose component will be displayed, instead of merely considering the position error like other scholars. The principle of the proposed method is to transform a time-dependent problem into a time-independent one with the help of the failure extreme points and endpoints, so as to enhance the operation efficiency under the premise of ensuring accuracy. Finally, the simulation results verify that the relative error of the envelope method is less than 6.0% compared with that of the Monte Carlo simulation method, and the computational efficiency is higher than that of the Monte Carlo method, which demonstrates that the envelope method has better comprehensive performance.
1. Introduction
Industrial robots are critical components of modern manufacturing because they offer low production costs, high stability, and high efficiency [1–3]. The 7-DOF redundant robot serves as a current research hotspot, famous for its tremendous flexibility in completing multiple tasks like avoiding singularity, avoiding obstacles, and optimizing pathways. Although the structure of multiaxis and multilink makes operations more convenient, they further magnify the uncertainty of the ontology, leading to the error accumulation problem becoming more prominent. Since errors universally exist in the manufacturing, assembly, and joint rotation processes [4,5], the actual pose of the end-effector may differ from the predetermined one where the precision is in perfect condition. Therefore, kinematic reliability, which is widely used in mechanical structures, is proposed as a method of uncertainty evaluation to quantify the impact of multiple factors on accuracy.
For the inverse kinematics of 7-DOF redundant robots, studies have shown that some algorithms can make the pose of the manipulator almost identical to the ideal state. Dereli and Köker [6] employed a group optimization technique called the firefly algorithm to narrow the error between the end-effector and the target point to mm. This did not, however, imply that robots will be able to achieve such precision in practical engineering owing to the influence of external factors that cannot be eliminated. In general, the link error of industrial robots ranges from 0.01 mm to 1 mm [7,8]. On the other hand, the rotation error, as well as other errors, cannot be overlooked. These variables, in contrast to the inverse kinematics algorithm, significantly disrupt the working condition of the manipulator. Furthermore, intelligent welding work performed by robots in conjunction with vision systems has been progressing apace in recent years [9], whose kinematic error can be effectively predicted based on the reliability theory. All of this suggests that kinematic robot reliability is an essential study area with several applications.
The kinematic reliability problem of robots can be divided into two categories: point reliability and time-dependent reliability [10–13]. The instantaneous failure rate at a point on the trajectory is reflected in point reliability regardless of whether motion is valid at other points, which commonly occurs during the measurement of the repeatability of positioning. Due to the low complexity and frequent occurrences of this problem, it has become the primary research type for scholars. The first-order reliability method (FORM) and the Monte Carlo simulation (MCS) method are two typically used calculation methods [14,15], with the latter serving as a benchmark for the new method in that it has higher accuracy based on a large number of sample simulations but higher calculation cost [16]. In addition, Zhang and Xu [17] used a combination of saddle point approximation (SPA) and dimension reduction (DR) to forecast the precision of a 6-DOF industrial robot, filling a gap in the nonlinear error literature. However, in engineering practice, the reliability of continuous trajectory is of higher application value since it can reveal the laws of robot terminal motion in global terms, which is why the time-dependent problem of kinematic reliability receives more attention.
To be more specific, it refers to the likelihood that an arbitrary discrete point on its terminal trajectory falls within a presupposed tolerance zone, the center of which is on the desired path [18–20]. Nevertheless, the true motion accuracy is difficult to estimate on account of the high correlation between discrete points at different times and the accumulation of failure rates over time. Regrettably, almost all time-dependent problem research objects are four-bar mechanisms with a single output and few variables in the last decade, where popular algorithms such as Rice’s formula, extreme value method, and envelope method [21–23] have been successfully implemented.
With the improvement of relevant theories over time, some algorithms are appropriately used for robot time-dependent reliability analysis after reasonable modification. Zhao et al. [24] introduced the first-passage method into the analysis of the 6-DOF robot’s pose error, showing that the method’s accuracy is fairly good when accounting for low relationships between discrete points, but the calculation area is limited to situations with a very low failure rate. Pendey and Zhang [25] effectively computed the fluctuation range around the intended motion, and the probability density function of the extreme value distribution points is derived in terms of the principle of maximum entropy. Whereas this method requires a large amount of sample data, its stability may be insufficient for widespread adoption. In order to solve the problem that the failure probability of the 3R planar robot is difficult to estimate, Vieira et al. [26] coupled the interval analysis method with the MCS method, taking advantage of both to formulate the law of plane motion, but it cannot guarantee its adaptability to spatial motion. Given the considerable limitations of the method mentioned above in the investigation of time-dependent reliability, a new method is indispensable to assess the robot from various perspectives.
Motivated by the abovementioned studies and the huge application potential of kinematic reliability, this article chooses the 7-DOF redundant robot as the research object. The main contributions are threefold: (1) compared with planar mechanisms or 6-DOF manipulators, the 7-DOF redundant robot involves more input variables and higher coupling in motion function, whose time-dependent reliability problem has not been studied yet. (2) The simulation analysis in this article, conducted by the enveloped method, including three positions and three attitude angles, is remarkably distinct from those that just focus on position errors to better display the components of the total error. (3) The error function of the comprehensive pose obeying nonlinear is approximated on the basis of the multidimensional normal integral.
The framework of this article is as follows. In Section 2, we elaborate on the robotic kinematics. Section 3 deduces error functions and kinematic reliability models. We propose the envelope method in Section 4. A simulation example is presented in Section 5, and the conclusions are summarized in Section 6.
2. Robot Kinematics
To begin with, the forward kinematics of the robot should be established to construct the error function, which describes the transfer relationship between connecting rod movements. A classic method, the improved DH parameter [27], identifies the internal relativity between two neighboring links through four parameters in the definite coordinate systems, as shown in Figure 1. The symbols , , , and of four parameters represent link length, link twist, link offset, and joint angle, respectively, listed in Table 1.

The homogeneous transformation matrix of adjacent frames is described aswhere is the abbreviation of the symbol and is the abbreviation of the symbol .
According to the progressive relationship between the connecting rods, the pose matrix of end-effector is given bywhere , , and are all unit vector projected onto the fixed coordinate frame. is a position vector in the fixed coordinate frame.
3. Kinematic Reliability Model
The kinematic reliability of a robot can be defined as the probability that the actual output trajectory of its end-effector in the specified interval falls within the allowable error range. In order to vividly describe the influence of uncertain factors on the accuracy of the trajectory, Figure 2 is used to depict the difference between the failed trajectory (red curve), the qualified trajectory (black curve), and the ideal trajectory (purple curve).

In the past, it has been discovered that the error of a robot is primarily caused by the geometric dimension of the connecting rod and the joint rotation [8]. Furthermore, link offset should be considered due to inevitable deviation from assembly and joint wear, which is easily overlooked. In engineering practice, however, the standard deviation of each variable is controlled within a small range, usually less than one-thousandth relative to the corresponding mean value. Based on the research results from the literature [11,17,28], the error distribution of each variable is assumed to be normal and mutually independent. Then, the error function of end-effector in three positions and three attitude angles, is described aswhere and are the actual output value and theoretical output value, respectively, and is the dimension of the vector .
When the magnitude of standard deviation of the random variable is much smaller than the counterpart of the mean value, the first-order Taylor formula [29] is an effective treatment for approximating the error function. As a result, the error function can be expanded at the mean value of , written aswhere , .
Remark 1. In this article, all variables follow normal distribution, so the correlation property can be used to further simplify the operation.
Property 1. If , there iswhere a and b are both real numbers.
Therefore, subject to the normal distribution should be transformed into subject to the standard normal distribution, given bySo equation (3) can be rewritten aswhere , , is the -th element of , and is the standard deviation of .
Since the joint angle changes with time , can be converted into a time-dependent function. If moving time is set in the interval , taking the position of the end-effector as an example, its kinematic reliability model can be defined asThe corresponding failure probability model can be expressed asSimilarly, the comprehensive reliability model of the robot end-effector isThe comprehensive pose failure model of the robot end-effector isAccording to the preceding derivation, if a trajectory is eligible in interval , the error of a point at any time must satisfy the condition that the error is within the tolerance sphere. Nonetheless, calculating the error at each point in reality is unaffordable. As a result, it is critical to develop a high-precision and high-efficiency algorithm to replace the MCS method for solving the kinematic reliability problem.
4. Envelope Method
The linchpin of the envelope method is to utilize failure boundary functions and converting reliability problem into a time-independent problem. and are limit state functions, which are tangent to hypersurfaces and , respectively.
According to this nature, there is
If , we have
Then, is expressed as follows:
We solve equation (14) to obtainwhere is the inverse function of with respect to .
So are as
There are also
After the above derivation about and , the kinematic reliability model of the robot becomes
Since it is very difficult to obtain the expression of the inverse function directly, and need to be simplified further.
Substituting equation (7) into equation (14), is rewritten aswhere is the abbreviation of .
Referring to the research results of literature [30], the gradient of and are collinear at some instants, as shown in Figure 3.

Hence, U is equal towhere is a constant.
We plug this into equation (19) to yield
Therefore, is given by
Equation (19) can be transformed into
On the basis of equation (23), the value of can be solved. On the other hand, the number of solutions may not be fixed due to different intricacies of trajectories. Ultimately, expansion point is represented by
In the same way, () can be obtained in line with whose expression is
And expansion point is described as the following expression:where () is the solution of equation (25).
Remark 2. Directly solving equations (19) and (25), both of which include about a thousand trigonometric functions, is difficult and takes an inordinate amount of time. An efficient way is to take points at appropriate intervals to find the closest solution among all samples.
Since the above derivation is reasonable, we use and instead of and to simplify the kinematic reliability model. There isBeside time instant and , start time and end time should be also added to improve calculation accuracy. Consequently, we recast the kinematic reliability model aswhere contains , and . is the sign function of , defined asCombining the previous descriptions, is similar to , which leads to the error distribution of being normal. In light of this property, the mean and covariance matrix can be derived to compute the reliability , whose formula iswhere is the cumulative distribution function (CDF) of multiple standard normal variables.
The formula of iswhere denotes the covariance matrix with elements, calculated by
Remark 3. It is worth noting that must be a positive definite matrix during the calculation process. This is to say, the is full-rank with a positive determinant value.
Otherwise, we must eliminate the point with the lowest failure rate one by one until the covariance matrix meets the above condition. The calculation formula for failure probability is written asImmediately after this process, the effective points are kept to reconstitute the new mean and matrix , namely,where .
Accordingly, the reliability of the robot end-effector is calculated byFinally, we summarize the flow of kinematic reliability, as shown in Figure 4.
When the reliability values for the specified tolerances are obtained, the performance of the robot can be evaluated according to the quality indicator of the industry. In general, product accuracy must meet at least Three-Sigma (reliability 99.7%) requirement for enterprises, but such a production situation is not clearly competitive due to a certain amount of profit applied to after-sales links and the low loyalty of customers. Therefore, the Six Sigma principle [31] (reliability 99.99966%) has gradually become a new quality measure since it was proposed in the 1990s. If the precision of the robot can meet such situation, it means that the product is almost impossible to defect.

5. Simulation Results
In this article, the simulation platform MATLAB 2019a is employed, and the CPU is an AMD Ryzen 7 3700X 8-Core. For example, to replicate the shipyard welding task, a 7-DOF redundant robot is chosen, with Table 2 displaying the corresponding variable distribution. Based on the robot kinematics literature, the cyclic coordinate descent (CCD) method [32] is used to calculate the corresponding joint angles of the robot, and we have and , which are the starting angle and the ending angle, respectively. Furthermore, the following motion function of the robot joint angles within 3 s is planned:where represents the value of the joint angle and represents the total time of executing the task.
Figure 5 depicts the motion trajectory of the robot end-effector in the Cartesian coordinate system from the starting point to the ending point.

In order to demonstrate the validity of the envelope method, the equivalent extreme value (EEV) method [33] and MCS are introduced for comparison, with the latter serving as a desired value benchmark with random trajectories simulated. For the MCS method, only 31 points on a trajectory with a 0.1 s interval are taken into account due to unaffordable computation of points at each time. If all of these points are within the allowable tolerance, this trajectory is regarded as valid. Otherwise, this course is unqualified. When is at 3 s, for example, failure samples and valid samples are displayed, as illustrated in Figure 6, with the comprehensive position error () of representative trajectories shown in Figure 7.


Figure 8 plots the results of three algorithms. It can be seen that the failure rate of three algorithms decreases with the increment of error. From a global perspective, however, the envelope method has higher agreement with the MCS method. Since the EEV is commonly applied in conditions where there is an extreme point, it is not suitable for the case of small fluctuations in trajectory. This is why the accuracy of EEV is unstable for kinematic reliability.

(a)

(b)

(c)
According to the trajectory given above, taking the direction as an example, the time is obtained by equation (25) as 0 s, 1.62 s, 3 s, and the covariance matrix can be further obtained aswhich is positive definite.
Therefore, reliability at different errors in the y direction can be calculated by equation (36). If , thenwhere .
The specific values at the error of 0.6 mm, 0.8 mm, and 0.1 mm are counted in Table 3. When the error is 0.6 mm in x-, y-, and z-directions, the relative errors of the envelope method are 5.8%, 4.1%, and 2.9% compared with the MCS method, and the counterparts of the EEV method are 43.8%, 24.9%, and 46.4%, respectively.
Figure 9 provides the difference in the failure rate of the robot end-effector between the envelope method and the MCS method. As the position error increases, the difference between them becomes smaller, with 90 percent of the data concentrated in the region. In conjunction with the preceding analysis, the envelope method is adequate for the kinematic reliability of the redundant robot in a specific direction.

Three algorithms are again introduced to compare their characteristics for the kinematic reliability of the attitude angles, as shown in Figure 10. The envelope method produces results that are slightly closer to those of the MCS method than to those of the EEV method. The reason for this phenomenon is that the deviation of a specific point, which has a large impact on the accuracy of the trajectory of attitude angles, is much larger than the deviation of other points on the trajectory.

(a)

(b)

(c)
According to the data in Table 4, which contains failure rate values for various errors, the envelope method is still superior to the EEV method. In addition, the error of the attitude angle can be ignored, which is too small relative to the other two angles. When the error is , the failure rates of the three attitude angles obtained by the envelope method are 4.4% and 5.5% relative to the MCS method, better than the result of the EEV. This further shows that the envelope method has high robustness.
On the other hand, we note that the EEV method performs better when calculating the reliability of the attitude angle rather than the position in light of Figures 8 and 10. The summary is important for two reasons: (1) when computing kinematic reliability of attitude angles, the correlation between points at various times is higher, resulting in the accuracy of a single point being similar to the total accuracy. (2) The error in attitude angles of the trajectory is highly volatile, with a certain point having a significant impact on the trajectory.
Figure 11 shows the difference in the failure rate of attitude angles between the envelope method and the MCS method from to . The curves in the figure show that the difference between the envelope method and the MCS method is less than and approaches 0 with the increment of error. This also proves that the envelope method has stronger stability and a wider range of applications.

Another key factor to consider when deciding whether or not a method may be implemented is its efficiency. The envelope method requires 5502 s to calculate, which is slower than the 2038 s from the EEV method but much faster than the 36050 s simulated by the MCS method. Although the calculation efficiency of the envelope method is not as good as that of the EEV method, its accuracy is higher. Therefore, it is more appropriate to use the envelope method to evaluate the kinematic reliability of the robot while taking into account the computation efficiency and accuracy. Taking advantage of the enveloped method, the tolerances of -, -, and -directions are 2 mm, 2.3 mm, and 2 mm, respectively. In the same way, the tolerances of -, -, and -attitude angles can be obtained as , , and , respectively.
The comprehensive error function of the robot does not usually follow a normal distribution, which is hard to analyze directly through the envelop method. As a result, in this article, the time expansion points in three directions re-approximated as those of the comprehensive error function since the error of three attitude angles is smaller to be ignored based on Tables 3 and 4. The reliability data collected are depicted in Figure 12, which is calculated approximately by I-PCM (Improved Product of Conditional Marginal) [34] for multidimensional normal integral. When the error is 0.6∼1.5 mm, the error of the envelope method relative to the MCS method is 7%∼68%, which is far inferior to 46%∼76% yielded by the EEV method, indicating the envelope method has better performance. Although the calculated result of the envelope method may not reach the ideal state, it obtains a more accurate value close to MCS, and the maximum difference between them is . On the whole, it still proves to be a good solution for the kinematic reliability of the robot.

6. Conclusion
The kinematic reliability of a robot is an important indicator for evaluating its performance, and it is inextricably linked to product quality. Although the MCS method is commonly applied to reliability analysis, it often comes with a huge computational cost, which is unbearable in reality. For a 7-DOF redundant robot, this work constructs the error model of the robot trajectory using probabilistic methods and determines the kinematic reliability of the trajectory using the envelope method. Then, the reliability value can be obtained by solving the equation and calculating the multidimensional normal integral, instead of a large number of stochastic simulations.
The simulation example shows that the envelope method is more accurate than the EEV method and closer to the MSC method, with relative errors between the envelope method and MCS of less than 6% in all orientations and attitude angles. Furthermore, the difference between the envelope method and MCS is no more than , implying that the envelope method is precise. On the other hand, in terms of efficiency, the envelope method, EEV method, and MCS method take 5502 s, 2038 s, and 36050 s, respectively, to calculate kinematic reliability. Considering the synthesizing capacity of the solution accuracy and calculation efficiency, the envelope method is the best method for evaluating the kinematic reliability of the robot.
Furthermore, we do not demonstrate whether this method can be used to calculate the reliability of other error models with nonnormal distributions. In the future, we plan to apply this method to implement other tasks, such as spraying and assembly performed by various robots, in order to fully exploit the potential of the proposed method.
Data Availability
The parameters of the aerial robot we used have been shown in the article. The data of simulation and experiments can be provided if necessary.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Authors’ Contributions
Jiahui Gu, Ziyi Li, Shaopeng Kang, and Rui Ma contributed equally to this work.
Acknowledgments
This work was partially supported by the National Natural Science Foundation of China (52005231), Social Development Science and Technology Support Project of Changzhou (Grant nos. CE20215050 and CE20209002), and Natural Science Research Foundation of Jiangsu Universities (Grant no. 20KJB580005).