Submit manuscript...
eISSN: 2574-8092

International Robotics & Automation Journal

Opinion Volume 3 Issue 2

The problem in accuracy compensation of industrial robot

Shuanglong Liu, Wenhe Liao, Wei Tian, Lin Zhang, Jiachen Jiao

Nanjing University of Aeronautics and Astronautics, China

Correspondence: Wei Tian, Nanjing University of Aeronautics and Astronautics, China, Tel 13851662331, Fax 84891836

Received: August 31, 2017 | Published: September 27, 2017

Citation: Liu S, Liao W, Tian W, et al. The problem in accuracy compensation of industrial robot. Int Rob Auto J. 2017;3(2):282-283. DOI: 10.15406/iratj.2017.03.00050

Download PDF

Abstract

Industrial robots are increasingly being used in automated production line with their high flexibility and low cost. However, the poor features of repeat positioning accuracy and absolute positioning accuracy of industrial robots have been a bottleneck restricting its application in higher accuracy situation. Many scholars have studied methods to improve the accuracy of robot. But there is congenital defect for industrial robots when their accuracy is further improved by compensation and few people pay attention to this. In this article, the backlash of industrial robot joint is discussed and its impact on accuracy compensation is analyzed. Finally, the method of eliminating the backlash and further improving the robot`s accuracy based on joint feedback is proposed.

Keywords: industrial robots, position accuracy, backlash, joint feedback

Introduction

From the published article, it seems that there is a 0.3mm accuracy limit which cannot be break no matter what model we used to compensate (for high-load 6R series industrial robots without external testing equipment). Veitschegger & Wu1 built the kinematic model of PUMA 560 and used the differential error transform to calibrate kinematic parameters; the worst case compensated error is 0.3mm (only in one direction). Nubiola2 also calibrated the parameters of robot and the result of positioning error is 0.696mm. Wei3 proposed the concept of error similarity and adopted the method of spatial interpolation to find the size of error. Result has shown that the maximum value of positioning error is 0.386mm. This is due to the existence of backlash errors in the joints, which are caused by the reverse movement of the joints and it can be expressed by the error of moving from different directions to the same position. In this paper, the grating scales are installed at the joints of the robot, by which the size of the backlashes are measured. Through the analysis of defects in parameter calibration and accuracy compensation, the method improving robot accuracy based on joint feedback is proposed. The preliminary experimental results show that the maximum absolute position error is 0.14mm.

The influences of backlash in accuracy compensation

The backlash is generated in the joints of the robot, which is related to the direction of movement of the robot. It can be regarded as an angle error of the joint angle. The position error of the robot in the Cartesian space is generated by the error of each joint, and it can be calculated by the parameter error model of the robot3 when the previous motion state of each joint is known. There are many models that describe the end-effector error of 6R serial robot and they can be separated into two categories: parameter calibration based on model and error model based on sampling. However, both two categories require sampling of robot and the backlash error will exists in sampling points. Figure 1 shows the backlash error of five random points in the Cartesian coordinate system. Those points move in the direction from p1 to p5, and then go back from p5 to p1. At the same time, the laser tracker records each point in different direction of movement and error is calculated by the equation of ΔP= Δ x 2 +Δ y 2 +Δ z 2 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqipy0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbmqaaeWacmGaaiaabmqaamaabaabaaGcbaqcfaOaeyiLdq Kaamiuaiabg2da9maakaaabaGaeyiLdqKaamiEamaaCaaabeqcfasa aiaaikdaaaqcfaOaey4kaSIaeyiLdqKaamyEamaaCaaajuaibeqaai aaikdaaaqcfaOaey4kaSIaeyiLdqKaamOEamaaCaaabeqcfasaaiaa ikdaaaaajuaGbeaaaaa@477F@ . The experiment was done 30 times and the error was averaged. In order to make the result contrast, the backlash error is also measured under the condition of the joints feedback. It`s clear that the value of backlash error is far more than the repeated accuracy of 6D serial robot (±0.06mm for KUKA KR210) in some points. The error of the sampling point not only includes the deterministic error caused by the connecting rod parameter but also the random error caused by the backlash and the repeated positioning accuracy of the robot. For the model based sampling, the error measured at one point is not the error when the robot moves to. This is because the robot could be reach one position from different direction. Once the movement direction of joint change, then the backlash error occur. As for the parameter calibration based on kinematic, the sampling point error contains a relatively large random part, and the calculated link parameter error is inaccurate.

Figure 1 Backlash error in Cartesian coordinate system.

The accuracy compensation technology based on joints feedback

This paper eliminates the backlash by installing grating at the first three joints. The backlash error at each joint can be nearly zero by establishing a PID model with the grating feedback. The relationship between the grating readings and the actual rotation angle becomes critical when it is known that raster feedback can accurately reach a joint angle. In fact, due to the existence of the robot geometric parameters error, the actual value of the joint rotation angle of the robot is not necessarily related to the grating reading, which is different in the case of other joint changes. At the same time, only three gratings are installed at robot joints, and it is impossible to calibrate the kinematic parameters of the robot under the condition of joint feedback. Referring to the principle of accuracy compensation, a method of error coupling is proposed in this paper. The spatial position error of the target point is corrected by only rotating the first three axes of the robot. Hence, based on the similarity principle of the space error,4 an estimation model of space error of robot is established. The error of the target point can be transformed into the angle error of the first three joints by the error coupling model. And finally, the error is corrected by the closed loop feedback of the first three joint. By means of Kriging interpolation,5,6 the position error estimation model is established:

e l ^ = w T e l , l=x,y,z MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D aebbnrfifHhDYfgiuL2zLjhasaacH8srps0lbbf9q8WrFfeuY=Hhbb f9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq =He9q8qqQ8frFve9Fve9Ff0dmeaabaqaceGacaGaaeWabaWaaeaaea aakeaajuaGdaqiaaqaaiaadwgadaWgaaqaaiaadYgaaeqaaaGaayPa daGaeyypa0Jaam4DamaaCaaabeqcfasaaiaadsfaaaqcfaOaamyzam aaBaaabaGaamiBaaqabaGaaiilaiaabccacaWGSbGaeyypa0JaamiE aiaacYcacaWG5bGaaiilaiaadQhaaaa@4B2D@  (1)

 Where e l ^ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D aebbnrfifHhDYfgiuL2zLjhasaacH8srps0lbbf9q8WrFfeuY=Hhbb f9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq =He9q8qqQ8frFve9Fve9Ff0dmeaabaqaceGacaGaaeqabaWaaeaaea aakeaajuaGdaqiaaqaaiaadwgadaWgaaqaaiaadYgaaeqaaaGaayPa daaaaa@3DD2@ is the estimation error of composition point; w R m MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D aebbnrfifHhDYfgiuL2zLjhasaacH8srps0lbbf9q8WrFfeuY=Hhbb f9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq =He9q8qqQ8frFve9Fve9Ff0dmeaabaqaceGacaGaaeWabaWaaeaaea aakeaajuaGcaWG3bGaeyicI4SaamOuamaaCaaajuaibeqaaiaad2ga aaaaaa@3FAF@ is the weighted vector between composition point and sampling points; e l R m MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D aebbnrfifHhDYfgiuL2zLjhasaacH8srps0lbbf9q8WrFfeuY=Hhbb f9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq =He9q8qqQ8frFve9Fve9Ff0dmeaabaqaceGacaGaaeWabaWaaeaaea aakeaajuaGcaWHLbWaaSbaaeaacaWGSbaabeaacqGHiiIZcaWGsbWa aWbaaeqajuaibaGaamyBaaaaaaa@40B3@  represents the vector of the errors of sampling points in the l MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D aebbnrfifHhDYfgiuL2zLjhasaacH8srps0lbbf9q8WrFfeuY=Hhbb f9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq =He9q8qqQ8frFve9Fve9Ff0dmeaabaqaceGacaGaaeWabaWaaeaaea aakeaajuaGcaWGSbaaaa@3C07@ direction.

 Weighted vector can be obtained by the following formula:

{ λ=( F T R 1 F)( F T R 1 rf) W= R 1 (rFλ) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D aebbnrfifHhDYfgiuL2zLjhasaacH8srps0lbbf9q8WrFfeuY=Hhbb f9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq =He9q8qqQ8frFve9Fve9Ff0dmeaabaqaceGacaGaaeWabaWaaeaaea aakeaajuaGdaGabaabaeqabaGaaq4Udiabg2da9iaacIcacaWGgbWa aWbaaeqajuaibaGaamivaaaajuaGcaWGsbWaaWbaaKqbGeqabaGaey OeI0IaaGymaaaaieGajuaGcaWFgbGaaiykaiaacIcacaWGgbWaaWba aeqajuaibaGaamivaaaajuaGcaWGsbWaaWbaaKqbGeqabaGaeyOeI0 IaaGymaaaajuaGcaWGYbGaeyOeI0IaamOzaiaacMcaaeaacaWGxbGa eyypa0JaamOuamaaCaaabeqcfasaaiabgkHiTiaaigdaaaqcfaOaai ikaiaadkhacqGHsislcaWGgbGaaq4UdiaacMcaaaGaay5Eaaaaaa@5B1D@  (2)

Where λ MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D aebbnrfifHhDYfgiuL2zLjhasaacH8srps0lbbf9q8WrFfeuY=Hhbb f9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq =He9q8qqQ8frFve9Fve9Ff0dmeaabaqaceGacaGaaeWabaWaaeaaea aakeaacaaH7oaaaa@3BCC@ is the Lagrange multiplier, f= [1  θ 1    θ n ] T MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D aebbnrfifHhDYfgiuL2zLjhasaacH8srps0lbbf9q8WrFfeuY=Hhbb f9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq =He9q8qqQ8frFve9Fve9Ff0dmeaabaqaceGacaGaaeWabaWaaeaaea aakeaajuaGcaWGMbGaeyypa0Jaai4waiaaigdacaqGGaGaaqiUdmaa BaaajuaibaGaaGymaaqabaqcfaOaaeiiaiabl+UimjaabccacaaH4o WaaSbaaKqbGeaacaWGUbaajuaGbeaacaGGDbWaaWbaaeqajuaibaGa amivaaaaaaa@4A6C@  ,and F R m MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D aebbnrfifHhDYfgiuL2zLjhasaacH8srps0lbbf9q8WrFfeuY=Hhbb f9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq =He9q8qqQ8frFve9Fve9Ff0dmeaabaqaceGacaGaaeWabaWaaeaaea aakeaajuaGcaWGgbGaeyicI4SaamOuamaaCaaabeqcfasaaiaad2ga aaaaaa@3F7E@ denote the angle matrix constructed by each point`s f MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D aebbnrfifHhDYfgiuL2zLjhasaacH8srps0lbbf9q8WrFfeuY=Hhbb f9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq =He9q8qqQ8frFve9Fve9Ff0dmeaabaqaceGacaGaaeWabaWaaeaaea aakeaajuaGcaWGMbaaaa@3C01@ ; r R m MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D aebbnrfifHhDYfgiuL2zLjhasaacH8srps0lbbf9q8WrFfeuY=Hhbb f9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq =He9q8qqQ8frFve9Fve9Ff0dmeaabaqaceGacaGaaeWabaWaaeaaea aakeaajuaGcaWGYbGaeyicI4SaamOuamaaCaaabeqcfasaaiaad2ga aaaaaa@3FAA@ is the correlation vector between composition point and sampling points; R R m×m MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D aebbnrfifHhDYfgiuL2zLjhasaacH8srps0lbbf9q8WrFfeuY=Hhbb f9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq =He9q8qqQ8frFve9Fve9Ff0dmeaabaqaceGacaGaaeWabaWaaeaaea aakeaajuaGcaWGsbGaeyicI4SaamOuamaaCaaabeqcfasaaiaad2ga cqGHxdaTcaWGTbaaaaaa@4293@ is the correlation vector between sampling points. The calculated method of correlation between two points based on Gaussian model:

r i,j =S(ξ, θ i , θ j )= k=1 n exp( ξ k | θ i,k θ j,k |) MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D aebbnrfifHhDYfgiuL2zLjhasaacH8srps0lbbf9q8WrFfeuY=Hhbb f9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq =He9q8qqQ8frFve9Fve9Ff0dmeaabaqaceGacaGaaeWabaWaaeaaea aakeaajuaGcaWGYbWaaSbaaKqbGeaacaWGPbGaaiilaiaadQgaaeqa aKqbakabg2da9iaadofacaGGOaqedmvETj2BSbacfaGae8NVdGNaai ilaiab=H7aXnaaBaaajuaibaGaamyAaaqcfayabaGaaiilaiab=H7a XnaaBaaajuaibaGaamOAaaqabaqcfaOaaiykaiabg2da9maarahaba GaciyzaiaacIhacaGGWbGaaiikaiabgkHiTiab=57a4naaBaaajuai baGaam4AaaqcfayabaWaaqWaaeaacqWF4oqCdaWgaaqaaKqbGiaadM gajuaGcaGGSaqcfaIaam4AaaqcfayabaGaeyOeI0Iae8hUde3aaSba aKqbGeaacaWGQbGaaiilaiaadUgaaeqaaaqcfaOaay5bSlaawIa7ai aacMcaaeaacaWGRbGaeyypa0JaaGymaaqaaiaad6gaaiabg+Givdaa aa@6DCF@  (3)

Where ξ k R n MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D aebbnrfifHhDYfgiuL2zLjhasaacH8srps0lbbf9q8WrFfeuY=Hhbb f9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq =He9q8qqQ8frFve9Fve9Ff0dmeaabaqaceGacaGaaeWabaWaaeaaea aakeaajuaGcaaH+oWaaSbaaKqbGeaacaWGRbaajuaGbeaacqGHiiIZ caWGsbWaaWbaaeqajuaibaGaamOBaaaaaaa@41C8@ denotes the parameter of Gaussian model, which can be obtained by maximum likelihood estimation.

The error coupling model is:

ΔX= [ ( P X ) T P X ] 1 ( P X ) T ΔP MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbmqaaeWacmGaaiaabeqaamaabaabaaGcbaqcfaOaeyiLdq Kaamiwaiabg2da9maadmaabaWaaeWaaeaadaWcaaqaaiabgkGi2kaa dcfaaeaacqGHciITcaWGybaaaaGaayjkaiaawMcaamaaCaaabeqcfa saaiaabsfaaaqcfa4aaSaaaeaacqGHciITcaWGqbaabaGaeyOaIyRa amiwaaaaaiaawUfacaGLDbaadaahaaqcfasabeaacqGHsislcaaIXa aaaKqbaoaabmaabaWaaSaaaeaacqGHciITcaWGqbaabaGaeyOaIyRa amiwaaaaaiaawIcacaGLPaaadaahaaqcfasabeaacaqGubaaaKqbak abfs5aejaadcfaaaa@54AC@  (4)

Where ΔX= [ Δ θ 1 Δ θ 2 Δ θ 3 ] T MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9 vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=x fr=xb9adbiqaaeGabiqaaiaabeqaamaabaabaaGcbaGaeyiLdqucfa Oaamiwaiabg2da9iaacUfafaqabeqadaaabaGaeyiLdqeccaGae8hU de3aaSbaaKqbGeaacaaIXaaabeaaaKqbagaacqGHuoarcqWF4oqCda Wgaaqcfasaaiaaikdaaeqaaaqcfayaaiabgs5aejab=H7aXnaaBaaa juaibaGaaG4maaqabaaaaKqbakaac2fadaahaaqcfasabeaacaWGub aaaaaa@4AE3@  and ΔP= [Δx,Δy,Δz] T MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D aebbnrfifHhDYfgiuL2zLjhasaacH8srps0lbbf9q8WrFfeuY=Hhbb f9v8qqaqFr0xc9pk0xbba9q8WqFfea0=yr0RYxir=Jbba9q8aq0=yq =He9q8qqQ8frFve9Fve9Ff0dmeaabaqaceGacaGaaeWabaWaaeaaea aakeaajuaGcqGHuoarcaWGqbGaeyypa0Jaai4waiabgs5aejaadIha caGGSaGaeyiLdqKaamyEaiaacYcacqGHuoarcaWG6bGaaiyxamaaCa aabeqcfasaaiaadsfaaaaaaa@49D0@ .

Result

Since the spatial position of the robot is mainly determined by the first three joints, the gratings are only mounted on the first three. The specification of the robot used in this paper is KUKA KR210 extra. Take a 100mm×100mm×100mm grid in the robot operation space, and randomly generate 30 sampling points and 15 test points in this grid. The results of the experiment of the robot are shown in Figure 2. And, the position errors of the robot are obtained under the three conditions no compensation, compensation without feedback and compensation based on joints feedback. As can be seen from the diagram, after eliminating the backlash of the joints by joint feedback, the maximum positioning error of the robot in the test space is 0.14mm, with an average value of 0.078mm. We also have detected that the backlash will lead to a positioning error of 0.1mm~0.25mm in Cartesian space.

Figure 2 Absolute positioning accuracy in three different compensation methods.

Conclusion

Because of the existence of joint backlash, it is impossible to further improve the absolute localization accuracy of the robot by conventional accuracy compensation methods. In this paper, a method of robot accuracy compensation based on joint feedback is proposed. From the preliminary experimental results, the absolute localization accuracy of robot is improved, which demonstrated that this method is effective. The follow-up should verify this method in a larger space and the joints of A4, A5, A6 should be mounted with gratings, since the backlashes of joints A4, A5, A6 still exist.

Acknowledgments

 No financial interest.

Conflict of interest

Author declares that there is none of the conflicts.

References

Creative Commons Attribution License

©2017 Liu, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.