Structural parameters identification for industrial robot using a hybrid algorithm

To improve the precision and reduce the movement uncertainty of the industrial robot, a novel hybrid optimization algorithm which combines adaptive genetic algorithm with simulated annealing algorithm is proposed in this article. First, for the sake of increasing the global exploring ability of relevant individuals, the adaptive crossover and mutation operator are used in the phase of adaptive genetic algorithm. If the population optimized by adaptive genetic algorithm is trapped in the local optimal area and simultaneously meets the transformation rule, then it is consequently optimized by simulated annealing to enhance the population diversity and hunt for a better solution so that the probability of finding the global optimal solution is greatly increased. Then, corresponding experiments based on single point repeatability are conducted to acquire data and identify the structural parameters of the industrial robot. Moreover, the single point repeatability test and length test are all implemented at the same time to verify the effectiveness of the proposed method. At last, the result reveals that the proposed method is effective to identify the real structural parameters of the industrial robot, thus enormously decreasing the single point repeatability and length deviation at the same time, which extremely increases the precision and decreases the movement uncertainty of the industrial robot.


Introduction
With the rapid development of manufacturing and industrial applications, the industrial robot 1 is the universally used apparatuses in various fields, which makes a profound and critical impact in our production and life. For example, it is not strange to find that the industrial robot is widely used in the fields of measurement, assembly, detection, transportation, and so on. 2 Composed of a base, two arms, three pairs of rotational joints, a probe, and some other relevant components, 2 one kind of industrial robot illustrated in Figure 1 is a special instrument which makes it able to work expediently, move easily, and display data conveniently.
Although this kind of industrial robot has many merits, but the coordinate of the probe tip modifies in a large scale if the structural parameters change in a small scale resulting from the fact that the coordinate of the probe tip is obtained School of Mechanical Engineering, Hubei University of Technology, Wuhan, China by a nonlinear transformation function from all the rotational angles of the joints to the probe tip. 2 Based on this fact and analysis, it is obviously deduced that the coordinate error of the probe tip is greatly and mainly impacted by the structural parameters of the industrial robot. 3 However, what should be noticed is that the real values of the structural parameters are not equal to the corresponding nominal ones in almost cases after the industrial robot has been utilized for a period of time as there exists the actual influential factors such as deviations of temperature and humidity, action of gravity, vibrations and interactions between various instruments, and so on. [3][4][5] Hence, affirmation can be made that the real coordinate of the probe tip is not commonly equal to the one which is calculated using the nominal structural parameters. Therefore, it is extremely required for us to improve the accuracy of the industrial robot in various fields, to which a great number of researchers devoted much works and energies.
Nowadays, there are two main methods to improve the accuracy of the industrial robot. 6,7 The first one is to use higher precise components such as high precise coders, sensors, arms, and so on, which in turn dramatically increases the manufacturing cost of the industrial robot. The second one is to identify the structural parameters and compensate them again using a little auxiliary instruments or algorithms without any increased cost of components. 7,8 But, because of the nonlinear transformation function between the coordinate of the probe tip and the structural parameters used in the industrial robot, 2-6 small deviations of the structural parameters may accumulate and lead to relative substantial coordinate errors of the probe tip. Thus, it is obvious that increasing only one or several components' precision is not capable of effectively decreasing the whole measurement errors in a good degree. As a result, it is essential and necessary for us to execute another method, that is, identify the real structural parameters by other approaches. For example, Fang and Dang 9 proposed a kinematic parameters identification method using quantumbehaved particle swarm optimization, which is able to simultaneously identify the open-chain series mechanism and parallel mechanism. Santolaria et al. 10 identified the structural parameters of the industrial robot by utilizing four laser tracers to establish a geometric multilateration system and obtained data from laser tracks. Liu et al. 11 combined a geometric identification with the parameters identification to propose a new two-step calibration method. He et al. 12 used a laser tracker for spraying robot calibration, where the precision of the robot was improved by 60%. An elasto-geometrical calibration method was used by Kamali et al. 13 in 2016, in which the geometric parameter errors and joint stiffness parameters were identified at the same time by using a laser tracker.
However, although the methods mentioned above give us a lot of selections and approaches to calibrate and compensate the structural parameters of the industrial robot, but they nearly have the same characteristic that almost all the methods previously mentioned rely on higher precise apparatus especially laser trackers 14,15 which are extreme expensive and affected easily by the external environment, resulting in a relative high cost and nondeterminacy. Moreover, by using lots of precise instruments, the calibration and compensation procedures of the methods are complicated and difficult for us to operate sometimes, which is their defect and bottleneck. To solve this problem and simplify the calibration procedure, a new method is proposed in this article based on single point repeatability by using a hybrid optimization algorithm of adaptive genetic algorithm (AGA) 16 and simulated annealing (SA) 17,18 algorithm. Without using higher precise instruments or complex calibration procedures in the new proposed method, only a conical hole is fixed on a base, to which the probe tip is inserted and the arms of the industrial robot are rotated to a great number of different configurations to generate the acquired data set. Then the data set is utilized by the hybrid algorithm to get the optimized solution, that is, the optimal structural parameters of the industrial robot. And then, the identified parameters of the industrial robot are used to recalculate the single point repeatability and length deviation and some remarkable conclusions are presented at last.
To complete the calibration process, four steps are commonly carried out step by step 19,20 : modeling, measurement, identification, and compensation. Modeling is to formulate the transformation function between the coordinate of the probe tip and the structural parameters of the industrial robot. As was presented by Roth et al., 19 a perfect kinematic model of the robot requires three conditions: (1) completeness, (2) continuousness, and (3) minimality. In recent years, it is known to us all that there are a lot of mathematical models used to calibrate the structural parameters of the industrial robot such as Denavit-Hartenberg (DH) model, 21,22 POE model, 23 CPC model, 24 MDH model, 25 S model, 26 and so on. However, what should be noted is that the DH model and MDH model are most widely utilized in robotic area on account that they have clear physical meaning and thus used conveniently. Therefore, without loss of generality, the MDH model is adopted and used in this article to formulate the transformation function between the coordinate of the probe tip and the structural parameters of the industrial robot.
The main contributions of this article can be briefly summarized as follows: 1. A new computation mechanism combining AGA with SA to formulate simulated annealing-based adaptive genetic algorithm (AGA-SA) is introduced in this article. Through this integration, the hybrid algorithm AGA-SA not only overcomes the weakness by SA that the AGA is easy to be trapped in the local optimal area but also strengthens the ability of local searching by AGA. That is to say, the AGA-SA is not only able to jump out of the local optimal area when it is trapped in it but also maintains strong local searching ability at the same time, which resulting in a higher probability of not only finding the preciser solution of the structural parameters but also increasing the precision and decreasing the movement uncertainty of the industrial robot. 2. The kinematic equation is mathematically analyzed in detail and conclusions are derived that it is impossible for us to precisely measure the real coordinates of the probe tip of the industrial robot, which reveals that it consequently nearly impossible for us to directly calculate the real structural parameters by just acquiring the joint angles and coordinate of the industrial robot. Hence, the algorithm AGA-SA and method of single point repeatability are naturally introduced and integrated to optimize the structural parameters and acquire their real ones. 3. Different from other methods, it is not necessary for us to set precise initial parameters of the structural parameters to be calibrated in the proposed calibration method, that is, the parameters are set arbitrarily as long as they are not too big or small, which simplifies the computation and specification process of the parameters.
The article is organized as follows. The second section formulates the mathematical model of the transformation function from the joint angles to the probe tip. The third section shows the relevant principles and analyses of the scheme. Then, the hybrid algorithm combing AGA and SA is introduced in the fourth section. Moreover, experiments including single point repeatability test and length deviation test are conducted in the fifth section. At last, remarkable conclusions are given in the sixth section.

Kinematic modeling
The aim of kinematic modeling is to construct the mathematical relationship between the coordinate of the probe tip and the structural parameters of the industrial robot. According to MDH model, 25 the coordinate of the probe tip can be computed by equation where q i is the joint angle, d i is the linkage length, a iÀ1 is the joint offset, a iÀ1 is the torsion angle, and l x , l y , and l z are the components of the probe length in x direction, y direction, and z direction, respectively. And, x, y, and z are the coordinates of the probe tip in x direction, y direction, and z direction, respectively.
It should be noted that there is a deviation between the zero position of the coder and that of the joint, so the joint angle q i is expressed by equation where q 0;i is the deviation of the zero position and q s;i is the rotational angle of the coder. Note that q 0;i is a constant and a structural parameter. The nominal structural parameters of the industrial robot used in this article are shown in Table 1.

Kinematic analysis
Set X ¼ ðq 0;1 ; q 0;2 ; q 0;3 ; q 0;4 ; q 0;5 ; q 0;6 ; d 1 ; d 2 ; d 3 ; d 4 ; d 5 ; d 6 ; a 0 ; a 1 ; a 2 ; a 3 ; a 4 ; a 5 ; a 0 ; a 1 ; a 2 ; a 3 ; a 4 ; a 5 ; l x ; l y ; l z Þ, then the coordinate ðx; y; zÞ of the probe tip can be rewritten as the structural parameter X and the readings of the coder by equation (3) The total differential with respect to every structural parameter of equation (3) is formulated by equation If the structural parameters of the industrial robot change in an extreme small scale, then the differentials of the structural parameters and the coordinates of the probe tip are referred as equal to their corresponding deviations, that is Â dq 0;1 ; dq 0;2 ; dq 0;3 ; dq 0;4 ; dq 0;5 ; dq 0;6 ; dd 1 ; dd 2 ; dd 3 ; dd 4 ; dd 5 ; dd 6 ; da 0 ; da 1 ; da 2 ; da 3 ; da 4 ; da 5 ; da 0 ; da 1 ; and ½dxdydz ¼ ½DxDyDz On the condition that the nominal structural parameters are given beforehand and 27 structural parameters need to be calibrated, if 9 groups of real coordinates of the probe tip and corresponding joint angles are obtained by some instruments, then there are 27 groups of equations derived by equations (4), (5), and (6), which can help to not only directly calculate the deviations of the structural parameters but also compensate the nominal structural parameters by these deviations in turn.

Analysis
Based on the analyses previously mentioned, the coordinates of the probe tip can be obtained and the structural parameters can be calibrated by equation (3). Specifically speaking, for a given X, equations can be formulated from equation (4) by using higher precise instruments to acquire N groups of different points k Pð k x; k y; k zÞðk ¼ 1; 2; . . . ; N Þ which are mathematically expressed by equation where k ¼ 1; 2; . . . ; N . Here, index k indicates the k groups of coordinates of the probe tip. By solving equation (7), the real structural parameters are obtained and can consequently be used as the new nominal parameters of the industrial robot. However, the fact that the general coordinates k Pð k x; k y; k zÞðk ¼ 1; 2; . . . ; N Þ in equation (7) are measured by higher precise instruments makes it difficult even impossible for us to implement the process in real fields. Therefore, it is extreme necessary for us to use other ways to settle this problem. One method to solve this difficulty and simplify this process is to use the single point repeatability calibration. For this process, a single point is measured for a great number of times with different poses for obtaining the corresponding joint angles and calibrating the real structural parameters, which is the main principle of this method.

Objective function
As mentioned above, after fixing the conical hole on the base, the probe tip is inserted to the conical hole still and the arms are rotated to N kinds of different configurations for obtaining N groups of different joint angles k q i ði ¼ 1; 2; :::; 6; k ¼ 1; 2; :::; N Þ. Suppose an arbitrarily given structural parameter X is known, then a series of points k Pð k x; k y; k zÞðk ¼ 1; 2; :::; N Þ can be computed by substituting k q i ði ¼ 1; 2; :::; 6; k ¼ 1; 2; :::; N Þ and X into equation (7) and the deviations between the measured coordinates and the real coordinates are thus calculated by equation where x, y, and z are the real coordinates of the probe tip in direction of x, y, and z, respectively. Note that it is not easy for us to measure the real coordinate ðx; y; zÞ of the probe tip especially in real fields, hence the coordinate ðx; y; zÞ is replaced by its mean value expressed by equation After given k E, then their average is However, the average deviation E cannot indicate the distribution regularity of this group of data, so the standard deviation s is introduced by equation Consequently, the objective function for X based on "3 sigma principle" is Let X Ã be the real structural parameter of the industrial robot, then Equation (12) is utilized as the objective function to search for the global optimal solution which is accepted as the real structural parameter of the industrial robot and compensates them in turn. Note that the closer X approaches to X Ã , the smaller FðX Þ is. Therefore, the main work thereafter is to obtain the different joint angles and use them to search for the global optimal solution by usage of optimization algorithm.

Adaptive genetic algorithm
Genetic algorithm (GA), 27 a natural bioinspired intelligent optimization algorithm simulating the biological mechanism of selection, crossover, and mutation, is a widely used evolutionary optimization method in various fields, which obeys the principle of "survival of the fitness" proposed by Charles Robert Darwin. In GA, the individuals' survival abilities are indicated by their own fitness values, that is, the individuals with larger fitness values have better survival abilities and thus have larger probabilities of executing the genetic operation and giving birth to child. Note that the individuals with larger fitness values have better survival abilities while the aim of this problem in this article is to search for the best solution with the minimal objective value, hence the fitness function takes the form of where f itnessðÁÞ is the fitness value of the corresponding variable, reallim is an extreme small real constant which ensures that the denominator is not equal to zero when FðX Þ is equal or extremely approaches to zero. However, although GA has a good global searching ability, it can be trapped in a local optimal area with a certain probability especially when the optimized problem is extreme nonlinear such as the problem in this article. Therefore, the AGA 16 is used in this article to increase the probability of hunting for the global optimum. In AGA, the probability of crossover and mutation are all self-adaptive dynamically changing by taking the form of and where k 1 ; k 2 ; k 3 ; k 4 1 are constants, f max , f , f 0 , f are the maximal fitness value, the mean fitness value, the larger fitness value of the two crossover individuals, and the mutation fitness value, respectively. In this article, k 1 ; k 2 ; k 3 ; k 4 are set to be 1, 0.5, 1, and 0.5, respectively. For clarification, the details about AGA are specified as follows.

I.
Selection. In selection phase, the individuals are selected by roulette wheel 28 . II. Crossover. Let X i ðtÞ and X j ðtÞ be any two arbitrary individuals in the tth generation, then their corresponding newly generated individuals in the ðt þ 1Þth generation after executed the operation of crossover are updated by equation and X j ðt þ 1Þ ¼ rX j ðtÞ þ ð1 À rÞX i ðtÞ where r is a constant given beforehand.

III. Mutation. Mutation operation is executed successively as the following steps:
ffi Generate the mutation position u ¼ randiðDÞ, where D is the dimension of the variable X and randi is a system function which can randomly selects an integer from integer array varying from 1 to D wherein this article D ¼ 27. ffl Define variable a ¼ X lowerlimitðuÞ, b ¼ X upperlimitðuÞ, where X lowerlimit and X upperlimit are the minimal limit vector and maximal limit vector of variable X, respectively. Let sðtÞ ¼ 1 À ð t G Þ b , where t is the current generation number, G and b are the maximal generation number and a constant, respectively, which are all given beforehand. Ð Let X i;u ðtÞ and X i;u ðt þ 1Þ be the values in the uth position of the tth and ðt þ 1Þth generation, respectively, of any arbitrary individual i; X i;u ðt þ 1Þ lower and X i;u ðt þ 1Þ upper be the lower value and upper value of X i;u ðt þ 1Þ, respectively, then X i;u ðt þ 1Þ lower ¼ X i;u ðtÞ À sðtÞðX i;u ðtÞ À aÞ (19) and X i;u ðt þ 1Þ upper ¼ X i;u ðtÞ þ sðtÞðb À X i;u ðtÞÞ (20) ð At last, the value X i;u ðt þ 1Þ in the uth position of X i ðt þ 1Þ is updated by equation where rand is a random function arbitrarily generating a number from 0 to 1.

SA algorithm and AGA-SA
SA 17,18 algorithm, which simulates the annealing process of the metal, has a good characteristic of escaping from the local optimal area when the best solution found so far is trapped in the local optimal region. Motivated by this trait, if the solution optimized by AGA is trapped in the local optimal region, then it is optimized by SA through a transformation rule to help the current best individual to escape from the local optimal region and wander widely in a broader space. Let t, X best ðtÞ and f itnessðX best ðtÞÞ be the current generation number, the best individual found so far until tth generation, and the fitness value of the best individual X best ðtÞ, respectively, then the transformation rule takes the form of jf itnessðX best ðtÞÞ À f itnessðX best ðt À KÞÞj e; where K is a positive integer, e and d are transformational thresholds. Note that K, e, and d are constants and set by ourselves beforehand. That is to say, when the solution X best ðtÞ optimized by AGA meets equation (22), then the optimization process enters the SA process, that is, the best individual X best ðtÞ optimized by AGA is optimized by SA furtherly. Note that if the best solution X best ðtÞ found so far is not the real parameter until AGA-SA runs a great number of generations and the best solution in the last generation is nearly kept still, then the whole optimization process cannot find the global optimum. Thus, a termination condition is set to finish the unnecessary redundant computation of AGA-SA if the optimized result meets the accuracy demand. In this article, the termination condition is set as where t is the generation number, G is the maximal generation number, t and h are thresholds which determine the ending of the whole optimization process. The relative parameters are set as follows: ffi Population number N ¼ 60; r ¼ 0: To make comparison, the AGA, SA, and AGA-SA are respectively utilized as the algorithm to get the structural parameters of the industrial robot and their corresponding obtained optimized structural parameters are respectively implanted into the corresponding controller to replace the corresponding old ones. Moreover, the industrial robot with the old structural parameters and the corresponding new structural parameters is reused again to conduct the experiments of single point repeatability and length deviation test and the relevant outcomes are respectively presented thereafter. The flow chart of AGA, SA, and AGA-SA for calibrating the structural parameters of the industrial robot are shown in Figures 2, 3, and 4, respectively. Figure 2. The flow chart of AGA. AGA: adaptive genetic algorithm.

Data acquisition
As is illustrated in Figure 5, after the conical hole and the base of the industrial robot are simultaneously fixed on the platform firstly, then the arms of the industrial robot are manipulated to make the probe tip inserted to the conical hole and fitted with it. Moreover, the arms are rotated again to N groups of different poses and the readings k q i ði ¼ 1; 2; . . . ; 6; k ¼ 1; 2; . . . ; N Þ of the encoder at each position are recorded from the man-machine interface illustrated in Figure 6 while maintaining the probe tip fitted with the conical hole during the whole process. The data set k q i ði ¼ 1; 2; . . . ; 6; k ¼ 1; 2; . . . ; N Þ is then acquired and utilized thereafter. Note that the top and bottom diameter of the conical hole are a little bigger and smaller than that of the probe tip, respectively. Therefore, the probe tip can be trapped and fitted in the conical hole to maintain the spatial position of the center of the probe tip unchanged when the arms are rotated to different poses.

Parameters identification
After obtaining N groups of joint angles, they are used as the data set by AGA, SA, and AGA-SA, respectively, to optimize the structural parameters. If a certain best individual X best ðtÞ satisfies equation (23) or equation (24), then it is acknowledged that the algorithm successfully finds the Figure 3. The flow chart of SA. SA: simulated annealing. optimal solution X best ðtÞ denoted by X optimum which is accepted as the real parameter of the industrial robot. Note that X optimum which is derived by relevant algorithm and test data is usually infinitely close but not equal to X Ã due to the existence of errors of real engineering problems.
Other critical problems are the range settings of the various variables of the industrial robot and their boundary treatments when the variables exceed the corresponding boundaries. In this article, because the values of the length parameters are nonnegative numbers and the nominal values of a iÀ1 ði ¼ 1; 2; 5; 6Þ, d 6 , and a iÀ1 ði ¼ 1; 3; 4Þ are all 0, so the deviation ranges of a iÀ1 ði ¼ 1; 2; 5; 6Þ and d 6 are set to be (0, 10); of a iÀ1 ði ¼ 3; 4Þ and d i ði ¼ 1; 2; :::; 5Þ are set to be +10; of a iÀ1 ði ¼ 1; 2; :::; 6Þ is set to be +10 ; of q 0;i ði ¼ 1; 2; :::; 6Þ is set to be +1 ; and of l x , l y , and l z are set to be +10, respectively.
When a certain variable y among the structural parameters surpasses the corresponding limit, then it is newly generated randomly by a generator, which is expressed by Criterion A illustrated as: Criterion A. If y exceeds the corresponding constraint y ¼ ðuplimit À lowerlimitÞ: Ã rand þ y end where uplimit and lowerlimit are the upper boundary and lower boundary of variable y, respectively.
Note that the closer the any two arbitrary individuals X i and X j are, the smaller of the gap between their fitness values are. That is to say, if jX i À X j j 2 < jX i À X k j 2 , then j f itnessðX i Þ À f itnessðX j Þj < jf itnessðX i Þ À f itnessðX k Þj, where j Á j 2 indicates 2-norm of the corresponding variable. Based on this fact, the term f itness max ðtÞ À f itness average ðtÞ is introduced in this article to indicate the individual density of the population in the tth generation, where f itness max ðtÞ and f itness average ðtÞ are the maximal fitness value and average fitness value in the tth generation, respectively. The changing curves of the objective function and f itness max ðtÞ À f itness average ðtÞ optimized by AGA-SA are shown in Figures 7 and 8, respectively. The calibrated structural parameters of the industrial robot optimized by

Verification of the calibrated parameters
Because it is difficult for us to measure the real structural parameters of the industrial robot, thus the single point repeatability test illustrated in Figure 5 and the length test illustrated in Figure 9 are simultaneously used in this phase to verify the effectiveness and correctness of the method used in this article. As is shown in Figure 5, the single point is tested for many times and repeatability results with respect to AGA, SA, and AGA-SA are shown in Figure 10.
As is shown in Figure 9, the two standard balls are attached to the fixture plate with a standard length L between the two ball centers O 1 and O 2 and the fixture plate is fixed by a vise on the platform. In this test, the arms of the industrial robot are rotated to make the ball of the probe tip attached to the standard balls in different positions to obtain different groups of joint angles, with which the length L between the two standard balls is computed through the transformation function. For the sake of     narrative, the error DL between the standard length L and the tested length L 0 is defined as Let L 0 1 and L 0 2 be the tested length using the nominal structural parameters and identified structural parameters of the industrial robot, respectively, then and For comparison, the standard length L is tested for many times and the corresponding error curves are shown in Figure 11. Note that the radiuses of the two standard balls are all R and the length between the two centers O 1 and O 2 of the corresponding standard balls is a fixed distance L, which is illustrated in Figure 9(b).

Discussion
After the experiments are completed and verifications are accomplished which are illustrated by the tables and figures, observations can be made and conclusions can be drawn consequently, which are divided into Part 1 and Part 2, respectively.  respectively; after optimized by SA are 0.1988 mm, À0.1940 mm, À0.0093 mm, and 0.1171 mm, respectively; and after optimized by AGA-SA are 0.0493 mm, À0.0499 mm, À0.0080 mm, and 0.0289 mm, respectively.
From the results presented above, analyses can be conducted and some conclusions can be acquired that the Max, Min, Mea, and Std of DL 2 optimized by AGA with respect to DL 1 are decreased by 74%, 64%, 84%, and 72%, respectively; optimized by SA with respect to DL 1 are decreased by 44%, 30%, 81%, and 42%, respectively; and optimized by AGA-SA with respect to DL 1 are decreased by 86%, 82%, 83%, and 86%, respectively.
Overall, from the observations and analysis presented above, it is evidently derived that: 1. It is nearly impossible for us to conduct the method of measuring the real coordinates of the probe tip of the industrial robot and directly computing the real structural parameters using the tested data in the real fields, thus a method that combines testing the rotational angles of the probe tip of the industrial robot and utilizing the hybrid algorithm AGA-SA is proposed to simplify the process and resolve the contradictions. 2. The Max, Min, Mea, and Std of the single point repeatability and length error, including their three components in the corresponding directions, are all decreased with a big scale after optimized by corresponding algorithms (AGA, SA, and AGA-SA). However, nearly all the Max, Min, Mea, and Std of the single point repeatability and length error are decreased with a greater scale after optimized by AGA-SA than those optimized by AGA or SA at the same time. 3. Through this analysis, it is shown that the precision and relevant errors of the industrial robot are respectively increased and decreased drastically with a bigger corresponding increasing and decreasing magnitude after optimized by AGA-SA than those optimized by AGA and SA, which indicates that this proposed method and algorithm are effective and more suitable for us to calibrate the real structural parameters of the industrial robot than AGA and SA simultaneously.

Conclusions
A kinematic model of the industrial robot is established based on MDH model firstly. Then, the kinematic model and error formula are analyzed and results are derived that: (1) it is extreme difficult even impossible for us to measure the real coordinates of the industrial robot and; (2) only using a simple conventional algorithm is easy to be trapped in the local optimal region and thus incapable of searching for the real structural parameters of the industrial robot. Moreover, based on these facts and to solve the problem presented previously, the algorithm of AGA and SA are combined to formulate a hybrid algorithm AGA-SA which simultaneously has a good local searing ability and global exploring capacity. And then, the experiments of the single point repeatability and the length test are conducted combining the hybrid algorithm to calibrate the structural parameters and verify the effectiveness of the proposed method. The result shows that the single point repeatability and the length deviation after calibrated by AGA-SA are all decreased substantially and with a bigger magnitude compared with those calibrated by AGA or SA, which suggests that this method is effective and more suitable for us to greatly increase the precision and decrease the movement uncertainty of the industrial robot than AGA and SA. The proposed modeling and identification method can be commonly expanded to the calibration of a serials of robots, where the identification method is required to increase the accuracy of the corresponding robots. This will be furtherly studied in our future work.

Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.