Research Article  Open Access
Xiaoyan Chen, Qiuju Zhang, Yilin Sun, "Evolutionary Robot Calibration and Nonlinear Compensation Methodology Based on GADNN and an Extra Compliance Error Model", Mathematical Problems in Engineering, vol. 2020, Article ID 3981081, 15 pages, 2020. https://doi.org/10.1155/2020/3981081
Evolutionary Robot Calibration and Nonlinear Compensation Methodology Based on GADNN and an Extra Compliance Error Model
Abstract
This study addresses the problem of nonlinear error predictive compensation to achieve high positioning accuracy for advanced industrial applications. An improved calibration method based on the generalisation performance evaluation is proposed to enhance the stability and accuracy of robot calibration. With the development of technology, a deep neural network (DNN) optimised by a genetic algorithm (GA) is applied to predict the nonlinear error of the calibrated robot. To address the change of external payload, an extra compliance error model is established with a linear piecewise method. A global compensation method combining the GADNN nonlinear regression prediction model and the compliance error model is then proposed to achieve the robot’s highprecision positioning performance under any external payload. Experimental results obtained on a Staubli RX160L robot with a FARO laser tracker are introduced to demonstrate the effectiveness and benefits of our proposed methodology. The enhanced positioning accuracy can reach 0.22 mm with 98% probability (i.e., the maximum positioning error in all test data).
1. Introduction
With the emergence and subsequent development of digital processing and intelligent manufacturing, the breadth and depth of industrial robotic applications (i.e., application fields and task requirements) are increasing [1]. To meet the highprecision requirements of complex tasks, the absolute positioning accuracy of robots is particularly important [2‒4]. Robot calibration and error compensation are two important ways to improve the absolute positioning accuracy of robots, which have been studied by several scholars [5‒8].
Robot calibration can be divided into robot kinematic calibration and nonkinematic calibration [9‒11]. The latter mainly includes the identification of structural robot deformation errors caused by the change of temperature [12], the link selfgravity of robots, and the change of external payload [13‒18], as well as the calibration of dynamic parameters [19‒22]. All the above calibration methods need the error model, which is based on the structural parameters of robots. Such a model contains a complex modelling process and identification steps. However, the accuracy of the error model is still insufficient for nonlinear errors. Hence, modelfree approaches have been developed for robot error compensation. In [23], a fuzzy interpolation compensation method is presented, and its effectiveness is confirmed by the simulation results. Zeng et al. [24] proposed an error compensation method with error similarity analysis based on the semivariogram function to improve the absolute positional accuracy of industrial robots. Other works [25, 26] used an inversedistance weight (IDW) interpolation model to reduce absolute position errors. The common disadvantage of the above interpolation methods is that the robot endeffector should be at a constant orientation relative to the reference frame. Undoubtedly, the interpolation compensation method has its disadvantages and limitations.
Guo et al. [9] presented a multilevel calibration technique for an industrial robot, which includes parallelogram structural error compensation, kinematic calibration, and nongeometric error compensation. In the last level of this method, contrary to the similar works mentioned above, the robot joint space (instead of the Cartesian space) is divided into a sequence of fanshaped cells, after which the inversedistance weight method is used to compensate the nongeometric errors of the robot. Although this method enlarges the space range of the robot compensation (divided by the motion range of the first three joints), it is still unable to compensate the error of the robot’s orientation and position because the inverse solution of the same position corresponds to multiple sets of joint values. To achieve global compensation (at any orientation and position), an artificial neural network (ANN) based nonlinear regression prediction method was proposed to facilitate learning ability and high adaptation [27, 28]. A radial basis function network (RBFN) algorithm was applied to interpolate the error workspace of the industrial robot to get continuous functions, and the results proved that the proposed scheme can improve the accuracy of the robot [27]. Nguyen et al. [28] proposed an effective solution for all robot error sources by combining the modelbased identification of the robot’s geometric errors via an extended Kalman filter and a compensation technique using the ANN. However, the ANNbased nongeometric correction method in [28] was based on a twolayer backpropagation neural network (BPNN) with one hidden layer and did not take into account the change of external payload.
As we have shown, modelbased kinematic calibration or nonkinematic calibration is essential for uncalibrated robots (i.e., the uncalibrated robots fresh out of a factory or the robots with numerous errors after being used for a certain period of time). Besides, the nonlinear errors including the intrinsic errors (i.e., elastic deformation of the link and joint caused by the link selfgravity of the robot, gear backlash, hysteresis, and friction) and the extrinsic alterations (i.e., the variable external payloads) of the robot cannot be eliminated by modelbased robot calibration. In the current study, the predictive compensation method for the nonlinear errors of robots, especially the residual errors of calibrated robots, is yet to be fully studied in the literature. Therefore, the goal of this work is the enhancement of calibration and error compensation techniques for achieving the robot’s highprecision positioning performance under any external payload. Three key problems that should be considered are as follows:(1)There is no unified uniform standard for the determination of the parameters in the calibration model, so it is of great significance to propose a parameter optimisation method for the calibration quality(2)Given the aforementioned nonlinear errors, a global compensation method should be presented to predict effectively and efficiently these positioning errors(3)Due to the large range of external payload variation of the heavy robot, the nonlinear deformation error caused by the change of the external payload should be solved
In this work, an improved calibration method, which takes into account the generalisation performance and robustness of geometric parameter correction, is introduced to enhance the essential positioning accuracy of robots. DNN has received substantial attention in both the signal processing field and the machine learning field with its strong regression capabilities [29‒31]. The design of DNN architecture must be optimised to make the DNN demonstrate the best predictive capacity. The critical factors that determine the performance of DNN prediction are as follows: the number of hidden layers, the number of nodes per layer, the training data, the training method, and the validation data in the training process. GA is an automatic optimisation method based on the principles of evolution postulated with the natural selection and the adaptation of individuals to the environment, which can offer an optimal solution for the DNN design [32]. Thus, a GADNN prediction model is proposed to achieve the optimal nonlinear error prediction compensation and to further enhance the positioning accuracy of the robot. Moreover, as the training data for the DNN prediction model are measured under a constant endeffector of the robot, the model does not consider the effect of the change of external payload on the positioning errors. According to the relation between force and torque, an extra compliance error model for predicting joint deformation caused by any external payload is established in this work. Thus, a multilevel accuracy improvement method based on the GADNN model and the extra compliance model can be obtained to realise robot optimisation calibration and global error compensation.
To address the aforementioned problems and the corresponding solutions, the rest of this paper is organized as follows: Section 2 proposes an improved calibration method for serial robots. Section 3 develops a GADNN regression model for nonlinear error prediction. To consider the effect of the change of external payload, an extra compliance error model based on torque division is described in Section 4. In Section 5, several experiments are implemented, and the results are presented for the Staubli RX160L robot. Finally, conclusions are drawn in Section 6.
2. Parameter Optimisation Method for the Calibration Model
In our preliminary work [11], the MDH model with five parameters is used for the robot kinematic modelling, which can describe all geometric errors continuously. The Staubli RX160L robot is a 6DOF serial robot with six rotational joints and can be regarded as a typical representative of industrial robots (see Figure 1(a)), and the kinematic model of this robot is shown in Figure 1(b), which defines the seven link frames associated with the robot.
(a)
(b)
Because the actual robot base frame (F_{0}) cannot be touched, the small transformation matrix ^{B}T_{0} is introduced to present the error vector ∆σ between the constructed robot base frame (F_{B}) and the actual robot base frame (F_{0}):
Then, the error equation is expressed as follows:
Hence, the robot kinematic position error model is established as follows:where E_{k} represents the position errors; H_{k} is the geometric identification matrix; λ_{k} denotes the geometric error parameters; and M_{θ}, M_{d}, M_{a}, M_{α}, M_{β}, and M_{σ} are the coefficient matrices corresponding to the error parameters.
By using the kinematic error model, all defined geometric parameters can be identified theoretically. However, the model is not the best choice to correct all geometric parameters in the process of robot calibration. Generalisation performance, also known as scalability, refers to the performance of the classifiers trained on a certain training set to adapt to the data outside the training set. In this work, a Kfold crossvalidation method is applied to evaluate the generalisation performance of the robot calibration model with different geometric error parameters. Firstly, the minimum number of parameters of the calibration model (i.e., several important joint angle parameters, setover parameters, and link length parameters) is determined. Taking the Staubli RX160L robot as an example, the minimum number is 9 (i.e., ∆θ_{2}–∆θ_{5}, ∆d_{3}, ∆d_{4}, and ∆a_{1}–∆a_{3}), and the maximum number is 24 (each link frame contains 4 parameters). Therefore, according to the limited range, different calibration models can be formed by the arbitrary combination of different geometric error parameters; from this, the number of calibration models can be calculated:where n is the difference value between the maximum and minimum number of geometric error parameters in the calibration model.
Then, the training data used for calibration model identification are divided into K parts, one of which is taken as the test set without repetition; the other K1 parts are used as the training set for the training calibration model. The mean square error (MSE) of the prediction accuracy of each training model can be calculated, and the final value EV_{K} obtained by averaging K times is the evaluation value of the generalisation performance. Hence, each calibration model should have a performance evaluation value from which the calibration model corresponding to the minimum value is the optimised calibration model of the robot:
It can be seen that the optimised calibration model is obtained by using these various calibration models with different geometric error parameters and the generalisation performance evaluation method. Moreover, the improved calibration method based on the parameter optimisation method can also be combined with engineering practice, thus increasing or reducing the effective constraints on the range of parameter correction so as to obtain the required calibration model. Hence, the improved calibration method proposed in this section can enhance the stability and accuracy of robot modelbased calibration. Note that, however, some errors are unable to be established by the robot kinematic position error model, which also affect the quality of the absolute positioning accuracy of robots. To tackle this problem, the nonlinear errors of the calibrated robot will be discussed later.
3. DNNBased Regression Prediction Method for the Nonlinear Errors
Given the complexity of the sources of the robot’s residual errors (i.e., elastic deformation of links and joints caused by the link selfgravity of the robot and the change of external payload, gear backlash, hysteresis, and friction) after modelbased calibration, these residual errors generally have strong nonlinearity. At present, most of the learning methods, such as regression, are shallow structure algorithms [8, 27, 28, 33]. These methods are limited in terms of their ability to represent complex functions in the case of limited samples and computational units, and their generalisation ability for strong nonlinear errors is restricted to a certain extent.
DNN can achieve complex function approximation by learning a deep nonlinear network structure (5, 6, or even 10 hidden layers are usually present) and show a strong capacity to learn the effective features of data from a small number of samples. Two types of joints of industrial robots exist, namely, rotational joints and prismatic joints, which correspond to the joint variables θ_{i} and d_{i}, respectively. In this work, all of the robot’s joint variables are regarded as feature objects, i.e., the input data set of the DNN. For a serial robot with six rotational joints (i.e., joints θ_{1}–θ_{6}), the input layer of the DNN has six neuron nodes.
To realise the prediction and compensation of the robot positioning error, the output of the DNN is the predicted position error with a threedimensional (3D) vector E_{r} (E_{r} = [∆x, ∆y, ∆z,]). Moreover, the original data acquisition used for DNN model training is vital in realising the exact regression prediction of the model, which plays a decisive role in the predictive capacity of the model. The measurement configurations utilised in this work for 6DOF series robots, in accordance with the distribution characteristics of the robot’s nonlinear errors and the observability index O_{1} [11], are shown in Figure 2. Using the optimised configurations with 46 groups of different orientations (in each configuration, a set of position points with uniform grid distribution in the whole workspace of the robot) can help achieve the robot’s overall acquisition of the nonlinear errors.
The MSE is selected as the loss function to evaluate the degree of inconsistency between the predicted value and the real value of the DNN regression prediction model, that is, the criterion to measure the predictive capacity of the model. As the DNN contains more hidden layers and adopts the backpropagation algorithm, the traditional activation functions (i.e., sigmoid function and tansigmoid function) can lead to gradient vanishing. To tackle this problem, the parametric rectified linear unit (PReLU) function [34] is applied in the first few layers of the DNN as shown in Figure 3(a), whereas the tansigmoid function is used in the latter layers of the DNN as shown in Figure 3(b).
(a)
(b)
Finally, the GA is applied to optimise the number of hidden layers, the number of nodes per layer, the learning rate, the training data, the training method (including the specific distribution of the activation functions), and the validation data in the training process for building the optimal DNN architecture. The loss function (MSE) of the DNN model is treated as the fitness function for the GA, and the size of the initial population of the GA is 100. Furthermore, the crossover operationmutation operation has 90% probability, and the mutation rate is 0.05. Hence, the GADNN nonlinear regression prediction method is developed as an effective solution for the highperformance compensation of the robot’s nonlinear errors, as shown in Figure 4.
It is important to note that the training data set for the GADNN model is always measured with a constant endeffector of the calibrated robot, and the nonlinear regression prediction model does not consider the effect of the change of external payload on the positioning error. To overcome this difficulty, the robot’s positioning errors under variable external payloads will be analysed in detail in the next section.
4. Compliance Error Model with the Variable External Payloads
According to the relation between force and torque, the formula of torque for each joint can be derived based on the robot’s kinematic model as follows:where J_{θ} represents the kinematic Jacobian matrix; equation (7) describes the transformation of the force between coordinate systems; ^{0}F_{m} is the force vector of the external payload; ^{0}R_{6} is the rotation matrix between the actual robot base frame (F_{0}) and the flange frame (F_{6}); ^{6}P_{m} is the vector between the centroid of the external payload and the origin of the flange frame (F_{6}); and E_{3} is a 3 × 3 identity matrix.
Ψ(P) is the skew symmetric matrix for a threedimensional vector , which can be expressed as follows:
To analyse the compliance error caused by the external payload, the robot joint is usually assumed to be a linear torsional spring [16]. However, Pashkevich et al. [13] demonstrated that the torquetorsion relation was nonlinear in the robot’s elastic behaviours under loading. Based on the material in [35], the relation between nonlinear torque and torsion can be defined as a thirdorder polynomial:
The thirdorder polynomial satisfies two conditions in equation (10), so the polynomial is monotonously increasing, i.e., for any torque, only one torsion value δθ_{i} is present. Then, combined with the torque obtained from equation (6), a linear piecewise method based on torque division is used to present the joint stiffness, and the simplified method of the two subsections is shown in Figure 5.
Based on equation (6), the range of the torque corresponding to each joint can be calculated under any external payload and the whole workspace of the robot. A simulation study on the Staubli RX160L robot is implemented to illustrate the linear piecewise method. Given that the nominal load of the robot is 14 kg, two external payload values of 0.4 kg and 14 kg are selected in this simulation experiment. Combined with the 46 configurations in Figure 2, 2,000 measurement points are selected in the whole workspace of the robot. The range of the torque corresponds to each joint (joints 2–6 of the Staubli RX160L robot, the axis of joint 1 is parallel to the direction of the force, so the torque of joint 1 can be ignored), as shown in Figure 6.
(a)
(b)
As can be observed from Figure 6, joint 6 (i.e., j6 in Figure 6) is subjected to a very small torque (approximate to zero), so its deformation effect can be neglected. The torque ranges of joints 2 and 3 are wider, and these torques are much greater than those of other joints. We can see that the fluctuation range of torque between joint 4 and joint 5 is similar. For general industrial robots, the stiffness of joint 4 is much larger than that of joint 5, so the stiffness of joint 4 can be set as a constant. For joints 2, 3, and 5, the corresponding stiffness can be divided into two or more piecewise values according to the magnitude of the torque.
Based on the linear piecewise method with the two subsections, the position errors of the origin of the tool frame (F_{T}) for 6DOF series robots can be given as follows:where ^{6}P_{T} is the vector between the origin of the robot tool frame (F_{T}) and the origin of the flange frame (F_{6}).
To obtain the best piecewise critical torque values (i.e.,,, and ) of the three joints, we use the positiondistance method based on torque difference to obtain the additional position errors after the change of external payload. The variation range of the external payload is determined according to the minimum load and the nominal load of the robot, the m external payload values of different sizes are selected in this range, and the corresponding positiondistance under each of the two loads can be obtained according to the torque difference of these two loads. Next, we adjust the critical values to minimise the residual error between the predictive positiondistance and the actual positiondistance. As such, the best stiffness value of each joint is obtained as follows:
Hence, the extra compliance error model (i.e., equation (13)) based on torque division is established to predict the additional position errors caused by the change of external payload. In summary, a global compensation method based on the combination of the GADNN nonlinear regression prediction model and the extra compliance error model is applied to correct the nonlinear errors of the calibrated robot for achieving highprecision positioning performance under any external payload. The global compensation model can be given as follows:where ∆ X_{(GA−DNN)} is the predictive position error by the GADNN nonlinear regression prediction model.
5. Experiments and Result Analysis
Figure 7 shows the experimental setup of this work, which consists of a FARO laser tracker with a resolution of 0.5 μm and an accuracy of 8 μm + 0.4 μm/m, a PC, and a 6DOF Staubli RX160L robot with a CS8 controller.
5.1. Performance of the Improved Calibration Method
To collect the robot measurement positions for the geometric parameter identification, 90 measurement positions are selected by the intelligent selection algorithm of optimal measurement poses [11]. In this experiment, the threefold crossvalidation method is applied to evaluate the generalisation performance of the robot calibration model with different geometric error parameters. Based on the performance evaluation values, an optimal calibration model for this robot can be obtained. The nominal value of MDH parameters and the identified geometric parameters of the optimal model are shown in Table 1.
 
“Nom.”: nominal value; “Iden.”: identified value; “—”: does not exist. 
To present the advantages of the proposed improved calibration method (Section 2), the results of the validation are compared with those obtained by the conventional kinematic calibration method (method presented in [28, 35]), as shown in Table 2. According to the calibration results, we can see that the number of parameters used in the optimal model and the conventional calibration model is 17 (i.e., ∆θ_{1}–∆θ_{5}, ∆d_{1}, ∆d_{3}–∆d_{6}, ∆a_{1}–∆a_{6}, ∆α_{2}, and ∆β_{2}) and 23 (i.e., ∆θ_{1}–∆θ_{6}, ∆d_{1}, ∆d_{3}–∆d_{6}, ∆a_{1}–∆a_{6}, ∆α_{1}–∆α_{5}, and ∆β_{2}), respectively.

Furthermore, 40 random positions are chosen as the test data in the overall workspace of the robot, and the enhancement results of position accuracy are shown in Figure 8. As can be seen, both the conventional calibration method and the proposed calibration method have evident effects on improving positioning accuracy. Notably, the residual position errors after calibration using the proposed calibration model are smaller than those after calibration using the conventional calibration model under the same measurement environment. The average residual position error after the improved calibration is reduced significantly to 0.29 mm from 4.51 mm (before calibration), and the maximum position error decreased substantially from 7.69 mm to 0.65 mm. Consequently, the proposed calibration method has higher identification stability and can offer a considerable improvement in the robot positioning accuracy.
(a)
(b)
(c)
5.2. Accuracy Enhancement by the GADNN Prediction Method
Due to the influence of nonlinear errors, the calibrated robot still has large residual errors which cannot meet the requirements in highprecision field applications. The Staubli RX160L robot has a general structure of a 6DOF serial robot, which is often used in highprecision processing, and its repeatability is about ±0.05mm. Hence, the repetitive accuracy is the limit value of the absolute positioning accuracy of the robot and also the target value of error compensation using our proposed approach.
Based on the 46 configurations of different orientations (Figure 2), 50 to 100 endpoints with approximately uniform spatial distribution in the overall workspace of the robot for each configuration are selected, and a total of 2,796 measurement positions are measured for DNN model training. In addition, 708 verification positions with 11 groups of arbitrary different orientations are also selected in the overall workspace and used as the test data. Consequently, a total of 3,504 endpoints are used in the experiment of nonlinear error compensation, as shown in Figure 9.
(a)
(b)
To better present the nonlinear characteristics of the error data, four sets of data are selected from 2,796 measurement positions, each with the same configuration and the same position in the Z direction, and the corresponding position errors (i.e., d in Figure 10) of the four groups of data are calculated. According to the X and Y coordinate data in these four sets of data, the corresponding position errors are displayed on the XY plane, as shown in Figure 10. We can see that the errors of the robot’s different configurations vary greatly in the approximately same position. Hence, the results illustrate the necessity of nonlinear regression prediction in achieving the residual error compensation of the calibrated robots.
(a)
(b)
(c)
(d)
To demonstrate the superiority of the DNN over the shallow neural network (SNN), the performance of the proposed GADNN model is compared with that of a 2layer BPNN model and another 5layer BPNN model. Gradient descent (GD) is one of the most popular algorithms with which to perform optimisation and by far the most common way to optimise neural networks. To optimise neural networks further, certain excellent gradient descent optimisation algorithms are proposed and applied to neural network training, i.e., momentum, Nesterov accelerated gradient (NAG), Adagrad, Adadelta, RMSprop, and adaptive moment estimation (Adam). Kingma and Ba [36] showed that Adam might be the overall best choice. Hence, the gradient descent algorithm and Adam algorithm are used in the comparative experiment.
In our approach, the GA is applied to select the most suitable gradient descent optimisation algorithm for training the DNN model. Based on the collected training data, the optimised GADNN model has a network architecture with 11 hidden layers, and the corresponding optimised node numbers are listed as follows: 10, 10, 13, 13, 16, 16, 15, 13, 11, 9, and 7. The results of the iteration, which we obtained by using the joint data and error data of the 2,796 measurement positions, are shown in Figure 11. Figure 11 shows that the prediction accuracy of the model greatly improves with the increase of the depth of the neural network, and the GADNN prediction method greatly accelerates the convergence speed of the algorithm. The advantages and necessity of the proposed method can be explained in terms of the iteration results.
The predictive compensation and validation results of the proposed GADNN nonlinear regression prediction method are compared with those obtained by the popular ANNbased nongeometric correction method from [28], the GRNN method from [33], and the RBFN method from [27]. These three methods are commonly used in nonlinear regression prediction. The results obtained from these four prediction methods are compared in Figure 12, which show that the RBFN method and the proposed GADNN method have better predictive ability than the other methods. Next, 708 test points with different positions and orientations are randomly selected in the robot workspace, and two prediction models established by the RBFN method and the proposed GADNN method are used to predict the position errors of the test data, respectively. The compensation results for the test data are shown in Figure 13; it can be seen that the proposed GADNN method has evident advantages. The average residual position error is reduced significantly from 0.52 mm (before compensation) to 0.07 mm using the proposed predictive model, and the positioning accuracy almost reaches the repetition accuracy of the robot.
5.3. Effectiveness of the Global Compensation Method
The 3,504 endpoints used in the above experiments are measured under the same endeffector (the mass is 3.6 kg), so there are no positioning errors caused by the change of external payload. To demonstrate the practical effectiveness and performance of the proposed global compensation method and to confirm the benefits of the method for the change of external payload, this section addresses the compliance error model of the Staubli RX160L robot with the change of external payload.
A measurement tool is designed for changing the external payload, which can fix different sizes and number of weight blocks and be used as the endeffector of the robot, as shown in Figure 14. In this experiment, six external payload values of different sizes (i.e., 0.4, 1.95, 3.95, 5.95, 7.95, and 9.95 kg) are selected, and 50 measurement positions are selected and measured under each load. Using the positiondistance method (Section 4), the best piecewise critical torque values (i.e., , , and ) of joints 2, 3, and 5 can be obtained, and the stiffness of each joint of the robot can be identified (Table 3). Then, the compliance error model caused by the change of external payload can be established.
(a)
(b)
 
“—”: does not exist. 
To verify the accuracy of the proposed global compensation method, two groups of different orientations in 708 verification positions, which include 168 test points, are selected. Using the proposed global compensation method, which combines the compliance error model and the GADNN nonlinear regression prediction model, we can predict and compensate the position errors of these test points with 5.95 and 9.95 kg of the external payloads, respectively.
Figure 15 depicts the residual position error of the test points after compensation. The positioning accuracy after the compensation using the GADNN prediction method is poor because this method does not consider the additional compliance errors caused by the change of the external payload. Furthermore, the robot’s positioning accuracy using the proposed global compensation method is improved with an average rate of 63% compared with that obtained using the GADNN model. According to the positioning accuracy of all test data, the enhanced positioning accuracy (i.e., the maximum positioning error) can reach 0.22 mm with 98% probability by using the proposed global compensation method.
(a)
(b)
6. Conclusions
In this paper, a multilevel accuracy improvement method is proposed, which includes the improved kinematic calibration, the GADNN nonlinear regression prediction, and the extra compliance error compensation. To enhance the stability and accuracy of robot calibration, a parameter optimisation method based on the Kfold crossvalidation and generalisation performance evaluation is applied. Given that the calibrated robot’s residual errors (i.e., elastic deformation of the links and joints caused by the robot’s link selfgravity and the change of the external payload, gear backlash, hysteresis, and friction) generally have strong nonlinearity, a DNN nonlinear regression prediction method based on the GA optimisation of the architecture (i.e., the number of hidden layers, the number of nodes per layer, the learning rate, the training data, the training method, and the validation data in the training process) is developed as an effective solution for achieving the highperformance compensation of the robot’s nonlinear errors. Furthermore, the robot positioning accuracy is further improved due to the use of the extra compliance error model to compensate for the compliance error caused by the change of the external payload.
The proposed theoretical methods have been validated via an experimental study involving a kinematic calibration and error compensation of a Staubli RX160L industrial robot. The experimental results using the GADNN prediction method are more accurate than the previous methods (i.e., methods from [27, 28, 33]), with the positioning accuracy of not more than 0.21 mm. The effectiveness and benefits of the global compensation method based on the GADNN prediction model and the extra compliance error model are demonstrated by the compensation results of the tests with different external payloads. The position accuracy after compensation using the proposed global compensation method improved with an average rate of 63% compared with that obtained using the GADNN method without considering the additional compliance error.
In the experimental study, data acquisition takes a long time in robot calibration and compensation, especially for the training of the DNN model (2,796 measurement points are used in this paper). Admittedly, the greater the amount of training data is, then the more accurate the DNN model becomes. Given that the large amount of data is needed for each robot for nonlinear compensation, the efficiency becomes relatively low. Consequently, we will focus on the application of the transfer learning method in this DNNbased prediction model, which only needs to use a small amount of data to share a similar DNN model for the same type of the robot.
Data Availability
No data were used to support this study.
Disclosure
Yilin Sun received his B.S. degree in School of Mechanical Engineering from Jiangnan University, Wuxi, China, in 2013. He is currently a Ph.D. candidate of Jiangnan University. His research is in the development of bionic grippers, email: 15261596813@163.com.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research was supported by the National Natural Science Foundation of China (no. 51575236) and Postgraduate Research and Practice Innovation Program of Jiangsu Province (no. KYCX18_1840).
References
 D.H. Kim, T. J. Y. Kim, X. Wang et al., “Smart machining process using machine learning: a review and perspective on machining industry,” International Journal of Precision Engineering and ManufacturingGreen Technology, vol. 5, no. 4, pp. 555–568, 2018. View at: Publisher Site  Google Scholar
 L. Wu, X. Yang, K. Chen, and H. Ren, “A minimal POEbased model for robotic kinematic calibration with only position measurements,” IEEE Transactions on Automation Science and Engineering, vol. 12, no. 2, pp. 758–763, 2015. View at: Publisher Site  Google Scholar
 Y. Qiao, Y. Chen, B. Chen, and J. Xie, “A novel calibration method for multirobots system utilizing calibration model without nominal kinematic parameters,” Precision Engineering, vol. 50, pp. 211–221, 2017. View at: Publisher Site  Google Scholar
 A. Olarra, D. Axinte, and G. Kortaberria, “Geometrical calibration and uncertainty estimation methodology for a novel selfpropelled miniature robotic machine tool,” Robotics and ComputerIntegrated Manufacturing, vol. 49, pp. 204–214, 2018. View at: Publisher Site  Google Scholar
 Y. Jiang, T. Li, L. Wang, and F. Chen, “Kinematic error modeling and identification of the overconstrained parallel kinematic machine,” Robotics and ComputerIntegrated Manufacturing, vol. 49, pp. 105–119, 2018. View at: Publisher Site  Google Scholar
 G. Palmieri, M.C. Palpacelli, L. Carbonari, and M. Callegari, “Visionbased kinematic calibration of a smallscale spherical parallel kinematic machine,” Robotics and ComputerIntegrated Manufacturing, vol. 49, pp. 162–169, 2018. View at: Publisher Site  Google Scholar
 L. Kong, G. Chen, Z. Zhang, and H. Wang, “Kinematic calibration and investigation of the influence of universal joint errors on accuracy improvement for a 3DOF parallel manipulator,” Robotics and ComputerIntegrated Manufacturing, vol. 49, pp. 388–397, 2018. View at: Publisher Site  Google Scholar
 D. Wang, Y. Bai, and J. Zhao, “Robot manipulator calibration using neural network and a camerabased measurement system,” Transactions of the Institute of Measurement and Control, vol. 34, no. 1, pp. 105–121, 2012. View at: Publisher Site  Google Scholar
 Y. Guo, S. Yin, Y. Ren, J. Zhu, S. Yang, and S. Ye, “A multilevel calibration technique for an industrial robot with parallelogram mechanism,” Precision Engineering, vol. 40, pp. 261–272, 2015. View at: Publisher Site  Google Scholar
 A. Joubair and I. A. Bonev, “Nonkinematic calibration of a sixaxis serial robot using planar constraints,” Precision Engineering, vol. 40, pp. 325–333, 2015. View at: Publisher Site  Google Scholar
 X. Chen, Q. Zhang, and Y. Sun, “Nonkinematic calibration of industrial robots using a rigidflexible coupling error model and a full pose measurement method,” Robotics and ComputerIntegrated Manufacturing, vol. 57, pp. 46–58, 2019. View at: Publisher Site  Google Scholar
 R. Li and Y. Zhao, “Dynamic error compensation for industrial robot based on thermal effect model,” Measurement, vol. 88, pp. 113–120, 2016. View at: Publisher Site  Google Scholar
 A. Pashkevich, A. Klimchik, and D. Chablat, “Enhanced stiffness modeling of manipulators with passive joints,” Mechanism and Machine Theory, vol. 46, no. 5, pp. 662–679, 2011. View at: Publisher Site  Google Scholar
 C. Dumas, S. Caro, M. Cherif, S. Garnier, and B. Furet, “Joint stiffness identification of industrial serial robots,” Robotica, vol. 30, no. 4, pp. 649–659, 2012. View at: Publisher Site  Google Scholar
 A. Klimchik, D. Chablat, and A. Pashkevich, “Stiffness modeling for perfect and nonperfect parallel manipulators under internal and external loadings,” Mechanism and Machine Theory, vol. 79, pp. 1–28, 2014. View at: Publisher Site  Google Scholar
 R. M. Colorado and J. M. Valenzuela, “Experimental parameter identification of flexible joint robot manipulators,” Robotica, vol. 36, no. 3, pp. 313–332, 2018. View at: Google Scholar
 P. Xu, B. Li, C.F. Cheung, and J.F. Zhang, “Stiffness modeling and optimization of a 3DOF parallel robot in a serialparallel polishing machine,” International Journal of Precision Engineering and Manufacturing, vol. 18, no. 4, pp. 497–507, 2017. View at: Publisher Site  Google Scholar
 L. Ma, P. Bazzoli, P. M. Sammons, R. G. Landers, and D. A. Bristow, “Modeling and calibration of highorder jointdependent kinematic errors for industrial robots,” Robotics and ComputerIntegrated Manufacturing, vol. 50, pp. 153–167, 2018. View at: Publisher Site  Google Scholar
 Z. Bingül and O. Karahan, “Dynamic identification of Staubli RX60 robot using PSO and LS methods,” Expert Systems with Applications, vol. 38, no. 4, pp. 4136–4149, 2011. View at: Publisher Site  Google Scholar
 C. Urrea and J. Pascal, “Parameter identification methods for real redundant manipulators,” Journal of Applied Research and Technology, vol. 15, no. 4, pp. 320–331, 2017. View at: Publisher Site  Google Scholar
 A. Bahloul, S. Tliba, and Y. Chitour, “Dynamic parameters identification of an industrial robot with and without payload,” IFACPapersOnLine, vol. 51, no. 15, pp. 443–448, 2018. View at: Publisher Site  Google Scholar
 C. Urrea and J. Pascal, “Design, simulation, comparison and evaluation of parameter identification methods for an industrial robot,” Computers & Electrical Engineering, vol. 67, pp. 791–806, 2018. View at: Publisher Site  Google Scholar
 Y. Bai, “On the comparison of modelbased and modeless robotic calibration based on a fuzzy interpolation method,” The International Journal of Advanced Manufacturing Technology, vol. 31, no. 1112, pp. 1243–1250, 2007. View at: Publisher Site  Google Scholar
 Y. Zeng, W. Tian, and W. Liao, “Positional error similarity analysis for error compensation of industrial robots,” Robotics and ComputerIntegrated Manufacturing, vol. 42, pp. 113–120, 2016. View at: Publisher Site  Google Scholar
 W. Tian, Y. Zeng, W. Zhou, and W. Liao, “Calibration of robotic drilling systems with a moving rail,” Chinese Journal of Aeronautics, vol. 27, no. 6, pp. 1598–1604, 2014. View at: Publisher Site  Google Scholar
 W. Tian, D. Mei, P. Li, Y. Zeng, P. Hong, and W. Zhou, “Determination of optimal samples for robot calibration based on error similarity,” Chinese Journal of Aeronautics, vol. 28, no. 3, pp. 946–953, 2015. View at: Publisher Site  Google Scholar
 J. H. Jang, S. H. Kim, and Y. K. Kwak, “Calibration of geometric and nongeometric errors of an industrial robot,” Robotica, vol. 19, no. 3, pp. 311–321, 2001. View at: Publisher Site  Google Scholar
 H.N. Nguyen, J. Zhou, and H.J. Kang, “A calibration method for enhancing robot accuracy through integration of an extended Kalman filter algorithm and an artificial neural network,” Neurocomputing, vol. 151, pp. 996–1005, 2015. View at: Publisher Site  Google Scholar
 B. Wu, M. Yang, K. Li et al., “A reverberationtimeaware DNN approach leveraging spatial information for microphone array dereverberation,” EURASIP Journal on Advances in Signal Processing, vol. 81, pp. 1–13, 2017. View at: Publisher Site  Google Scholar
 H. B. Huang, X. R. Huang, J. H. Wu, M. L. Yang, and W. P. Ding, “Novel method for identifying and diagnosing electric vehicle shock absorber squeak noise based on a DNN,” Mechanical Systems and Signal Processing, vol. 124, pp. 439–458, 2019. View at: Publisher Site  Google Scholar
 H. B. Huang, J. H. Wu, X. R. Huang, M. L. Yang, and W. P. Ding, “The development of a deep neural network and its application to evaluating the interior sound quality of pure electric vehicles,” Mechanical Systems and Signal Processing, vol. 120, pp. 98–116, 2019. View at: Publisher Site  Google Scholar
 E. M. Ozuna, D.L. Flores, E. G. López, D. Cervantes, and P. Juárez, “Model based on GA and DNN for prediction of mRNASmad7 expression regulated by miRNAs in breast cancer,” Theoretical Biology and Medical Modelling, vol. 15, no. 1, pp. 1–12, 2018. View at: Publisher Site  Google Scholar
 B. N. Panda, M. V. A. R. Bahubalendruni, and B. B. Biswal, “A general regression neural network approach for the evaluation of compressive strength of FDM prototypes,” Neural Computing and Applications, vol. 26, no. 5, pp. 1129–1136, 2015. View at: Publisher Site  Google Scholar
 K. He, X. Zhang, S. Ren, and J. Sun, “Delving deep into rectifiers: surpassing humanlevel performance on imageNet classification,” in Proceedings of IEEE International Conference on Computer Vision, pp. 1026–1034, Santiago, Chile, December 2015. View at: Google Scholar
 J. Zhou and H.J. Kang, “A hybrid leastsquares genetic algorithm–based algorithm for simultaneous identification of geometric and compliance errors in industrial robots,” Advances in Mechanical Engineering, vol. 7, no. 6, pp. 1–12, 2015. View at: Publisher Site  Google Scholar
 D. P. Kingma and J. L. Ba, “Adam: a method for stochastic optimization,” in Proceedings of International Conference on Learning Representations, pp. 1–13, San Diego, CA, USA, 2015. View at: Google Scholar
Copyright
Copyright © 2020 Xiaoyan Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.